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

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


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


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


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


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


== <math>\textstyle m^3 h</math>と<math>\textstyle n^4</math> ==
==<math>\textstyle m^3 h</math>と<math>\textstyle n^4</math>==


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


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


2つの状態1と2をとる事の出来る系を考え、それぞれの状態にある確率を<math>\textstyle p1</math>と<math>\textstyle p2</math> とする。<math>\textstyle p1</math>と<math>\textstyle p2</math>は時刻<math>\textstyle t</math>の関数であり、<math>\textstyle p1(t)</math>と<math>\textstyle p2(t)</math>と表わされる。<math>\textstyle p1(t)</math>と<math>\textstyle p2(t)</math>は確率であるから、  
\par 2つの状態1と2をとる事の出来る系を考え、それぞれの状態にある確率を<math>\textstyle p1</math>と<math>\textstyle p2</math> とする。<math>\textstyle p1</math>と<math>\textstyle p2</math>は時刻<math>\textstyle t</math>の関数であり、<math>\textstyle p1(t)</math>と<math>\textstyle p2(t)</math>と表わされる。<math>\textstyle p1(t)</math>と<math>\textstyle p2(t)</math>は確率であるから、


:<span class="texhtml">''p''1(''t'') + ''p''2(''t'') = 1</span>
:<math> p1(t) + p2(t) = 1 </math>


の関係にある。いま状態1から状態2へ移っていく単位時間での割合(遷移率)を<math>\textstyle \alpha</math>とし、状態2から状態1への遷移率を<math>\textstyle \beta</math>とする。 <math>\textstyle p1(t)</math>と<math>\textstyle p2(t)</math>の時間的経過を表わす微分方程式は、  
の関係にある。いま状態1から状態2へ移っていく単位時間での割合(遷移率)を<math>\textstyle \alpha</math>とし、状態2から状態1への遷移率を<math>\textstyle \beta</math>とする。 <math>\textstyle p1(t)</math>と<math>\textstyle p2(t)</math>の時間的経過を表わす微分方程式は、


:<math> \frac{dp1(t)}{dt} = -\alpha p1(t) + \beta p2(t)</math>
:<math> \frac{dp1(t)}{dt} = -\alpha p1(t) + \beta p2(t)</math>
27行目: 25行目:
:<math> \frac{dp2(t)}{dt} = \alpha p1(t) - \beta p2(t)</math>
:<math> \frac{dp2(t)}{dt} = \alpha p1(t) - \beta p2(t)</math>


\par と表される。<math>\textstyle \alpha</math>と<math>\textstyle \beta</math>が定数であるとして、定常状態になれば、  
\par と表される。<math>\textstyle \alpha</math>と<math>\textstyle \beta</math>が定数であるとして、定常状態になれば、


:<math> \frac{dp1(\infty)}{dt} = -\alpha p1(\infty) + \beta p2(\infty) = 0</math>
:<math> \frac{dp1(\infty)}{dt} = -\alpha p1(\infty) + \beta p2(\infty) = 0</math>
33行目: 31行目:
:<math> \frac{dp2(t)}{dt} = \alpha p1(\infty) - \beta p2(\infty) = 0</math>
:<math> \frac{dp2(t)}{dt} = \alpha p1(\infty) - \beta p2(\infty) = 0</math>


<math>\textstyle p1(\infty) + p2(\infty) = 1</math> であるから、  
<math>\textstyle p1(\infty) + p2(\infty) = 1</math> であるから、


:<math>p1(\infty) = \frac{\beta}{\alpha+\beta}</math>
:<math>p1(\infty) = \frac{\beta}{\alpha+\beta}</math>
39行目: 37行目:
:<math>p2(\infty) = \frac{\alpha}{\alpha+\beta}</math>
:<math>p2(\infty) = \frac{\alpha}{\alpha+\beta}</math>


となる。また微分方程式を解析的に解くと、  
となる。また微分方程式を解析的に解くと、


:<math>p1(t) = \left(p1(0)-\frac{\beta}{\alpha+\beta}\right) e^{-(\alpha+\beta)t} + \frac{\beta}{\alpha+\beta} </math>
:<math>p1(t) = \left(p1(0)-\frac{\beta}{\alpha+\beta}\right) e^{-(\alpha+\beta)t} + \frac{\beta}{\alpha+\beta} </math>
45行目: 43行目:
:<math>p2(t) = \left(p2(0)-\frac{\alpha}{\alpha+\beta}\right) e^{-(\alpha+\beta)t} + \frac{\alpha}{\alpha+\beta} </math>
:<math>p2(t) = \left(p2(0)-\frac{\alpha}{\alpha+\beta}\right) e^{-(\alpha+\beta)t} + \frac{\alpha}{\alpha+\beta} </math>


となる。これらの式は、<math>\textstyle p1(t)</math>と<math>\textstyle p2(t)</math>はそれぞれ指数関数的に<math>\textstyle p1(\infty)</math>と<math>\textstyle p2(\infty)</math>に近づいていき、その時定数<math>\textstyle \tau</math>は<math>\textstyle 1/(\alpha+\beta)</math>であること、およびこれらの値<math>\textstyle p1(\infty)</math>、<math>\textstyle p2(\infty)</math>、<math>\textstyle \tau</math>は、初期値<math>\textstyle p1(0)</math>、<math>\textstyle p2(0)</math>には依存しないことを示している。さらに、  
となる。これらの式は、<math>\textstyle p1(t)</math>と<math>\textstyle p2(t)</math>はそれぞれ指数関数的に<math>\textstyle p1(\infty)</math>と<math>\textstyle p2(\infty)</math>に近づいていき、その時定数<math>\textstyle \tau</math>は<math>\textstyle 1/(\alpha+\beta)</math>であること、およびこれらの値<math>\textstyle p1(\infty)</math>、<math>\textstyle p2(\infty)</math>、<math>\textstyle \tau</math>は、初期値<math>\textstyle p1(0)</math>、<math>\textstyle p2(0)</math>には依存しないことを示している。さらに、


\par  
\par


:<math>q1(t) = p1(t) - \frac{\beta}{\alpha+\beta} </math>
:<math>q1(t) = p1(t) - \frac{\beta}{\alpha+\beta} </math>
53行目: 51行目:
:<math>q2(t) = p2(t) - \frac{\alpha}{\alpha+\beta} </math>
:<math>q2(t) = p2(t) - \frac{\alpha}{\alpha+\beta} </math>


とすると、  
とすると、


:<span class="texhtml">''q''1(''t'') = ''q''1(0)''e''<sup> − (α + β)''t''</sup></span>
:<math>q1(t)= q1(0)e^{-(\alpha+\beta)t}</math>


:<span class="texhtml">''q''2(''t'') = ''q''2(0)''e''<sup> − (α + β)''t''</sup></span>
:<math>q2(t)= q2(0)e^{-(\alpha+\beta)t}</math>


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


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


Hodgkin-Huxley以前に、電気生理学の実験が行われていなかったわけではない。電流と電位変化に関する研究は、かなり多く行われていた。しかしながら、細胞にはいろいろなイオンチャネルを通して電流が流れるため、細胞の電位<math>\textstyle v</math>と外部から流す電流<math>\textstyle I_{ext}</math>の間の関係は、  
Hodgkin-Huxley以前に、電気生理学の実験が行われていなかったわけではない。電流と電位変化に関する研究は、かなり多く行われていた。しかしながら、細胞にはいろいろなイオンチャネルを通して電流が流れるため、細胞の電位<math>\textstyle v</math>と外部から流す電流<math>\textstyle I_{ext}</math>の間の関係は、


:<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>


となり、実験データの解釈は単純ではない。電位をコントロールして行う実験方法であるvoltage clamp(電位固定法)は、1940年代にアメリカの生物物理学者Kenneth Cole (1900 - 1984)らにより開発された。HodgkinとHuxleyはこのvoltage-clampを巧みに利用して大きな成果を得る事が出来たと言える。上記の式で<math>\textstyle v</math>が一定となるように外部電流を<math>\textstyle I_{clamp}</math>を流すと、左辺は<math>\textstyle 0</math>となるため、  
となり、実験データの解釈は単純ではない。電位をコントロールして行う実験方法であるvoltage clamp(電位固定法)は、1940年代にアメリカの生物物理学者Kenneth Cole (1900 - 1984)らにより開発された。HodgkinとHuxleyはこのvoltage-clampを巧みに利用して大きな成果を得る事が出来たと言える。上記の式で<math>\textstyle v</math>が一定となるように外部電流を<math>\textstyle I_{clamp}</math>を流すと、左辺は<math>\textstyle 0</math>となるため、


:<span class="texhtml">
:<math>I_{clamp} = \sum_X G_{X}(v-E_X)</math>
</span>


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


:<math>I_{clamp} = G_{A}(v-E_A)</math>


<br>
となる。ここで<math>\textstyle I_{clamp}</math>は実験の測定値、<math>\textstyle v</math>は実験の設定値、<math>\textstyle E_A</math>は実験条件で定まる定数なので、イオンチャネル<math>\textstyle A</math>のコンダクタンス<math>\textstyle G_A</math>を、
 
{|
|- style="text-align: center;"
| ''I''<sub>''c''''l''''a''''m''''p''</sub> =
| <span style="font-size: x-large; font-family: serif;">∑</span>
| ''G''<sub>''X''</sub>(''v'' − ''E''<sub>''X''</sub>)
|- style="text-align: center; vertical-align: top;"
|
| ''X''
|
|}
 
<br>
 
という関係が得られる。もし溶液の組成を工夫しチャネルのブロッカーなどを用いて、イオンチャネル<math>\textstyle A</math>を流れる電流が測れたとすると、
 
:<span class="texhtml">''I''<sub>''c''''l''''a''''m''''p''</sub> = ''G''<sub>''A''</sub>(''v'' − ''E''<sub>''A''</sub>)</span>
 
となる。ここで<math>\textstyle I_{clamp}</math>は実験の測定値、<math>\textstyle v</math>は実験の設定値、<math>\textstyle E_A</math>は実験条件で定まる定数なので、イオンチャネル<math>\textstyle A</math>のコンダクタンス<math>\textstyle G_A</math>を、  


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


と算出できることになる。  
と算出できることになる。


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


Single-channel recording  
Single-channel recording


Markovモデル  
Markovモデル


Fractalモデルとの論争  
Fractalモデルとの論争


== 現在におけるHHモデル ==
==現在におけるHHモデル==


== References ==
==References==


<references/>
<references/>
66

回編集