「ミカエリス・メンテン式」の版間の差分

出典: フリー百科事典『ウィキペディア(Wikipedia)』
削除された内容 追加された内容
en:Michaelis–Menten kinetics から一部和訳
タグ: サイズの大幅な増減 曖昧さ回避ページへのリンク
2行目: 2行目:
'''ミカエリス・メンテン式'''(ミカエリス・メンテンしき、{{lang-en-short|Michaelis–Menten equation}})とは、[[酵素反応速度論|酵素の反応速度論]]に大きな業績を残した[[レオノール・ミカエリス]]と[[モード・メンテン|モード・レオノーラ・メンテン]]にちなんだ、[[酵素]]の反応速度''v'' に関する式で、
'''ミカエリス・メンテン式'''(ミカエリス・メンテンしき、{{lang-en-short|Michaelis–Menten equation}})とは、[[酵素反応速度論|酵素の反応速度論]]に大きな業績を残した[[レオノール・ミカエリス]]と[[モード・メンテン|モード・レオノーラ・メンテン]]にちなんだ、[[酵素]]の反応速度''v'' に関する式で、
:<math>v \equiv \frac{d[\mathrm{P}]}{dt} = \frac{V_{\mathrm{max}}[\mathrm{S}]}{K_{\mathrm{m}} + [\mathrm{S}]}</math>
:<math>v \equiv \frac{d[\mathrm{P}]}{dt} = \frac{V_{\mathrm{max}}[\mathrm{S}]}{K_{\mathrm{m}} + [\mathrm{S}]}</math>
で表される。ここで、[P]は反応産物の濃度、[S]は[[基質]]濃度、''V''<sub>max</sub> は基質濃度が無限大のときの反応速度である。また、''K<sub>m</sub>'' は'''ミカエリス・メンテン定数'''と言い、''v'' = ''V''<sub>max</sub> /2(最大速度の半分の速度)を与える基質濃度を表す。この式をもとにしたモデルを'''ミカエリス・メンテン動力学''' ({{en|Michaelis–Menten kinetics}})という。
で表される。ここで、[P]は反応産物の濃度、[S]は[[基質]]濃度、''V''<sub>max</sub> は基質濃度が無限大のときの反応速度である。また、''K<sub>m</sub>'' は'''ミカエリス・メンテン定数'''と言い、''v'' = ''V''<sub>max</sub> /2(最大速度の半分の速度)を与える基質濃度を表す。この式をもとにしたモデルを'''ミカエリス・メンテン動力学'''という。


この式により、反応速度''v'' は
この式により、反応速度''v'' は
8行目: 8行目:
* 基質濃度が高い([S] ≫ ''K<sub>m</sub>'' )ときはその濃度に無関係に最大速度 ''V''<sub>max</sub> に収束
* 基質濃度が高い([S] ≫ ''K<sub>m</sub>'' )ときはその濃度に無関係に最大速度 ''V''<sub>max</sub> に収束
となることが分かる<ref>{{cite|和書 |author=大竹伝雄 |title=化学工学Ⅲ |publisher=岩波書店 |year=1978 |isbn=4-00-021103-X |page=16}}</ref>。
となることが分かる<ref>{{cite|和書 |author=大竹伝雄 |title=化学工学Ⅲ |publisher=岩波書店 |year=1978 |isbn=4-00-021103-X |page=16}}</ref>。

== 解析解 ==
ミカエリス・メンテン式は左辺に基質 S の濃度の[[時間微分]]を持つ[[微分方程式]]である。この方程式は[[変数分離]]形にして積分することで解くことができるが、解は[[初等関数]]で表現できる形にはならず、[[ランベルトのW関数]] ''W'' を用いて以下のように表される。

:<math>[\mathrm{S}](t) = K_{\mathrm{m}} W \left( \frac{1}{K_{\mathrm{m}}} \exp \left( -\frac{V_{\mathrm{max}} t + C}{K_{\mathrm{m}}} \right) \right)</math>

ここで ''C'' は[[積分定数]]である。

== 実験によるパラメータの決定 ==
2つのパラメータ、''V''<sub>max</sub> とミカエリス・メンテン定数''K<sub>m</sub>'' の値は実験から、[[非線形最小二乗法]]によって求めることができる。

歴史的には、[[ラインウィーバー=バークプロット|ラインウィーバー・バークのプロット]]<ref>{{cite|和書|author=A. G. Whittaker|author2=A. R. Mount|author3=M. R. Heal|coauthor=中村亘男訳|title=キーノート化学シリーズ 物理化学キーノート|publisher=シュプリンガー・フェアラーク東京|year=2002|isbn=4-431-70956-8}}</ref>、[[ヘインズ=ウルフプロット]]、[[イーディー=ホフステー図]]といった普通の[[最小二乗法]]を適用できる方法や[[コーニッシュボーデンの直接的直線プロット]]<ref>{{cite journal|和書|title=原典からの酵素反応速度論|author=知名秀泰、岡田豊|journal=生物工学|volume=82|issue=1|year=2014|pages=20–25|url=https://www.sbj.or.jp/wp-content/uploads/file/sbj/9201/9201_yomoyama.pdf}}</ref>などの方法が使われたが、今日ではこれらの方法を使うべきではないとされる。


== 導出 ==
== 導出 ==
=== 迅速平衡法による導出 ===
=== 迅速平衡法による導出 ===
酵素(以下 E)が基質(以下 S)と結合して酵素基質複合体(以下 ES)を形成、ES E S に戻るか反応生成物(以下 P)を生成する一連の反応機構を以下のように仮定する。
酵素(以下E)が基質(以下S)と結合して酵素基質複合体(以下ES)を形成、ESがEとSに戻るか反応生成物(以下P)を生成する一連の反応機構を以下のように仮定する。


:<math>\mathrm{E} + \mathrm{S} \, \overset{k_{+1}}{\underset{k_{-1}} \rightleftharpoons} \, \mathrm{ES} \, \overset{k_{+2}} {\longrightarrow} \, \mathrm{E} + \mathrm{P}</math>
:<math>\mathrm{E} + \mathrm{S} \, \overset{k_{+1}}{\underset{k_{-1}} \rightleftharpoons} \, \mathrm{ES} \, \overset{k_{+2}} {\longrightarrow} \, \mathrm{E} + \mathrm{P}</math>


この反応は <math>\mathrm{E}+\mathrm{S} \Leftrightarrow \mathrm{ES}</math> と <math>\mathrm{ES} \rarr \mathrm{E}+\mathrm{P}</math> の 2 つの反応過程からできている。後者の反応を律速段階と仮定し[[反応速度定数]]を ''k''<sub>+2</sub> と設定する。
この反応は <math>\mathrm{E}+\mathrm{S} \Leftrightarrow \mathrm{ES}</math> と <math>\mathrm{ES} \rarr \mathrm{E}+\mathrm{P}</math> の2つの反応過程からできている。後者の反応を律速段階と仮定し[[反応速度定数]]を ''k''<sub>+2</sub>と設定する。
<math>\mathrm{E}+\mathrm{S} \Leftrightarrow \mathrm{ES}</math> の反応は迅速に[[化学平衡]]に達していると仮定し、[[解離定数]]を''K''<sub>s</sub> と設定する。
<math>\mathrm{E}+\mathrm{S} \Leftrightarrow \mathrm{ES}</math> の反応は迅速に[[化学平衡]]に達していると仮定し、[[解離定数]]を''K''<sub>s</sub> と設定する。
:<math>K_\mathrm{s} = \frac{[\mathrm{E}][\mathrm{S}]}{[\mathrm{ES}]} \cdots (1)</math>
:<math>K_\mathrm{s} = \frac{[\mathrm{E}][\mathrm{S}]}{[\mathrm{ES}]} \cdots (1)</math>


仮定されている反応系に存在する酵素種は、基質と結合していない酵素 E と、基質 S と結合した酵素 ES 2 種類のみなので、全酵素濃度 &#x5B;E&#x5D;<sub>0</sub> は両者の濃度の和に等しい。
仮定されている反応系に存在する酵素種は、基質と結合していない酵素Eと、基質Sと結合した酵素ESの2種類のみなので、全酵素濃度 &#x5B;E&#x5D;<sub>0</sub> は両者の濃度の和に等しい。
:<math>[\mathrm{E}] + [\mathrm{ES}] = [\mathrm{E}]_0 \cdots (2)</math>
:<math>[\mathrm{E}] + [\mathrm{ES}] = [\mathrm{E}]_0 \cdots (2)</math>


40行目: 28行目:
:<math>v = \frac{d[\mathrm{P}]}{dt} = k_{+2} [\mathrm{ES}] \cdots (4)</math>
:<math>v = \frac{d[\mathrm{P}]}{dt} = k_{+2} [\mathrm{ES}] \cdots (4)</math>


(3)を(4)に代入して、
(3) (4) に代入して、
:<math>v = \frac{ k_{+2} [\mathrm{E}]_0 [\mathrm{S}]}{K_\mathrm{s} + [\mathrm{S}]} \cdots (5)</math>
:<math>v = \frac{ k_{+2} [\mathrm{E}]_0 [\mathrm{S}]}{K_\mathrm{s} + [\mathrm{S}]} \cdots (5)</math>


(4)式から反応速度 ''v'' は [ES] に比例することがわかるが、[ES] の最大値は(2)式より [E]<sub>0</sub> である。したがって反応速度 ''v'' の最大値 ''V''<sub>max</sub> は次式となる。
(4) 式から反応速度''v''は [ES] に比例することがわかるが、[ES] の最大値は (2) 式より [E]<sub>0</sub>である。したがって反応速度''v''の最大値''V''<sub>max</sub>は次式となる。
:<math>V_\mathrm{max} = k_{+2} [\mathrm{E}]_0 \cdots (6)</math>
:<math>V_\mathrm{max} = k_{+2} [\mathrm{E}]_0 \cdots (6)</math>


(5), (6)より、
(5), (6) より、
:<math>v = \frac{V_\mathrm{max}[\mathrm{S}]}{K_\mathrm{s} + [\mathrm{S}]}</math>
:<math>v = \frac{V_\mathrm{max}[\mathrm{S}]}{K_\mathrm{s} + [\mathrm{S}]}</math>


=== 定常状態法による導出 ===
=== 定常状態法による導出 ===
系の別の解析がイギリスの植物学者{{仮リンク|ジョージ・エドワード・ブリッグズ|en|George Edward Briggs}}とイギリスの遺伝学者[[J・B・S・ホールデン]]によって1925年に行われた<ref name="briggs25">
迅速平衡法では <math>\mathrm{E}+\mathrm{S} \Leftrightarrow \mathrm{ES}</math> が迅速に平衡に達すると仮定されているため、<math>\mathrm{ES} \rarr \mathrm{E}+\mathrm{P}</math> の速度定数が <math>\mathrm{E}+\mathrm{S} \Leftrightarrow \mathrm{ES}</math> の速度定数よりもはるかに小さい反応にしか成り立たない。定常状態法によって求めることで一般の反応でも同様の式が成り立つことが証明される。
{{cite journal
| last1 = Briggs
| first1 = G.E.
| last2 = Haldane
| first2 = J.B.S.
| year = 1925
| title = A note on the kinematics of enzyme action
| journal = Biochem J
| volume = 19
| pages = 338–339
| doi =10.1042/bj0190338
| pmid = 16743508
| issue = 2
| pmc = 1259181
}}
</ref><ref name=LaidlerPCBio>{{cite book |last1=Laidler |first1=Keith J. |title=Physical Chemistry with Biological Applications |date=1978 |publisher=Benjamin/Cummings |isbn=0-8053-5680-0 |pages=428–430}}</ref>。この迅速平衡法では <math>\mathrm{E}+\mathrm{S} \Leftrightarrow \mathrm{ES}</math> が迅速に平衡に達すると仮定されているため、<math>\mathrm{ES} \rarr \mathrm{E}+\mathrm{P}</math> の速度定数が <math>\mathrm{E}+\mathrm{S} \Leftrightarrow \mathrm{ES}</math> の速度定数よりもはるかに小さい反応にしか成り立たない。定常状態法によって求めることで一般の反応でも同様の式が成り立つことが証明される。


反応機構は同様で、<math>\mathrm{E}+\mathrm{S} \Leftrightarrow \mathrm{ES}</math> について右向きの速度定数を''k''<sub>+1</sub> 、左向きの速度定数を''k''<sub>-1</sub> とする。
反応機構は同様で、<math>\mathrm{E}+\mathrm{S} \Leftrightarrow \mathrm{ES}</math> について右向きの速度定数を''k''<sub>+1</sub> 、左向きの速度定数を''k''<sub>-1</sub> とする。
59行目: 63行目:
\frac{d[\mathrm{ES}]}{dt} &= k_{+1}[\mathrm{E}][\mathrm{S}] - (k_{-1} + k_{+2})[\mathrm{ES}] = 0 \cdots (2)\end{align}</math>
\frac{d[\mathrm{ES}]}{dt} &= k_{+1}[\mathrm{E}][\mathrm{S}] - (k_{-1} + k_{+2})[\mathrm{ES}] = 0 \cdots (2)\end{align}</math>


この反応機構では E ES しか酵素種が存在しないので
この反応機構ではEとESしか酵素種が存在しないので
:<math>[\mathrm{E}] + [\mathrm{ES}] = [\mathrm{E}]_0 \cdots (3)</math>
:<math>[\mathrm{E}] + [\mathrm{ES}] = [\mathrm{E}]_0 \cdots (3)</math>


反応産物は ES より''k''<sub>+2</sub> の速度で生成されるので
反応産物はESより''k''<sub>+2</sub>の速度で生成されるので
:<math>v = \frac{d[\mathrm{P}]}{dt} = k_{+2}[\mathrm{ES}] \cdots (4)</math>
:<math>v = \frac{d[\mathrm{P}]}{dt} = k_{+2}[\mathrm{ES}] \cdots (4)</math>


(1)または(2)式と(3)式を連立方程式とみなして [ES] を求めると
(1) または (2) 式と (3) 式を連立方程式とみなして [ES] を求めると
:<math>[\mathrm{ES}] = \frac{k_{+1}[\mathrm{E}]_0[\mathrm{S}]}{k_{-1} + k_{+2} + k_{+1}[\mathrm{S}]} \cdots (5)</math>
:<math>[\mathrm{ES}] = \frac{k_{+1}[\mathrm{E}]_0[\mathrm{S}]}{k_{-1} + k_{+2} + k_{+1}[\mathrm{S}]} \cdots (5)</math>


(5)式を(4)式に代入して速度''v'' を得た後、分子分母を''k''<sub>+1</sub> で割る。
(5) 式を (4) 式に代入して速度''v'' を得た後、分子分母を''k''<sub>+1</sub> で割る。
:<math>v = \frac{k_{+2} [\mathrm{E}]_0 [\mathrm{S}] }{ (k_{-1} + k_{+2})/k_{+1} + [\mathrm{S}]} \cdots (6)</math>
:<math>v = \frac{k_{+2} [\mathrm{E}]_0 [\mathrm{S}] }{ (k_{-1} + k_{+2})/k_{+1} + [\mathrm{S}]} \cdots (6)</math>


速度パラメーターとして
速度パラメーターとして
:<math>V_{\mathrm{max}} = k_{+2} [\mathrm{E}]_0, \quad K_{\mathrm{m}} = \frac{k_{-1} + k_{+2}}{k_{+1}}</math>
:<math>V_{\mathrm{max}} = k_{+2} [\mathrm{E}]_0, \quad K_{\mathrm{m}} = \frac{k_{-1} + k_{+2}}{k_{+1}}</math>
と定義すれば、(6)式は
と定義すれば、(6) 式は
:<math>v = \frac{V_{\mathrm{max}}[\mathrm{S}]}{K_{\mathrm{m}} + [\mathrm{S}]}</math>
:<math>v = \frac{V_{\mathrm{max}}[\mathrm{S}]}{K_{\mathrm{m}} + [\mathrm{S}]}</math>
となる。
となる。

===仮定と限界===
導出における第一段階は、自由[[拡散]]を頼りにしている{{仮リンク|質量作用の法則|en|Law of mass action}}を適用することである。しかし、高濃度の[[タンパク質]]が存在する生細胞の環境では、[[細胞質]]は自由に流れる液体というよりも粘性のあるゲルのように振る舞うことが多く、拡散による分子の動きが制限され、反応速度が変化する<ref name="zhou08">
{{Cite journal
| last1 = Zhou
| first1 = H.X.
| last2 = Rivas
| first2 = G.
| last3 = Minton
| first3 = A.P.
| title = Macromolecular crowding and confinement: biochemical, biophysical, and potential physiological consequences
| journal = Annu Rev Biophys
| volume = 37
| issue =1
| pages = 375–97
| year = 2008
| doi = 10.1146/annurev.biophys.37.032807.125817
| pmid = 18573087
| pmc = 2826134
}}
</ref>。質量作用の法則は不均質な環境では有効であるが<ref name="grima07">
{{Cite journal
| last1 = Grima
| first1 = R.
| last2 = Schnell
| first2 = S.
| title = A systematic investigation of the rate laws valid in intracellular environments
| journal = Biophys Chem
| volume = 124
| issue = 1
| pages = 1–10
|date=Oct 2006
| doi = 10.1016/j.bpc.2006.04.019
| pmid = 16781049
}}
</ref>、細胞質を[[フラクタル]]としてモデル化する方が、その限定された移動性の動力学を表現するのに適している<ref name="schnell04">
{{Cite journal
| last1 = Schnell
| first1 = S.
| last2 = Turner
| first2 = T.E.
| title = Reaction kinetics in intracellular environments with macromolecular crowding: simulations and rate laws
| journal = Prog Biophys Mol Biol
| volume = 85
| issue = 2–3
| pages = 235–60
| year = 2004
| doi = 10.1016/j.pbiomolbio.2004.01.012
| pmid = 15142746
| citeseerx = 10.1.1.117.1997
}}
</ref>。

上記の2つのアプローチによって予測される結果の反応速度は類似しているが、唯一の違いは、迅速平衡近似では定数を<math>K_d</math>と定義するのに対し、準定常近似では<math>K_\mathrm{M}</math>を使用することである。しかしばがら、それぞれのアプローチは異なる仮定に基づいている。ミカエリス・メンテンの平衡解析は、基質が生成物の形成よりもはるかに速い時間スケールで平衡に達する場合、より正確には
:<math>\varepsilon_d = \frac{k_\mathrm{cat}}{k_r} \ll 1</math>
の時に妥当である<ref name="keener08">
{{cite book
| last1 = Keener
| first1 = J.
| last2 = Sneyd
| first2 = J.
| title = Mathematical Physiology: I: Cellular Physiology
| edition = 2
| year = 2008
| publisher = Springer
| isbn = 978-0-387-75846-6
}}
</ref>。

対照的にブリッグズ・ホールデン準定常状態解析はもし
:<math chem>\varepsilon_m = \frac\ce{[E]_0}{[\ce S]_0 + K_\ce{M}} \ll 1</math>
ならば妥当である<ref name="segel89"/><ref name="murray02">
{{cite book
| last1 = Murray
| first1 = J.D.
| title = Mathematical Biology: I. An Introduction
| edition = 3
| year = 2002
| publisher = Springer
| isbn = 978-0-387-95223-9
}}
</ref>。

したがって、酵素濃度が基質濃度または<math>K_\mathrm{M}</math>、あるいはその両方よりもはるかに低い場合に成立する。

ミカエリス・メンテン解析、ブリッグズ・ホールデン解析のいずれにおいても、近似の質は<math>\varepsilon\,\!</math>が小さくなるにつれて向上する。しかし、モデル構築の際には、その前提条件を無視してミカエリス・メンテン速度論が用いられることが多い<ref name="keener08" />。

重要なことは、不可逆性は扱いやすい解析解を得るために必要な単純化であるが、一般的な場合では生成物の形成は実際には不可逆的ではないということである。酵素反応はより正確には次のように記述される。

:<chem> E{} + S <=>[\mathit{k_{f_1}}][\mathit{k_{r_1}}] ES <=>[\mathit{k_{f_2}}][\mathit{k_{r_2}}] E{} + P</chem>

一般的に、不可逆性の仮定は、以下のいずれかが真である状況で良い仮定である。

: 1. 基質の濃度が生成物の濃度よりも非常に大きい。
::<chem> [S] \gg [P].</chem>
これは、標準的な[[in vitro]]のアッセイ条件でも真であるし、[[in vivo]]の多くの生物学的反応、特に生成物が後続の反応によって継続的に除去される場合にも真である。

: 2. 反応で放出されるエネルギーが非常に大きい、つまり
::<math>\Delta{G} \ll 0</math>

この2つの条件が成立しない場合(すなわち、反応が低エネルギーであり、生成物のプールがかなり存在する場合)、ミカエリス・メンテン方程式は破綻し、酵素の生物学を理解するためには、順反応と逆反応をあらわに考慮した、より複雑なモデリングアプローチが必要となる。


== 阻害がある場合のミカエリス・メンテン式 ==
== 阻害がある場合のミカエリス・メンテン式 ==
84行目: 189行目:
* [[不競合阻害]](不競争阻害)
* [[不競合阻害]](不競争阻害)
* [[混合阻害]]
* [[混合阻害]]
などの種類がある。前述の両逆数プロットを使うとこれらを見分けることができる。
などの種類がある。[[ラインウィーバー=バークプロット|両逆数プロット]]を使うとこれらを見分けることができる。


=== 基質阻害 ===
=== 基質阻害 ===
93行目: 198行目:
:<math>K_\mathrm{ss} = \frac{[\mathrm{ES}][\mathrm{S}]}{[\mathrm{ESS}]}</math>
:<math>K_\mathrm{ss} = \frac{[\mathrm{ES}][\mathrm{S}]}{[\mathrm{ESS}]}</math>


含まれる酵素種は E, ES ESS の3種類なので、全酵素濃度は
含まれる酵素種はEESとESSの3種類なので、全酵素濃度は
:<math>[\mathrm{E}]_0 = [\mathrm{E}] + [\mathrm{ES}] + [\mathrm{ESS}]</math>
:<math>[\mathrm{E}]_0 = [\mathrm{E}] + [\mathrm{ES}] + [\mathrm{ESS}]</math>


103行目: 208行目:


=== 競合阻害 ===
=== 競合阻害 ===
競合阻害(競争阻害、拮抗阻害とも)とは、基質と[[阻害剤]](以下、I)が酵素の同じ活性中心に結合する場合に起こる阻害のこと。
競合阻害(競争阻害、拮抗阻害とも)とは、基質と[[阻害剤]]以下、I)が酵素の同じ活性中心に結合する場合に起こる阻害のこと。


反応機構は次式の通り。
反応機構は次式の通り。
112行目: 217行目:
:<math>K_i = \frac{[\mathrm{E}][\mathrm{I}]}{[\mathrm{EI}]}</math>
:<math>K_i = \frac{[\mathrm{E}][\mathrm{I}]}{[\mathrm{EI}]}</math>


含まれる酵素種は E, ES EI の3種類なので、全酵素濃度は
含まれる酵素種はEESとEIの3種類なので、全酵素濃度は
:<math>[\mathrm{E}]_0 = [\mathrm{E}] + [\mathrm{ES}] + [\mathrm{EI}]</math>
:<math>[\mathrm{E}]_0 = [\mathrm{E}] + [\mathrm{ES}] + [\mathrm{EI}]</math>


132行目: 237行目:
K_i &= \frac{[\mathrm{E}][\mathrm{I}]}{[\mathrm{EI}]} = \frac{[\mathrm{ES}][\mathrm{I}]}{[\mathrm{ESI}]}
K_i &= \frac{[\mathrm{E}][\mathrm{I}]}{[\mathrm{EI}]} = \frac{[\mathrm{ES}][\mathrm{I}]}{[\mathrm{ESI}]}
\end{align}</math>
\end{align}</math>
反応機構の中にある酵素種は E, ES, EI ESI の4種類なので、全酵素濃度は
反応機構の中にある酵素種はEESEIとESIの4種類なので、全酵素濃度は
:<math>[\mathrm{E}]_0 = [\mathrm{E}] + [\mathrm{ES}] + [\mathrm{EI}] + [\mathrm{ESI}]</math>
:<math>[\mathrm{E}]_0 = [\mathrm{E}] + [\mathrm{ES}] + [\mathrm{EI}] + [\mathrm{ESI}]</math>


141行目: 246行目:
:<math>v = \frac{\frac{V_\mathrm{max}}{1 + [\mathrm{I}]/K_i} [\mathrm{S}]}{K_s + [\mathrm{S}]}</math>
:<math>v = \frac{\frac{V_\mathrm{max}}{1 + [\mathrm{I}]/K_i} [\mathrm{S}]}{K_s + [\mathrm{S}]}</math>


=== [[不競合阻害]] ===
=== 不競合阻害 ===
基質と阻害剤が酵素の異なる部位に結合するが、阻害剤は遊離の酵素には結合できず、酵素基質複合体のみに結合できる場合を不競合阻害(不競争阻害、不拮抗阻害とも)という。
基質と阻害剤が酵素の異なる部位に結合するが、阻害剤は遊離の酵素には結合できず、酵素基質複合体のみに結合できる場合を不競合阻害(不競争阻害、不拮抗阻害とも)という。


151行目: 256行目:
:<math>K_i = \frac{[\mathrm{ES}][\mathrm{I}]}{[\mathrm{ESI}]}</math>
:<math>K_i = \frac{[\mathrm{ES}][\mathrm{I}]}{[\mathrm{ESI}]}</math>


反応機構の中にある酵素種は E, ES, ESI の3種類なので、全酵素濃度は
反応機構の中にある酵素種はEESESIの3種類なので、全酵素濃度は
:<math>[\mathrm{E}]_0 = [\mathrm{E}] + [\mathrm{ES}] + [\mathrm{ESI}]</math>
:<math>[\mathrm{E}]_0 = [\mathrm{E}] + [\mathrm{ES}] + [\mathrm{ESI}]</math>


162行目: 267行目:
=== 混合型阻害 ===
=== 混合型阻害 ===
{{節スタブ}}
{{節スタブ}}

== 定数の決定 ==
定数<math>V_\max</math>および<math>K_\mathrm{M}</math>を決定する典型的な方法は、基質濃度 <math>[S]</math> を変化させて酵素アッセイを行い、初期反応速度 <math>v_0</math>を測定することである。ここでいう「初期」とは、比較的短い時間後に反応速度を測定することを意味し、その間、酵素-基質複合体は形成されているが、基質濃度はほぼ一定であり、平衡または準定常近似が有効であると仮定している<ref name="segel89" />。反応速度を濃度に対してプロットし、ミカエリス-メンテン方程式の[[非線形回帰]]を用いることで、パラメータを得ることができる<ref name="leskovac03">
{{Cite book
| last1 = Leskovac
| first1 = V.
| title = Comprehensive enzyme kinetics
| year = 2003
| publisher = Kluwer Academic/Plenum Pub.
| location = New York
| isbn = 978-0-306-46712-7
}}
</ref>。

非線形回帰を行うための計算機が利用できるようになる前は、方程式の線形化を伴うグラフ的な方法が用いられていた。[[イーディー=ホフステー図]]、[[ヘインズ=ウルフプロット]]、[[ラインウィーバー=バークプロット]]などが提案されており、これらの中ではヘインズ=ウルフプロットが最も正確である<ref name="leskovac03"/>。しかし、これら3つの方法は、可視化のためには有用であるが、データの誤差構造を歪めてしまい、非線形回帰には劣る<ref name="greco79">
{{cite journal
| last1 = Greco
| first1 = W.R.
| last2 = Hakala
| first2 = M.T.
| title = Evaluation of methods for estimating the dissociation constant of tight binding enzyme inhibitors
| journal = J Biol Chem
| volume = 254
| issue = 23
| pages = 12104–12109
| year = 1979
| doi = 10.1016/S0021-9258(19)86435-9
| pmid = 500698
| doi-access = free
}}
</ref>。 <math>v_0</math>について同様の[[誤差]]<math>dv_0</math>を仮定した場合、逆数表現では、<math>1/v_0</math>について<math>dv_0/v_0^2</math>の誤差が生じる({{仮リンク|不確かさの伝播|en|Propagation of uncertainty}})。<math>dv_0</math>の値の適切な見積りなしには、線形化を避けるべきである。また、[[最小二乗法]]を用いた[[回帰分析]]は、誤差が[[正規分布]]していることを前提としているが、これは<math>v_0</math>の値を変換した後は妥当ではない。にもかかわらず、現代の文献にはまだこれらの線形変換手法の使用が見られる<ref name="hayakawa06">
{{cite journal
| last1 = Hayakawa
| first1 = K.
| last2 = Guo
| first2 = L.
| last3 = Terentyeva
| first3 = E.A.
| last4 = Li
| first4 = X.K.
| last5 = Kimura
| first5 = H.
| last6 = Hirano
| first6 = M.
| last7 = Yoshikawa
| first7 = K.
| last8 = Nagamine
| first8 = T.
| last9 = Katsumata
| first9 = N.
| last10 = Ogata
| first10 = Tsutomu
| last11 = Tanaka
| first11 = Toshiaki
| title = Determination of specific activities and kinetic constants of biotinidase and lipoamidase in LEW rat and Lactobacillus casei (Shirota)
| journal = J Chromatogr B
| volume = 844
| issue = 2
| pages = 240–50
| year = 2006
| doi = 10.1016/j.jchromb.2006.07.006
| pmid = 16876490
| display-authors = 8
}}
</ref>。

その他、歴史的にはコーニッシュボーデンの直接的直線プロット<ref>{{cite journal|和書|title=原典からの酵素反応速度論|author=知名秀泰、岡田豊|journal=生物工学|volume=82|issue=1|year=2014|pages=20–25|url=https://www.sbj.or.jp/wp-content/uploads/file/sbj/9201/9201_yomoyama.pdf}}</ref>が知られている。

1997年、{{仮リンク|サンティアゴ・シュネル|en|Santiago Schnell}}とClaudio Mendozaは、ミカエリス・メンテン反応速度論の経時的反応速度解析のために、[[ランベルトのW関数]]の解に基づいた閉形式解を提案した<ref>{{cite journal |last1=Schnell |first1=S. |last2=Mendoza |first2=C. |title=A closed form solution for time-dependent enzyme kinetics |journal=Journal of Theoretical Biology |volume=187 |issue= 2|year=1997 |pages=207–212 |doi=10.1006/jtbi.1997.0425 |bibcode=1997JThBi.187..207S }}</ref>。

== 解析解 ==
ミカエリス・メンテン式は左辺に基質 S の濃度の[[時間微分]]を持つ[[微分方程式]]である。この方程式は[[変数分離]]形にして積分することで解くことができるが、解は[[初等関数]]で表現できる形にはならず、[[ランベルトのW関数]] ''W'' を用いて以下のように表される。
:<math chem>\frac{[\ce S]}{K_\mathrm{M}} = W(F(t)) </math>
上式において''W''はランベルトのW関数、
:<math chem>F(t) = \frac{[\ce S]_0}{K_\mathrm{M}} \exp\!\left(\frac{[\ce S]_0}{K_\mathrm{M}} - \frac{V_\max}{K_\mathrm{M}}\,t \right)</math>
である。

上記の方程式は経時データから<math>V_\max</math>と<math>K_\mathrm{M}</math>を推定するために使うことができる<ref>{{cite journal |first1=C. T. |last1=Goudar |first2=J. R. |last2=Sonnad |first3=R. G. |last3=Duggleby |year=1999 |title=Parameter estimation using a direct solution of the integrated Michaelis–Menten equation |journal=Biochimica et Biophysica Acta (BBA) - Protein Structure and Molecular Enzymology |volume=1429 |issue= 2|pages=377–383 |doi= 10.1016/s0167-4838(98)00247-7 |pmid=9989222}}</ref><ref>{{cite journal |first1=C. T. |last1=Goudar |first2=S. K. |last2=Harris |first3=M. J. |last3=McInerney |first4=J. M. |last4=Suflita |year=2004 |title=Progress curve analysis for enzyme and microbial kinetic reactions using explicit solutions based on the Lambert ''W'' function |journal=Journal of Microbiological Methods |volume=59 |issue=3 |pages=317–326 |doi=10.1016/j.mimet.2004.06.013 |pmid=15488275}}</ref>。


== 参考文献 ==
== 参考文献 ==

2021年10月15日 (金) 01:31時点における版

ミカエリス・メンテン式のプロット

ミカエリス・メンテン式(ミカエリス・メンテンしき、: Michaelis–Menten equation)とは、酵素の反応速度論に大きな業績を残したレオノール・ミカエリスモード・レオノーラ・メンテンにちなんだ、酵素の反応速度v に関する式で、

で表される。ここで、[P]は反応産物の濃度、[S]は基質濃度、Vmax は基質濃度が無限大のときの反応速度である。また、Kmミカエリス・メンテン定数と言い、v = Vmax /2(最大速度の半分の速度)を与える基質濃度を表す。この式をもとにしたモデルをミカエリス・メンテン動力学という。

この式により、反応速度v

  • 基質濃度が低い([S] ≪ Km )ときはその濃度に比例
  • 基質濃度が高い([S] ≫ Km )ときはその濃度に無関係に最大速度 Vmax に収束

となることが分かる[1]

導出

迅速平衡法による導出

酵素(以下E)が基質(以下S)と結合して酵素基質複合体(以下ES)を形成、ESがEとSに戻るか反応生成物(以下P)を生成する一連の反応機構を以下のように仮定する。

この反応は の2つの反応過程からできている。後者の反応を律速段階と仮定し反応速度定数k+2と設定する。 の反応は迅速に化学平衡に達していると仮定し、解離定数Ks と設定する。

仮定されている反応系に存在する酵素種は、基質と結合していない酵素Eと、基質Sと結合した酵素ESの2種類のみなので、全酵素濃度 [E]0 は両者の濃度の和に等しい。

[ES] を未知数として (1), (2) の連立方程式を解くと、

最初に仮定した反応機構では単位時間当たりに産生される反応産物Pの量は酵素基質複合体ESと速度定数 k +2 の積で与えられる。

(3) を (4) に代入して、

(4) 式から反応速度vは [ES] に比例することがわかるが、[ES] の最大値は (2) 式より [E]0である。したがって反応速度vの最大値Vmaxは次式となる。

(5), (6) より、

定常状態法による導出

系の別の解析がイギリスの植物学者ジョージ・エドワード・ブリッグズ英語版とイギリスの遺伝学者J・B・S・ホールデンによって1925年に行われた[2][3]。この迅速平衡法では が迅速に平衡に達すると仮定されているため、 の速度定数が の速度定数よりもはるかに小さい反応にしか成り立たない。定常状態法によって求めることで一般の反応でも同様の式が成り立つことが証明される。

反応機構は同様で、 について右向きの速度定数をk+1 、左向きの速度定数をk-1 とする。

定常状態では各酵素種の経時的濃度変化はないので、

この反応機構ではEとESしか酵素種が存在しないので

反応産物はESよりk+2の速度で生成されるので

(1) または (2) 式と (3) 式を連立方程式とみなして [ES] を求めると

(5) 式を (4) 式に代入して速度v を得た後、分子分母をk+1 で割る。

速度パラメーターとして

と定義すれば、(6) 式は

となる。

仮定と限界

導出における第一段階は、自由拡散を頼りにしている質量作用の法則を適用することである。しかし、高濃度のタンパク質が存在する生細胞の環境では、細胞質は自由に流れる液体というよりも粘性のあるゲルのように振る舞うことが多く、拡散による分子の動きが制限され、反応速度が変化する[4]。質量作用の法則は不均質な環境では有効であるが[5]、細胞質をフラクタルとしてモデル化する方が、その限定された移動性の動力学を表現するのに適している[6]

上記の2つのアプローチによって予測される結果の反応速度は類似しているが、唯一の違いは、迅速平衡近似では定数をと定義するのに対し、準定常近似ではを使用することである。しかしばがら、それぞれのアプローチは異なる仮定に基づいている。ミカエリス・メンテンの平衡解析は、基質が生成物の形成よりもはるかに速い時間スケールで平衡に達する場合、より正確には

の時に妥当である[7]

対照的にブリッグズ・ホールデン準定常状態解析はもし

ならば妥当である[8][9]

したがって、酵素濃度が基質濃度または、あるいはその両方よりもはるかに低い場合に成立する。

ミカエリス・メンテン解析、ブリッグズ・ホールデン解析のいずれにおいても、近似の質はが小さくなるにつれて向上する。しかし、モデル構築の際には、その前提条件を無視してミカエリス・メンテン速度論が用いられることが多い[7]

重要なことは、不可逆性は扱いやすい解析解を得るために必要な単純化であるが、一般的な場合では生成物の形成は実際には不可逆的ではないということである。酵素反応はより正確には次のように記述される。

一般的に、不可逆性の仮定は、以下のいずれかが真である状況で良い仮定である。

1. 基質の濃度が生成物の濃度よりも非常に大きい。

これは、標準的なin vitroのアッセイ条件でも真であるし、in vivoの多くの生物学的反応、特に生成物が後続の反応によって継続的に除去される場合にも真である。

2. 反応で放出されるエネルギーが非常に大きい、つまり

この2つの条件が成立しない場合(すなわち、反応が低エネルギーであり、生成物のプールがかなり存在する場合)、ミカエリス・メンテン方程式は破綻し、酵素の生物学を理解するためには、順反応と逆反応をあらわに考慮した、より複雑なモデリングアプローチが必要となる。

阻害がある場合のミカエリス・メンテン式

阻害とは何らかの理由で反応が遅くなることで、酵素反応の阻害には

などの種類がある。両逆数プロットを使うとこれらを見分けることができる。

基質阻害

次の図のように酵素基質複合体がさらに基質と結合して不活性となる場合を言う。

この場合の解離定数は

含まれる酵素種はE、ESとESSの3種類なので、全酵素濃度は

その他の値は迅速平衡法でミカエリス・メンテン式を求めた時のままで連立方程式を作り [ES] を求めると、

これを v = k+2 [ES] に代入しパラメーターVmax に変えると、

競合阻害

競合阻害(競争阻害、拮抗阻害とも)とは、基質と阻害剤(以下、I)が酵素の同じ活性中心に結合する場合に起こる阻害のこと。

反応機構は次式の通り。

基質と阻害剤の解離定数は

含まれる酵素種はE、ESとEIの3種類なので、全酵素濃度は

他の値は迅速平衡法で求めた時と同じで連立方程式を作り [ES] を求めると、

これを v = k+2 [ES] に代入しパラメーターVmax に変えると、

非競合阻害

基質と阻害剤が酵素の異なる部位に結合し、両者が互いに他の結合に影響を及ぼさない場合を非競合阻害(非競争阻害、非拮抗阻害とも)という。

反応機構は次の図の通り。

解離定数は、基質が遊離の酵素に結合するときも酵素阻害剤複合体に結合するときも同じで、阻害剤の結合定数も同様なので、

反応機構の中にある酵素種はE、ES、EIとESIの4種類なので、全酵素濃度は

この3つの式で連立方程式を作り、[ES] について求めると、

これを v = k+2 [ES] に代入しパラメーターVmax に変えると、

不競合阻害

基質と阻害剤が酵素の異なる部位に結合するが、阻害剤は遊離の酵素には結合できず、酵素基質複合体のみに結合できる場合を不競合阻害(不競争阻害、不拮抗阻害とも)という。

反応機構は次の図の通り。

解離定数は、

反応機構の中にある酵素種はE、ES、ESIの3種類なので、全酵素濃度は

この3つの式から連立方程式を作り、[ES] について求めると、

これを v = k+2 [ES] に代入しパラメーターVmax に変えると、

混合型阻害

定数の決定

定数およびを決定する典型的な方法は、基質濃度 を変化させて酵素アッセイを行い、初期反応速度 を測定することである。ここでいう「初期」とは、比較的短い時間後に反応速度を測定することを意味し、その間、酵素-基質複合体は形成されているが、基質濃度はほぼ一定であり、平衡または準定常近似が有効であると仮定している[8]。反応速度を濃度に対してプロットし、ミカエリス-メンテン方程式の非線形回帰を用いることで、パラメータを得ることができる[10]

非線形回帰を行うための計算機が利用できるようになる前は、方程式の線形化を伴うグラフ的な方法が用いられていた。イーディー=ホフステー図ヘインズ=ウルフプロットラインウィーバー=バークプロットなどが提案されており、これらの中ではヘインズ=ウルフプロットが最も正確である[10]。しかし、これら3つの方法は、可視化のためには有用であるが、データの誤差構造を歪めてしまい、非線形回帰には劣る[11]について同様の誤差を仮定した場合、逆数表現では、についての誤差が生じる(不確かさの伝播英語版)。の値の適切な見積りなしには、線形化を避けるべきである。また、最小二乗法を用いた回帰分析は、誤差が正規分布していることを前提としているが、これはの値を変換した後は妥当ではない。にもかかわらず、現代の文献にはまだこれらの線形変換手法の使用が見られる[12]

その他、歴史的にはコーニッシュボーデンの直接的直線プロット[13]が知られている。

1997年、サンティアゴ・シュネル英語版とClaudio Mendozaは、ミカエリス・メンテン反応速度論の経時的反応速度解析のために、ランベルトのW関数の解に基づいた閉形式解を提案した[14]

解析解

ミカエリス・メンテン式は左辺に基質 S の濃度の時間微分を持つ微分方程式である。この方程式は変数分離形にして積分することで解くことができるが、解は初等関数で表現できる形にはならず、ランベルトのW関数 W を用いて以下のように表される。

上式においてWはランベルトのW関数、

である。

上記の方程式は経時データからを推定するために使うことができる[15][16]

参考文献

  1. ^ 大竹伝雄『化学工学Ⅲ』岩波書店、1978年、16頁。ISBN 4-00-021103-X 
  2. ^ Briggs, G.E.; Haldane, J.B.S. (1925). “A note on the kinematics of enzyme action”. Biochem J 19 (2): 338–339. doi:10.1042/bj0190338. PMC 1259181. PMID 16743508. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1259181/. 
  3. ^ Laidler, Keith J. (1978). Physical Chemistry with Biological Applications. Benjamin/Cummings. pp. 428–430. ISBN 0-8053-5680-0 
  4. ^ Zhou, H.X.; Rivas, G.; Minton, A.P. (2008). “Macromolecular crowding and confinement: biochemical, biophysical, and potential physiological consequences”. Annu Rev Biophys 37 (1): 375–97. doi:10.1146/annurev.biophys.37.032807.125817. PMC 2826134. PMID 18573087. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2826134/. 
  5. ^ Grima, R.; Schnell, S. (Oct 2006). “A systematic investigation of the rate laws valid in intracellular environments”. Biophys Chem 124 (1): 1–10. doi:10.1016/j.bpc.2006.04.019. PMID 16781049. 
  6. ^ Schnell, S.; Turner, T.E. (2004). “Reaction kinetics in intracellular environments with macromolecular crowding: simulations and rate laws”. Prog Biophys Mol Biol 85 (2–3): 235–60. doi:10.1016/j.pbiomolbio.2004.01.012. PMID 15142746. 
  7. ^ a b Keener, J.; Sneyd, J. (2008). Mathematical Physiology: I: Cellular Physiology (2 ed.). Springer. ISBN 978-0-387-75846-6 
  8. ^ a b 引用エラー: 無効な <ref> タグです。「segel89」という名前の注釈に対するテキストが指定されていません
  9. ^ Murray, J.D. (2002). Mathematical Biology: I. An Introduction (3 ed.). Springer. ISBN 978-0-387-95223-9 
  10. ^ a b Leskovac, V. (2003). Comprehensive enzyme kinetics. New York: Kluwer Academic/Plenum Pub.. ISBN 978-0-306-46712-7 
  11. ^ Greco, W.R.; Hakala, M.T. (1979). “Evaluation of methods for estimating the dissociation constant of tight binding enzyme inhibitors”. J Biol Chem 254 (23): 12104–12109. doi:10.1016/S0021-9258(19)86435-9. PMID 500698. 
  12. ^ Hayakawa, K.; Guo, L.; Terentyeva, E.A.; Li, X.K.; Kimura, H.; Hirano, M.; Yoshikawa, K.; Nagamine, T. et al. (2006). “Determination of specific activities and kinetic constants of biotinidase and lipoamidase in LEW rat and Lactobacillus casei (Shirota)”. J Chromatogr B 844 (2): 240–50. doi:10.1016/j.jchromb.2006.07.006. PMID 16876490. 
  13. ^ 知名秀泰、岡田豊「原典からの酵素反応速度論」『生物工学』第82巻第1号、2014年、20–25頁。 
  14. ^ Schnell, S.; Mendoza, C. (1997). “A closed form solution for time-dependent enzyme kinetics”. Journal of Theoretical Biology 187 (2): 207–212. Bibcode1997JThBi.187..207S. doi:10.1006/jtbi.1997.0425. 
  15. ^ Goudar, C. T.; Sonnad, J. R.; Duggleby, R. G. (1999). “Parameter estimation using a direct solution of the integrated Michaelis–Menten equation”. Biochimica et Biophysica Acta (BBA) - Protein Structure and Molecular Enzymology 1429 (2): 377–383. doi:10.1016/s0167-4838(98)00247-7. PMID 9989222. 
  16. ^ Goudar, C. T.; Harris, S. K.; McInerney, M. J.; Suflita, J. M. (2004). “Progress curve analysis for enzyme and microbial kinetic reactions using explicit solutions based on the Lambert W function”. Journal of Microbiological Methods 59 (3): 317–326. doi:10.1016/j.mimet.2004.06.013. PMID 15488275. 

関連項目

外部リンク