Positive feedback, bistability, and hysteresis

Slides:



Advertisements
Similar presentations
シミュレーション演習 G. 総合演習 ( Mathematica 演 習) システム創成情報工学科 テキスト作成: 藤尾 光彦 講義担当: 尾下 真樹.
Advertisements

Division of Process Control & Process Systems Engineering Department of Chemical Engineering, Kyoto University
情報基礎実習 I (第6回) 木曜4・5限 担当:北川 晃. Stream クラスを用いたファイルの接続 … Dim インスタンス名 As New IO.StreamReader( _ “ ファイルの絶対パス ”, _ System.Text.Encoding.Default) … s = インスタンス名.
統計学 第3回 西山. 第2回のまとめ 確率分布=決まっている分布の 形 期待値とは平均計算 平均=合計 ÷ 個数から卒業! 平均=割合 × 値の合計 同じ平均値でも 同じ分散や標準偏差でも.
陰関数定理と比較静学 モデルの連立方程式体系で表されるとき パラメータが変化したとき 如何に変数が変化するか 至るところに出てくる.
円線図とは 回路の何らかの特性を複素平面上の円で表したもの 例えば、ZLの変化に応じてZinが変化する様子 Zin ZL
Ruth Onn, Alfred Bruckstein (Int J Comp Vision 1990)
第3回:ボールを上下に動かそう! (オブジェクトの移動、一次元)
Fortran と有限差分法の 入門の入門の…
・力のモーメント ・角運動量 ・力のモーメントと角運動量の関係
伝達事項 皆さんに数学と物理の全国統一テストを受けても らいましたが、この時の試験をまた受けていただ きます。
プログラミング入門 電卓番外編 ~エクセルで関数表示~.
コリオリ力の復習資料 見延 庄士郎(海洋気候物理学研究室)
ゲーム理論・ゲーム理論Ⅰ (第8回) 第5章 不完全競争市場の応用
(Fri) Astrophysics Laboratory MATSUO Kei
4.3 連立1次方程式   Ax = b   (23) と書くことができる。
3 一次関数 1章 一次関数とグラフ §3 一次関数の式を求めること          (3時間).
スペクトル法による数値計算の原理 -一次元線形・非線形移流問題の場合-
C言語 配列 2016年 吉田研究室.
基礎プログラミング (第五回) 担当者: 伊藤誠 (量子多体物理研究室) 内容: 1. 先週のおさらいと続き (実習)
大阪工業大学 情報科学部 情報システム学科 宇宙物理研究室 B 木村悠哉
次に 円筒座標系で、 速度ベクトルと加速度ベクトルを 求める.
ブロック線図によるシミュレーション ブロック線図の作成と編集 ブロック線図の保存と読込み ブロック線図の印刷 グラフの印刷
第3章 重回帰分析 ー 計量経済学 ー.
第3章 重回帰分析 ー 計量経済学 ー.
多変数関数の積分(6/3~24) 重積分(2重積分) 第6章(§5は除く) 重積分の定義 「連続関数は積分可能」
実習4 2次元テーブルの利用 フローチャートの作成.
4.2 連立非線形方程式 (1)繰返し法による方法
不等式の表す領域  直線で分けられた領域.
透視投影(中心射影)とは  ○ 3次元空間上の点を2次元平面へ投影する方法の一つ  ○ 投影方法   1.投影中心を定義する   2.投影平面を定義する
流体のラグランジアンカオスとカオス混合 1.ラグランジアンカオス 定常流や時間周期流のような層流の下での流体の微小部分のカオス的運動
第3回:ボールを上下に動かそう! (オブジェクトの移動、一次元)
重力レンズ効果を想定した回転する ブラックホールの周りの粒子の軌道
(ラプラス変換の復習) 教科書には相当する章はない
シミュレーション演習 狙い 初日 イントロ 2日目 プローブと信号 3日目 実験結果からの数理モデル作成
生物機能工学基礎実験 2.ナイロン66の合成・糖の性質 から 木村 悟隆
動力学(Dynamics) 運動方程式のまとめ 2008.6.17
Lorenz modelにおける 挙動とそのカオス性
勉強会 「紙と鉛筆で学ぶシステム生物学の数理」
システムモデルと伝達関数 1. インパルス応答と伝達関数 キーワード : 伝達関数、インパルス応答、 ステップ応答、ランプ応答
アーランの即時式モデル.
実習アウトライン Matlabを用いた生化学反応シミュレーション
形式言語とオートマトン Formal Languages and Automata 第4日目
形式言語とオートマトン Formal Languages and Automata 第4日目
奈良女子大集中講義 バイオインフォマティクス (9) 相互作用推定
シミュレーション演習 狙い 初日 イントロ 2日目 プローブと信号 3日目 実験結果からの数理モデル作成
磁気リコネクション (Craig-Henton解の安定性) ~シミュレーションサマースクール@千葉大より~
独立成分分析 5 アルゴリズムの安定性と効率 2007/10/24   名雪 勲.
6. ラプラス変換.
今後の予定 4日目 10月22日(木) 班編成の確認 講義(2章の続き,3章) 5日目 10月29日(木) 小テスト 4日目までの内容
基本システムのボード線図 ボード線図による基本システムの同定
デジタル画像とC言語.
Data Clustering: A Review
9.通信路符号化手法1 (誤り検出と誤り訂正の原理)
シミュレーション演習 狙い 初日 イントロ 2日目 プローブと信号 3日目 実験結果からの数理モデル作成
プログラムの基本構造と 構造化チャート(PAD)
4. システムの安定性.
「データ学習アルゴリズム」 第3章 複雑な学習モデル 報告者 佐々木 稔 2003年6月25日 3.1 関数近似モデル
計測工学 計測工学8 最小二乗法3 計測工学の8回目です。 最小二乗法を簡単な一時関数以外の関数に適用する方法を学びます。
データ解析 静岡大学工学部 安藤和敏
22・3 積分形速度式 ◎ 速度式: 微分方程式 ⇒ 濃度を時間の関数として得るためには積分が必要
【第六講義】非線形微分方程式.
卒論中間発表 2001/12/21 赤道の波動力学の基礎 北海道大学理学部 地球科学科 4年 山田 由貴子.
モデルの微分による非線形モデルの解釈 明治大学 理工学部 応用化学科 データ化学工学研究室 金子 弘昌.
パターン認識特論 カーネル主成分分析 和田俊和.
コンパイラ 2012年10月11日
3 一次関数 1章 一次関数とグラフ §4 方程式とグラフ         (3時間).
* Ehime University, Japan
信号データの変数代入と変数参照 フィードバック制御系の定常特性 フィードバック制御系の感度特性
コンピュータの高速化により, 即座に計算できるようになってきたが, 手法的にはコンピュータ出現以前に考え出された 方法が数多く使われている。
Presentation transcript:

Positive feedback, bistability, and hysteresis 柚木克之 (YUGI, Katsuyuki) Kuroda Lab., The University of Tokyo

実習課題:大腸菌lacオペロンの双安定性 Ozbudakモデルを動かしてヒ ステリシスを確認する 大腸菌は糖の濃度を記憶する ヒステリシスには次の2つの要 素が必要なことを確認 ポジティブ・フィードバック 協調作用 ヒステリシスが起きる理由を 力学系の言葉で説明 Ozbudak et al. (2004) Nature

Ozbudakらは2変数の式で ポジティブ・フィードバックをモデル化した TMGout (Xout) + TMGin (x) LacY-GFP (y) LacI Ozbudak et al. (2004) Nature の モデルを簡略化したもの TMG (methyl-β-D-thiogalactopyranoside) LacYによって取り込まれ、LacIを阻害する

Ozbudakモデルの変数と定数の意味 Variables Constants x: TMGin α, β: Rate constants, y: LacY-GFP α, β: Rate constants, r :Tightness of LacI repression

微分方程式と生化学反応系の対応 x, yの分解 TMGout (Xout) + TMGin (x) LacY-GFP (y) LacI

微分方程式と生化学反応系の対応 x, yの分解 TMGout (Xout) + TMGin (x) LacY-GFP (y) LacI

微分方程式と生化学反応系の対応 x, yの分解 TMGout (Xout) + TMGin (x) LacY-GFP (y) LacI

微分方程式と生化学反応系の対応 x, yの分解 TMGout (Xout) + TMGin (x) LacY-GFP (y) LacI

演習1:まずはOzbudak modelを動かそう

演習1(MATLAB): Ozbudak modelを動かす 初期値(変数) : x(t=0) = 0, y(t=0) = 0 パラメータ(定数): α = 0.1, β = 10, ρ = 25, Xout = 5

ファイルをDLして空欄を埋めなさい シミュレーション モデル function ozbudak( ) time = 0:0.1:30; x_init = 0; y_init = 0; alpha = 0.1; beta = 10; rho = 25; x_out = 5; s0 = [ x_init , y_init ]; % Initial values param = [ alpha , beta , rho , x_out ]; % Constants [t,time_course] = ode15s(@(t,s) ODE(t,s,param),time,s0); plot(t,time_course) end function dsdt = ODE(t,s,param) x = s(1); y = s(2); alpha = param(1); beta = param(2); rho = param(3); x_out = param(4); dsdt(1,:) = dsdt(2,:) = シミュレーション モデル

Dimensionless concenration うまくいくとこうなるはず x (TMGin) y (LacY-GFP) time Dimensionless concenration

演習2:Ozbudak modelがヒステリシスを 示すことを確認する

t = 60 時点で、TMGout (Xout) を 0 にする 実験的には、培地から集菌して wash Xout の初期値によって、ヒステリシスが生じない場合と生じる場合 がある ヒステリシスあり Xout≧8 Xout y (LacY-GFP) ヒステリシスなし Xout≦7 time time

演習2: 空欄を埋めなさい function ozbudak_hysteresis_timeseries( ) (略) [t,time_course] = ode15s(@(t,s) ODE(t,s,param),time,s0); plot(t,time_course(:,2)) % y のみプロットする hold on % 重ね描きを許可する命令 x_init = % time_course の終端値を x_init, y_init にセット y_init = s0 = [ x_init , y_init ]; x_out = % TMGout の値を 0 にする param = [ alpha , beta , rho , x_out ]; time = 60:0.1:120; plot(t,time_course(:,2)) end function dsdt = ODE(t,s,param) (略)

出力例 Xout (t=0) = 5 Xout (t=0) = 10 y (LacY-GFP) y (LacY-GFP) ヒステリシスなし ヒステリシスあり time time

演習3:ヒステリシスが起きるTMGの濃度を図示する

時系列を dose response curve に描き直す y (LacY-GFP) y (LacY-GFPend) Xout(TMGout) Ozbudak et al. (2004) Nature Xout X軸: TMGout と Y軸: LacY-GFPの最終値をプロット

演習3(MATLAB): ヒステリシスのシミュレーション 右図の実験に相当するシミ ュレーションを行いなさい 右図下段は次のような実験 TMG = X μM で培養 集菌 TMG = (X + ΔX) μM で培 養 以後この繰り返し 集菌の際、xはリセットされ るが、yは保存される y (LacY-GFPend) Xout(TMGout) Ozbudak et al. (2004) Nature

空欄を埋めなさい function ozbudak_fig2b( ) time = 0:0.1:30; x_init = 0; y_init = 0; alpha = 0.1; beta = 10; rho = 25; x_out = 0; param = [ alpha , beta , rho , x_out ]; % Constants figure(1) hold on while x_out <= 20; % このwhile文は Fig.2b 下段に相当 s0 = [ x_init , y_init ]; % Initial values [t,time_course] = ode15s(@(t,s) ODE(t,s,param),time,s0); y_end = % yの時系列の最終値(集菌時の値)を y_end に代入 plot( x_out , y_end , 'ko'); x_out = x_out + 1; % y_endを、次のシミュレーションの初期値にする end while x_out >= 0 % このwhileループには Fig.2b 上段に相当する内容を書く % ここには上の while ループと似た内容を書く。ただし % プロットする際、点は赤の十字(‘r+’)とする。 function dsdt = ODE(t,s,param) (以下略)

○はTMG徐々に増加、+はTMG減少 y (LacY-GFPend) Xout(TMGout)

演習4:ポジティブ・フィードバックを切ると ヒステリシスは起きなくなる

ポジティブ・フィードバックを遮断する x, yの分解 TMGout (Xout) + TMGin (x) LacY-GFP (y) LacI

ポジティブ・フィードバックを遮断する x, yの分解 TMGout (Xout) TMGin (x) GFP (y) LacI

ポジティブ・フィードバックを遮断すると ヒステリシスを示さなくなる Xout増加、減少のどちらの場合も同じS字カーブを描く

演習5:協調作用がなくなると ヒステリシスは起きなくなる

協調作用がない ⇔ Hill係数が1 x, yの分解 TMGout (Xout) + TMGin (x) LacY-GFP (y) LacI

Hill係数とヒステリシスの関係を調べてみよう ポジティブ・フィードバックあり Hill係数を振ってみる

演習5:Hill係数を1にしてみると ポジティブ・フィードバックあり、Hill係数=1

結論:ヒステリシスには次の2因子が必要である ポジティブ・フィードバック 協調作用 高いHill係数

なぜヒステリシスが起きるのか? ポジティブ・フィードバック、協調作用が必要な理 由を知りたい 力学系の道具を使って説明できる 相平面 ヌルクライン 固定点 (発展) ベクトル場 (発展)

道具その1:相平面 y(LacY-GFP) x(TMG) 時間変化する2変数を縦軸・横軸 にとった平面 ここでは x(TMG)-y(LacY-GFP) 平面のこと y(LacY-GFP) x(TMG)

道具その2:ヌルクライン ẋ = 0 ẏ = 0 y(LacY-GFP) x(TMG) ẋ = 0, ẏ = 0 となる点の集まりから なる曲線 ヌルクラインの交点は「固定点」と 呼ばれる ẋ = 0 ẏ = 0 y(LacY-GFP) x(TMG)

固定点は定常状態のことである ヌルクラインの交点3つ 解軌跡(右図では黒線)が 固定点に到達してそこにと どまる すなわち「固定点」 TMG, LacY-GFPともに 変化しないので定常状態 解軌跡(右図では黒線)が 固定点に到達してそこにと どまる

固定点は定常状態のことである 時系列 相平面 時系列が定常状態に達したら、それは相平面上では解軌跡が固定点に到 達したことに対応する

演習6: 相平面とヌルクライン Ozbudakモデルの解軌跡を相平面に描きなさい(MATLAB) 解軌跡 = 時間が経つにつれて点( x , y )が2次元平面上でどのように動くか ヌルクラインの式を紙と鉛筆で求め、MATLABでグラフ化しなさい

空欄を埋めなさい function ozbudak_phaseplane( ) (略。初期値、パラメータはx_out以外、演習1と同じものとする。) x_out = % ヒステリシスを起こす値と起こさない値を試すこと function dsdt = ODE(t,s,param) (略)    function plot_phase_plane( time_course , alpha , beta , rho )    x_time_course = time_course(:,1); y_time_course = time_course(:,2);    figure(1); hold on;    plot( , , 'k'); % 解の軌跡を描く    xlabel('TMG') ylabel('LacY-GFP')    y = 0:0.01:0.8;           ; % dx/dt = 0     plot( x , y , 'r' ); % ‘r’はredのr     x = 0:0.1:10;    y = x.^2 ./ ; % dy/dt = 0 plot( , , ‘b’ ); hold off; end

解答出力例 (黒線が解軌跡) ヒステリシスなし ヒステリシスあり

解答出力例 (黒線が解軌跡) ヒステリシスなし ヒステリシスあり

ポジティブ・フィードバックなし y(LacY-GFP) x0 x(TMG)

協調作用なし y(LacY-GFP) x0 x(TMG)

まとめ Ozbudakモデルはヒステリシスを再現する ヒステリシスには次の2つの要素が必要である 大腸菌は糖の濃度を記憶する ヒステリシスには次の2つの要素が必要である ポジティブ・フィードバック 協調作用 ヒステリシスが起きる理由は力学系の言葉で説明できる

発展課題

力学系の道具を使ってみよう 相平面 ヌルクライン 固定点 ベクトル場

演習6: 固定点 Ozbudakモデルに3つある固定点の座標を求めなさい 固定点 = ヌルクラインの交点

空欄を埋めなさい function S = ozbudak_fixed_point( ) syms x y; S = solve('-x + 0.1 * 5 + 10 * y = 0', ... ‘ ', 'x', 'y'); % 固定点の座標を表示 end

道具その4:ベクトル場 点 ( x , y ) に大きさ ( ẋ , ẏ ) のベクトルを 描いたもの 「( ẋ , ẏ ) のベクトル」 はその地点の風の 強さと向きに相当

ベクトル場を描くと、系内の「風」の向きと強さがわかる 初期値を相平面のどこ か一点に置く 解軌跡は矢印の方向 をなぞる 最後は固定点に到達

相平面上にベクトル場を描く 点 ( x , y ) に大きさ ( ẋ , ẏ ) のベクトルを描く 描き方の基本 ヌルクラインの性質を利用した描き方 x のヌルクライン上では ẋ = 0 すなわちベクトルは垂直 ベクトルの向きは? 例えば ẋ > 0 となる領域を求める場合、 ẋ = f( x , y ) > 0 とおいて式変形 ヌルクラインを境にどちら側の領域が ẋ > 0 となるかがわかる 点 ( x , y ) に大きさ ( ẋ , ẏ ) のベクトルを描く

ヌルクラインとベクトル場 ẋ = 0 ẋ = 0 のヌルクラインの式は y(LacY-GFP) x(TMG) まず、 ẋ > 0 となる領域を探す。 ẋ = 0 y(LacY-GFP) より、 ẋ > 0 となる領域は である。 ẋ = 0 のヌルクラインの式は であるから、 ẋ = 0 を x(TMG) 境として上側が ẋ > 0、下側が ẋ < 0 である。

ヌルクラインとベクトル場 ẋ = 0 ẋ > 0 は右向きのベクトル ẋ < 0 は左向きのベクトル y(LacY-GFP) 「ẋ = 0 を境として下側が ẋ < 0、上側 が ẋ > 0 である」 ẋ = 0 y(LacY-GFP) ẋ > 0 は右向きのベクトル ẋ < 0 は左向きのベクトル x(TMG)

MATLABでベクトル場を描く方法 相平面をメッシュ状に区切る [ x , y ] = meshgrid( 0.01:0.2:2 , 0.01:0.2:2 ); 各点における時間変化を計算 DX = -x + a*y + x.^2 .* y;   % 右辺の方程式はあくまでも一例 DY = b - a*y - x.^2 .* y; % 課題とは関係ない 点( x , y )に大きさ( DX , DY )の矢印を描く quiver( x , y , DX , DY );

演習: ベクトル場を描く(MATLAB) MATLABを使ってOzbudakモデルのベクトル場を描きなさい。

空欄を埋めなさい function ozbudak_vector_field() (略) function vector_field( alpha , beta , rho , x_out ) [x,y] = meshgrid( ); DX = DY = figure(1) quiver( ); end

正解出力例

Further Readings Steven Strogatz, “Nonlinear dynamics and chaos”, Perseus Books Publishing, 生物・化学の例が多く比較的読みやすい非線形力学系教科書 Ozbudak, E.M. et al., “Multistability in the lactose utilization network of Escherichia coli”, Nature 427(6976):737-40, 2004. 大腸菌 lac オペロンのような古典的な系で、ヒステリシスが起きることを示した Bhalla, U.S. and Iyengar, R., “Emergent properties of networks of biological signaling pathways”, Science 283(5400):381-7. シグナル伝達経路に生じる双安定性が記憶装置の基盤である可能性を指摘 Tanaka, R.J. et al. “Skin barrier homeostasis in atopic dermatitis: feedback regulation of kallikrein activity”, PLoS ONE 6(5):e19895, 2011. アトピー性皮膚炎を起こした外部刺激を取り除いても炎症がおさまらない理由として双安 定性の可能性を指摘