Team:Paris/Transduction modeling Deterministic and Stochastic Simulations
From 2009.igem.org
(→Stochastic Simulations) |
Christophe.R (Talk | contribs) (→Deterministic and Stochastic Simulations) |
||
(22 intermediate revisions not shown) | |||
Line 5: | Line 5: | ||
- | |||
- | + | ||
- | + | ||
- | + | ==Deterministic and Stochastic Simulations== | |
- | + | ||
<html> | <html> | ||
+ | <style type="text/css"> | ||
+ | #left-side { | ||
+ | position: absolute; | ||
+ | height: 23px; | ||
+ | width: 30px; | ||
+ | top: 0px; | ||
+ | left: 120px; | ||
+ | margin-top:10px; | ||
+ | padding-top: 7px; | ||
+ | background: url(https://static.igem.org/mediawiki/2009/1/1b/Left_menu_pari.png); | ||
+ | z-index:4; | ||
+ | } | ||
+ | |||
+ | #middle-side { | ||
+ | height: 25px; | ||
+ | width: 400px; | ||
+ | position: absolute; | ||
+ | top: 0px; | ||
+ | left: 130px; | ||
+ | margin-top:10px; | ||
+ | padding-top: 5px; | ||
+ | background: #dadada; | ||
+ | z-index:5; | ||
+ | } | ||
+ | |||
+ | #right-side { | ||
+ | position: absolute; | ||
+ | height: 23px; | ||
+ | width: 30px; | ||
+ | margin-top:10px; | ||
+ | padding-top: 7px; | ||
+ | top: 0px; | ||
+ | left: 520px; | ||
+ | background: url(https://static.igem.org/mediawiki/2009/4/40/Right_menu_paris.png); | ||
+ | z-index:4; | ||
+ | } | ||
+ | #left-side2 { | ||
+ | position: absolute; | ||
+ | height: 23px; | ||
+ | width: 30px; | ||
+ | top: 35px; | ||
+ | left: 30px; | ||
+ | margin-top:10px; | ||
+ | padding-top: 7px; | ||
+ | background: url(https://static.igem.org/mediawiki/2009/4/40/Left_menu_paris2.png); | ||
+ | z-index:4; | ||
+ | } | ||
+ | |||
+ | #middle-side2 { | ||
+ | height: 25px; | ||
+ | width: 600px; | ||
+ | position: absolute; | ||
+ | top: 35px; | ||
+ | left: 40px; | ||
+ | margin-top:10px; | ||
+ | padding-top: 5px; | ||
+ | background: white; | ||
+ | z-index:5; | ||
+ | } | ||
+ | |||
+ | #right-side2 { | ||
+ | position: absolute; | ||
+ | height: 23px; | ||
+ | width: 30px; | ||
+ | margin-top:10px; | ||
+ | padding-top: 7px; | ||
+ | top: 35px; | ||
+ | left: 630px; | ||
+ | background: url(https://static.igem.org/mediawiki/2009/f/f8/Right_menu_paris2.png); | ||
+ | z-index:4; | ||
+ | } | ||
+ | |||
+ | a.menu_sub { | ||
+ | padding-left: 7px; | ||
+ | padding-right: 7px; | ||
+ | } | ||
+ | |||
+ | a.menu_sub_active { | ||
+ | padding-left: 7px; | ||
+ | padding-right: 7px; | ||
+ | color:#b0310e; | ||
+ | font-weight:bold; | ||
+ | } | ||
+ | </style> | ||
+ | <div id="left-side"></div> | ||
+ | <div id="middle-side"><center> | ||
+ | <a class="menu_sub"href="https://2009.igem.org/Team:Paris/DryLab#bottom"> Main </a>| | ||
+ | <a class="menu_sub"href="https://2009.igem.org/Team:Paris/Production_modeling#bottom"> Delay model</a>| | ||
+ | <a class="menu_sub"href="https://2009.igem.org/Team:Paris/Production_modeling2#bottom"> Vesicle model</a>| | ||
+ | <a class="menu_sub_active"href="https://2009.igem.org/Team:Paris/Transduction_modeling#bottom"> Fec simulation</a> | ||
+ | </center> | ||
</div> | </div> | ||
- | <div id=" | + | <div id="right-side"></div> |
+ | <div id="left-side2"></div> | ||
+ | <div id="middle-side2"><center> | ||
+ | <a class="menu_sub"href="https://2009.igem.org/Team:Paris/Transduction_modeling#bottom"> Introduction </a>| | ||
+ | <a class="menu_sub"href="https://2009.igem.org/Team:Paris/Transduction_modeling_The Fec Operon as used in our system : chemical equations and kinetics#bottom"> Fec Operon</a>| | ||
+ | <a class="menu_sub_active"href="https://2009.igem.org/Team:Paris/Transduction_modeling_Deterministic and Stochastic Simulations#bottom"> Deterministic and stochastic simulations</a>| | ||
+ | <a class="menu_sub"href="https://2009.igem.org/Team:Paris/Transduction_modeling_Results and discussion#bottom"> Results and conclusion</a> | ||
+ | </center> | ||
</div> | </div> | ||
- | <div id=" | + | <div id="right-side2"></div> |
</html> | </html> | ||
- | |||
- | |||
- | |||
To study our model, we decided to run two type of simulations. | To study our model, we decided to run two type of simulations. | ||
First of all, we have started with deterministic simulation in order to see whether our description of the chemical cascade reflected our understanding and expected behaviour of the reception system. | First of all, we have started with deterministic simulation in order to see whether our description of the chemical cascade reflected our understanding and expected behaviour of the reception system. | ||
- | Then, to check whether our system was | + | Then, to check whether our system was sensitive to molecular noise, we decided to run exact stochastic simulations based on the Gillespie method. |
+ | The whole model was built on the Simbiology toolbox of Matlab ; chemical equations as presented [https://2009.igem.org/Team:Paris/Transduction_modeling_The_Fec_Operon_as_used_in_our_system_:_chemical_equations_and_kinetics#Overview here] were written and a the resulting diagram of interaction between species was ploted : | ||
- | + | [[Image:Schéma sim bio.tiff|500px|center]] | |
- | + | Simulations were performed either in the case of complexation between FecA and FecR or with no complexation between the two species. In each case, different calues for parameter kf_2 were chosen. | |
+ | |||
+ | ===Deterministic Model=== | ||
+ | |||
+ | We used classical algorithms to solve the system of differential equations and get deterministic simulations ; in each case, according to Cauchy Lipschitz, there exist AN UNIQUE solution for a given set of initial conditions. The results of the deterministic simulation are all gathered here : | ||
===Stochastic Simulations=== | ===Stochastic Simulations=== | ||
+ | |||
+ | Several runs of stochastic simulations were run for the different values of k and the different initial amounts of FecA. We used a Gillespie algorithm to compute the amounts of molecules at each step ; the simulations were run with Matlab Simbiology toolbox. | ||
+ | |||
+ | Below, two graphes issued from our simuations : in the first one, the signal is transducted correctly and the receptor is activated while in the second one, the initial amount of FecA is unable to activate the positive retro action loop. | ||
[[Image:Stochastic.png|600px|center]] | [[Image:Stochastic.png|600px|center]] | ||
+ | |||
+ | |||
+ | [[Image:Stochstic inactivation.jpg|600px|center]] | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | <html> | ||
+ | <a href="https://2009.igem.org/Team:Paris/Transduction_modeling_Deterministic_and_Stochastic_Simulations#bottom"><img style="width:40px; height:40px;" src="https://static.igem.org/mediawiki/2009/1/10/Paris_Up.png"/></a> | ||
+ | </html> | ||
+ | |||
+ | |||
+ | {{Template:Paris2009_guided|Transduction_modeling_The_Fec_Operon_as_used_in_our_system_:_chemical_equations_and_kinetics#bottom|Transduction_modeling_Results_and_discussion#bottom}} |
Latest revision as of 20:39, 21 October 2009
iGEM > Paris > Reception > Modeling
Deterministic and Stochastic Simulations
To study our model, we decided to run two type of simulations.
First of all, we have started with deterministic simulation in order to see whether our description of the chemical cascade reflected our understanding and expected behaviour of the reception system.
Then, to check whether our system was sensitive to molecular noise, we decided to run exact stochastic simulations based on the Gillespie method.
The whole model was built on the Simbiology toolbox of Matlab ; chemical equations as presented here were written and a the resulting diagram of interaction between species was ploted :
Simulations were performed either in the case of complexation between FecA and FecR or with no complexation between the two species. In each case, different calues for parameter kf_2 were chosen.
Deterministic Model
We used classical algorithms to solve the system of differential equations and get deterministic simulations ; in each case, according to Cauchy Lipschitz, there exist AN UNIQUE solution for a given set of initial conditions. The results of the deterministic simulation are all gathered here :
Stochastic Simulations
Several runs of stochastic simulations were run for the different values of k and the different initial amounts of FecA. We used a Gillespie algorithm to compute the amounts of molecules at each step ; the simulations were run with Matlab Simbiology toolbox.
Below, two graphes issued from our simuations : in the first one, the signal is transducted correctly and the receptor is activated while in the second one, the initial amount of FecA is unable to activate the positive retro action loop.