Tokyo/Formulation/1.toggle model

From 2007.igem.org

< Tokyo/Formulation(Difference between revisions)
 
(38 intermediate revisions not shown)
Line 1: Line 1:
-
<中身>
+
__NOTOC__
-
(相平面解析)
+
<br>[[Tokyo/Works|Works top]]&nbsp;&nbsp;&nbsp;0.[[Tokyo/Works/Hybrid promoter|Hybrid promoter]]&nbsp;&nbsp;&nbsp;'''1.[[Tokyo/Works/Formulation |Formulation]]'''&nbsp;&nbsp;&nbsp;2.[[Tokyo/Works/Assay |Assay1]]&nbsp;&nbsp;&nbsp;3.[[Tokyo/Works/Simulation |Simulation]]&nbsp;&nbsp;&nbsp;4.[[Tokyo/Works/Assay2 |Assay2]]&nbsp;&nbsp;&nbsp;5.[[Tokyo/Works/Future works |Future works]]
-
パラメータによって相平面が変わり,平衡点が1つのときと3つのときがある.
+
<br><br>[[Tokyo/Formulation/1.toggle model |Step1]]&nbsp;&nbsp;&nbsp;[[Tokyo/Formulation/2.toggle model with hybrid promoter |Step2]]&nbsp;&nbsp;&nbsp;[[Tokyo/Formulation/3.AHL-experssing model|Step3]]  
-
平衡点が3つのときbistableになり,A状態,B状態ができる.
+
<br>
-
(シミュレーション)
+
== 1.Single cell model:mutual inhibition ==
-
初期値を安定点付近から始めるとちゃんとそのまま安定してるよ
+
First, the ordinary differential equations (ODEs) of the toggle switch were derived as
 +
<br>[[Image:expression1-1.jpg|200px|left|thumb|Ex1-1 ]]  [[Image:parameter1-1.jpg|200px|]]
-
<br>[[Image:式-1.JPG|200px|]]
+
<br>These equations were normalized as follows:
-
<br>定常状態では時間変化が0.つまり,右辺=0を解けばよい.すると,
+
<br>[[Image:expression1-2.jpg|200px|none|thumb|Ex1-2 ]]
-
<br>[[Image:Siki2.jpg|200px|]]
+
<br>In the steady state, time derivatives are zero:
-
<br>となり,グラフにすると下のようになる.この線のことを一般的にnullclineと言い,nullclineの交点が
+
<br>[[Image:expression1-3.jpg|80px|none|thumb|Ex1-3]]
-
平衡点となる.
+
-
<br>[[Image:toggle1.jpg|200px|]] [[Image:toggle2.jpg|200px|]]
+
<br>As a result, the nullclines of this system were derived as
-
<br>平衡点には安定な平衡点と不安定な平衡点があり,今回は黒丸が安定な平衡点で白丸が不安定な平衡点である.
+
<br>[[Image:Siki2.jpg|200px|none|thumb|Ex1-4]]
-
<br>実際にシミュレーションした結果と比較すると明らかである.
+
-
<br>左端のグラフが左上の安定点付近を初期値として計算した結果.真ん中のグラフが中心の不安定点付近を初期値として計算した結果.右端のグラフが,右下の安定点付近を初期値として計算した結果.
+
-
<br>[[Image:toggle3.jpg|200px|]] [[Image:toggle4.JPG|200px|]] [[Image:toggle5.JPG|200px|]]
+
<br>which indicate the nullclines of the system shown in Fig 1.1.A-C. Where about parameters, we use three sets of parameters.
 +
 
 +
<br>  A)the maximum expression rate of repressor A and repressor B is balanced,and hill coefficient of both A and B is three.
 +
<br>[[Image:parameter1-2.jpg|150px|center|thumb|Table1.A]]
 +
<br>  B)the maximum expression rate of repressor A and repressor B is equal,and hill coefficient of A is one.
 +
<br>[[Image:parameter1-3.JPG|150px|center|thumb|Table1.B]]
 +
<br>  C)the maximum expression rate of repressor A and repressor B is not balanced,and hill coefficient of both A and B is three.
 +
<br>[[Image:parameter1-4.JPG|150px|center|thumb|Table1.C]]
 +
 
 +
[[Image:toggle1.jpg|260px|left|thumb|Figure 1.1.A]]
 +
[[Image:toggle2.jpg|270px|left|thumb|Figure 1.1.B]]
 +
[[Image:Toggle1-4.jpg|270px|none|thumb|Figure 1.1.C]]
 +
<br>First,we carried out kinetic simulations in the condition of Fig 1.1.A. The results are shown in Fig 1.2.A-C.
 +
 
 +
[[Image:toggle3.jpg|200px|left|thumb|Figure 1.2.A  (Ra(0),Rb(0))=(0.0,2.5)]] [[Image:toggle4.JPG|200px|left|thumb|Figure 1.2.B  (Ra(0),Rb(0))=(2.5,0.0)]]
 +
[[Image:toggle5.JPG|200px|left|thumb|Figure 1.2.C  (Ra(0),Rb(0))=(1.5,1.3)]]
 +
[[Image:toggle1-1.jpg|200px|none|thumb|Figure 3  bistable]]
 +
 
 +
 
 +
<br>Fig 1.2.A-C indicate that when the initial condition is (Ra,Rb)=(0.0,2.5), which is near the stable equilibrium point B, the values of Ra and Rb go to stable equilibrium point B, and when the initial condition is (Ra,Rb)=(2.5,0.0), which is near the stable equilibrium point A, the values of Ra and Rb go to the stable equilibrium point B.
 +
 
 +
<br>Next, when the number of stable equilibrium point is one(Fig 1.1.B), the result of simulation are shown in Fig 1.4.A-C.
 +
 
 +
[[Image:toggle6.JPG|200px|left|thumb|Figure 1.4.A  (Ra(0),Rb(0))=(0.0,2.5)]] [[Image:toggle7.JPG|200px|left|thumb|Figure 1.4.B  (Ra(0),Rb(0))=(2.5,0.0)]]
 +
[[Image:toggle8.JPG|200px|left|thumb|Figure 1.4.C  (Ra(0),Rb(0))=(1.5,1.3)]]
 +
[[Image:toggle1-2.jpg|220px|none|thumb|Figure 1.5  monostable]]
 +
 
 +
<br><br>
 +
<br>
 +
<br>Fig 4.A-C indicate that the value of Ra and Rb go to the stable equilibrium point A regardless of an initial value in case of monostable state.
 +
 
 +
<br>'''As a result, taking two stable states needs the phaseplane of two stable equilibrium points and Hill coefficients was very important.'''
 +
 
 +
== ==
 +
[[Tokyo/Formulation/1.toggle model|Step.1]] >> [[Tokyo/Formulation/2.toggle model with hybrid promoter|Step.2]]

Latest revision as of 02:45, 27 October 2007


Works top   0.Hybrid promoter   1.Formulation   2.Assay1   3.Simulation   4.Assay2   5.Future works

Step1   Step2   Step3  

1.Single cell model:mutual inhibition

First, the ordinary differential equations (ODEs) of the toggle switch were derived as


Ex1-1
Parameter1-1.jpg


These equations were normalized as follows:


Ex1-2


In the steady state, time derivatives are zero:


Ex1-3


As a result, the nullclines of this system were derived as


Ex1-4


which indicate the nullclines of the system shown in Fig 1.1.A-C. Where about parameters, we use three sets of parameters.


  A)the maximum expression rate of repressor A and repressor B is balanced,and hill coefficient of both A and B is three.


Table1.A


  B)the maximum expression rate of repressor A and repressor B is equal,and hill coefficient of A is one.


Table1.B


  C)the maximum expression rate of repressor A and repressor B is not balanced,and hill coefficient of both A and B is three.


Table1.C
Figure 1.1.A
Figure 1.1.B
Figure 1.1.C


First,we carried out kinetic simulations in the condition of Fig 1.1.A. The results are shown in Fig 1.2.A-C.

Figure 1.2.A (Ra(0),Rb(0))=(0.0,2.5)
Figure 1.2.B (Ra(0),Rb(0))=(2.5,0.0)
Figure 1.2.C (Ra(0),Rb(0))=(1.5,1.3)
Figure 3 bistable



Fig 1.2.A-C indicate that when the initial condition is (Ra,Rb)=(0.0,2.5), which is near the stable equilibrium point B, the values of Ra and Rb go to stable equilibrium point B, and when the initial condition is (Ra,Rb)=(2.5,0.0), which is near the stable equilibrium point A, the values of Ra and Rb go to the stable equilibrium point B.


Next, when the number of stable equilibrium point is one(Fig 1.1.B), the result of simulation are shown in Fig 1.4.A-C.

Figure 1.4.A (Ra(0),Rb(0))=(0.0,2.5)
Figure 1.4.B (Ra(0),Rb(0))=(2.5,0.0)
Figure 1.4.C (Ra(0),Rb(0))=(1.5,1.3)
Figure 1.5 monostable





Fig 4.A-C indicate that the value of Ra and Rb go to the stable equilibrium point A regardless of an initial value in case of monostable state.


As a result, taking two stable states needs the phaseplane of two stable equilibrium points and Hill coefficients was very important.

Step.1 >> Step.2