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⁴

(未完成)

電位変化

(未完成)

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

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

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

以下は数式的な説明。

---

2つの状態1と2をとる事の出来る系を考え、それぞれの状態にある確率をp₁とp₂とする。p₁とp₂は時刻tの関数であり、p₁(t)とp2(t)と表わすとする。p₁(t)とp₂(t)は確率であるから、

$p1(t)+p2(t)=1\,$

の関係が常に成り立つ。いま状態1から状態2へ移っていく単位時間での割合(遷移率)をαとし、状態2から状態1への遷移率をβとする。 p₁(t)とp₂(t)の時間的経過を表わす微分方程式は、

${\frac {dp_{1}(t)}{dt}}=-\alpha p_{1}(t)+\beta p_{2}(t)$

${\frac {dp_{2}(t)}{dt}}=\alpha p_{1}(t)-\beta p_{2}(t)$

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

${\frac {dp_{1}(\infty )}{dt}}=-\alpha p_{1}(\infty )+\beta p_{2}(\infty )=0\,$

${\frac {dp_{2}(\infty )}{dt}}=\alpha p_{1}(\infty )-\beta p_{2}(\infty )=0\,$

ここで、

$p_{1}(\infty )+p_{2}(\infty )=1\,$

であるから、

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

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

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

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

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

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

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

さらに、

$q_{1}(t)=p_{1}(t)-{\frac {\beta }{\alpha +\beta }}\,$

$q_{2}(t)=p_{2}(t)-{\frac {\alpha }{\alpha +\beta }}\,$

と表すとすると、

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

$q_{2}(t)=q_{2}(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)$

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

未完成