Team:Paris/Transduction modeling Deterministic and Stochastic Simulations

From 2009.igem.org

Revision as of 19:19, 15 October 2009 by CChabbert (Talk | contribs)

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 sensible 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 :

Schéma sim bio.tiff


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 ;

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.

Stochastic.png


Stochstic inactivation.jpg