Team:BNU-China/modeling.html

From 2014.igem.org

(Difference between revisions)
 
(25 intermediate revisions not shown)
Line 27: Line 27:
<div class="clear"></div>
<div class="clear"></div>
-
<br/><br/>
+
<div class="hr2"></div>
<h2>Abstract</h2>
<h2>Abstract</h2>
-
<p>Analyzing dynamics of E.coli’s chemotaxis and estimating effectiveness of the "Prometheus" require high accuracy, and thus hard to realize with experiments. Here, we construct a virtual peanut root-E.coli system (by c++ and OpenGL), storing the 3-D environmental information by k-d tree. We quantify the process in which “Promethei”, our engineered bacteria, carrying Mo, move towards peanut roots. We at last estimate around 20-40% of bacteria could work in reality, and that the best “memory” time be around 40 min (See here). All codes (in c++) of the project are available*.</p>
+
<p>Analyzing dynamics of <i>E.coli</i>’s chemotaxis and estimating effectiveness of the "Prometheus" require high accuracy, and thus hard to realize with experiments. Here, we construct a virtual peanut root-<i>E.coli</i> system (by c++ and OpenGL), storing the 3-D environmental information by k-d tree. We quantify the process in which “Promethei”, our engineered bacteria, carrying Mo, move towards peanut roots. We at last estimate around 20-40% of bacteria could work in reality, and that the best “memory” time be around 40 min (<a href="#res">See here</a>). All codes (in c++) of the project are available*.</p>
-
<p class="fig">*The size of the code is too big, please contract Li Si-yao (201211211043@mail.bnu.edu.cn) to get it.</p>
+
<p class="fig">*The size of the code is too big, please contact Li Si-yao (201211211043@mail.bnu.edu.cn) to get it.</p>
<br/>
<br/>
-
<a title="Fig. 1 The Motion Graph of the screenshots while running the Simulation Program. The purple points are E.coli tending towards the root surface." href="https://static.igem.org/mediawiki/2014/e/e9/Bnu_root_chemotacix_gif.gif" rel="prettyPhoto"> <span class="overlay zoom" style="display: none;"></span><img style="margin-left:200px; width:50%" src="https://static.igem.org/mediawiki/2014/e/e9/Bnu_root_chemotacix_gif.gif" /></a>
+
<a title="Fig. 1 The Motion Graph of the screenshots while running the Simulation Program. The purple points are <i>E.coli</i> tending towards the root surface." href="https://static.igem.org/mediawiki/2014/e/e9/Bnu_root_chemotacix_gif.gif" rel="prettyPhoto"> <span class="overlay zoom" style="display: none;"></span><img style="margin-left:200px; width:50%" src="https://static.igem.org/mediawiki/2014/e/e9/Bnu_root_chemotacix_gif.gif" /></a>
<br/>
<br/>
-
<p class="fig" style="margin-left:30px; width:86%">Fig. 1 The Motion Graph of the screenshots while running the Simulation Program. The purple points are E.coli tending towards the root surface.</p>
+
<p class="fig" style="margin-left:30px; width:86%">Fig. 1 The Motion Graph of the screenshots while running the Simulation Program. The purple points are <i>E.coli</i> tending towards the root surface.</p>
-
<hr>
+
 
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
 +
<div class="hr2"></div>
<h2 id="eco"><i>E.coli</i> model</h2>
<h2 id="eco"><i>E.coli</i> model</h2>
<h3>Mechanism</h3>
<h3>Mechanism</h3>
-
<p>We divide Escherichia coli’s behavior into two parts: “running” straight lines, and “tumbling” for reorientation[1]. When attractant concentration uniformed, E.coli can be described as doing “random walk”[4], and the running time fits a normal distribution. With spatial concentration gradients of chemical attractance,when bacteria sense a higher attractant concentration than before,the mean of the statistical distribution of running time increases. Otherwise, the time keeps the same as that in the uniform environment, no matter how high or low the concentration is[4].</p>
+
<p>We divide <i>Escherichia coli</i>’s behavior into two parts: “running” straight lines, and “tumbling” for reorientation[1]. When attractant concentration uniformed, <i>E.coli</i> can be described as doing “random walk”[4], and the running time fits a normal distribution. With spatial concentration gradients of chemical attractance,when bacteria sense a higher attractant concentration than before,the mean of the statistical distribution of running time increases. Otherwise, the time keeps the same as that in the uniform environment, no matter how high or low the concentration is[4].</p>
<div class="clear"></div>
<div class="clear"></div>
<br/>
<br/>
<h4>Running</h4>
<h4>Running</h4>
-
<p>The running of E.Coli can be descibed as a straight line. The distribution of running duration time fits a normal curve, with mean <strong>τ<sub>run0</sub></strong> and standard division <strong>σ<sub>run</sub></strong>.</p>
+
<p>The running of <i>E.coli</i> can be descibed as a straight line. The distribution of running duration time fits a normal curve, with mean <strong>τ<sub>run0</sub></strong> and standard division <strong>σ<sub>run</sub></strong>.</p>
<p>Brown and Berg[4] suggested that the mean <strong>τ</strong> has a functional relationship, shown as follows, with the current attractant (in their case, glutamate) concentration and the change of concentration over time. (See formula (1),(2)).
<p>Brown and Berg[4] suggested that the mean <strong>τ</strong> has a functional relationship, shown as follows, with the current attractant (in their case, glutamate) concentration and the change of concentration over time. (See formula (1),(2)).
Line 65: Line 66:
<a title="Fig. 3 the result from experiment." href="https://static.igem.org/mediawiki/2014/4/43/Bnu_modeling2.png" rel="prettyPhoto"> <span class="overlay zoom" style="display: none;"></span><img class="left" width=48% src="https://static.igem.org/mediawiki/2014/4/43/Bnu_modeling2.png"></a>
<a title="Fig. 3 the result from experiment." href="https://static.igem.org/mediawiki/2014/4/43/Bnu_modeling2.png" rel="prettyPhoto"> <span class="overlay zoom" style="display: none;"></span><img class="left" width=48% src="https://static.igem.org/mediawiki/2014/4/43/Bnu_modeling2.png"></a>
-
<a title="Fig. 4 Fitting the experiment by <strong>a</strong>( <strong>K<sub>D</sub>C</strong> / (<strong>K<sub>D</sub></strong> + <strong>C</strong>)<sup>2</sup>) + <strong>d</strong>, where <strong>a</strong> = 389, <strong>d</strong> = 260 and <strong>K<sub>D</sub></strong> is 3.5mM. This work is made by c++ and OpenGL. The code to get <strong>K<sub>D</sub></strong> and <strong>a</strong> and to make sketch is here." href="https://static.igem.org/mediawiki/2014/3/3f/Bnu_modeling3.png" rel="prettyPhoto"> <span class="overlay zoom" style="display: none;"></span><img class="right" width=48% src="https://static.igem.org/mediawiki/2014/3/3f/Bnu_modeling3.png"></a>
+
<a title="Fig. 4 Fitting the experiment by <strong>a</strong>( <strong>K<sub>D</sub>C</strong> / (<strong>K<sub>D</sub></strong> + <strong>C</strong>)<sup>2</sup>) + <strong>d</strong>, where <strong>a</strong> = 389, <strong>d</strong> = 260 and <strong>K<sub>D</sub></strong> is 3.5mM. This work is made by c++ and OpenGL. The code to get <strong>K<sub>D</sub></strong> and <strong>a</strong> and to make sketch is here." href="https://static.igem.org/mediawiki/2014/3/3f/Bnu_modeling3.png" rel="prettyPhoto"> <span class="overlay zoom" style="display: none;"></span><img class="right" width=48% src="https://static.igem.org/mediawiki/2014/3/3f/Bnu_modeling3.png"></a>
<p class="fig" style="margin-left:100px; width:180px; margin-top:260px">Fig. 3 the result from experiment.</p>
<p class="fig" style="margin-left:100px; width:180px; margin-top:260px">Fig. 3 the result from experiment.</p>
-
<p class="fig" style="width:305px; margin-left:500px; margin-top:-40px">Fig. 4 Fitting the experiment by <strong>a</strong>( <strong>K<sub>D</sub>C</strong> / (<strong>K<sub>D</sub></strong> + <strong>C</strong>)<sup>2</sup>) + <strong>d</strong>, where <strong>a</strong> = 389, <strong>d</strong> = 260 and <strong>K<sub>D</sub></strong> is 3.5mM. This work is made by c++ and OpenGL. The code to get <strong>K<sub>D</sub></strong> and <strong>a</strong> and to make sketch is here.</p>
+
<p class="fig" style="width:305px; margin-left:500px; margin-top:-40px">Fig. 4 Fitting the experiment by <strong>a</strong>( <strong>K<sub>D</sub>C</strong> / (<strong>K<sub>D</sub></strong> + <strong>C</strong>)<sup>2</sup>) + <strong>d</strong>, where <strong>a</strong> = 389, <strong>d</strong> = 260 and <strong>K<sub>D</sub></strong> is 3.5mM. This work is <br/>made by c++ and OpenGL. The code to get <strong>K<sub>D</sub></strong> and <strong>a</strong> and to make sketch is available. Please contact us.</p>
<div class="clear"></div>
<div class="clear"></div>
Line 76: Line 77:
<p>In conclusion, modeling running is based on following statement:</p>
<p>In conclusion, modeling running is based on following statement:</p>
<ul>
<ul>
-
<li>The velocity during a run is constant (always equal to the mean)*
+
<li>The velocity during a run is constant (always equal to the mean);
-
<li>E.coli move straightly during a run, and neither die nor divide in the process. The effect of quorum sensing and nutritional factors are eliminated.
+
<li><i>E.coli</i> move straightly during a run, and neither die nor divide in the process. The effect of quorum sensing and nutritional factors are eliminated;
-
<li>The E.coli’s running time t<sub>run</sub> is a random number from N(<τ<sub>run</sub>>,σ<sub>run</sub>)
+
<li>The <i>E.coli</i>’s running time t<sub>run</sub> is a random number from N(<τ<sub>run</sub>>,σ<sub>run</sub>);
-
<li>The mean duration time <τ<sub>run</sub>> can be derived from formula (1) and (2), given the value of C and dC/dt.
+
<li>The mean duration time <τ<sub>run</sub>> can be derived from formula (1) and (2), given the value of C and dC/dt;
<li>The specific duration_time of a single bacterium changes in the same pace with the changes of the statistical distribution’s mean (relying on C and dC/dt). (Fig. 5)
<li>The specific duration_time of a single bacterium changes in the same pace with the changes of the statistical distribution’s mean (relying on C and dC/dt). (Fig. 5)
</ul>
</ul>
Line 91: Line 92:
<h4>Tumbling</h4>
<h4>Tumbling</h4>
<p><i>E. coli</i> reorient between two runs.In modeling this process, we assume:</p>
<p><i>E. coli</i> reorient between two runs.In modeling this process, we assume:</p>
-
<p>E. coli do not change positions as they are tumbling;</p>
+
<ul>
-
<p>the direction after reorientation is randomly chosen.</p>
+
<li><i>E. coli</i> do not change positions as they are tumbling;
 +
<li>the direction after reorientation is randomly chosen.
 +
</ul>
<a title="Fig. 6 Schematic drawing of a tumble, from [2], page 45" href="https://static.igem.org/mediawiki/2014/9/93/Bnu_modeling4.png" rel="prettyPhoto"> <span class="overlay zoom" style="display: none;"></span><img class="center" style="width:50%" src="https://static.igem.org/mediawiki/2014/9/93/Bnu_modeling4.png"></a>
<a title="Fig. 6 Schematic drawing of a tumble, from [2], page 45" href="https://static.igem.org/mediawiki/2014/9/93/Bnu_modeling4.png" rel="prettyPhoto"> <span class="overlay zoom" style="display: none;"></span><img class="center" style="width:50%" src="https://static.igem.org/mediawiki/2014/9/93/Bnu_modeling4.png"></a>
<p align="center" class="fig" style="width:320px; margin-left:250px">Fig. 6 Schematic drawing of a tumble, from [2], page 45</p>
<p align="center" class="fig" style="width:320px; margin-left:250px">Fig. 6 Schematic drawing of a tumble, from [2], page 45</p>
Line 102: Line 105:
<br/>
<br/>
-
<p class="fig">*The statement is made to calculate the process more easily. We suppose it is the <sub>run</sub> duration time, not the difference of velocity, is the main variant of the whole chemotaxis system.</p>
+
<p class="fig">*The statement is made to calculate the process more easily. We suppose it is the run duration time, not the difference of velocity, is the main variant of the whole chemotaxis system.</p>
<h3>Coefficients</h3>
<h3>Coefficients</h3>
-
<p>Velocity = 100μm</p>
+
<p>velocity = 100μm</p>
-
<p><sub>run</sub> duration time mean in uniform environment = 1.3</p>
+
<p>run duration time mean in uniform environment = 1.3</p>
-
<p><sub>run</sub> duration standard division = 3</p>
+
<p>run duration standard division = 3</p>
<p> mean tumble duration time = 0.14 sec</p>
<p> mean tumble duration time = 0.14 sec</p>
<p>alpha (See equation (1)) = 660 sec</p>
<p>alpha (See equation (1)) = 660 sec</p>
-
<p>KD = 3500 μmol/L</p>
+
<p><strong>K<sub>D</sub></strong> = 3500 μmol/L</p>
-
<br/><br/>
+
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
 +
<div class="hr2"></div>
<h2 id="mod">Modeling Root</h2>
<h2 id="mod">Modeling Root</h2>
<h3>Constructing Virtual Root from Laser Scanning data</h3>
<h3>Constructing Virtual Root from Laser Scanning data</h3>
<p>We use laser scanner to get the root surface topological information. The scanner is Konica Minolta Vivid910, with lens of middle focal length, offered by the College of Information Science and Technology of Beijing Normal University. We obtained 142735 vertices and 253404 faces from it (Fig. 8).</p>
<p>We use laser scanner to get the root surface topological information. The scanner is Konica Minolta Vivid910, with lens of middle focal length, offered by the College of Information Science and Technology of Beijing Normal University. We obtained 142735 vertices and 253404 faces from it (Fig. 8).</p>
<a title="Fig. 8 The virtual peanut from the scanning, scanned by Konica Minolta Vivid910, with lens of middle focal length, offered by the College of Information Science and Technology of Beijing Normal University." href="https://static.igem.org/mediawiki/2014/2/20/Bnu_modeling7.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img class="center" style="width:40%" src="https://static.igem.org/mediawiki/2014/2/20/Bnu_modeling7.png"></a>
<a title="Fig. 8 The virtual peanut from the scanning, scanned by Konica Minolta Vivid910, with lens of middle focal length, offered by the College of Information Science and Technology of Beijing Normal University." href="https://static.igem.org/mediawiki/2014/2/20/Bnu_modeling7.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img class="center" style="width:40%" src="https://static.igem.org/mediawiki/2014/2/20/Bnu_modeling7.png"></a>
-
<p class="fig" >Fig. 8 The virtual peanut from the scanning, scanned by Konica Minolta Vivid910, with lens of middle focal length, offered by the College of Information Science and Technology of Beijing Normal University.</p>
+
<p class="fig" >Fig. 8 The virtual peanut from the scanning, scanned by Konica Minolta Vivid910, with lens of middle focal length, offered by the College of Information Science <br/>and Technology of Beijing Normal University.</p>
<h3>Building concentration gradient</h3>
<h3>Building concentration gradient</h3>
<p>The concentration gradient of the root’s exudate can be regarded as stable in the soil. Exudate concentration decreases as the distance from the root surface increases. (Fig. 9.) </p>
<p>The concentration gradient of the root’s exudate can be regarded as stable in the soil. Exudate concentration decreases as the distance from the root surface increases. (Fig. 9.) </p>
-
<a title="Fig. 9 Sketch of the expected process of chemotaxis. Root exudate’s concentration gradient are shown in different colors. Green points are imaginary E.coli <sub>run</sub>ning up the gradient towards the root." href="https://static.igem.org/mediawiki/2014/1/11/Bnu_modeling8.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img class="center" style="width:40%" src="https://static.igem.org/mediawiki/2014/1/11/Bnu_modeling8.png"></a>
+
<a title="Fig. 9 Sketch of the expected process of chemotaxis. Root exudate’s concentration gradient are shown in different colors. Green points are imaginary <i>E.coli</i> <sub>run</sub>ning up the gradient towards the root." href="https://static.igem.org/mediawiki/2014/1/11/Bnu_modeling8.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img class="center" style="width:40%" src="https://static.igem.org/mediawiki/2014/1/11/Bnu_modeling8.png"></a>
-
<p class="fig">Fig. 9 Sketch of the expected process of chemotaxis. Root exudate’s concentration gradient are shown in different colors. Green points are imaginary E.coli <sub>run</sub>ning up the gradient towards the root.</p>
+
<p class="fig">Fig. 9 Sketch of the expected process of chemotaxis. Root exudate’s concentration gradient are shown in different colors. Green points are imaginary <i>E.coli</i> <br/>running up the gradient towards the root.</p>
-
<p>The concentration value of points is calculated from [8]. Despite limited study directly discussing the spatial distribution patterns of peanut’s root exudates, we can still describe, as [8] proposed for maize and wheat, the spatial distribution of succinate as Y = A<sub>1</sub> * X<sup>-B</sup> (Fig.11), where Y represents the concentration of succinate at the distance X from the peanut root surface, when X >= 1mm. </p>
+
<p>The concentration value of points is calculated from [8]. Despite limited study directly discussing the spatial distribution patterns of peanut’s root exudates, we can still describe, as [8] proposed for maize and wheat, the spatial distribution of succinate as Y = A * X<sup>-B</sup> (Fig.11), where Y represents the concentration of succinate at the distance X from the peanut root surface, when X >= 1mm. </p>
<a title="Fig. 10 Concentration gradient of maize and wheat proposed in [8]. Y-axis is the percentage of C14 activity, reflecting the concentration. " href="https://static.igem.org/mediawiki/2014/3/3f/Bnu_modeling9.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img class="left" style="width:40%" src="https://static.igem.org/mediawiki/2014/3/3f/Bnu_modeling9.png"></a>
<a title="Fig. 10 Concentration gradient of maize and wheat proposed in [8]. Y-axis is the percentage of C14 activity, reflecting the concentration. " href="https://static.igem.org/mediawiki/2014/3/3f/Bnu_modeling9.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img class="left" style="width:40%" src="https://static.igem.org/mediawiki/2014/3/3f/Bnu_modeling9.png"></a>
-
<p class="fig" style="float:left; width:50%; margin-top:300px;">Fig. 10 Concentration gradient of maize and wheat proposed in [8]. Y-axis is the percentage of C14 activity, reflecting the concentration. </p>
+
<p class="fig" style="float:left; width:50%; margin-top:300px;">Fig. 10 Concentration gradient of maize and wheat proposed in [8]. Y-axis is the <br/>percentage of C14 activity, reflecting the concentration. </p>
<div class="clear"></div>
<div class="clear"></div>
<br/>
<br/>
Line 133: Line 136:
<p>We derive the concentration at the bacterium’s vincinity from the formula (3) described above.</p>
<p>We derive the concentration at the bacterium’s vincinity from the formula (3) described above.</p>
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
 +
<div class="hr2"></div>
<h2 id="pro">Programming</h2>
<h2 id="pro">Programming</h2>
<p>The program simulates movement per second, herein attractant is succinate. The pseudo-code reads as follows:</p>
<p>The program simulates movement per second, herein attractant is succinate. The pseudo-code reads as follows:</p>
-
<p><strong>move</strong> (a bunch of E.coli , the concentration gradient)</p>
+
<p><strong>move</strong> (a bunch of <i>E.coli</i> , the concentration gradient)</p>
-
<p>each E.coli has <strong>[state, timer, position, direction, duration_time]</strong></p>
+
<p>each <i>E.coli</i> has <strong>[state, timer, position, direction, duration_time]</strong></p>
<a title="" href="https://static.igem.org/mediawiki/2014/f/fe/Bnu_modeling09.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img  style="width:70%" src="https://static.igem.org/mediawiki/2014/f/fe/Bnu_modeling09.png"></a>
<a title="" href="https://static.igem.org/mediawiki/2014/f/fe/Bnu_modeling09.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img  style="width:70%" src="https://static.igem.org/mediawiki/2014/f/fe/Bnu_modeling09.png"></a>
<p class="fig">
<p class="fig">
Line 143: Line 147:
*** This step is by calculating duration_time (this moment) - timer, and if the result > 0, then it should go on this behavior.
*** This step is by calculating duration_time (this moment) - timer, and if the result > 0, then it should go on this behavior.
</p>
</p>
-
<br/><br/>
+
 
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
 +
<div class="hr2"></div>
<h2 id="res">Results and analysis</h2>
<h2 id="res">Results and analysis</h2>
<h3>Assessment of “Prometheus”</h3>
<h3>Assessment of “Prometheus”</h3>
-
<p>In each <sub>run</sub> estimating the effectiveness of the whole system, the program simulates behaviors of 1000 bacteria over 1 hour (3600 sec). Fig. 11 encompasses 4 different orders of initial succinate concentration on root surface. Compared with the null control (0 M), the engineered E. coli can be successfully recruited around root. Under the concentration of 10<sup>-2</sup> M, the number of bacteria within 1mm around root surface doubled that of null control.</p>
+
<p>In each run estimating the effectiveness of the whole system, the program simulates behaviors of 1000 bacteria over 1 hour (3600 sec). Fig. 11 encompasses 4 different orders of initial succinate concentration on root surface. Compared with the null control (0 M), the engineered <i>E. coli</i> can be successfully recruited around root. Under the concentration of 10<sup>-2</sup> M, the number of bacteria within 1mm around root surface doubled that of null control.</p>
<a title="Fig. 11 Concentration gradient of maize and wheat proposed in [8]. Y-axis is the percentage of C14 activity, reflecting the concentration. " href="https://static.igem.org/mediawiki/2014/d/d2/Bnu_modeling11.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img style="width:80%; margin-left:100px; " src="https://static.igem.org/mediawiki/2014/d/d2/Bnu_modeling11.png"></a>
<a title="Fig. 11 Concentration gradient of maize and wheat proposed in [8]. Y-axis is the percentage of C14 activity, reflecting the concentration. " href="https://static.igem.org/mediawiki/2014/d/d2/Bnu_modeling11.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img style="width:80%; margin-left:100px; " src="https://static.igem.org/mediawiki/2014/d/d2/Bnu_modeling11.png"></a>
-
<p class="fig" style="width:100%;">Fig. 11 The react of E. coli to different concentration of succinate. The numbers of bacteria within 1mm of root surface (y-axis) represents the bacteria successfully attracted and utilized overtime. We assumed the total number of bacteria being 1000. Under the model, the peak performance of “Promethei” are expected between 10<sup>-2</sup> and 10<sup>-3</sup>. Our simulation fits the experiments and the theoretical curve K<sub>D</sub>C / (K<sub>D</sub> + C)<sup>2</sup> well. </p>
+
<p class="fig" style="width:95%;">Fig. 11 The react of <i>E. coli</i> to different concentration of succinate. The numbers of bacteria within 1mm of root surface (y-axis) represents the bacteria successfully <br/>attracted and utilized overtime. We assumed the total number of bacteria being 1000. Under the model, the peak performance of “Promethei” are expected <br/>between 10<sup>-2</sup> and 10<sup>-3</sup>. Our simulation fits the experiments and the theoretical curve K<sub>D</sub>C / (K<sub>D</sub> + C)<sup>2</sup> well. </p>
-
<p>The results amazingly fit the 96-well plate experiment and the theoretical curve (Fig. 12). The <sub>run</sub>ning of the program has the similar mechanism with the experiments done in a laboratory. That can prove our codes are correct from the other aspect.</p>
+
<p>The results amazingly fit the 96-well plate experiment and the theoretical curve (Fig. 12). The running of the program has the similar mechanism with the experiments done in a laboratory. That can prove our codes are correct from the other aspect.</p>
<br/>
<br/>
Line 163: Line 168:
<p class="fig" style="width:2%;float:left; margin-left:450px">(c) </p>
<p class="fig" style="width:2%;float:left; margin-left:450px">(c) </p>
<div class="clear"></div>
<div class="clear"></div>
-
<p class="fig">Fig.12 The contract of the program result in computer (a), the experiment result (b) ( get to know more click here) and the theoretical curve of K<sub>D</sub>C / (K<sub>D</sub>+C)<sup>2</sup> (c).</p>
+
<p class="fig">Fig.12 The contract of the program result in computer (a), the experiment result (b) ( <a href="https://2014.igem.org/Team:BNU-China/Chemotaxis.html" target="_blank">get to know more click here</a>) and the theoretical curve of K<sub>D</sub>C / (K<sub>D</sub>+C)<sup>2</sup> (c).</p>
<p>Considering the overall amount of succinate secreted via root by a plant could reach 13 mg [9], we estimate the concentration on root surface be between 0.1-10 mM. We conservatively reckon that around 20-40% bacteria could eventually “transport” cargos to the root (cut a half for unconsidered conditions in soil). Thus about 60-80% would lead to indirect use in the process, raising the issue of potential Mo pollution.</p>
<p>Considering the overall amount of succinate secreted via root by a plant could reach 13 mg [9], we estimate the concentration on root surface be between 0.1-10 mM. We conservatively reckon that around 20-40% bacteria could eventually “transport” cargos to the root (cut a half for unconsidered conditions in soil). Thus about 60-80% would lead to indirect use in the process, raising the issue of potential Mo pollution.</p>
<h3>Estimation of “Memory” Time</h3>
<h3>Estimation of “Memory” Time</h3>
-
<p>We further estimate the best time interval to set the “memory” time, the time most appropriate for E. coli to suicide. Such time is determined as immediately after most of bacteria has fulfilled their task. With results in Fig. 2, the “memory” time is estimated around 40 minutes. </p>
+
<p>We further estimate the best time interval to set the “memory” time, the time most appropriate for <i>E. coli</i> to suicide. Such time is determined as immediately after most of bacteria has fulfilled their task. With results in Fig. 2, the “memory” time is estimated around 40 minutes. </p>
 +
 
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
 +
<div class="hr2"></div>
<h2 id="rev">Review and forecast</h2>
<h2 id="rev">Review and forecast</h2>
<p>We think our model applies to the situation where bacteria swim freely without casualty. However, for the model has many statements, we could expect nuances from the real complex world. The differences in condition are listed in Table 1.</p>
<p>We think our model applies to the situation where bacteria swim freely without casualty. However, for the model has many statements, we could expect nuances from the real complex world. The differences in condition are listed in Table 1.</p>
<a title="Table 1 The difference between conditions in the real world and those assumed in our model. Shown in red are factors that could significantly alter the eventual result. Viscid layers around root could potentially boost the system’s performance. " href="https://static.igem.org/mediawiki/2014/d/dc/Bnu_modeling13.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img style="width:60%; margin-left:150px;" src="https://static.igem.org/mediawiki/2014/d/dc/Bnu_modeling13.png"></a>
<a title="Table 1 The difference between conditions in the real world and those assumed in our model. Shown in red are factors that could significantly alter the eventual result. Viscid layers around root could potentially boost the system’s performance. " href="https://static.igem.org/mediawiki/2014/d/dc/Bnu_modeling13.png" rel="prettyPhoto"> <span class="overlay zoom" ></span><img style="width:60%; margin-left:150px;" src="https://static.igem.org/mediawiki/2014/d/dc/Bnu_modeling13.png"></a>
-
<p class="fig" style="width:100%">Table 1 The difference between conditions in the real world and those assumed in our model. Shown in red are factors that could significantly alter the eventual result. Viscid layers around root could potentially boost the system’s performance. </p>
+
<p class="fig" style="width:95%">Table 1 The difference between conditions in the real world and those assumed in our model. Shown in red are factors that could significantly alter the eventual result. <br/>Viscid layers around root could potentially boost the system’s performance. </p>
-
<p>We can probably improve the performance by triggering suicide with a threshold concentration, instead of time. But since we lack corresponding element for such function, and with certain biosafety concerns, we stick to time trigger.</p><p>
+
<p>We can probably improve the performance by triggering suicide with a threshold concentration, instead of time. But since we lack corresponding element for such function, and with certain biosafety concerns (to prevent gene pollution), we stick to time trigger.</p><p>
Still, the modeling, albeit simple, provides us with an otherwise infeasible evaluation of the whole project. The Prometheus has certain potential to work as we wished.</p>
Still, the modeling, albeit simple, provides us with an otherwise infeasible evaluation of the whole project. The Prometheus has certain potential to work as we wished.</p>
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
<a href="#top"><p class="fig" style="margin-left:790px">Back to Top</p></a>
 +
<div class="hr2"></div>
<h2 id="ref">Reference</h2>
<h2 id="ref">Reference</h2>
-
<p>[1]Chemotaxis in Escherichia coli analyed by Three-dimensional Tracking, Nature 239:500-504, by Howard C. Berg & Douglas A. Brown</p>
+
<p>[1]Chemotaxis in <i>Escherichia coli</i> analyed by Three-dimensional Tracking, Nature 239:500-504, by Howard C. Berg & Douglas A. Brown</p>
-
<p>[2]E. coli in Motion, Springer, by Howard C. Berg</p>
+
<p>[2]<i>E. coli</i> in Motion, Springer, by Howard C. Berg</p>
<p>[3]The Gradient-Sensing Mechanism in Bacterial Chemotaxis, Proc. Natl. Acad. Sci. USA 69:2509-2512, by Macnab, R. M., D. E. Koshland</p>
<p>[3]The Gradient-Sensing Mechanism in Bacterial Chemotaxis, Proc. Natl. Acad. Sci. USA 69:2509-2512, by Macnab, R. M., D. E. Koshland</p>
-
<p>[4]Temporal Stimulation of Chemotaxis in Escherichia coli, Proc. Natl. Acad. Sci. USA 71:1388-1392, by Brown, D. A., and H. C. Berg</p>
+
<p>[4]Temporal Stimulation of Chemotaxis in <i>Escherichia coli</i>, Proc. Natl. Acad. Sci. USA 71:1388-1392, by Brown, D. A., and H. C. Berg</p>
<p>[5]Directional persistence of chemotactic bacteria in a traveling concentration wave, PNAS, by J. Saragosti, V. Calvez, N. Bournaves, B. Perthame, A. Buguin, and P. Silberzan</p>
<p>[5]Directional persistence of chemotactic bacteria in a traveling concentration wave, PNAS, by J. Saragosti, V. Calvez, N. Bournaves, B. Perthame, A. Buguin, and P. Silberzan</p>
-
<p>[6]Modeling E.coli Tumbles by Rotational Diffusion Implications for Chemotaxis, PLOS one, by Jonathan Saragosti, Pascal SilBerzran, Axel Buguin</p>
+
<p>[6]Modeling <i>E.coli</i> Tumbles by Rotational Diffusion Implications for Chemotaxis, PLOS one, by Jonathan Saragosti, Pascal SilBerzran, Axel Buguin</p>
<p>[7]Algorithms 4th edition, Pearson, Robert Sedgwick and Kane Wayne
<p>[7]Algorithms 4th edition, Pearson, Robert Sedgwick and Kane Wayne

Latest revision as of 03:18, 18 October 2014

1 2

Modeling

Abstract

Analyzing dynamics of E.coli’s chemotaxis and estimating effectiveness of the "Prometheus" require high accuracy, and thus hard to realize with experiments. Here, we construct a virtual peanut root-E.coli system (by c++ and OpenGL), storing the 3-D environmental information by k-d tree. We quantify the process in which “Promethei”, our engineered bacteria, carrying Mo, move towards peanut roots. We at last estimate around 20-40% of bacteria could work in reality, and that the best “memory” time be around 40 min (See here). All codes (in c++) of the project are available*.

*The size of the code is too big, please contact Li Si-yao (201211211043@mail.bnu.edu.cn) to get it.



Fig. 1 The Motion Graph of the screenshots while running the Simulation Program. The purple points are E.coli tending towards the root surface.

Back to Top

E.coli model

Mechanism

We divide Escherichia coli’s behavior into two parts: “running” straight lines, and “tumbling” for reorientation[1]. When attractant concentration uniformed, E.coli can be described as doing “random walk”[4], and the running time fits a normal distribution. With spatial concentration gradients of chemical attractance,when bacteria sense a higher attractant concentration than before,the mean of the statistical distribution of running time increases. Otherwise, the time keeps the same as that in the uniform environment, no matter how high or low the concentration is[4].


Running

The running of E.coli can be descibed as a straight line. The distribution of running duration time fits a normal curve, with mean τrun0 and standard division σrun.

Brown and Berg[4] suggested that the mean τ has a functional relationship, shown as follows, with the current attractant (in their case, glutamate) concentration and the change of concentration over time. (See formula (1),(2)).

KD is the dissociation constant of the complex formed by glutamate and its receptor. Pb is the fraction of receptors (protein) bound with the attractant.

In their model, the mean velocity during a run, however, is independent from C and dC/dt [5,6].

Brown and Berg found the curve KDC/(KD + C)2 fits the chemotaxis sensitivity assays well (Fig. 2).

Obtain KD

Brown and Berg found the curve KDC/(KD + C)2 fits the chemotaxis sensitivity assays well (Fig. 2).


Fig. 2


We cannot determine the precise dissociation constants KD by experiment. However, the trend revealed in out experiments (Fig. 3) fits with the model well. With data from McfR-Succinate experiment, we set the KD value of succinate-receptor complex as 3.5 mM by Least Square method.

Fig. 3 the result from experiment.

Fig. 4 Fitting the experiment by a( KDC / (KD + C)2) + d, where a = 389, d = 260 and KD is 3.5mM. This work is
made by c++ and OpenGL. The code to get KD and a and to make sketch is available. Please contact us.

Given a concrete constant number of α in function (1), we would be able to calculate the value of τrun. Series of experiments[4] convinced us to* set α as 660 sec.

*We have no time to repeat the experiments with succinate.

In conclusion, modeling running is based on following statement:

  • The velocity during a run is constant (always equal to the mean);
  • E.coli move straightly during a run, and neither die nor divide in the process. The effect of quorum sensing and nutritional factors are eliminated;
  • The E.coli’s running time trun is a random number from N(<τrun>,σrun);
  • The mean duration time <τrun> can be derived from formula (1) and (2), given the value of C and dC/dt;
  • The specific duration_time of a single bacterium changes in the same pace with the changes of the statistical distribution’s mean (relying on C and dC/dt). (Fig. 5)

Fig. 5

For a single bacterium, its surroundings decides the statistical distribution of its own running time. In other words, the discrepancy doesn’t change:

where trun is the specified starting and stopping time for a single bacterium. τrun2 is the current mean duration time and τrun1 is the mean duration time at last second; both are derived from the function (1) and (2).


Tumbling

E. coli reorient between two runs.In modeling this process, we assume:

  • E. coli do not change positions as they are tumbling;
  • the direction after reorientation is randomly chosen.

Fig. 6 Schematic drawing of a tumble, from [2], page 45


The angle dθ is not strictly random [1,5,6] (its distribution shown in Fig. 7). However, the distribution is not formularily clear and would require too much effort to build a random number fitting this distribution, thus we roughly regard it as a random process.

Fig. 7 the distribution of change in direction ( A figure from [1] ).


*The statement is made to calculate the process more easily. We suppose it is the run duration time, not the difference of velocity, is the main variant of the whole chemotaxis system.

Coefficients

velocity = 100μm

run duration time mean in uniform environment = 1.3

run duration standard division = 3

mean tumble duration time = 0.14 sec

alpha (See equation (1)) = 660 sec

KD = 3500 μmol/L

Back to Top

Modeling Root

Constructing Virtual Root from Laser Scanning data

We use laser scanner to get the root surface topological information. The scanner is Konica Minolta Vivid910, with lens of middle focal length, offered by the College of Information Science and Technology of Beijing Normal University. We obtained 142735 vertices and 253404 faces from it (Fig. 8).

Fig. 8 The virtual peanut from the scanning, scanned by Konica Minolta Vivid910, with lens of middle focal length, offered by the College of Information Science
and Technology of Beijing Normal University.

Building concentration gradient

The concentration gradient of the root’s exudate can be regarded as stable in the soil. Exudate concentration decreases as the distance from the root surface increases. (Fig. 9.)

Fig. 9 Sketch of the expected process of chemotaxis. Root exudate’s concentration gradient are shown in different colors. Green points are imaginary E.coli
running up the gradient towards the root.

The concentration value of points is calculated from [8]. Despite limited study directly discussing the spatial distribution patterns of peanut’s root exudates, we can still describe, as [8] proposed for maize and wheat, the spatial distribution of succinate as Y = A * X-B (Fig.11), where Y represents the concentration of succinate at the distance X from the peanut root surface, when X >= 1mm.

Fig. 10 Concentration gradient of maize and wheat proposed in [8]. Y-axis is the
percentage of C14 activity, reflecting the concentration.


As the relationship is applicable only when x >= 1 mm. We estimate the percentage decrease linearly from 0 to 1 mm, and the Y value is 70% at x = 1 mm. Thus, the function of the concentration C on the distance to the root surface x is

C = 100 - 0.3 * x if x < 100 and

C = 70 * x^(-1.2) if x > 100.         (3)

Using k-d tree to store points

K-dimension tree[7] is a smart and widely-used data structure to store information of points in k-dimension space. It is very fast, typically requiring O(logN) time, to find the nearest point with a given position. Here we implement k-d tree to construct all 10,000 points.

We derive the concentration at the bacterium’s vincinity from the formula (3) described above.

Back to Top

Programming

The program simulates movement per second, herein attractant is succinate. The pseudo-code reads as follows:

move (a bunch of E.coli , the concentration gradient)

each E.coli has [state, timer, position, direction, duration_time]

*The points of concentration gradient are stored in a 3-d tree It will expounded clearly in the part below.
** The mean value is calculate by the formula(1), (2) and (4).
*** This step is by calculating duration_time (this moment) - timer, and if the result > 0, then it should go on this behavior.

Back to Top

Results and analysis

Assessment of “Prometheus”

In each run estimating the effectiveness of the whole system, the program simulates behaviors of 1000 bacteria over 1 hour (3600 sec). Fig. 11 encompasses 4 different orders of initial succinate concentration on root surface. Compared with the null control (0 M), the engineered E. coli can be successfully recruited around root. Under the concentration of 10-2 M, the number of bacteria within 1mm around root surface doubled that of null control.

Fig. 11 The react of E. coli to different concentration of succinate. The numbers of bacteria within 1mm of root surface (y-axis) represents the bacteria successfully
attracted and utilized overtime. We assumed the total number of bacteria being 1000. Under the model, the peak performance of “Promethei” are expected
between 10-2 and 10-3. Our simulation fits the experiments and the theoretical curve KDC / (KD + C)2 well.

The results amazingly fit the 96-well plate experiment and the theoretical curve (Fig. 12). The running of the program has the similar mechanism with the experiments done in a laboratory. That can prove our codes are correct from the other aspect.


(a)


(b)

(c)

Fig.12 The contract of the program result in computer (a), the experiment result (b) ( get to know more click here) and the theoretical curve of KDC / (KD+C)2 (c).

Considering the overall amount of succinate secreted via root by a plant could reach 13 mg [9], we estimate the concentration on root surface be between 0.1-10 mM. We conservatively reckon that around 20-40% bacteria could eventually “transport” cargos to the root (cut a half for unconsidered conditions in soil). Thus about 60-80% would lead to indirect use in the process, raising the issue of potential Mo pollution.

Estimation of “Memory” Time

We further estimate the best time interval to set the “memory” time, the time most appropriate for E. coli to suicide. Such time is determined as immediately after most of bacteria has fulfilled their task. With results in Fig. 2, the “memory” time is estimated around 40 minutes.

Back to Top

Review and forecast

We think our model applies to the situation where bacteria swim freely without casualty. However, for the model has many statements, we could expect nuances from the real complex world. The differences in condition are listed in Table 1.

Table 1 The difference between conditions in the real world and those assumed in our model. Shown in red are factors that could significantly alter the eventual result.
Viscid layers around root could potentially boost the system’s performance.

We can probably improve the performance by triggering suicide with a threshold concentration, instead of time. But since we lack corresponding element for such function, and with certain biosafety concerns (to prevent gene pollution), we stick to time trigger.

Still, the modeling, albeit simple, provides us with an otherwise infeasible evaluation of the whole project. The Prometheus has certain potential to work as we wished.

Back to Top

Reference

[1]Chemotaxis in Escherichia coli analyed by Three-dimensional Tracking, Nature 239:500-504, by Howard C. Berg & Douglas A. Brown

[2]E. coli in Motion, Springer, by Howard C. Berg

[3]The Gradient-Sensing Mechanism in Bacterial Chemotaxis, Proc. Natl. Acad. Sci. USA 69:2509-2512, by Macnab, R. M., D. E. Koshland

[4]Temporal Stimulation of Chemotaxis in Escherichia coli, Proc. Natl. Acad. Sci. USA 71:1388-1392, by Brown, D. A., and H. C. Berg

[5]Directional persistence of chemotactic bacteria in a traveling concentration wave, PNAS, by J. Saragosti, V. Calvez, N. Bournaves, B. Perthame, A. Buguin, and P. Silberzan

[6]Modeling E.coli Tumbles by Rotational Diffusion Implications for Chemotaxis, PLOS one, by Jonathan Saragosti, Pascal SilBerzran, Axel Buguin

[7]Algorithms 4th edition, Pearson, Robert Sedgwick and Kane Wayne

[8]Spatial distribution of root exudates of five plant species as assessed by labeling, J. Plant Nutr. Soil Sci. 2006, 169,360-362 , by Daniela Sauer, Yakov Kuzyakov, and Karl Stahr

[9] Mineral Nutrition of Higher Plants. China Agricultural University Press. 2008.



Prev Page Next Page