NSM-031

チャネルを流れる電流 — 開状態確率から電流時間経過へ（式28・29・30）

作成日: 2026-06-03 / §5 電流変換 / Mathematics in Neuroscience プレゼン準備

このノートの中心命題:
マルコフ連鎖モデルで得られた状態確率 $p(t)$ に、コンダクタンスベクトル $\mathbf{v}$ を乗じて駆動力 $(V - V_\mathrm{rev})$ とチャネル数 $N$ でスケールすれば、観測電流 $I(t)$ が得られる。スペクトル分解 $e^{Qt} = \sum_i A_i e^{\lambda_i t}$ を使うと、$I(t)$ は定常電流 $I(\infty)$ と $k-1$ 個の指数減衰成分の和として書ける（式28・29・30）。

1. 設定の整理

チャネルは $k$ 個の離散状態を持つ。先頭 $k_A$ 個が開状態 (Open)、残りの $k - k_A$ 個が閉状態 (Closed)である。

記号	意味	次元 / 単位
$k$	全状態数（開 + 閉）	スカラー（整数）
$k_A$	開状態数（添字 $1, 2, \ldots, k_A$）	スカラー（整数）
$Q$	Q 行列（遷移速度行列）	$k \times k$（$\mathrm{s}^{-1}$）
$p_i(t)$	時刻 $t$ に状態 $i$ にいる確率	スカラー（無次元）
$\mathbf{p}(t)$	確率行ベクトル $[p_1(t), \ldots, p_k(t)]$	$1 \times k$
$\gamma_i$	開状態 $i$ のコンダクタンス	pS（ピコジーメンス）
$\mathbf{v}$	コンダクタンスベクトル	$k \times 1$（pS）
$V$	膜電位	mV
$V_\mathrm{rev}$	反転電位	mV
$N$	チャネル数	スカラー（整数）
$\lambda_i$	Q 行列の第 $i$ 固有値（$\lambda_1 = 0$）	$\mathrm{s}^{-1}$
$A_i$	第 $i$ 固有値の射影行列	$k \times k$
$\tau_i$	第 $i$ 指数成分の時定数 $= 1/\|\lambda_i\|$	ms または s

図 1. 状態遷移の模式図。開状態（四角）はコンダクタンス $\gamma_i > 0$ を持ち電流を流す。閉状態（丸）は $\gamma = 0$。

2. コンダクタンスと駆動力

開状態 $i$ のチャネルを通じて流れる電流は、オームの法則の生物版として次のように書ける。

$$I_i = \gamma_i \cdot (V - V_\mathrm{rev})$$

$V_\mathrm{rev}$ は反転電位（reversal potential）と呼ばれ、電流がちょうどゼロになる膜電位である。 $(V - V_\mathrm{rev})$ を駆動力（driving force）と呼ぶ。駆動力が正なら電流は一方向に流れ、負なら逆向きに流れる。

図 2. 電流–電圧関係（I-V 曲線）。傾きがコンダクタンス $\gamma_i$、ゼロ交差点が $V_\mathrm{rev}$。$(V - V_\mathrm{rev}) > 0$ で内向き電流（負の符号規約によっては逆）。

3. 1チャネルの期待電流

チャネルが確率的に各開状態を占めるとき、時刻 $t$ における単一チャネルの期待電流は、全開状態にわたって「確率 × コンダクタンス」を足した値に駆動力を掛ける。

$$E[I_\text{single}(t)] = (V - V_\mathrm{rev}) \sum_{i=1}^{k_A} \gamma_i \, p_i(t)$$

閉状態のコンダクタンスはゼロなので、閉状態の項は和に現れない。 $N$ 個の独立なチャネルがあれば、全体の期待電流はこれを $N$ 倍する（大数の法則による）。

コンダクタンスベクトル $\mathbf{v}$
$$\mathbf{v} = \underbrace{[\gamma_1, \gamma_2, \ldots, \gamma_{k_A}]^\top}_{k_A \text{ 個（開状態）}}, \underbrace{[0, 0, \ldots, 0]^\top}_{k - k_A \text{ 個（閉状態）}}$$ $\mathbf{v}$ は「開状態だけを拾うフィルター」として機能し、$\mathbf{p}(t) \mathbf{v}$ が期待コンダクタンスを与える。

4. 式(28) — $I(t)$ の指数和表現

$Q$ 行列のスペクトル分解（NSM-028 参照）を使うと、確率ベクトルは $$\mathbf{p}(t) = \mathbf{p}(0) e^{Qt} = \mathbf{p}(0) \sum_{i=1}^{k} A_i e^{\lambda_i t}$$ と書ける。$Q$ 行列は各行の和がゼロという制約から必ず $\lambda_1 = 0$ の固有値を持つ（NSM-021 参照）。この $\lambda_1 = 0$ 成分が定常電流 $I(\infty)$ を生み、残りの $k-1$ 個（$i = 2, \ldots, k$）が時間的に減衰する。

したがって電流の時間経過は次の式(28)で表される。

(28) $$I(t) = I(\infty) + \sum_{i=2}^{k} b_i \exp\!\left(-\frac{t}{\tau_i}\right)$$ ここで $\tau_i = 1/|\lambda_i|$

$i = 2$ から始まる理由: $\lambda_1 = 0$ に対応する項は時定数が無限大（$\tau_1 = \infty$）であり、時間によって変化しない定数項 $I(\infty)$ に吸収されるからである。

図 3. $I(t)$ の時間経過と各指数成分の分解。複数の指数成分が重なって $I(\infty)$ に収束する。

5. 式(29) — 係数 $b_i$ の行列表現

各指数成分の係数 $b_i$ は次の式(29)で与えられる。

(29) $$b_i = N (V - V_\mathrm{rev}) \, \mathbf{p}(0) \, A_i \, \mathbf{v}$$

各因子の意味を順に確認する。

因子	型	意味
$N$	スカラー	チャネル数。全チャネルの合計電流にするための倍率。
$(V - V_\mathrm{rev})$	スカラー（mV）	駆動力。電圧固定実験では定数として外に出せる。
$\mathbf{p}(0)$	$1 \times k$ 行ベクトル	初期確率分布。どの状態からスタートするかを指定。
$A_i$	$k \times k$ 行列	固有値 $\lambda_i$ に対応する射影行列（スペクトル分解 $e^{Qt} = \sum_i A_i e^{\lambda_i t}$ の係数）。
$\mathbf{v}$	$k \times 1$ 列ベクトル	コンダクタンスベクトル。開状態のみ $\gamma_i > 0$、閉状態は 0。

$\mathbf{v}$ は「開状態を拾うフィルター」
$\mathbf{v}$ の先頭 $k_A$ 成分のみが非ゼロである。そのため行列積 $(\cdots) \mathbf{v}$ は、行列の先頭 $k_A$ 列だけを拾い出し、閉状態に対応する列を自動的にゼロにする。

図 4. $b_i$ の行列積の次元フロー。$(1 \times k) \times (k \times k) \times (k \times 1) = $ スカラー。

6. 式(30) — 係数 $b_i$ の成分展開

式(29)の行列積を成分で書き下すと式(30)が得られる。

(30) $$b_i = N (V - V_\mathrm{rev}) \sum_{r=1}^{k} \sum_{j=1}^{k_A} p_r(0) \, \gamma_j \, a_{rj}^{(i)}$$

各記号の意味:

$a_{rj}^{(i)}$: 射影行列 $A_i$ の $(r, j)$ 成分
$r$: 初期状態の添字（$1$ から $k$ まで、開状態も閉状態も含む）
$j$: 開状態の添字（$1$ から $k_A$ まで、コンダクタンスが非ゼロ）
$p_r(0)$: 時刻 $0$ に状態 $r$ にいる確率
$\gamma_j$: 開状態 $j$ のコンダクタンス

式(30)の直感
二重和は「全ての初期状態 $r$ から全ての開状態 $j$ への寄与の組み合わせ」を足している。 $p_r(0)$ は「最初の確率の重み」、$a_{rj}^{(i)}$ は「固有モード $i$ が状態 $r$ から状態 $j$ をどれだけ経由するか」の指標、$\gamma_j$ は「状態 $j$ がどれだけ電流を流すか」の大きさである。これらの積を全ての $(r, j)$ ペアについて足すことで、固有モード $i$ が電流全体に与える寄与 $b_i$ が決まる。

式(29)と式(30)の対応関係

式(29)を成分表記で展開すると次のようになり、式(30)と一致することがわかる。

$$\mathbf{p}(0) A_i \mathbf{v} = \sum_{r=1}^{k} p_r(0) [A_i \mathbf{v}]_r = \sum_{r=1}^{k} p_r(0) \sum_{j=1}^{k} a_{rj}^{(i)} v_j$$

$v_j = \gamma_j$（$j \leq k_A$）、$v_j = 0$（$j > k_A$）なので、第二の和は $j = 1, \ldots, k_A$ のみ残り、式(30)の二重和に一致する。

7. $N$ チャネルと単一チャネル平均の等価性

$N$ 個の独立かつ同一のチャネルが存在するとき、観測される総電流は次のように書ける。

$$I_\text{total}(t) = N \cdot E[I_\text{single}(t)] = N (V - V_\mathrm{rev}) \mathbf{p}(0) e^{Qt} \mathbf{v}$$

これは「単一チャネルのステップ実験を $N$ 回繰り返して平均した電流時間経過」とも等価である。理由は大数の法則により、$N$ 回試行の標本平均が期待値に収束するからである。

実験との対応
マクロパッチ（多数チャネルを含むパッチ）での電圧ステップ実験は、単一チャネルの確率的なふるまいを大数の法則で平均した結果を測定している。式(28)の指数和は、その巨視的な電流波形の数学的記述である。

8. まとめ — $p(t)$ から $I(t)$ への変換の全体像

図 5. $p(t)$ から $I(t)$ への変換フロー。3ステップで式(28)の指数和が導かれる。

核心の流れ:
(1) Q 行列のスペクトル分解 → $p(t) = \sum_i A_i e^{\lambda_i t}$ を成分で持つ
(2) コンダクタンスベクトル $\mathbf{v}$ で開状態を抽出 → 各固有モードの「電流への寄与率」を計算
(3) 駆動力 $(V - V_\mathrm{rev})$ とチャネル数 $N$ でスケール → $b_i$ 決定
(4) $\lambda_1 = 0$ 成分は $I(\infty)$ に、$\lambda_i \neq 0$ 成分は指数減衰 $b_i e^{-t/\tau_i}$ に対応

確信度ラベルと参考文献

🟢 式(28)(29)(30) の定式化 — Colquhoun & Hawkes (1977, 1981) の原典に基づく標準的な枠組み
🟢 $\lambda_1 = 0$ が定常成分に対応すること — Q 行列の性質（行和ゼロ）から厳密に導かれる（NSM-021）
🟢 $\mathbf{v}$ の構成（開状態のみ非ゼロ）— Colquhoun & Hawkes の符号規約に準拠
🟡 「$N$ 倍 = $N$ 回平均」の等価性 — 大数の法則による（数学的には i.i.d. 仮定が必要。実際にすべてのチャネルが同一の Q 行列を持つかは仮定）

参考文献

Colquhoun, D., & Hawkes, A. G. (1977). Relaxation and fluctuations of membrane currents that flow through drug-operated channels. Proceedings of the Royal Society B, 199(1135), 231–262.
Colquhoun, D., & Hawkes, A. G. (1981). On the stochastic properties of single ion channels. Proceedings of the Royal Society B, 211(1183), 205–235.
Colquhoun, D., & Hawkes, A. G. (1995). A Q-matrix cookbook. In Single-Channel Recording (2nd ed.), Ch. 20. Plenum Press.

[NSM-011] 行列指数関数 $e^{Qt}$ と Q 行列 — $e^{Qt}$ の対角化とスペクトル分解の計算法
[NSM-021] Kolmogorov 前進方程式と平衡分布 $\pi$ — $\lambda_1 = 0$ が生じる理由と定常分布 $\pi$ の意味
[NSM-028] 固有値の計算 — Q 行列から時定数へ — $\lambda_i$ と $\tau_i = 1/|\lambda_i|$ の求め方
[NSM-003] 観測 dwell time のスペクトル分解 — 係数 $b_i$ の符号と dwell time 分布のピーク条件

NSM-031 / 作成: 2026-06-03 / Mathematics in Neuroscience プレゼン準備