「Hodgkin-Huxley方程式」の版間の差分

編集の要約なし
編集の要約なし
編集の要約なし
 
(5人の利用者による、間の58版が非表示)
1行目: 1行目:
Hodgkin-Huxley Equations
<div align="right"> 
<font size="+1">[http://researchmap.jp/read0012307 井本 敬二]</font><br>
''自然科学研究機構生理学研究所 生体情報研究系 岡崎国立共同研究機構 ''<br>
DOI:<selfdoi /> 原稿受付日:2012年12月7日 原稿完成日:2018年1月4日<br>
担当編集委員:[http://researchmap.jp/michisukeyuzaki 柚崎 通介](慶應義塾大学 医学部生理学)<br>
</div>


== 概略  ==
英語名:Hodgkin-Huxley Equations


Alan Lloyd Hodgkin (1914-1998)とAndrew Fielding Huxley (1917- )は、ともにイギリスの電気生理学者である。イカの巨大軸索における活動電位の発生と伝搬を測定し、その解析から現在の電気生理学の基礎となる概念を生み出すとともに、興奮性細胞(神経細胞、心筋、骨格筋)の電気現象を定量的に扱う道を開いた。HodgkinとHuxleyは、電気生理学の基礎を築いた功績により、同じく電気生理学者のJohn Carew Ecclesとともに、1963年のノーベル医学・生理学賞を受賞している。
{{box|text= [[wj:アラン・ロイド・ホジキン|Alan Lloyd Hodgkin]]と[[wj:アンドリュー・フィールディング・ハクスリー|Andrew Fielding Huxley]]は[[wj:イカ|イカ]]の[[巨大軸索]]の[[活動電位]]を、[[ナトリウムチャネル|Na<sup>+</sup>チャネル]]、[[カリウムチャネル|K<sup>+</sup>チャネル]]の開閉を実験的に測定し、開閉の[[wj:非線形|非線形]]な動態を[[wj:微分方程式|微分方程式]]を含む数式で表わす事に成功した。これをHodgkin-Huxley方程式という。比較的少ない数のパラメータで神経[[軸索]]の[[活動電位]]の発生と伝播を示す事に成功した。チャネルの開閉特性を比較的少ないパラメータでかなり正確に表すことができる点で、今でも高く評価されている。現在広く行われている[[興奮性細胞]]([[神経細胞]]、[[wj:心筋|心筋細胞]]、[[wj:骨格筋|骨格筋]]など)の電位シミュレーションでは、Hodgkin-Huxley方程式や(Hodgkin-Huxleyモデル)、それに類似のモデルが用いられる。
}}


HodgkinとHuxleyの業績の意義は次のように要約できる。
== Hodgkin-Huxley方程式とは ==
 Alan Lloyd Hodgkin (1914-1998)とAndrew Fielding Huxley (1917-2012)は、ともに[[wj:イギリス|イギリス]]の電気生理学者である。イカの巨大軸索における活動電位の発生と伝播を測定し、その解析から現在の電気生理学の基礎となる概念を生み出すとともに、興奮性細胞(神経細胞、心筋細胞、骨格筋細胞など)の電気現象を定量的に扱う道を開いた<ref><pubmed> 12991237 </pubmed></ref><ref>Journal of Physiologyは、Hodgkin &amp; Huxley (1952)論文の60周年を記念して、2012年5月にオンライン版の特別号を出版している。Hogkin &amp; Huxleyおよび関係する論文は、このサイトからリンクされている。</ref> 。HodgkinとHuxleyは、電気生理学の基礎を築いた功績により、同じく電気生理学者の[[wj:ジョン・C・エックルス|John Carew Eccles]]と3人で、1963年の[[wj:ノーベル医学・生理学賞|ノーベル医学・生理学賞]]を受賞している。


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


== ''m''<sup>3</sup>''h''と''n''<sup>4</sup>  ==
== ''m''<sup>3</sup>''h''と''n''<sup>4</sup>  ==


(未完成)  
 HodgkinとHuxleyは、[[#.E9.9B.BB.E4.BD.8D.E5.9B.BA.E5.AE.9A.E6.B3.95:_.E5.9F.BA.E7.A4.8E.E3.81.A8.E3.81.AA.E3.81.A3.E3.81.9F.E6.8A.80.E8.A1.93|電位固定法]](voltage-clamp)を用いて活動電位に伴うNa<sup>+</sup>とK<sup>+</sup>のコンダクタンスの変化を定量的に解析し、Na<sup>+</sup>とK<sup>+</sup>には別々の通り道があることを示した。そしてNa<sup>+</sup>とK<sup>+</sup>のコンダクタンスが[[wj:電位|電位]]に依存的な[[ゲート]](gate)により開閉されると考えた。
 
*Na<sup>+</sup>チャネルは3つの活性化ゲート''m''と不活性化ゲート''h''により開閉される。
*K<sup>+</sup>チャネルは4つの活性化ゲート''n''により開閉される。
 
''m''、''h''、''n''は、ゲートが開いている割合を示し、0と1の間で時間的に変化する変数であり、単純な[[#Two-state_model:_.E5.9F.BA.E7.A4.8E.E7.9A.84.E3.81.AA.E8.80.83.E3.81.88.E6.96.B9|Two-stateモデル]]に従うと仮定されている。''m''と''n''は、静止時に閉じており(すなわち0の値をとる)、脱分極した時に開く(すべてが開けば1となる)。一方、''h''は静止時に開いており脱分極時に閉じる。''m''と''n''ではなく、''m''<sup>3</sup>および''n''<sup>4</sup>としたのは、主に電流の立ち上がりの形をよく再現するためである。


== 電位変化  ==
 [[wj:電流|電流]]はコンダクタンスと電圧に比例する(''I'' = ''GV''; Ohmの法則)。電圧の大きさは、膜電位(v)から細胞膜内外のイオン濃度差による電位(平衡電位; E<sub>X</sub>)を補正しなくてはならない(v-E<sub>X</sub>)。 従って、 Na<sup>+</sup>とK<sup>+</sup>により担われる電流''I''<sub>Na</sub>と''I''<sub>K</sub>は、Na<sup>+</sup>とK<sup>+</sup>の最大コンダクタンスをそれぞれ ''G''<sup>max</sup><sub>Na</sub>、''G''<sup>max</sup><sub>K</sub> 、平衡電位を''E''<sub>Na</sub>、''E''<sub>K</sub>とすると、&nbsp;


(未完成)  
::<math>I_{Na} = G^{max}_{Na} m^3 h (v-E_{Na})\, </math>


== Two-state model: 基礎的な考え方*  ==
::<math>I_{K} = G^{max}_{K} n^4 (v-E_{K})\, </math>


状態が2つの系(例えばOpenとClosed)で、他の状態に移る率が一定の場合、次の性質がある。<br>
という式で表される。また膜電位の変化に対してオームの法則に従うコンダクタンスとしてリーク(leak)チャネルを想定した。これは、


*指数関数的に変化する
::<math>I_{leak} = G_{leak}(v-E_{leak})\, </math>
*近づく値、変化の速さは、初期条件に依存しない<br>


以下は数式的な説明。
と表され、イカ巨大軸索では主に[[wj:塩化物|クロライド]]イオン(Cl<sup>-</sup>)を通し静止電位保持に主要な役割を果たす。


----
----


2つの状態1と2をとる事の出来る系を考え、それぞれの状態にある確率を''p''<sub>1</sub>と''p''<sub>2</sub>とする。''p''<sub>1</sub>と''p''<sub>2</sub>は時刻''t''の関数であり、''p''<sub>1</sub>(''t'')と''p''2(''t'')と表わすとする。''p''<sub>1</sub>(''t'')と''p''<sub>2</sub>(''t'')は確率であるから、
 ''m''''h''''n''は[[#Two-state_model:_.E5.9F.BA.E7.A4.8E.E7.9A.84.E3.81.AA.E8.80.83.E3.81.88.E6.96.B9|Two-stateモデル]]モデルに従う値である。 開く速度定数αと閉じる速度定数βはいずれも電位に依存する。 HodgkinとHuxleyは''m''''h''''n''のそれぞれについていろいろな電位での αとβの値を実験的に測定し、一連のイカ巨大軸索を用いた実験値を便宜的に次の数式で表した。これらの式は何らかの理論から導きだされたものではない。なお注意すべき点として、以下の式で用いられている変数''v''は、膜電位''V''<sub>''m''</sub>そのものではなく静止電位を基準とした電位を意味している。またこれらの値は 6.3℃での速度定数であり、他の温度ではQ<sub>10</sub>=3で補正される。
 
::<math>\alpha_m = \frac{0.1(-v+25)}{\exp\left(\frac{-v+25}{10}\right)-1}</math>
 
::<math>\beta_m = 4\exp\left(\frac{-v}{18}\right)</math>
 
::<math>\alpha_h = 0.07\exp\left(\frac{-v}{20}\right)</math>
 
::<math>\beta_h = \frac{1}{\exp\left(\frac{-v+30}{10}\right)+1}</math>
 
::<math>\alpha_n = \frac{0.01(-v+10)}{\exp\left(\frac{-v+10}{10}\right)-1}</math>


::<math>p1(t) + p2(t) = 1\, </math>
::<math>\beta_n = 0.125\exp\left(\frac{-v}{80}\right)</math>


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


::<math> \frac{dp_1(t)}{dt} = -\alpha p_1(t) + \beta p_2(t)</math>  
 電位固定実験から得られた上記の式を用いると、電位の変化に応じて変化するNa<sup>+</sup>チャネルとK<sup>+</sup>チャネルのコンダクタンスを計算することができる。 一方、細胞の電位は、電流の入出により変化する。もし電位の初期条件と電流の時間経過がわかっていれば、電位を計算する事が出来る。これらを組み合わせることにより、細胞の電位の変化を計算できる。この関係は、現在広く行われている興奮性細胞の電位のシミュレーションの基本である。
::<math> \frac{dp_2(t)}{dt} = \alpha p_1(t) - \beta p_2(t)</math>


と表される。αとβが定数であるとして、定常状態になれば、
 以下は数式的な説明。


::<math> \frac{dp_1(\infty)}{dt} = -\alpha p_1(\infty) + \beta p_2(\infty) = 0\, </math>
----
::<math> \frac{dp_2(\infty)}{dt} = \alpha p_1(\infty) - \beta p_2(\infty) = 0\, </math>


ここで、
 細胞を[[wj:コンデンサ|キャパシタ]]([[wj:コンデンサ|コンデンサ]])と考える。細胞外液の電位を0とする。膜電位を''v''、蓄えられる電荷を''Q''、[[wj:静電容量|静電容量]](capacitance)をCとすると、


::<math>p_1(\infty) + p_2(\infty) = 1\, </math>
::<math>Q = Cv\, </math>


であるから、
の関係が成り立つ。両辺を時間''t''で微分すると、
 
::<math>\frac{dQ}{dt} = C\frac{dv}{dt}\, </math>


::<math>p_1(\infty) = \frac{\beta}{\alpha+\beta}\, </math>
となる。一方、流れる電流''I''はすなわち電荷量''Q''であるから、
::<math>p_2(\infty) = \frac{\alpha}{\alpha+\beta}\, </math>


となる。また微分方程式は解析的に解けて、
::<math>I = - \frac{dQ}{dt}\, </math>


::<math>p_1(t) = \left(p_1(0)-\frac{\beta}{\alpha+\beta}\right) e^{-(\alpha+\beta)t} + \frac{\beta}{\alpha+\beta}\, </math>
ここで右辺にマイナス符号が付くのは、細胞に流れ込む電流を負とする習わしによる。
::<math>p_2(t) = \left(p_2(0)-\frac{\alpha}{\alpha+\beta}\right) e^{-(\alpha+\beta)t} + \frac{\alpha}{\alpha+\beta}\, </math>


となる。 これらの式は次のことを示している。
電流''I''は、チャネルを流れる電流の和であるから、


*''p''<sub>1</sub>(''t'')と''p''<sub>2</sub>(''t'')はそれぞれ指数関数的に''p''<sub>1</sub>(∞)と''p''<sub>2</sub>(∞)に近づいていく
::<math>I = I_{Na} + I_{K} + I_{leak}\, </math>
*その時定数(time constant)τは1/(α+β)である
*これらの値''p''<sub>1</sub>(∞)、''p''<sub>2</sub>(∞)、τは、初期値''p''<sub>1</sub>(0)、''p''<sub>2</sub>(0)に依存しない。


さらに、
すなわち、


::<math>q_1(t) = p_1(t) - \frac{\beta}{\alpha+\beta}\, </math>
::<math>I = G^{max}_{Na} m^3 h (v-E_{Na}) + G^{max}_{K} n^4 (v-E_{K}) + G_{leak}(v-E_{leak})\, </math>
::<math>q_2(t) = p_2(t) - \frac{\alpha}{\alpha+\beta}\, </math>


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


::<math> q_1(t) = q_1(0)e^{-(\alpha + \beta)}\, </math>
::<math>\frac{dv}{dt} = -\frac{1}{C}\left(G^{max}_{Na} m^3 h (v-E_{Na}) + G^{max}_{K} n^4 (v-E_{K}) + G_{leak}(v-E_{leak})\right)\, </math>
::<math> q_2(t) = q_2(0)e^{-(\alpha + \beta)}\, </math>


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


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


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


以下は数式的な説明。<br>
 以下は数式的な説明。


----
----


外部より電流''I''<sub>ext</sub>を流した場合、電位の変化は、次の式で示される。  
 外部より電流''I''<sub>ext</sub>を流した場合、電位の変化は、次の式で示される。  


::<math>\frac{dv}{dt} = -\frac{1}{C}\left(\sum_X G_{X}(v-E_X) - I_{ext}\right)</math>
::<math>\frac{dv}{dt} = -\frac{1}{C}\left(\sum_X G_{X}(v-E_X) - I_{ext}\right)</math>


この式から、''I''<sub>ext</sub>と電位との関係を理解する事は難しい。しかし''v''が一定となるような外部電流''I''<sub>clamp</sub>を流すと、左辺は0となるため、  
''X''はいろいろな種類のチャネルを示している。
 
 この式から、''I''<sub>ext</sub>と電位との関係を理解する事は難しい。しかし''v''が一定となるような外部電流''I''<sub>clamp</sub>を流すと、左辺は0となるため、  


::<math> I_{clamp} = \sum_X G_X (v - E_X)\, </math>
::<math> I_{clamp} = \sum_X G_X (v - E_X)\, </math>


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


::<math>I_{clamp} = G_A (v - E_A)\, </math>
::<math>I_{clamp} = G_A (v - E_A)\, </math>


となる。これはOhmの法則である。ここで''I''<sub>clamp</sub>は実験の測定値、''v''は実験の設定値、''E''<sub>A</sub>は実験条件で定まる定数なので、  
となる。これは[[wj:オームの法則|Ohmの法則]]である。ここで''I''<sub>clamp</sub>は実験の測定値、''v''は実験の設定値、''E''<sub>A</sub>は実験条件で定まる定数なので、  


::<math>G_{A} = \frac{I_{clamp}}{v-E_A}\, </math>
::<math>G_{A} = \frac{I_{clamp}}{v-E_A}\, </math>


の関係式を用いて、実験データよりイオンチャネル''A''のコンダクタンス''G''<sub>A</sub>を算出できることになる。  
の関係式を用いて、実験データよりイオンチャネル''A''のコンダクタンス''G''<sub>A</sub>を算出できることになる。
 
== Two-state model: 基礎的な考え方  ==
 
 OpenとClosedの2つの状態がある系で、他の状態に移る率が一定の場合、次の性質がある。<br>
 
*[[wj:指数関数|指数関数]]的に変化し、一定の値に近づいていく。
*近づく値、変化の速さは、初期条件に依存しない。<br>
 
 以下は数式的な説明。
 
----
 
 OpenとClosedの2つの状態がある系を考え、Openの状態にある確率を''p''とする。Closedの状態にある確立は, 1-''p''となる。''p''は時刻''t''の関数であり、''p''(''t'')と表わすとする。
 
 いま状態Closedから状態Openへ移っていく単位時間での割合(速度定数、rate constant)をαとし、状態Openから状態Closedへの速度定数をβとする。 ''p''(''t'')の時間的経過を表わす微分方程式は、
 
::<math> \frac{dp(t)}{dt} = \alpha (1-p(t)) - \beta p(t)\, </math>
 
と表される。αとβが定数であるとして、''t'' =∞で定常状態になれば、
 
::<math> \frac{dp(\infty)}{dt} = \alpha (1-p(\infty)) - \beta p(\infty) = 0\, </math>
 
であるから、
 
::<math>p(\infty) = \frac{\alpha}{\alpha+\beta}\, </math>
 
となる。また微分方程式は解析的に解けて、
 
::<math>p(t) = \left(p(0)-\frac{\alpha}{\alpha+\beta}\right) e^{-(\alpha+\beta)t} + \frac{\alpha}{\alpha+\beta}\, </math>
 
となる。 これらの式は次のことを示している。
 
*''p''(''t'')は指数関数的に''p''(∞)に近づいていく。
*その時定数(time constant)τは1/(α+β)である。
*これらの値''p''(∞)、τは、初期値''p''(0)に依存しない。
 
 さらに、
 
::<math>q(t) = p(t) - \frac{\alpha}{\alpha+\beta}\,  </math>
 
と表すとすると、
 
::<math> q(t) = q(0)e^{-(\alpha + \beta)t}\, </math>
 
とより単純な形式となる。この関係は微分方程式の数値計算でよく用いられる。
 
 なお、このTwo-stateモデルは''n''-stateモデルに拡張することが可能であり、拡張された式はMarkov過程のシミュレーション等に用いられている。


== HHモデルに対する批判 ==
==批判 ==


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


== 現在におけるHHモデル  ==
#ゲート電流(gating current): チャネルタンパクの動きを反映すると考えられるゲート電流の電位依存性は、電流の電位依存性よりも過分極側にずれており、チャネルが最終的に開く過程は電位依存的ではないと推測された。<ref><pubmed>4700900</pubmed></ref> &nbsp;
#Single-channel recording: Na<sup>+</sup>チャネルの開口時間は、Hodgkin-Huxleyモデルで予測されるよりも短く、脱分極から遅れてチャネルが開く場合が観察された。 <ref><pubmed>6316158</pubmed></ref>
#Markov過程モデル:&nbsp;Hodgkin-Huxleyモデルよりもより複雑な過程を示すことの出来るMarkov過程モデルの方が、より詳細なチャネルの性質や、変異による性質の変化を表すことが出来た。 <ref><pubmed>6094703</pubmed></ref> <ref><pubmed>10448858</pubmed></ref>


(未完成)
== 現在におけるHodgkin-Huxleyモデル  ==


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


HH方程式の数値計算は、ノートパソコンでも十分行う事が出来る。
== Hodgkin-Huxley方程式を使ってみる  ==


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


== References  ==
==関連項目==
* [[積分発火モデル]]


未完成
== 参考文献  ==


<references />
<references />