Team:UT-Tokyo/Data/Modeling/Model02

From 2011.igem.org

(Difference between revisions)
(Created page with "{{:Team:UT-Tokyo/Templates/BeginContent|fullpagename=Team:UT-Tokyo/Data/Modeling/Model02|subpagename=Model2: A simple model for E.coli chemotaxis}} =Aim= We derived a simple rel...")
 
(18 intermediate revisions not shown)
Line 1: Line 1:
-
{{:Team:UT-Tokyo/Templates/BeginContent|fullpagename=Team:UT-Tokyo/Data/Modeling/Model02|subpagename=Model2: A simple model for E.coli chemotaxis}}
+
{{:Team:UT-Tokyo/Templates/BeginContent|fullpagename=Team:UT-Tokyo/Data/Modeling/Model02|subpagename=Model02}}
 +
=[[Team:UT-Tokyo/Data/Modeling|Modeling]]/Model2: A simple model for ''E.coli'' chemotaxis=
-
=Aim=
+
=[[Team:UT-Tokyo/Data/Modeling/Model02/applet|Interactive demo]]=
-
We derived a simple relationship between Asp concentration and E.coli chemotaxis so that
+
-
we were able to simulate macroscopic ¥(n~10^8) E.coli colony chemotacitc behavior.
+
-
=Background=
+
==Aim==
-
We considered internal condition of E.coli (receptor methylation etc) by SPECS model<html><sup class="ref">[1]</sup></html>.
+
We tried to derive a simple relationship between Asp concentration and ''E.coli'' chemotactic behavior so that
-
SPECS is a relatively accurate taxis model, and is suitable for simulation of approximately ten thousand E.coli for approximately 10 minutes.
+
it become possible to simulate macroscopic (n~10<html><sup>8</sup></html>) ''E.coli''-colony chemotactic behavior in the entire simulation.
-
This model, however, take too much time to reproduce our system, because the system includes much more E.coli(about one hundred of million) and takes a relatively long time (a few days).
+
-
So, by adapting SPECS model to our system, we have devised a new model that provides less calculation amount and enough accuracy.
+
-
Concretely, we divided chemotactic motion of E.coli groups into two factors, that is, parallel translation and diffusion. Then we derived approximation equations that figure out the two respective values from Asp concentration and its gradient.
+
-
=Method=
+
==Background==
-
In this study, we used SPECS method to simulate E.coli chemotactic motion.
+
SPECS <html><sup class="ref">[1]</sup></html> is one of the models simulating ''E.coli'' chemotactic behavior, and it considers internal condition of ''E.coli'' (receptor methylation etc).  It is suitable for semi-macroscopic simulation which has around ten thousand ''E.coli'' and the timescale is within a hour.
-
Here we introduce the brief overview of SPECS.
+
For our system, however, it takes too much time to replicate entire system, because the system includes much more ''E.coli''(about one hundred of million) and the time scale is relatively long (a few days).
 +
Therefore we devised a new model based on SPECS which requires less calculation amount.
-
We processed Monte Carlo simulation for 20000 E.colis which behavior was based on SPECS in some different Asp concentration and its gradient.  The gradient we used were exponential ¥[eqn 1].
+
==Method==
-
We approximated E.colis movement as the combination of parallel translation and diffusion.
+
We approximated the motion of ''E.coli'' groups doing a process of chemotaxis into two factors, that is, parallel translation and diffusion.
-
We recorded the time course of average position <x> and root-mean-square <<html>x<sup>2</sup></html>> and then
+
-
calculated average velocity ¥[eqn 2] and the diffusion constant ¥[eqn 3].
+
-
For these results we derived the relationship between:
+
[[File:utt_m2_figure.PNG]]
-
* Asp concentration and average velocity of E.coli
+
-
* Asp concentration and diffusion constant of E.coli
+
-
* gradient of Asp concentration and average velocity of E.coli
+
-
* gradient of Asp concentration and diffusion constant of E.coli
+
-
Finally we fitted approximated curve to these data.
+
 +
Let 'V' represent the parallel translation velocity, 'D' represent the diffusion coefficient of ''E.coli'' groups and 'A' represent L-Asp concentration.
 +
''E.coli'' can detect A and the gradient ∇A.
 +
They have an internal circuit to translate these signals into motions of their flagella and try to move towards regions with more Asp.
 +
So we thought V and D could be represented as a function of A and ∇A.
 +
And then we used the SPECS model to derive these functions.
-
=Result=
+
Here we introduce SPECS briefly.
-
¥[fig1. relationship between Asp conc. and avg. velocity]
+
An ''E.coli'' belongs to one of two mobile states: "Run" and "Tumble".
-
¥[fig2. relationship between gradient of Asp conc. and avg. velocity]
+
Run state indicates that ''E.coli'' goes straight until the state changes.
-
¥[fig3. relationship between Asp conc. and diffusion coef.]
+
Tumble state indicates that ''E.coli'' is changing its direction and doesn't move around.
-
¥[fig4. relationship between gradient of Asp conc. and diffusion coef.]
+
''E.coli'' changes its state from "Run" to "Tumble" with probability P<html><sub>rt</sub> </html> [1/sec].
-
fitted curveは次の式で与えられる。
+
From "Tumble" to "Run" with P<html><sub>tr</sub></html> [1/sec].
-
¥[eqn4](tanh(0.66*(log(A)+11.35))+1.0)*0.055*grad(log(A))
+
These probabilities depend on the ''E.coli'' internal state.
-
¥[eqn5](tanh(0.65*(log(A)+11.5))+1.0)*(tanh(0.71*(grad(log(A))-2.05)+0.9)*0.0037+0.0047
+
The internal state is represented by two variables: "a" and "m".
-
¥φ means Asp concentration.
+
"a" indicates kinase activity and "m" represents the receptor methylation level of the cell.
 +
These vary depending on L-Asp concentration.
 +
In addition, SPECS takes into account Brownian fluctuation.
 +
Using this model, SPECS tries to replicate the internal circuit of ''E.coli'' which translates an Asp concentration distribution into a mobile state.
 +
 
 +
We processed a Monte Carlo simulation for 20000 ''E.coli'' behaving according to SPECS when placed in fields containing different concentrations and gradients of Asp.  The shape of the Asp gradient used was exponential according to the following relationship:
 +
 
 +
[[File:utt_m2_eqn1.png]]
 +
 
 +
where A0 represents the concentration of Asp at x = 0 and x0 a parameter that represents the slope of the gradient.
 +
 
 +
We approximated ''E.coli'' movement by a combination of parallel translation and diffusion.
 +
We recorded the average position <x> and root-mean-square <x<html><sup>2</sup></html>> over time and then
 +
calculated the average velocity:
 +
 
 +
[[File:utt_m2_eqn2.png]]
 +
 
 +
and the diffusion constant:
 +
 
 +
[[File:utt_m2_eqn3.png]]
 +
 
 +
We derived the relationship between
 +
* Asp concentration and average velocity of ''E.coli''
 +
* Asp concentration and diffusion constant of ''E.coli''
 +
* gradient of Asp concentration and average velocity of ''E.coli''
 +
* gradient of Asp concentration and diffusion constant of ''E.coli''
 +
Finally we fitted approximated curves (represented as functions of A and ∇A) to these data so that we could derive functions we needed in the entire system simulation.
 +
 
 +
==Result==
 +
 
 +
{{:Team:UT-Tokyo/Templates/Image|file=utt_m2_fig1.png|caption=Figure 1. relationship between Asp conc. and avg. velocity}}
 +
{{:Team:UT-Tokyo/Templates/Image|file=utt_m2_fig2.png|caption=Figure 2. relationship between gradient of Asp conc. and avg. velocity}}
 +
{{:Team:UT-Tokyo/Templates/Image|file=utt_m2_fig3.png|caption=Figure 3. relationship between Asp conc. and diffusion coef}}
 +
{{:Team:UT-Tokyo/Templates/Image|file=utt_m2_fig4.png|caption=Figure 4. relationship between gradient of Asp conc. and diffusion coef}}
 +
 
 +
The fitted curves are represented by the following equations.
 +
 
 +
[[File:utt_m2_eqn4.png]]
 +
 
 +
[[File:utt_m2_eqn5.png]]
<html>
<html>
Line 43: Line 77:
<div id="references">
<div id="references">
<ul>
<ul>
-
   <li id="ref_1">[1] Jiang L, Ouyang Q, Tu Y (2010) Quantitative Modeling of Escherichia coli Chemotactic Motion in Environments Varying in Space and Time. PLoS Comput Biol 6(4): e1000735. doi:10.1371/journal.pcbi.1000735</li>
+
   <li id="ref_1">[1] Jiang L, Ouyang Q, Tu Y (2010) Quantitative Modeling of "Escherichia coli" Chemotactic Motion in Environments Varying in Space and Time. PLoS Comput Biol 6(4): e1000735. doi:10.1371/journal.pcbi.1000735</li>
</ul>
</ul>
</div>
</div>
</html>
</html>
{{:Team:UT-Tokyo/Templates/EndContent}}
{{:Team:UT-Tokyo/Templates/EndContent}}

Latest revision as of 02:08, 6 October 2011

Modeling/Model2: A simple model for E.coli chemotaxis

Interactive demo

Aim

We tried to derive a simple relationship between Asp concentration and E.coli chemotactic behavior so that it become possible to simulate macroscopic (n~108) E.coli-colony chemotactic behavior in the entire simulation.

Background

SPECS [1] is one of the models simulating E.coli chemotactic behavior, and it considers internal condition of E.coli (receptor methylation etc). It is suitable for semi-macroscopic simulation which has around ten thousand E.coli and the timescale is within a hour. For our system, however, it takes too much time to replicate entire system, because the system includes much more E.coli(about one hundred of million) and the time scale is relatively long (a few days). Therefore we devised a new model based on SPECS which requires less calculation amount.

Method

We approximated the motion of E.coli groups doing a process of chemotaxis into two factors, that is, parallel translation and diffusion.

Utt m2 figure.PNG

Let 'V' represent the parallel translation velocity, 'D' represent the diffusion coefficient of E.coli groups and 'A' represent L-Asp concentration. E.coli can detect A and the gradient ∇A. They have an internal circuit to translate these signals into motions of their flagella and try to move towards regions with more Asp. So we thought V and D could be represented as a function of A and ∇A. And then we used the SPECS model to derive these functions.

Here we introduce SPECS briefly. An E.coli belongs to one of two mobile states: "Run" and "Tumble". Run state indicates that E.coli goes straight until the state changes. Tumble state indicates that E.coli is changing its direction and doesn't move around. E.coli changes its state from "Run" to "Tumble" with probability Prt [1/sec]. From "Tumble" to "Run" with Ptr [1/sec]. These probabilities depend on the E.coli internal state. The internal state is represented by two variables: "a" and "m". "a" indicates kinase activity and "m" represents the receptor methylation level of the cell. These vary depending on L-Asp concentration. In addition, SPECS takes into account Brownian fluctuation. Using this model, SPECS tries to replicate the internal circuit of E.coli which translates an Asp concentration distribution into a mobile state.

We processed a Monte Carlo simulation for 20000 E.coli behaving according to SPECS when placed in fields containing different concentrations and gradients of Asp. The shape of the Asp gradient used was exponential according to the following relationship:

Utt m2 eqn1.png

where A0 represents the concentration of Asp at x = 0 and x0 a parameter that represents the slope of the gradient.

We approximated E.coli movement by a combination of parallel translation and diffusion. We recorded the average position <x> and root-mean-square <x2> over time and then calculated the average velocity:

Utt m2 eqn2.png

and the diffusion constant:

Utt m2 eqn3.png

We derived the relationship between

  • Asp concentration and average velocity of E.coli
  • Asp concentration and diffusion constant of E.coli
  • gradient of Asp concentration and average velocity of E.coli
  • gradient of Asp concentration and diffusion constant of E.coli

Finally we fitted approximated curves (represented as functions of A and ∇A) to these data so that we could derive functions we needed in the entire system simulation.

Result

Figure 1. relationship between Asp conc. and avg. velocity
Figure 1. relationship between Asp conc. and avg. velocity

Figure 2. relationship between gradient of Asp conc. and avg. velocity
Figure 2. relationship between gradient of Asp conc. and avg. velocity

Figure 3. relationship between Asp conc. and diffusion coef
Figure 3. relationship between Asp conc. and diffusion coef

Figure 4. relationship between gradient of Asp conc. and diffusion coef
Figure 4. relationship between gradient of Asp conc. and diffusion coef

The fitted curves are represented by the following equations.

Utt m2 eqn4.png

Utt m2 eqn5.png

References

  • [1] Jiang L, Ouyang Q, Tu Y (2010) Quantitative Modeling of "Escherichia coli" Chemotactic Motion in Environments Varying in Space and Time. PLoS Comput Biol 6(4): e1000735. doi:10.1371/journal.pcbi.1000735