Revision as of 20:46, 26 September 2012

iGEM Grenoble 2012

iGEM 2012
main page

The goal of this part was only to check the speed of the diffusion. We use the same model as the 2011 IGEM team of Grenoble, which is based on the Bangalore 2007 IGEM team model. We first have to add the diffusion terms in the equations thus we get:

Because of the temporal and special derivatives, we could not use a classic matlab solver to solve this set of equations. The approximation we used consisted in dividing the space into a grid:

It is the finite difference method. We thus get:

with h=l_x/N, where l_x is the length of the grid and N is the number of points of discretization along x. By using the same approximation on y_i, and assuming that l_y=l_x, and that we have the same number of points of discretization, we get:

In addition, we also use the finite difference method on the time scale. Thus, we get:

where p=∆t/Nt where ∆t is the interval of time we work in and N_t the number of points of our time discretization.

With these approximations, we could solve the equations. Here are the scripts we used for the simulations.

First visual simulation

Here we give a video, just in the goal to illustrate the diffusion, you can find the scripts here. We model the diffusion in a petri dish. We put a drop of cAMP in the middle of the dish, and we can check how it diffuses through the bacteria. What we can see on the video is the evolution of adenylate cyclase. We give the evolution of the system during 280 minutes. Here is the colors scale, representing the intensity of the signal, and the link to the video:

Quorum sensing video

Temporal evolution

Here we give the temporal evolution of Ca at the drop of cAMP.

Then we give the temporal evolution of Ca at 0,5 cm of the drop of cAMP:

We notice that it takes 100 minutes more than at the drop.

Conclusion: If we use a test tube to make the detection, as soon as one bacterium is turned on (as we can get a signal), we will have to wait 100 minutes to have the whole test tube green.

@@ Line 4: / Line 4: @@
 <body id="Modeling">
 <div id="cadre">
-<!--<section>
-<center>
-<a href="https://2012.igem.org/wiki/index.php?title=Team:Grenoble/Modeling/Amplification/ODE"><img src="https://static.igem.org/mediawiki/2012/e/ef/ODE.png" alt="" /></a>
-&nbsp;&nbsp;&nbsp;&nbsp;
-<a href="https://2012.igem.org/wiki/index.php?title=Team:Grenoble/Modeling/Amplification/Sensitivity"><img src="https://static.igem.org/mediawiki/2012/a/a4/Sensitivity_and_parameters.png" alt="" /></a>
-&nbsp;&nbsp;&nbsp;&nbsp;
-<a href="https://2012.igem.org/wiki/index.php?title=Team:Grenoble/Modeling/Amplification/Quorum"><img src="https://static.igem.org/mediawiki/2012/6/65/Quorum_Sensing.png" alt="" /></a>
-&nbsp;&nbsp;&nbsp;&nbsp;
-<a href="https://2012.igem.org/wiki/index.php?title=Team:Grenoble/Modeling/Amplification/Stochastic"><img src="https://static.igem.org/mediawiki/2012/a/ad/Stochastic_analysis.png" alt="" /></a>
-</center>
-</section>-->
 <section>
 The goal of this part was only to check the speed of the
@@ Line 23: / Line 13: @@
 </br>
 </br>
-Because of the temporal and special derivatives, we couldn’t use a classic matlab solver to solve this set of equations. The approximation we used consisted in dividing the space into a grid:
+Because of the temporal and special derivatives, we could not use a classic matlab solver to solve this set of equations. The approximation we used consisted in dividing the space into a grid:
 </br>
 </br>