Team:WHU-China/Project/Modeling

  This is our gene circuit designed to equip our E.coli with the ability to change color with time.
  The key players are LsrR and LsrK. The alternation of their expression results in the oscillating concentration of the reporter—GFP, and thus yield the periodically changing color.
  Theoretically, at the very beginning, LsrR is expressed and inhibits LsrK. After some time, LsrR inhibits itself and its concentration drops. Therefore LsrK has a chance to express and climbs due to its self enlarging. However, LsrK induces the reexpression of LsrR at the same time. In this way, their concentrations wave alternatively. C1’s function is to delay the oscillation and stabilize the system.

  In the part of waving time, there is a two nodes system .Using the method of the mathematic theory of the ODE and the stochastic process.,we get our models.In the models, X1 denotes LsrK, X2 denotes LsrR.
  First we generate a ODE to describe the system:

  the solution diagram of the function is:

  Taking the fluctuations arising from the molecules into account, we need to resort to stochastic simulations. We use here the Gillespie algorithm to simulate a stochastic version of the model.Suppose the system size Ω=1000,then we can get the result of simulatio.

  The hofp Bifurcation diagram of the parameter a1 a2 a3: