マクロ電流 vs 単一チャネル滞在時間ヒストグラム
両者は何が違うのか — そして、なぜ我々は単一チャネルを使うのか
結論を先に:
マクロ電流 I(t) と 単一チャネル滞在時間ヒストグラム f(t) は、形は違うが 同じ Q 行列の固有値 $\lambda_i$ を反映している。
教科書はマクロ電流で書くが、状態数推定には滞在時間ヒストグラムのほうが圧倒的に情報量が多い。
1. 3つの「似て非なるもの」を区別する
| 名前 | 横軸 | 縦軸 | データ源 |
|---|
| A. 単一チャネル電流トレース | 時間 (秒) | 電流 (pA) | 1個のチャネル |
| B. マクロ電流トレース (教科書) | 時間 (秒) | 電流 (nA, pA) | 多数のチャネルの和 |
| C. 滞在時間ヒストグラム (我々の主役) | 継続時間 (ms) | 頻度 | 1個のチャネルの統計 |
A・B は「時間軸上での電流の動き」を見る。C は「どの長さの滞在時間が、どれだけよく起きるか」を集計したもの — 横軸が時刻ではなく、滞在の長さ。
2. A: 単一チャネル電流トレース
図1. 単一チャネル電流トレース。開 (-2 pA) と閉 (0 pA) を不規則に行き来する方形波
3. B: マクロ電流トレース (教科書でよく見る形)
図2. マクロ電流。1000個などの大集団の合計。これが教科書でよく見る曲線
3.1 物理的に何が起きているか
例えば 1000 個の AChR があるとき:
- $t=0$ (神経伝達物質到着): 全部が閉状態
- $t=1\,\text{ms}$: 30% が開いた
- $t=3\,\text{ms}$: 70% が開いた (peak)
- $t=10\,\text{ms}$: 50% が開いている
- $t=50\,\text{ms}$: 5% しか開いていない
電流 $I(t) = (\text{開いているチャネル数}) \times (\text{1チャネル電流})$
つまり開状態にいるチャネルの割合 $p_O(t)$ の時間変化がマクロ電流の波形になる。
4. C: 滞在時間ヒストグラム (我々の主役)
電流トレースから、各 open/closed の継続時間を取り出して集計する:
図3. 電流トレースから滞在時間ヒストグラムを作る流れ。Step 2 の横軸は「時刻」ではなく「滞在の長さ (ms)」
5. なぜ両者が同じ Q行列を見ているのか
5.1 マクロ電流の数式
$p_O(t)$ (時刻 $t$ で開状態にいる確率) は Kolmogorov 方程式の解として:
$p_O(t) = \pi_O + \sum_i a_i\, e^{\lambda_i t}$
マクロ電流は $I(t) = N \cdot i_{\text{single}} \cdot p_O(t)$ なので、形は指数の和。$\lambda_i$ は Q行列の固有値。
5.2 滞在時間ヒストグラムの数式
$f(t) = \sum_i b_i\, |\lambda_i|\, e^{\lambda_i t}$
同じ $\lambda_i$ — Q行列の固有値が出てくる。
結論: 両者は同じ Q行列の固有値 $\lambda_i$ を共有する。
- マクロ電流: 時間軸上で観察される「緩和の時定数 $\tau_i = 1/|\lambda_i|$」
- 滞在時間ヒストグラム: 継続時間軸上で観察される「指数成分」
形は違うが、情報源は同一。
6. では、なぜ単一チャネルを使うのか
| 観点 | マクロ電流 (B) | 滞在時間ヒストグラム (C) |
|---|
| 測定の容易さ | 簡単 (whole-cell patch) | 難しい (single channel patch) |
| 得られる情報 | 平均的な挙動のみ | 個別イベントの完全な分布 |
| 状態数推定 | 難しい (時定数の分離が難しい) | 得意 (mixture vs Erlang が一目) |
| 並列 vs 直列の区別 | 困難 | peak の有無で直接見える |
| ノイズへの強さ | 1本のフィットに依存 | 多数イベントで統計的に強い |
7. 教科書が B を使う理由
- 歴史: Hodgkin-Huxley (1952) はマクロ電流から始まった。Single channel patch (Neher & Sakmann) は1976年・Nobel 1991。
- 測定の容易さ: 大きな電流が出るのでオシロで見える
- 薬理学: 薬の効果は IV曲線・時定数で議論される
- 直感性: 「ON時定数 / OFF時定数」が学習しやすい
8. プレゼンへの組み込み方
発表中、§2 (single channel recording) のあたりで以下のメッセージを挿入すると、聴衆 (数学者) も既存知識 (教科書のマクロ電流) と接続できる:
「教科書ではマクロ電流の指数減衰でこの話が出ます。しかし状態数を当てるには、単一チャネルの滞在時間ヒストグラムを使うほうが情報量が圧倒的に多い。だからこの発表はそこから出発します。」
9. まとめ — 1分でわかる対比
| マクロ電流 | 滞在時間ヒストグラム |
|---|
| 横軸 | 時間 (時刻) | 継続時間 (長さ) |
| 縦軸 | 電流 | 頻度 |
| 形 | peak + 指数減衰 | 単調減少 or peak付き |
| 情報源 | 同じ Q 行列の固有値 $\lambda_i$ |
| 状態数推定 | 難しい | 得意 |
作成日: 2026-05-18 / プロジェクト: Mathematics in Neuroscience プレゼン準備 / §3 補足