Team:Paris/Analysis/Math+Sim

From 2008.igem.org

Revision as of 10:31, 26 October 2008 by Felipe (Talk | contribs)


Mathematical analysis and simulations


FIFO

  • In order to go further into the details of mathematical analysis of the FIFO, let us remind you the set of equations used and the corresponding network (note that you can have an explanation on the way they have been built on the Model Construction page):


Subsystem1.jpg
FliA dynamics.jpg
CFP.jpg
YFP.jpg
RFP.jpg

where CFP, YFP, and RFP will be denoted below as respectively Z1,Z2 and Z3. We have implemented this system using Matlab (see the corresponding codes)









  • Then, in order to obtain simulations of the expected FIFO behaviour, we set ideal conditions, by imposing flhDC as a step function, equal to one then to zero.


FlhDC Test FIFO cresc.jpg
FlhDC Test FIFO decresc.jpg
  • We saw during the overall description of the core system that without FliA, the FIFO would presumably not work. We then simulated a first system, where [FliA] stays to zero value to confirm that qualitative conclusion :


Essai without fliA.jpg

Indeed, we may see that there is a LIFO behavior rather than the FIFO we expect...

  • Then, we simulated the entire system, to check if we had
    • the lasting burst due to FliA (more important for Z3 than for Z2, and more important for Z2 than Z1) in the increasing phase.
    • the effect of fliA which maintained the concentrations to their maximum (more important for Z3 than for Z2, and more important for Z2 than Z1) in the decreasing phase. We have obtained the following curves, which are in perfect agreement with the theoretical definition of a 'FIFO' process :


Essai with fliA.jpg

  • FliA enables the curves to cross, and adds a delay on the genes that are most affected, with gives a better observability of the FIFO behavior.

Oscillations

  • In addition to the previous FIFO system, we will add a negative feedback on FlhDC promoter to observe oscillations. Thus, we obtained the following network and the corresponding equations :
Oscillations1.jpg
Eqn flhDC.jpg
FliA dynamics.jpg
CFP.jpg
YFP.jpg
Eqn EnvZ-RFP.jpg

Parameters' changes

Since the core system has difficulties to produce oscillations, we perform a parameter search and we study the impact of parameters' changes in the core system. We can brain-storm about the possible modifications on the core system. The following table sumarize our main findings:

Description System Evaluation Observations
Core system
Core system mini.png
Orgsad.gif Non oscillating system
Removing the activation of envZ via FlhDC Core system delta flhDC envZ dotted mini.png Orgstar.gif Modest damped oscillations
Removing the auto-activation of flia Core system delta flia flia.png Orgstar.gifOrgstar.gif Damped oscillations
Removing both: the activation of envZ via FlhDC and the auto-activation of flia Core system delta flhDC envZ delta flia flia.png Orgstar.gifOrgstar.gif Damped oscillations with greater amplitud
Enhancing the inhibition of FlhDC via envZ; Removing both: the activation of envZ via FlhDC and the auto-activation of flia Core system envZ plus.png Orgstar.gifOrgstar.gifOrgstar.gif Conserved oscillations
Synchronized oscillations should be searched taking the last system as a starting point.