lunelukkio.com / 公開エリア › 神経科学

イオンチャネルの状態遷移

開口・閉口の確率的切り替わりを数学でとらえる — Q行列への道筋

Section 1

イオンチャネルとは

イオンチャネルは細胞膜に埋め込まれたタンパク質の孔（チャネル）であり、 Na+・K+・Ca2+ などのイオンを選択的に通過させる。 神経細胞が活動電位を発生させるしくみの根幹を担っており、 筋収縮・心拍リズム・感覚受容にも深く関わる。

チャネルの特徴的な性質は、確率的な開閉である。 ある瞬間には開いている（Open）が、次の瞬間には閉じている（Closed）—— この切り替わりはランダムに起こり、古典的な化学反応速度論だけでは扱いにくい。 そこで登場するのが確率過程の数学、とくに連続時間マルコフ連鎖である。

細胞外 細胞内 Closed イオン通過不可 β（開口速度） α（閉口速度） Open イオンが流れる

図1. イオンチャネルの Closed（閉）と Open（開）の模式図。閉口速度 α と開口速度 β でランダムに切り替わる。

Section 2

2状態モデル — Open ⇌ Closed

最も単純なモデルは、チャネルが「Open（O）」か「Closed（C）」の 2つの状態しか取らないと仮定するモデルである。

状態は連続時間でランダムに切り替わる。キーとなる量は遷移速度定数（rate constant）である。

遷移速度定数意味
Open → Closed α（単位: s-1） Open 状態から微小時間 dt の間に Closed へ移る確率 ≈ α · dt
Closed → Open β（単位: s-1） Closed 状態から微小時間 dt の間に Open へ移る確率 ≈ β · dt
α/β の慣用について: 本稿では Colquhoun & Hawkes の慣用に従い、α を閉口速度（O→C）、β を開口速度（C→O）と定義している。Hodgkin-Huxley モデルの慣用（α が活性化方向 C→O）とは逆向きなので注意。
Open O Closed C α β

図2. 2状態モデルの状態遷移図。α と β は遷移速度定数（s-1）。

マルコフ性（無記憶性）: 次の状態への遷移確率は、現在の状態のみに依存し、過去の履歴とは無関係。 たとえば「Open に 10 ms いた」か「2 ms しかいない」かにかかわらず、 次の dt 秒間に Closed へ移る確率は等しく α · dt。 これが連続時間マルコフ連鎖の核心的な仮定である。

Open 状態の滞在時間は指数分布に従う

マルコフ性と速度定数 α から、Open 状態に留まる時間 T は指数分布 T ∼ Exp(α) に従うことが示せる。平均滞在時間は 1/α [秒]。 同様に Closed の平均滞在時間は 1/β [秒]。 これは単一チャネル記録の実験で直接測定可能な量である。

Section 3

複数状態への拡張とQ行列の導入

実際のイオンチャネルは 2 状態よりも複雑な構造を持つ。 たとえば Closed 状態にも複数のサブ状態（C₁, C₂ …）が存在したり、 不活性化状態（Inactivated）を経由したりする。

C₂ Closed 2 k₂₁ k₁₂ C₁ Closed 1 β α O Open δ γ I Inactivated

図3. 多状態モデルの例（C₂ ⇌ C₁ ⇌ Open ⇌ Inactivated）。各矢印に速度定数が対応する。

速度定数の行列表現 — Q行列

状態が n 個になると、すべての遷移速度定数を管理するのに行列表現が有効である。 Q行列（生成行列、generator matrix）は n × n の行列で、 成分 qij が「状態 i から状態 j への遷移速度定数」を表す。

Q行列の3つの性質:
  • 非対角成分 qij ≥ 0（i ≠ j）— 別状態への遷移速度
  • 対角成分 qii = −∑j≠i qij ≤ 0 — 外へ出る総速度（符号反転）
  • 各行の和 = 0 — 確率の保存（チャネルは消えも現れもしない）

2状態モデル（O ⇌ C）の Q行列は次の形になる:

Q = [ ] −α α β −β ← Open 行 ← Closed 行

図4. 2状態モデルのQ行列。対角成分は負（外へ出る総速度）、非対角は正（i から j への遷移速度）。

Q行列があれば、「時刻 0 に状態 i にいたとき、時刻 t に状態 j にいる確率」を すべての状態ペアについて一括で計算できる。この計算に使われる手法が行列指数関数 eQt である。

行列指数関数の数学的基礎へ →