Team:Paris/Analysis/Math+Sim
From 2008.igem.org
(→Oscillations) |
(→Oscillations) |
||
(33 intermediate revisions not shown) | |||
Line 1: | Line 1: | ||
{{Paris/Menu}} | {{Paris/Menu}} | ||
- | + | {{Paris/Header|Mathematical Analysis and Simulations}} | |
- | + | {{Paris/Section_contents_analysis}} | |
- | + | ||
= FIFO = | = FIFO = | ||
- | * In order to go further into the details of mathematical analysis of the FIFO, let us remind you the | + | * In order to go further into the details of mathematical analysis of the FIFO, let us remind you the genetic network that has been designed to produce a First In - First Out behaviour and the result of the corresponding model : |
<br> | <br> | ||
- | [[Image: | + | |
- | [[Image:FliA_dynamics.jpg| | + | [[Image:FIFO.png|400px|left]][[Image:essai_with_fliAbis.jpg|500px|center]] |
- | [[Image:CFP.jpg| | + | |
- | [[Image:YFP.jpg| | + | Here is the set of equations used for the numerical simulations, with the parameters found in bibliographical studies (note that you can have an explanation on the way they have been built on the [[Team:Paris/Network_analysis_and_design/Core_system/Model_construction|Model Construction]] page): |
- | [[Image:RFP.jpg| | + | |
+ | [[Image:FliA_dynamics.jpg|center]] | ||
+ | [[Image:CFP.jpg|center]] | ||
+ | [[Image:YFP.jpg|center]] | ||
+ | [[Image:RFP.jpg|center]] | ||
where CFP, YFP, and RFP will be denoted below as respectively Z1,Z2 and Z3. | where CFP, YFP, and RFP will be denoted below as respectively Z1,Z2 and Z3. | ||
We have implemented this system using Matlab (see the [[Team:Paris/Modeling/BOB/Simulations/Codes#FIFO|corresponding codes]]) | We have implemented this system using Matlab (see the [[Team:Paris/Modeling/BOB/Simulations/Codes#FIFO|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. | * 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. | ||
Line 31: | Line 24: | ||
[[Image:FlhDC_Test_FIFO_decresc.jpg|center]] | [[Image:FlhDC_Test_FIFO_decresc.jpg|center]] | ||
- | * | + | * We then simulated a first system, where [FliA] stays to zero value to confirm the qualitative conclusion that '''without FliA''' the FIFO would presumably not work : |
<br> | <br> | ||
- | [[Image: | + | [[Image:essai_without_fliAbis.jpg|600px|center]]<br> |
Indeed, we may see that there is a LIFO behavior rather than the FIFO we expect... | Indeed, we may see that there is a LIFO behavior rather than the FIFO we expect... | ||
- | * | + | For a qualitative understanding of what is going on with and without FliA : |
+ | {{Paris/Toggle|Qualitative reasoning about the FIFO|Team:Paris/FIFO}} | ||
+ | |||
+ | * The main differences that appear between the entire system (that means '''with FliA''') and the sparse on ('''without FliA'''), are the following : | ||
** 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 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 | + | ** 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 as previously shown. |
- | + | ||
- | + | ||
* 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. | * 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 = | = 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 (notice that we decide to model in a single arrowthe process through which that EnvZ repress FlhDC): | * 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 (notice that we decide to model in a single arrowthe process through which that EnvZ repress FlhDC): | ||
- | [[Image: | + | [[Image:Reseauoscillfinal.png|400px|center]] |
<br> | <br> | ||
[[Image:eqn_flhDC.jpg|center]] | [[Image:eqn_flhDC.jpg|center]] | ||
Line 61: | Line 56: | ||
<br> | <br> | ||
Nonetheless, there are two questions we could focus on: | Nonetheless, there are two questions we could focus on: | ||
- | * Is it possible to prove theoretically that it will not produce oscillations ? For detailled answers : | + | * Is it possible to prove theoretically that it will not produce oscillations ? For detailled answers : |
- | * Does this result depend on the integration methods used ? For more detailled considerations : {{Paris/Toggle|Impact of integration methods|Team:Paris/ | + | {{Paris/Toggle|mathematical analysis|Team:Paris/Network_analysis_and_design/Core_system/Mathematical_analysis_and_simulations/Mathematical_analysis}} |
- | + | * Does this result depend on the integration methods used ? For more detailled considerations : | |
+ | |||
+ | {{Paris/Toggle|Impact of integration methods|Team:Paris/integration}} | ||
+ | |||
Finally, the logical continuation of the process will be to try to improve the biological system in order to find oscillations. This will be described in the part below, as well as in the section [[Team:Paris/Network_analysis_and_design/System_improvements|System Improvements]]. | Finally, the logical continuation of the process will be to try to improve the biological system in order to find oscillations. This will be described in the part below, as well as in the section [[Team:Paris/Network_analysis_and_design/System_improvements|System Improvements]]. | ||
- | = System analysis = | + | = System interactions analysis = |
Simulations of the core system displayed above reveal that it does not exhibit an oscillatory behavior. In this section we use the model of the core system to try to figure out the contribution of some key characteristics of the network topology on the dynamics of the system. This analysis is done by successively simulating altered forms of the system. | Simulations of the core system displayed above reveal that it does not exhibit an oscillatory behavior. In this section we use the model of the core system to try to figure out the contribution of some key characteristics of the network topology on the dynamics of the system. This analysis is done by successively simulating altered forms of the system. | ||
Line 99: | Line 97: | ||
| style=" text-align:center;" |[[Image:Core5.png|280px|center]] | | style=" text-align:center;" |[[Image:Core5.png|280px|center]] | ||
| style=" text-align:center;" |[[Image:Orgstar.gif]] | | style=" text-align:center;" |[[Image:Orgstar.gif]] | ||
- | | Non oscillating system. A | + | | Non oscillating system. This variant consist in decreasing the threshold value of EnvZ above which flhDC is inhibited. A lower threshold slightly improves the behavior of the system. This value of this parameter is unknown and can thus be used as control to obtain oscillations. |
|- style="background: #dddddd;" | |- style="background: #dddddd;" | ||
| style="background: #D4E2EF; text-align:center;"|[[Image:Core_system4.png|center|150px]]Combining three modifications all together. | | style="background: #D4E2EF; text-align:center;"|[[Image:Core_system4.png|center|150px]]Combining three modifications all together. | ||
Line 108: | Line 106: | ||
- | The core system cannot provide oscillations. | + | The core system cannot provide oscillations. Even with the alternatives examined above it can only show damped oscillations. Moreover, despite its negative effects on oscillations, activation of envZ by FlhDC should not be removed as it is required for the FIFO behavior of the system. |
- | + | ||
- | + | ||
- | + | ||
- | + | ||
+ | In the [[Team:Paris/Network_analysis_and_design/System_improvements/Description|next section]], we examine significant topology improvements of the core system that are potentially better for oscillations and that preserve FIFO required interactions. | ||
- | |||
- | |||
- | + | {{Paris/Navig|Team:Paris/Analysis}} |
Latest revision as of 03:18, 30 October 2008
Mathematical Analysis and Simulations
Other pages:
FIFO
Here is the set of equations used for the numerical simulations, with the parameters found in bibliographical studies (note that you can have an explanation on the way they have been built on the Model Construction page): 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)
Indeed, we may see that there is a LIFO behavior rather than the FIFO we expect... For a qualitative understanding of what is going on with and without FliA : ↓ Qualitative reasoning about the FIFO ↑
Oscillations
If you want more details about the way we have built them, you can go to the model construction page. If one tries to draw simulations of this system, he will obtain this kind of results : In a nutshell : It does not work !
↓ mathematical analysis ↑
↓ Impact of integration methods ↑
System interactions analysisSimulations of the core system displayed above reveal that it does not exhibit an oscillatory behavior. In this section we use the model of the core system to try to figure out the contribution of some key characteristics of the network topology on the dynamics of the system. This analysis is done by successively simulating altered forms of the system. The following table lists the different forms of the system considered and their respective outcome :
In the next section, we examine significant topology improvements of the core system that are potentially better for oscillations and that preserve FIFO required interactions.
|