Aritalab:Lecture/Biochem/Modularity/NF

From Metabolomics.JP
< Aritalab:Lecture | Biochem | Modularity(Difference between revisions)
Jump to: navigation, search
m
m (MAPK カスケードと負のフィードバック)
 
(4 intermediate revisions by one user not shown)
Line 1: Line 1:
 +
__FORCETOC__
 +
 
==MAPK カスケードと負のフィードバック==
 
==MAPK カスケードと負のフィードバック==
 
ここでは Cell Designer ソフトウェアを用いて、MAPK カスケードのシミュレーションを実感することを目標としています。
 
ここでは Cell Designer ソフトウェアを用いて、MAPK カスケードのシミュレーションを実感することを目標としています。
Cell Designer ソフトウェアのサンプルファイルの中から、 MAPK を選択してください。
+
Cell Designer ソフトウェアのデータベースBioModelsから、 MAPK を検索してBIOMD0000000010を選択してください。
 
{|
 
{|
 
|valign="top"|
 
|valign="top"|
 
Mitogen-activated Protein Kinase (MAPK) は全身の細胞で発現し、細胞外のシグナルを核内に伝える代表的なシグナル伝達経路です。
 
Mitogen-activated Protein Kinase (MAPK) は全身の細胞で発現し、細胞外のシグナルを核内に伝える代表的なシグナル伝達経路です。
MAPK とは総称で、代表的なものにERK1/2などがあります。ここでは MAPK が MAPキナーゼキナーゼ (MAPKK) により 2 段階のリン酸化を受け、MAPKK はMAPKKキナーゼにより 2 段階のリン酸化を受けるモデルになっています。最終的にできるMAPK二リン酸は、MKKK を抑制しています。
+
MAPK とは総称で、代表的なものにERK1/2があります。ここでは MAPK が MAPキナーゼキナーゼ (MAPKK) により 2 段階のリン酸化を受け、MAPKK はMAPKKキナーゼにより 2 段階のリン酸化を受けるモデルになっています。最終的にできるMAPK二リン酸は、MKKK を抑制しています。
左図の右下から左上に戻ってくる経路 (J0という反応を┴型で抑制している) が、抑制効果を表しています。
+
左図の右下から左上に戻ってくる経路 (抑制する反応を┴型で表現) が、抑制効果を表しています。
 
このサイクルには、一つだけ「抑制」の効果があるので、負のフィードバックループになっています。
 
このサイクルには、一つだけ「抑制」の効果があるので、負のフィードバックループになっています。
 
|valign="top"|
 
|valign="top"|
Line 13: Line 15:
 
{|
 
{|
 
|valign="top"|
 
|valign="top"|
とりあえず経路全体をシミュレーションしてみます。8種のタンパク質が複雑な動きをすることがわかります。
+
BioModelsからダウンロードしたモデルが上のモデルと同じかどうかは一見わかりにくいですが、エディタでタンパク質の位置関係を整理すると、同じモデルであることがわかります。
 +
* MAPK  ... Erk2
 +
* MKK  ... Mek1
 +
* MKKK ... Mos
 +
という対応関係です。最後のMAPKによるフィードバック部分は
 +
MKKK / ((1 + pow(MAPK_PP / Ki, n)) * (K1 + MKKK))
 +
というヒル式になっています。
 +
|valign="top"|
 +
[[File:0000000010.png|400px|MAPK pathway by Cell Designer]]
 +
|}
 +
 
 +
{|
 +
|valign="top"|
 +
とりあえず経路全体をシミュレーションしてみます(Simulationパネル)。8種のタンパク質が複雑な動きをすることがわかります。
  
初期値をみると MKKK, MKK, MAPK がそれぞれ 90, 280, 280 になっています。(3種に色分けされた行方向のグループについて、それぞれ一番左の要素です。) それ以外の初期値はすべて 10 になっていて、それぞれタンパク質のリン酸化状態を表現しています。MKK, MAPK は2段階のリン酸化がおこなわれるのでPがはいる丸印がボックスに2箇所ずつ用意されています。
+
初期値をみると Mos, Mek1, Erk2 がそれぞれ 90, 280, 280 になっています。(リン酸化されていない初期状態。) シミュレーション結果のグラフでは、黄色、赤、緑色の線がそれぞれ 90, 280, 280 の値からスタートしているのがわかります。
 +
残りの初期値はすべて 10 になっていて、全てリン酸化されたタンパク質を表現しています。Mek1, Erk2 は2段階のリン酸化がおこなわれます。
 
|valign="top"|
 
|valign="top"|
 
[[File:Lecture-Biochem-Modularity-MAPKsim.png|300px|MKKK and MKK]]
 
[[File:Lecture-Biochem-Modularity-MAPKsim.png|300px|MKKK and MKK]]
Line 22: Line 38:
 
{|
 
{|
 
|valign="top"|
 
|valign="top"|
ここで一番上の行(黄色)に注目します、つまりMKKKのリン酸化だけを考えてみましょう。一番上の行(黄色)だけを残して残りを消去し、シミュレーションをおこなってみます。
+
ここでフィードバックに当たる部分のヒル式を MKKK / (K1 + MKKK) 形に修正し、フィードバックをなくします。ループがないため、MosとMos-P(リン酸化された状態)が平衡状態、つまり量の変動がない状態に落ち着きます。
フィードバックループがないため、MKKKとMKKK-P(リン酸化された状態)が平衡状態、つまり量の変動がない状態に落ち着くはずです。
+
  
左図で黄色が MKKK-P, 赤色が MKKK です。全ての MKKK がリン酸化されてそのまま動かないことがわかります。
+
左図で黄色が Mos-P, 赤色が Mos です。全ての Mos がリン酸化されてそのまま動かないことがわかります。
MKKKとMKKK-Pはそれぞれ 90, 10 が初期値だったので合計した 100 が最終的な MKKK-P の量になっています。
+
Mos とMos-Pはそれぞれ 90, 10 が初期値だったので合計した 100 が最終的な Mos-P の量になっています。
 
|valign="top"|
 
|valign="top"|
 
[[File:Lecture-Biochem-Modularity-MKKKonly.png|300px|MAPK pathway by Cell Designer]]
 
[[File:Lecture-Biochem-Modularity-MKKKonly.png|300px|MAPK pathway by Cell Designer]]
Line 33: Line 48:
 
{|
 
{|
 
|valign="top"|
 
|valign="top"|
つぎにMKKKとMKKだけに注目したネットワークを作って動きを見てみます。ネットワーク図でいうと、黄色と緑の要素だけ残してシミュレーションをおこないます。フィードバックループがないため、やはり量の変動がない状態に落ち着くはずです。
+
つぎにシミュレーションパネルのチェックボックスを利用して、Mos, Mek1だけの表示にしてみます。フィードバックループがないため、やはり量の変動がない状態に落ち着くはずです。
 
|valign="top"|
 
|valign="top"|
 
[[File:Lecture-Biochem-Modularity-MKKKandMKK.png|300px|MKKK and MKK]]
 
[[File:Lecture-Biochem-Modularity-MKKKandMKK.png|300px|MKKK and MKK]]
Line 39: Line 54:
 
{|
 
{|
 
|valign="top"|
 
|valign="top"|
上の図にあった黄色と赤のライン (MKKK) はこの図において黄色と緑になっています。一つ前のシミュレーション結果とはy軸方向のスケーリングが異なることに注意します。MKKの三つが左から順にそれぞれ赤、水色、青です。赤と青に注目してください。MKKKの動きと同じく、相補的な挙動をみせています。MKK, MKK-P, MKK-PPの初期値はそれぞれ 280, 10, 10 だったので合計した 300 に近い値が MKK-PP (青)の最終値です。 また MKK-P (水色)は初めに量が増えるものの、最終的には小さい値で平衡を保つことになります。
+
Mosとそのリン酸化体、Mek1とそのリン酸化体はいずれも相補的な挙動をみせています。二段階リン酸化を経る場合、一段階目のリン酸化体は量が大きく増えない点にも注意しましょう。
 
|valign="top"|
 
|valign="top"|
 
[[File:Lecture-Biochem-Modularity-MKKKandMKKsim.png|300px|MKKK and MKK]]
 
[[File:Lecture-Biochem-Modularity-MKKKandMKKsim.png|300px|MKKK and MKK]]
Line 45: Line 60:
 
{|
 
{|
 
|valign="top"|
 
|valign="top"|
ここまでの知見はMKKK, MKK, MAPK とカスケードが3段階になっても同じ動きのはずです。(何段階になっても同じ。)フィードバックがない状態では、全くリン酸化されていないMAPK と 二回リン酸化された MAPK の関係は、初期状態から量が逆転しておわるだけと予想されます。フィードバックがないだけのネットワークをシミュレーションしてみます。
+
ここまでの知見はカスケードが何段階になっても同じです。フィードバックがない状態では、それぞれのタンパク質がリン酸化されておわります。
 
|valign="top"|
 
|valign="top"|
 
[[File:Lecture-Biochem-Modularity-MKKKandMKKandMAPK.png|300px|MKKK and MKK and MAPK]]
 
[[File:Lecture-Biochem-Modularity-MKKKandMKKandMAPK.png|300px|MKKK and MKK and MAPK]]
Line 52: Line 67:
 
|valign="top"|
 
|valign="top"|
 
フィードバックがない状態で8種類のタンパク質を観測すると
 
フィードバックがない状態で8種類のタンパク質を観測すると
: MKKK: 黄色と水色
+
: Mos: 黄色と水色
: MKK: 赤とマゼンタ
+
: Mek1: 赤とマゼンタ
: MAPK: 緑と薄いグレー
+
: Erk2: 緑と薄いグレー
という3組の動きがどれも似ていることがわかります。(ちょっとわかりにくいですが色を確認してください。)
+
という3組の動きがどれも似て、いずれも相補的に動いていることがわかります。(ちょっとわかりにくいですが色を確認してください。)
量が変化する時間をみると、その順序は MKKK &rarr; MKK &rarr; MAPK です。MAPKの薄いグレーだけは増加のスピードが速くマゼンタを追い抜いて量が増えますが、増加し始めるタイミングはマゼンタより遅いことに注目してください。
+
量が変化しはじめる時間をみると、その順序は Mos &rarr; Mek1 &rarr; Erk2 です。Mek2の薄いグレーだけは増加のスピードが速くマゼンタを追い抜いて量が増えますが、増加し始めるタイミングはマゼンタより遅いことに注目してください。
 
|valign="top"|
 
|valign="top"|
 
[[File:Lecture-Biochem-Modularity-MKKKandMKKandMAPKsim.png|300px|MKKK and MKK and MAPK]]
 
[[File:Lecture-Biochem-Modularity-MKKKandMKKandMAPKsim.png|300px|MKKK and MKK and MAPK]]
Line 63: Line 78:
 
{|
 
{|
 
|valign="top"|
 
|valign="top"|
ここで最初のネットワークに戻して、フィードバックを考慮した状態でシミュレーションをやってみましょう。
+
ここで最初のネットワークに戻って、フィードバックを考慮した状態でシミュレーションをやってみましょう。
 
黄色と水色、赤とマゼンタ、緑と薄いグレーが相補的な動きを見せるところは一緒です。
 
黄色と水色、赤とマゼンタ、緑と薄いグレーが相補的な動きを見せるところは一緒です。
ネットワーク図における MKKK, MKK, MAPK の各グループをモジュールと考え、代表タンパク質 (MKKK-P, MKK-PP, MAPK-PP) だけを表示させてみます。
+
ネットワーク図における各タンパク質をモジュールと考え、代表タンパク質 (Mos-P, Mek1-PP, Erk2-PP) だけを表示させてみます。
 
この3者がネットワークの挙動を示す鍵となっています。
 
この3者がネットワークの挙動を示す鍵となっています。
薄いグレー色が上がってくるにつれて量が増えていた水色 (MKKK-P) が減ってくることに注意してください。
+
薄いグレー色が上がってくるにつれて量が増えていた水色 (Mos-P) が減ってくることに注意してください。
水色が減るということは、これと相補的な挙動をする黄色 (MKKK) が増えるということで、シミュレーションの初期に近い状態に戻っていきます。
+
水色が減るということは、これと相補的な挙動をする黄色 (Mos) が増えるということで、シミュレーションの初期に近い状態に戻っていきます。
 
|
 
|
 
|valign="top"|
 
|valign="top"|
Line 75: Line 90:
  
 
ここまで読んできた人は、最初のシミュレーション結果を複雑で難しいとは思わなくなります。なぜなら、3つのモジュールに分けて考えればよく、水色、マゼンタ、薄いグレーの挙動だけ把握すればあとは補えることがわかったからです。
 
ここまで読んできた人は、最初のシミュレーション結果を複雑で難しいとは思わなくなります。なぜなら、3つのモジュールに分けて考えればよく、水色、マゼンタ、薄いグレーの挙動だけ把握すればあとは補えることがわかったからです。
下の図で左が「難しく見える」シミュレーション結果です。このうち、3モジュールに注目すると真ん中のようになります。結局、3つのタンパク質グループが順番に量を増加させ、順番に量を減少させるだけです。
+
下の図で左が「難しく見える」シミュレーション結果です。このうち、3モジュールに注目すると真ん中のようになります。結局、3つのタンパク質グループ(ネットワーク図における黄色、黄緑、青の3種)が順番に量を増加させ、順番に量を減少させるだけです。
 
ですから、時間軸を長くとると右図のようになり、振動していることがわかります。
 
ですから、時間軸を長くとると右図のようになり、振動していることがわかります。
  

Latest revision as of 04:35, 7 June 2019


Contents

[edit] MAPK カスケードと負のフィードバック

ここでは Cell Designer ソフトウェアを用いて、MAPK カスケードのシミュレーションを実感することを目標としています。 Cell Designer ソフトウェアのデータベースBioModelsから、 MAPK を検索してBIOMD0000000010を選択してください。

Mitogen-activated Protein Kinase (MAPK) は全身の細胞で発現し、細胞外のシグナルを核内に伝える代表的なシグナル伝達経路です。 MAPK とは総称で、代表的なものにERK1/2があります。ここでは MAPK が MAPキナーゼキナーゼ (MAPKK) により 2 段階のリン酸化を受け、MAPKK はMAPKKキナーゼにより 2 段階のリン酸化を受けるモデルになっています。最終的にできるMAPK二リン酸は、MKKK を抑制しています。 左図の右下から左上に戻ってくる経路 (抑制する反応を┴型で表現) が、抑制効果を表しています。 このサイクルには、一つだけ「抑制」の効果があるので、負のフィードバックループになっています。

MAPK pathway by Cell Designer

BioModelsからダウンロードしたモデルが上のモデルと同じかどうかは一見わかりにくいですが、エディタでタンパク質の位置関係を整理すると、同じモデルであることがわかります。

  • MAPK ... Erk2
  • MKK ... Mek1
  • MKKK ... Mos

という対応関係です。最後のMAPKによるフィードバック部分は

MKKK / ((1 + pow(MAPK_PP / Ki, n)) * (K1 + MKKK))

というヒル式になっています。

MAPK pathway by Cell Designer

とりあえず経路全体をシミュレーションしてみます(Simulationパネル)。8種のタンパク質が複雑な動きをすることがわかります。

初期値をみると Mos, Mek1, Erk2 がそれぞれ 90, 280, 280 になっています。(リン酸化されていない初期状態。) シミュレーション結果のグラフでは、黄色、赤、緑色の線がそれぞれ 90, 280, 280 の値からスタートしているのがわかります。 残りの初期値はすべて 10 になっていて、全てリン酸化されたタンパク質を表現しています。Mek1, Erk2 は2段階のリン酸化がおこなわれます。

MKKK and MKK

ここでフィードバックに当たる部分のヒル式を MKKK / (K1 + MKKK) 形に修正し、フィードバックをなくします。ループがないため、MosとMos-P(リン酸化された状態)が平衡状態、つまり量の変動がない状態に落ち着きます。

左図で黄色が Mos-P, 赤色が Mos です。全ての Mos がリン酸化されてそのまま動かないことがわかります。 Mos とMos-Pはそれぞれ 90, 10 が初期値だったので合計した 100 が最終的な Mos-P の量になっています。

MAPK pathway by Cell Designer

つぎにシミュレーションパネルのチェックボックスを利用して、Mos, Mek1だけの表示にしてみます。フィードバックループがないため、やはり量の変動がない状態に落ち着くはずです。

MKKK and MKK

Mosとそのリン酸化体、Mek1とそのリン酸化体はいずれも相補的な挙動をみせています。二段階リン酸化を経る場合、一段階目のリン酸化体は量が大きく増えない点にも注意しましょう。

MKKK and MKK

ここまでの知見はカスケードが何段階になっても同じです。フィードバックがない状態では、それぞれのタンパク質がリン酸化されておわります。

MKKK and MKK and MAPK

フィードバックがない状態で8種類のタンパク質を観測すると

Mos: 黄色と水色
Mek1: 赤とマゼンタ
Erk2: 緑と薄いグレー

という3組の動きがどれも似て、いずれも相補的に動いていることがわかります。(ちょっとわかりにくいですが色を確認してください。) 量が変化しはじめる時間をみると、その順序は Mos → Mek1 → Erk2 です。Mek2の薄いグレーだけは増加のスピードが速くマゼンタを追い抜いて量が増えますが、増加し始めるタイミングはマゼンタより遅いことに注目してください。

MKKK and MKK and MAPK

ここで最初のネットワークに戻って、フィードバックを考慮した状態でシミュレーションをやってみましょう。 黄色と水色、赤とマゼンタ、緑と薄いグレーが相補的な動きを見せるところは一緒です。 ネットワーク図における各タンパク質をモジュールと考え、代表タンパク質 (Mos-P, Mek1-PP, Erk2-PP) だけを表示させてみます。 この3者がネットワークの挙動を示す鍵となっています。 薄いグレー色が上がってくるにつれて量が増えていた水色 (Mos-P) が減ってくることに注意してください。 水色が減るということは、これと相補的な挙動をする黄色 (Mos) が増えるということで、シミュレーションの初期に近い状態に戻っていきます。

MKKK-P and MKK-PP and MAPK-PP only

ここまで読んできた人は、最初のシミュレーション結果を複雑で難しいとは思わなくなります。なぜなら、3つのモジュールに分けて考えればよく、水色、マゼンタ、薄いグレーの挙動だけ把握すればあとは補えることがわかったからです。 下の図で左が「難しく見える」シミュレーション結果です。このうち、3モジュールに注目すると真ん中のようになります。結局、3つのタンパク質グループ(ネットワーク図における黄色、黄緑、青の3種)が順番に量を増加させ、順番に量を減少させるだけです。 ですから、時間軸を長くとると右図のようになり、振動していることがわかります。

MKKK and MKK and MAPK MKKK and MKK and MAPK MKKK and MKK and MAPK

[edit] まとめ

負のフィードバックは、振動現象を実現する。
Personal tools
Namespaces

Variants
Actions
Navigation
metabolites
Toolbox