Team:Montreal/Modeling

From 2008.igem.org

(Difference between revisions)

Revision as of 20:47, 1 August 2008

Home	The Team	The Project	Parts	Modeling	Notebook	Links	Sponsors

Modeling

Theorists: Vincent Quenneville-Bélair and Alexandra Ortan

Background

Using coupled differential equations, we are modeling the repressilator, which is a network of three genes, whose product proteins are repressing each other's growth. This cycle is taking place in each of a colony of cells, who communicate amongst themselves by exchanging an autoinducer molecule. The model attempts to take into account a sparse, heterogeneous distribution of cells with depletion of the autoinducer molecule and leakage.

For this purpose, the model is currently being coded up in Mathematica. The simulations, based on that continuous model, are generated using xCellerator and NDelayDSolve for different cell configurations (fig. 1); the results obtained are graphs of the concentration of molecules in the system versus time (fig. 2). That is, our interest lies in the phase difference between clusters of cells. As of now, a low number of cells is being used for testing, but a higher one will be reached later.

Journal

The Mathematica notebook used to model the repressillator incorporates xCellerator to create the differential equations representing the chemical reactions at play. The system was then solved using NDSolve, a standard numerical solver in Mathematica -- which is, however, unable to deal with delay equations. Hence, the equations were not evaluated at retarded times at first.

The first solution found was to add null cells in the network. The nulls cells, or pseudo cells, are just grid point where no internal cell reaction occurs, but only diffusion to the next grid point. However, the results depended to heavily on the number of such null cells and, as we did not have a mean of estimating what would be realistic, did not seem reliable enough.

We are now replacing NDSolve with NDelayDSolve to solve the system and we are therefore able to evaluated our equations at retarded times. We are not using null cells anymore. The value of the delays corresponds to the time the molecules of AI diffuse the intercellular distance: the reasoning being that each cells effectively "sees" the concentration of another cell at an earlier time -- the time the molecules take to diffuse between them.

For a given distance between the two clusters of cells (fig. 1), we plot the concentration of AI in each cells versus time. From there, we extract two informations: synchronization time and phase difference (fig. 3-6).

Figure 1: Although the cells are shown using a Voronoi graph, the cells are represented internally by points.	Figure 2: Simulation using Mathematica of the Repressilator Network with four cells shown in fig. 1. The plot shows the concentration of AI against time (in arbitrary units).	Figure 3: Time taken to have the two clusters synchronized (error of 5) versus distance (in arbitrary units).	Figure 4: Phase difference between the clusters with step 1.
Figure 5: Phase difference between the clusters with step 2.	Figure 6: Phase difference between the clusters with step 10.

@@ Line 35: / Line 35: @@
 Image:Synch-speed.png|Figure 3: Time taken to have the two clusters synchronized (error of 5) versus distance (in arbitrary units).
 Image:Phase-vs-distance-1.png|Figure 4: Phase difference between the clusters with step 1.
-Image:Phase-vs-distance-5.png|Figure 5: Phase difference between the clusters with step 5.
+Image:Phase-vs-distance-2.png|Figure 5: Phase difference between the clusters with step 2.
 Image:Phase-vs-distance-10.png|Figure 6: Phase difference between the clusters with step 10.
 </gallery>
 |}