Tokyo/Formulation/1.toggle model
From 2007.igem.org
Line 2: | Line 2: | ||
First,the ordinary differential equations(ODEs) of the toggle switch were derived as | First,the ordinary differential equations(ODEs) of the toggle switch were derived as | ||
- | <br>[[Image:expression1-1.jpg|200px|]] [[Image:parameter1-1.jpg|200px|]] | + | <br>[[Image:expression1-1.jpg|200px|left|thumb|Ex1-1]] [[Image:parameter1-1.jpg|200px|]] |
<br>These equations were normalized as follows: | <br>These equations were normalized as follows: | ||
- | <br>[[Image:expression1-2.jpg|200px|]] | + | <br>[[Image:expression1-2.jpg|200px|none|thumb|Ex1-2]] |
<br>In the steady state,time derivatives are zero: | <br>In the steady state,time derivatives are zero: | ||
- | <br>[[Image: | + | <br>[[Image:expression1-3.jpg|80px|none|thumb|Ex1-3]] |
<br>As a result,the nullclines of this system were derived as | <br>As a result,the nullclines of this system were derived as | ||
- | <br>[[Image:Siki2.jpg|200px|]] | + | <br>[[Image:Siki2.jpg|200px|none|thumb|Ex1-4]] |
- | <br>which indicate the nullclines of the system shown in | + | <br>which indicate the nullclines of the system shown in Fig 1.A-C.c.Where about parameters,we use three sets of parameters. |
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | <br> | + | <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|Tble1.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.A]] | |
+ | [[Image:toggle2.jpg|270px|left|thumb|Figure 1.B]] | ||
+ | [[Image:Toggle1-4.jpg|270px|none|thumb|Figure 1.C]] | ||
+ | <br>First,we carried out kinetic simulations in the condition of Fig1.A. The results are shown in Fig2.A-C. | ||
- | + | [[Image:toggle3.jpg|200px|left|thumb|Figure 2.A (Ra(0),Rb(0))=(0.0,2.5)]] [[Image:toggle4.JPG|200px|left|thumb|Figure 2.B (Ra(0),Rb(0))=(2.5,0.0)]] | |
+ | [[Image:toggle5.JPG|200px|left|thumb|Figure 2.C (Ra(0),Rb(0))=(1.5,1.3)]] | ||
+ | [[Image:toggle1-1.jpg|200px|none|thumb|Figure 3 two stable]] | ||
- | |||
- | <br> | + | <br>Fig 2.A-C indicate that when the initial condition is (Ra,Rb)=(0.0,2.5),which is near stable equilibrium point B, the values of Ra and Rb goes to stable equilibrium point B and when the initial condition is (Ra,Rb)=(2.5,0.0),which is near stable equilibrium point A, the values of Ra and Rb goes to stable equilibrium point B. |
- | <br> | + | <br>Next,when the number of stable equilibrium point is one, the result of simulation are shown in Fig 1.B. |
<br>Fig● indicate that when the initial condition is (Ra,Rb)=(0.0,2.5) and when the initial condition is (Ra,Rb)=(2.5,0.0), the values of Ra and Rb goes to stable equilibrium point B. | <br>Fig● indicate that when the initial condition is (Ra,Rb)=(0.0,2.5) and when the initial condition is (Ra,Rb)=(2.5,0.0), the values of Ra and Rb goes to stable equilibrium point B. | ||
<br>'''As a result,taking two stable status need the phaseplane of two stable equilibrium points and we have to set proper parameters.''' | <br>'''As a result,taking two stable status need the phaseplane of two stable equilibrium points and we have to set proper parameters.''' |
Revision as of 04:36, 24 October 2007
1.toggle model
First,the ordinary differential equations(ODEs) of the toggle switch were derived as
These equations were normalized as follows:
In the steady state,time derivatives are zero:
As a result,the nullclines of this system were derived as
which indicate the nullclines of the system shown in Fig 1.A-C.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.
B)the maximum expression rate of repressor A and repressor B is equal,and hill coefficient of A is one.
C)the maximum expression rate of repressor A and repressor B is not balanced,and hill coefficient of both A and B is three.
First,we carried out kinetic simulations in the condition of Fig1.A. The results are shown in Fig2.A-C.
Fig 2.A-C indicate that when the initial condition is (Ra,Rb)=(0.0,2.5),which is near stable equilibrium point B, the values of Ra and Rb goes to stable equilibrium point B and when the initial condition is (Ra,Rb)=(2.5,0.0),which is near stable equilibrium point A, the values of Ra and Rb goes to stable equilibrium point B.
Next,when the number of stable equilibrium point is one, the result of simulation are shown in Fig 1.B.
Fig● indicate that when the initial condition is (Ra,Rb)=(0.0,2.5) and when the initial condition is (Ra,Rb)=(2.5,0.0), the values of Ra and Rb goes to stable equilibrium point B.
As a result,taking two stable status need the phaseplane of two stable equilibrium points and we have to set proper parameters.