Hodgkin-Huxley方程式

Hodgkin-Huxley Equations

概略

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

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

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

m³hとn⁴

HodgkinとHuxleyは、電位固定法（voltage-clamp)を用いて活動電位に伴うNa⁺とK⁺のコンダクタンス（通りやすさ、抵抗の逆数)の変化を定量的に解析し、Na⁺とK⁺には別々の通り道があることを示した。そしてNa⁺とK⁺のコンダクタンスが電位に依存的なゲート（gate）により開閉されると考えた。

• Na⁺チャネルは3つの活性化ゲートmと不活性化ゲートhにより開閉される。
• K⁺チャネルは4つの活性化ゲートnにより開閉される。

m、h、nは、ゲートが開いている割合を示す値で、単純なTwo-state modelに従う。mとnは、静止時に閉じており脱分極した時に開く。一方、hは静止時に開き脱分極時に閉じる。mとnではなく、m³およびn⁴としたのは、主に電流の立ち上がりの形をよく再現するためである。

電流はコンダクタンスと電圧に比例する（I = GV; Ohmの法則）。電圧の大きさは、細胞膜内外のイオン濃度差による電位（平衡電位）を補正しなくてはならない（v-E_X）。 従って、 Na⁺とK⁺により担われる電流I_NaとI_Kは、Na⁺とK⁺の最大コンダクタンスをそれぞれ G^max_Na、G^max_K 、平衡電位をE_Na、E_Kとすると、 

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

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

という式で表される。また主にクロライドイオン（Cl^-）を通し静止電位保持に主要な役割を果たすリーク(leak)チャネルは、コンダクタンスは電位に依存しないとして、

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

と表される。

---

m、h、nはTwo-stateモデルに従う値である。 開く速度定数αと閉じる速度定数βはいずれも電位に依存する。 HodgkinとHuxleyはm、h、nのそれぞれについていろいろな電位での αとβの値を実験的に測定し、一連のイカ巨大軸索を用いた実験値を便宜的に次の数式で表した。これらの式は何らかの理論から導きだされたものではない。なお注意すべき点として、以下の式で用いられている変数vは、膜電位V_mそのものではなく静止電位を基準とした電位を意味している。

$\alpha _{m}={\frac {0.1(-v+25)}{\exp \left({\frac {-v+25}{10}}\right)-1}}$

$\beta _{m}=4\exp \left({\frac {-v}{18}}\right)$

$\alpha _{h}=0.07\exp \left({\frac {-v}{20}}\right)$

$\beta _{h}={\frac {1}{\exp \left({\frac {-v+30}{10}}\right)+1}}$

$\alpha _{n}={\frac {0.01(-v+10)}{\exp \left({\frac {-v+10}{10}}\right)-1}}$

$\beta _{n}=0.125\exp \left({\frac {-v}{80}}\right)$

電位変化

電位固定実験から得られた上記の式を用いると、電位の変化に応じて変化するNa⁺チャネルとK⁺チャネルのコンダクタンスを計算することができる。 一方、細胞の電位は、電流の入出により変化する。もし電位の初期条件と電流の時間経過がわかっていれば、電位を計算する事が出来る。 これらを組み合わせることにより、細胞の電位の変化を計算できる。この関係は、現在広く行われている興奮性細胞の電位のシミュレーションの基本である。

以下は数式的な説明。

---

細胞をキャパシタ（コンデンサ）と考える。細胞外液の電位を0とする。細胞膜の電位をv、蓄えられる電荷をQ、静電容量(capacitance)をCとすると、

$Q=Cv\,$

の関係が成り立つ。両辺を時間tで微分すると、

${\frac {dQ}{dt}}=C{\frac {dv}{dt}}\,$

となる。一方、流れる電流Iはすなわち電荷量Qであるから、

$I=-{\frac {dQ}{dt}}\,$

ここで右辺にマイナス符号が付くのは、細胞に流れ込む電流を負とする習わしによる。

電流Iは、チャネルを流れる電流の和であるから、

$I=I_{Na}+I_{K}+I_{leak}\,$

すなわち、

$I=G_{Na}^{max}m^{3}h(v-E_{Na})+G_{K}^{max}n^{3}(v-E_{K})+G_{leak}(v-E_{leak})\,$

と表され、電流Iと電位Vの関係を示す微分方程式に合わせると、

${\frac {dv}{dt}}=-{\frac {1}{C}}\left(G_{Na}^{max}m^{3}h(v-E_{Na})+G_{K}^{max}n^{3}(v-E_{K})+G_{leak}(v-E_{leak})\right)\,$

となる。

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

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)$

Xはいろいろな種類のチャネルを示している。

この式から、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を算出できることになる。

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 )t}\,$

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

なお、このTwo-stateモデルはn-stateモデルに拡張することが可能であり、拡張された式はMarkov過程のシミュレーション等に用いられている。

HHモデルに対する批判

HHモデルは、比較的少ない数のパラメータで神経軸索の活動電位の発生と伝播を示す事に成功した。しかしその後、イオンチャネルの存在が明らかになり、いろいろな測定が可能になって来ると、HHモデルでは説明できない事が見つかって来た。

1. ゲート電流(gating current): チャネルタンパクの動きを反映すると考えられるゲート電流の電位依存性は、電流の電位依存性よりも過分極側にずれており、チャネルが最終的に開く過程は電位依存的ではないと推測された。 
2. Single-channel recordin: Na⁺チャネルの開口時間は、HHモデルで予測されるよりも短く、脱分極から遅れてチャネルが開く場合が観察された。
3. Markov過程モデル: HHモデルよりもより複雑な過程を示すことの出来るMarkov過程モデルの方が、より詳細なチャネルの性質や、変異による性質の変化を表すことが出来た。
   

現在におけるHHモデル

HHモデルは、チャネルの開閉特性を比較的少ないパラメータでかなり正確に表すことができる点で、今でも高く評価されている。現在広く行われている興奮性細胞（神経細胞、心筋細胞、骨格筋細胞）の電位シミュレーションでは、通常、HHモデルもしくはそれに類似のモデルが用いられる。Markov過程モデルの方がチャネル分子の動きと関係付けることが容易でありが、HHモデルの方が計算量の面で効率的である。

HH方程式を使ってみる

HH方程式の数値計算には、ノートPCレベルのコンピュータがあれば十分である。基本的には常微分方程式の数値積分なので、C/C++やFortranなどの通常の言語でプログラムを作成すればよい。しかしそのような計算を目的に開発されて来たNEURONシミュレータを用いると、電位固定でのチャネルの性質や、チャネルを細胞に組み込んだ時の電位の変化などを比較的容易に検証する事が出来る。