Hodgkin-Huxley方程式

Hodgkin-Huxley Equations

概略

Alan Lloyd Hodgkin (1914-1998)とAndrew Fielding Huxley (1917- )は、ともにイギリスの電気生理学者である。イカの巨大軸索における活動電位の発生と伝搬を測定し、その解析から現在の電気生理学の基礎となる概念を生み出すとともに、興奮性細胞(神経細胞、心筋、骨格筋)の電気現象を定量的に扱う道を開いた。HodgkinとHuxleyは、電気生理学の基礎を築いた功績により、同じく電気生理学者のJohn Carew Ecclesとともに、1963年のノーベル医学・生理学賞を受賞している。

HodgkinとHuxleyの業績の意義は次のように要約できる。

1. 活動電位発生時に、ナトリウムイオン(Na⁺)とカリウムイオン(K⁺)が、細胞膜の別々の通路を通ることを示した。この発見はイオンチャネルの存在を予測するものであり、その後のイオンチャネル研究の源となった。なお当時の論文では、イオンチャネル・チャネルという用語は用いられておらず、コンダクタンスという用語が使用されている。
2. Na⁺チャネル、K⁺チャネルが開閉する非線形な動態を微分方程式を含む数式で表した。これらの式はまとめてHodgkin-Huxley方程式と呼ばれる。
3. Na⁺チャネル、K⁺チャネルおよびleakチャネルを示す数式を組み合わせ、活動電位の発生・伝播を数値的に再現した。現在行われている興奮性細胞の電位シミュレーションは、要素が増えるなどして複雑になっているが基本は変わらない。

m³hとn⁴

活動電位の発生時に、Na⁺、K⁺、Cl^-が個別に膜を透過することを見いだしたHodgkinとHuxleyは、voltage clamp法を用いてそれぞれのイオンの通りやすさ(コンダクタンス、抵抗の逆数)を測定した。
K⁺チャネルの場合、細胞膜を脱分極させるとK⁺チャネルは開いていき定常状態(steady state)に達する。 実験結果より、定常状態のコンダクタンスおよび定常状態に達するまでのコンダクタンス変化の速度定数が、電位によって一定であることが見いだされた。 2つの状態(OpenとClosed)があり状態間移行の速度定数が一定な系は、簡単な微分方程式を用いて表すことができる(two-stateモデル)。この考え方を適用したところ、4つの独立したゲートがあり、4つすべてが開いた時に電流が流れる、とすると実験データに合致することが示された。K⁺電流は、次の式で表される。

$I_{K}=G_{K}n^{4}(v-E_{K})\,$

G_Kは最大コンダクタンス、nはゲートが開いている確率、vは電位、E_KはK⁺の平行電位。

Na⁺チャネルの場合は、脱分極するとチャネルは開くが、不活性化により閉じていく。開くゲートmだけでなく閉じるゲートhを考える事により、不活性化を説明する事が出来る。実験データより、3つの活性化ゲートと1つの不活性化ゲートが想定された。

$I_{Na}=G_{Na}m^{3}h(v-E_{Na})\,$

G_Naは最大コンダクタンス、mは活性化ゲートが開いている確率、hは不活性化ゲートが開いている確率、E_NaはNa⁺の平行電位。

電位に依存しないリークチャネルは、

$I_{l}=G_{l}(v-E_{l})\,$

G_lは最大コンダクタンス、E_lはリーク電流の平行電位。

(未完成)

電位変化

(未完成)

Two-state model: 基礎的な考え方*

OpenとClosedの2つの状態がある系で、他の状態に移る率が一定の場合、次の性質がある。

• 指数関数的に変化する
• 近づく値、変化の速さは、初期条件に依存しない
  

以下は数式的な説明。

---

OpenとClosedの2つの状態がある系を考え、Openの状態にある確率をpとする。Closedの状態にある確立は, 1-pとなる。pは時刻tの関数であり、p(t)と表わすとする。

いま状態Closedから状態Openへ移っていく単位時間での割合(速度定数、rate constant)をαとし、状態Openから状態Closedへの速度定数をβとする。 p(t)の時間的経過を表わす微分方程式は、

${\frac {dp(t)}{dt}}=\alpha (1-p(t))-\beta p(t)\,$

と表される。αとβが定数であるとして、t =∞で定常状態になれば、

${\frac {dp(\infty )}{dt}}=\alpha (1-p(\infty ))-\beta p(\infty )=0\,$

であるから、

$p(\infty )={\frac {\alpha }{\alpha +\beta }}\,$

となる。また微分方程式は解析的に解けて、

$p(t)=\left(p(0)-{\frac {\alpha }{\alpha +\beta }}\right)e^{-(\alpha +\beta )t}+{\frac {\alpha }{\alpha +\beta }}\,$

となる。 これらの式は次のことを示している。

• p(t)は指数関数的にp(∞)に近づいていく。
• その時定数(time constant)τは1/(α+β)である。
• これらの値p、τは、初期値p(0)に依存しない。

さらに、

$q(t)=p(t)-{\frac {\alpha }{\alpha +\beta }}\,$

と表すとすると、

$q(t)=q(0)e^{-(\alpha +\beta )}\,$

とより単純な形式となる。この関係は微分方程式の数値計算でよく用いられる。

電位固定法: 基礎となった技術*

Hodgkin-Huxley以前に、電気生理学の実験が行われていなかったわけではない。電流と電位変化に関する研究はかなり多く行われていた。しかしながら、細胞にはいろいろなイオンチャネルを通して電流が流れるため、細胞の電位vと外部から流す電流I_extの間の関係は、単純ではない。そこでHodgkinとHuxleyは、 voltage clamp(電位固定法)を用いて、コンダクタンスの変化を測定して解析した。 voltage clampは1940年代にアメリカの生物物理学者Kenneth Cole (1900 - 1984)らにより開発された。

以下は数式的な説明。

---

外部より電流I_extを流した場合、電位の変化は、次の式で示される。

${\frac {dv}{dt}}=-{\frac {1}{C}}\left(\sum _{X}G_{X}(v-E_{X})-I_{ext}\right)$

この式から、I_extと電位との関係を理解する事は難しい。しかしvが一定となるような外部電流I_clampを流すと、左辺は0となるため、

$I_{clamp}=\sum _{X}G_{X}(v-E_{X})\,$

という関係が得られる。もし溶液の組成を工夫しさらにチャネルのブロッカーなどを用いて、イオンチャネルAを流れる電流を単離して測れたとすると、

$I_{clamp}=G_{A}(v-E_{A})\,$

となる。これはOhmの法則である。ここでI_clampは実験の測定値、vは実験の設定値、E_Aは実験条件で定まる定数なので、

$G_{A}={\frac {I_{clamp}}{v-E_{A}}}\,$

の関係式を用いて、実験データよりイオンチャネルAのコンダクタンスG_Aを算出できることになる。

HHモデルに対する批判

1. ゲート電流
2. Single-channel recording
3. Markovモデル
4. Fractalモデルとの論争

現在におけるHHモデル

(未完成)

HH方程式を使ってみる

HH方程式の数値計算は、ノートパソコンでも十分行う事が出来る。

References

未完成