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

編集の要約なし
編集の要約なし
編集の要約なし
1行目: 1行目:
Hodgkin-Huxley Equations  
英語名:Hodgkin-Huxley Equations  


== 概略  ==
== 概略  ==


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


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


#活動電位発生時に、ナトリウムイオン(Na<sup>+</sup>)とカリウムイオン(K<sup>+</sup>)が、脱分極により開く細胞膜の別々の通路を通ることを示した。この発見はイオンチャネルの存在を予測するものであり、その後のイオンチャネル研究の源となった。なお当時の論文では、イオンチャネル・チャネルという用語は用いられておらず、コンダクタンスという用語が使用されている。  
#活動電位発生時に、ナトリウムイオン(Na<sup>+</sup>)とカリウムイオン(K<sup>+</sup>)が、脱分極により開く細胞膜の別々の通路を通ることを示した。この発見はイオンチャネルの存在を予測するものであり、その後のイオンチャネル研究の源となった。なお当時の論文では、イオンチャネル・チャネルという用語は用いられておらず、コンダクタンスという用語が使用されている。  
13行目: 13行目:
== ''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>のコンダクタンスが電位に依存的なゲート(gate)により開閉されると考えた。  
 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>のコンダクタンスが電位に依存的なゲート(gate)により開閉されると考えた。  


*Na<sup>+</sup>チャネルは3つの活性化ゲート''m''と不活性化ゲート''h''により開閉される。  
*Na<sup>+</sup>チャネルは3つの活性化ゲート''m''と不活性化ゲート''h''により開閉される。  
20行目: 20行目:
''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モデル]]に従うと仮定されている。''m''と''n''は、静止時に閉じており(すなわち0の値をとる)、脱分極した時に開く(すべてが開けば1となる)。一方、''h''は静止時に開いており脱分極時に閉じる。''m''と''n''ではなく、''m''<sup>3</sup>および''n''<sup>4</sup>としたのは、主に電流の立ち上がりの形をよく再現するためである。  
''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モデル]]に従うと仮定されている。''m''と''n''は、静止時に閉じており(すなわち0の値をとる)、脱分極した時に開く(すべてが開けば1となる)。一方、''h''は静止時に開いており脱分極時に閉じる。''m''と''n''ではなく、''m''<sup>3</sup>および''n''<sup>4</sup>としたのは、主に電流の立ち上がりの形をよく再現するためである。  


電流はコンダクタンスと電圧に比例する(''I'' = ''GV''; Ohmの法則)。電圧の大きさは、細胞膜内外のイオン濃度差による電位(平衡電位)を補正しなくてはならない(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;  
 電流はコンダクタンスと電圧に比例する(''I'' = ''GV''; Ohmの法則)。電圧の大きさは、細胞膜内外のイオン濃度差による電位(平衡電位)を補正しなくてはならない(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>  
::<math>I_{Na} = G^{max}_{Na} m^3 h (v-E_{Na})\, </math>  
33行目: 33行目:
----
----


''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>そのものではなく静止電位を基準とした電位を意味している。  
 ''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>そのものではなく静止電位を基準とした電位を意味している。  
 
<br>


::<math>\alpha_m = \frac{0.1(-v+25)}{\exp\left(\frac{-v+25}{10}\right)-1}</math>  
::<math>\alpha_m = \frac{0.1(-v+25)}{\exp\left(\frac{-v+25}{10}\right)-1}</math>  
45行目: 43行目:
::<math>\alpha_n = \frac{0.01(-v+10)}{\exp\left(\frac{-v+10}{10}\right)-1}</math>  
::<math>\alpha_n = \frac{0.01(-v+10)}{\exp\left(\frac{-v+10}{10}\right)-1}</math>  
::<math>\beta_n = 0.125\exp\left(\frac{-v}{80}\right)</math>
::<math>\beta_n = 0.125\exp\left(\frac{-v}{80}\right)</math>
<br>


== 電位変化  ==
== 電位変化  ==


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


以下は数式的な説明。
 以下は数式的な説明。


----
----


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


::<math>Q = Cv\, </math>
::<math>Q = Cv\, </math>
82行目: 78行目:
::<math>\frac{dv}{dt} = -\frac{1}{C}\left(G^{max}_{Na} m^3 h (v-E_{Na}) + G^{max}_{K} n^3 (v-E_{K}) + G_{leak}(v-E_{leak})\right)\, </math>
::<math>\frac{dv}{dt} = -\frac{1}{C}\left(G^{max}_{Na} m^3 h (v-E_{Na}) + G^{max}_{K} n^3 (v-E_{K}) + G_{leak}(v-E_{leak})\right)\, </math>


となる。 <br>
となる。  
 
<br>
 
<br>


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


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


この式から、''I''<sub>ext</sub>と電位との関係を理解する事は難しい。しかし''v''が一定となるような外部電流''I''<sub>clamp</sub>を流すと、左辺は0となるため、  
 この式から、''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>
115行目: 107行目:


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


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


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


*指数関数的に変化し、一定の値に近づいていく。  
*指数関数的に変化し、一定の値に近づいていく。  
*近づく値、変化の速さは、初期条件に依存しない。<br>
*近づく値、変化の速さは、初期条件に依存しない。<br>


以下は数式的な説明。
 以下は数式的な説明。


----
----


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


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


::<math> \frac{dp(t)}{dt} = \alpha (1-p(t)) - \beta p(t)\, </math>
::<math> \frac{dp(t)}{dt} = \alpha (1-p(t)) - \beta p(t)\, </math>
153行目: 143行目:
*これらの値''p''(∞)、τは、初期値''p''(0)に依存しない。
*これらの値''p''(∞)、τは、初期値''p''(0)に依存しない。


さらに、
 さらに、


::<math>q(t) = p(t) - \frac{\alpha}{\alpha+\beta}\,  </math>
::<math>q(t) = p(t) - \frac{\alpha}{\alpha+\beta}\,  </math>
163行目: 153行目:
とより単純な形式となる。この関係は微分方程式の数値計算でよく用いられる。  
とより単純な形式となる。この関係は微分方程式の数値計算でよく用いられる。  


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


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


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


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


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


<br>
== HH方程式を使ってみる  ==


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


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


<references /> <br>
(執筆者:井本敬二 担当編集委員:柚崎通介)