Team:UPO-Sevilla/Modeling

From 2010.igem.org

(Difference between revisions)
 
(6 intermediate revisions not shown)
Line 12: Line 12:
<html>
<html>
<div class=contentBC>
<div class=contentBC>
 +
<h1>Models</h1>
<h1>Models</h1>
-
 
The members of the dry lab are simulating the different components of the full system. Three main components can be identified:
The members of the dry lab are simulating the different components of the full system. Three main components can be identified:
Line 23: Line 23:
</ol>
</ol>
-
 
+
<p>For the first two points, we have developed a simulation using Java. The explanation of the simulation can be found in the Chemotaxis link below. Moreover, you can <b>download the bacterial crowding simulation</b> <a href="https://static.igem.org/mediawiki/2010/3/3f/BacterialCrowdingSimulation.zip">here</a>
-
 
+
-
<h2>Chemoattractant Diffusion</h2>
+
-
 
+
-
<p>
+
-
The basic equations for the diffusion of the chemoattractant in the medium are the <a href="http://en.wikipedia.org/wiki/Fick%27s_laws_of_diffusion" target="_blank">Fick laws of diffusion</a>, which govern the variation of the concentration of a substance within a medium.
+
</p>
</p>
<p>
<p>
-
The flux <b>J</b> (that is, the amount of substance that flows through a given surface per unit of time mol m<sup>-2</sup>s<sup>-1</sup>) is given by:
+
The circuits and transduction signals inside the cell is simulated by using Matlab's Simbiology Toolbox. The model and results can be found in the Signaling link below.
</p>
</p>
-
 
+
<div class="center">
-
<center>
+
                  <a href="https://2010.igem.org/Team:UPO-Sevilla/Modeling/Chemotaxis"> <img class="subBanner" src="https://static.igem.org/mediawiki/2010/2/24/BacterialCrowdingChemotaxis.png" alt="Chemotaxis" /></a>
-
<img src="https://static.igem.org/mediawiki/2010/6/6f/UPO-Model-Eq1.png" width="150 "alt="First Fick Law"/>
+
                  <a href="https://2010.igem.org/Team:UPO-Sevilla/Modeling/Signaling"><img class="subBanner" src="https://static.igem.org/mediawiki/2010/5/59/BacterialCrowdingSignaling.png" alt="Signaling" /></a>
-
</center>
+
</div>  
-
 
+
-
 
+
-
 
+
-
<p>
+
-
where &phi; is the concentration (mol m<sup>-3</sup>) in a given point. <i>D</i> is a constant called the diffusion coefficient, and that depends on the medium .
+
-
</p>
+
-
<p>
+
-
Basically, the equation states that the is directed towards places with lower concentration (thus the minus sign). If the concentration is constant in the space (&nabla;&phi;=0) there is no flux.
+
-
</p>
+
-
 
+
-
<p>
+
-
If the flux is known, it is possible to determine the amount of substance that goes through a small surface <b>S</b> and a small amount of time <i>dt</i></p>
+
-
 
+
-
<center>
+
-
<img src="https://static.igem.org/mediawiki/2010/f/f8/UPO-Model-Eq2.png" width="170 "alt="First Fick Law"/>
+
-
</center>
+
-
 
+
-
<p>
+
-
In order to simulate the diffusion, we define the environment and discretize it in very small cells. Each cell determines a given volume <i>V</i>, and has a surface <b>S</b>. At a given time instant, the cell has an amount of substance <i>c</i> (and then a concentration <i>c</i>\<i>V</i>).
+
-
</p>
+
-
 
+
-
<p>
+
-
If the cells and time step &Delta;<i>t</i> are small, we can consider that the gradient of concentration can be approximated though the differences in concentration between a cell <i>i</i> and 4 (or 8) neighbors <i>j</i>. Thus:
+
-
</p>
+
-
 
+
-
<center>
+
-
<img src="https://static.igem.org/mediawiki/2010/e/ea/UPO-Model-Eq3.png" width="190" alt="First Fick Law"/>
+
-
</center>
+
-
 
+
-
<p>
+
-
and then, the amount of substance that diffusses from <i>i</i> to <i>j</i>:
+
-
</p>
+
-
 
+
-
<center>
+
-
<img src="https://static.igem.org/mediawiki/2010/6/6f/UPO-Model-Eq4.png" width="190" alt="First Fick Law"/>
+
-
</center>
+
-
 
+
-
<p>
+
-
The following figure illustrates the basic elements of the simulation. The flux <b>J</b> between cells is computed by the difference of concentrations. Then, this flux is used to compute the amount of substance that will flow to the neighbour cell. The amount is proportional to the flux and the common surface between cells.
+
-
</p>
+
-
 
+
-
<center>
+
-
<img src="https://static.igem.org/mediawiki/2010/9/97/UPO-Model-diffusion.png" width="300" alt="Basic parameters of the diffusion simulation."/>
+
-
</center>
+
-
 
+
-
<h2> Bacteria motion</h2>
+
-
 
+
-
<p>
+
-
The main actuator of <i>E. Coli</i> is a flagellar motor that can rotate clock and counterclockwise. Through a set of transmembrane receptors proteins, <i>E. Coli</i> is able to detect chemoattractants. Moreover, this detection influences the motion of the flagellar motor <a href=#topp>[Topp and Gallivan, 2007]</a>.
+
-
</p>
+
-
 
+
-
<p>
+
-
E. Coli has two main motion modes, which we will name:
+
-
 
+
-
<ol>
+
-
<li> Random Walk </li>
+
-
<li> Gradient climbing </li>
+
-
</ol>
+
-
</p>
+
-
 
+
-
<p><b>Random walk mode</b></p>
+
-
 
+
-
<p>
+
-
When no gradient of chemoattractant is present, <i>E. Coli</i> is in random walk mode. In this case, the bacteria performs smooth runs followed by tumbles.
+
-
</p>
+
-
 
+
-
<p>
+
-
Mathematically, we will model this as a Brownian motion:
+
-
</p>
+
-
 
+
-
 
+
-
<center>
+
-
<img src="https://static.igem.org/mediawiki/2010/b/bc/UPO-Model-Eq5.png" width="300 "alt="First Fick Law"/>
+
-
</center>
+
-
 
+
-
<p>
+
-
where <b>x</b> is the position of the bacteria and <b>v</b> is the velocity. This (vector) velocity is randomly sampled from a normal distribution of zero mean and a certain covariance matrix that models the potential .
+
-
</p>
+
-
 
+
-
 
+
-
<center>
+
-
<img src="https://static.igem.org/mediawiki/2010/b/bd/UPO-Model-Eq6.png" width="150 "alt="First Fick Law"/>
+
-
</center>
+
-
 
+
-
<p>
+
-
The bacteria will move for a time &Delta; <i>t</i><sub>rm</sub> with constant velocity <b>v</b>(<i>t</i>). At next time instant, a new velocity is (randomly selected) selected 
+
-
</p>
+
-
 
+
-
 
+
-
<b>Gradient climbing behaviour</b>
+
-
 
+
-
<p>
+
-
When a positive difference of concentration (a gradient) on the chemoattractant is detected, the bacteria enters into a new mode that we will call gradient climbing.
+
-
</p>
+
-
 
+
-
<p>
+
-
In this mode, the flagellar motor tends to move counterclockwise; as a result, the smooth runs last for more time, and the tumbling frequency decreases.
+
-
</p>
+
-
 
+
-
In order to model that, we use the same model than above, but a larger tumbling interval (smaller tumbling frequency) &Delta; <i>t</i><sub>gc</sub>;
+
-
 
+
-
<p id="topp">
+
-
<small>[Topp and Gallivan, 2007] Topp, S. and Gallivan, J. (2007). Guiding Bacteria with Small Molecules and RNA. J. Am. Chem. Soc., 129:6807–6811.</small>
+
-
</p>
+
-
 
+
-
<h2>Circuits and devices</h2>
+
-
 
+
-
 
+
-
<center>
+
-
<a href="https://2010.igem.org/Image:IGEMUPO-model-small.png">
+
-
<img src="https://static.igem.org/mediawiki/2010/8/8a/IGEMUPO-model-small.png" width="700" alt="Simbiology model"/>
+
-
</a>
+
-
</center>
+
-
 
+
-
 
+
-
<h1>Modeling Tools</h1>
+
</div>
</div>
</html>
</html>
 +
{{:Team:UPO-Sevilla/footer}}
{{:Team:UPO-Sevilla/footer}}
</div>
</div>

Latest revision as of 20:09, 26 October 2010

Models

The members of the dry lab are simulating the different components of the full system. Three main components can be identified:
  1. The difussion of the chemoattractant through the medium.
  2. The motion of the bacterias through the medium due to the gradient on the chemoattractant concentration.
  3. The circuits and devices for the chemoattractant generation within the bacteria.

For the first two points, we have developed a simulation using Java. The explanation of the simulation can be found in the Chemotaxis link below. Moreover, you can download the bacterial crowding simulation here

The circuits and transduction signals inside the cell is simulated by using Matlab's Simbiology Toolbox. The model and results can be found in the Signaling link below.


Footer