Team:Paris/Modeling/BOB/Simulations

From 2008.igem.org

(Difference between revisions)
(Oscillations)
Line 49: Line 49:
** Here, with n=5, we begin to see some oscillations at the begining...
** Here, with n=5, we begin to see some oscillations at the begining...
[[Image:essai_simul_n=5_2.jpg|500px|center]]
[[Image:essai_simul_n=5_2.jpg|500px|center]]
-
... and we may see that by zooming, there are some subdued oscillations.[[Image:essai_simul_n=5_zoom.jpg|500px|center]]** We then tried to increase drastically the value of n, to 300. Even though this had no biological reality, this met our simulating purpose
+
... and we may see that by zooming, there are some subdued oscillations.[[Image:essai_simul_n=5_zoom.jpg|500px|center]]
 +
** We then tried to increase drastically the value of n, to 300. Even though this had no biological reality, this met our simulating purpose
[[Image:essai_simul_n=300_2.jpg|500px|center]]
[[Image:essai_simul_n=300_2.jpg|500px|center]]

Revision as of 14:02, 2 September 2008

(Under Construction)

Simulations and Mathematical analysis

Oscillations

  • We wanted to see if, from a mathematical point of view, it was possible for the "short" system presented above to hover.
Oscillations1.jpg
  • Here are the equations we took into account :
FlhDC dynamics simul 3.jpg
FliA dynamics simul 3.jpg
TetR simul 3.jpg
  • The equations are normalized (thus the degradation term set to 1), as well as the parameters :
Parameter Table
Parameter Normalized Value
βFlhDC 1
θFlhDC 0.2222
n 2
βFliA 0.1429
β'FliA 0.8581
βLasI 0.2222
β'LasI 0.7778
  • At first, we quickly run simulations with different values of the parameters of the Hill function. Here are plot examples :
    • With n=2, we see no oscillations
      Essai simul 2.jpg
    • Here, with n=5, we begin to see some oscillations at the begining...
Essai simul n=5 2.jpg
... and we may see that by zooming, there are some subdued oscillations.
Essai simul n=5 zoom.jpg
    • We then tried to increase drastically the value of n, to 300. Even though this had no biological reality, this met our simulating purpose
Essai simul n=300 2.jpg

FIFO

  • The goal here is to present the results of the simulations we made concerning the FIFO part of the system.

Here is the system we implementated using Matlab (see the corresponding codes)

Subsystem1.jpg

and the corresponding equations (for more detailed information see our establishment of the model).

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 wanted to see if our predictions were accurate or not. We then solved the equations, forcing the behavior of FhlDC. In a first step, we imposed a constant production term of 1. Then, at a certain time, we set this production term to zero :
FlhDC Test FIFO cresc.jpg
FlhDC Test FIFO decresc.jpg

In fact we assumed that this behavior for FlhDC was acceptable regarding its estimated behavior in the whole system.

  • We saw previously that without FliA, the FIFO would presumably not work. We then simulated a first system, where [FliA] stays to zero value.
Essai without fliA.jpg

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.
Essai with fliA.jpg

  • We observe on these plots that the behavior is quite the one we expected, and that the FIFO is realized. 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.

Synchronisation


simul plsrs cells
Etude maths