Team:Peking/Modeling/Luminesensor/Simulation

From 2012.igem.org

Revision as of 02:05, 23 September 2012 by Niasw (Talk | contribs)

ODE model

According to the previous network and ODE model, we listed all the differential equations (detail here) and simulated this system with MATLAB. We specifically watched three expressions to understand the mechanism for Luminesensor.

Expression Description Remark
rb = 1 - [DL]/[DT] This indicates the repressing degree. DT indicates the total specific binding sites, while the DL indicates the free ones among DT.
rd = 2[LA2X]/[LT] This indicates the dimerizing degree. LT indicates the total Luminesensor molecules, and LA2X indicates all dimered Luminesensor molecules, i.e. LA2 + LA2DL.
ra =
( [LAX] + 2[LA2X] ) /[LT]
This indicates the activating degree. LAX indicates all monomer Luminesensor molecules, i.e. LA + LADL.

The simulation result is shown below:

Simulation Result

Fig 2. ODE Simulation Result of the prototype Luminesensor.

From the figure above, we discovered that the activation and decay of Luminesensor are the pioneers of progress, and the activating rate is the most completely switched variable as lighting varies. The promoter sequences in the DNA are repressed even though the Luminesensor has not all dimered.

Stochastic Simulation

In order to check the working stability of Luminesensor, we simulated this reaction network with a stochastic model. By estimating the volume of a cell, we converted the concentration of a component into the number of molecules by 1 n mol/L : 1. The result are shown below:

Simulation Result

Fig 3. Stochastic Simulation Result of the prototype Luminesensor.

According to the figure above, the noise did not influence this system. Thus, the Luminesensor is expected to work theoretically. Besides, the average value of stochastic simulation is coupled with the result of ODE model, which in turn proves the self-consistency of our ODE model.

  • Totop Totop