Team:Peking/Modeling/Ring/Simulation
From 2012.igem.org
Line 10: | Line 10: | ||
<h3 id="title1">ODE Model</h3> | <h3 id="title1">ODE Model</h3> | ||
<p> | <p> | ||
- | According to the previous circuit and ODE model, we listed all the differential equations <!--(<a href="/Team:Peking/Modeling/Appendix/ODE">detail here</a>)--> and simulated this system | + | According to the previous circuit and ODE model, we listed all the differential equations <!--(<a href="/Team:Peking/Modeling/Appendix/ODE">detail here</a>)--> and simulated this system in MATLAB with equations listed as below: |
</p> | </p> | ||
<div class="floatC"> | <div class="floatC"> | ||
Line 20: | Line 20: | ||
</div> | </div> | ||
<p> | <p> | ||
- | + | We applied the inverse square law to describe the light intensity distribution on the plate according to different radius, with a central intensity <i>I<sub>0</sub></i> in a region of <i>r=1mm</i>. Here, parameters are: | |
</p> | </p> | ||
<div class="floatC"> | <div class="floatC"> | ||
Line 46: | Line 46: | ||
</tr><tr> | </tr><tr> | ||
<td>r<sub>L</sub></td><td>2.31x10<sup>-2</sup></td><td>min<sup>-1</sup></td><td>LacI and LacIM1 dissociation rate constant</td><td><a href="#ref1" title="Subhayu Basu et al.(2005), A synthetic multicellular system for programmed pattern formation. <i>Nature</i>, vol.434: 1130: 1134">[1]</a></td> | <td>r<sub>L</sub></td><td>2.31x10<sup>-2</sup></td><td>min<sup>-1</sup></td><td>LacI and LacIM1 dissociation rate constant</td><td><a href="#ref1" title="Subhayu Basu et al.(2005), A synthetic multicellular system for programmed pattern formation. <i>Nature</i>, vol.434: 1130: 1134">[1]</a></td> | ||
- | |||
- | |||
</tr><tr> | </tr><tr> | ||
<td>I<sub>0</sub></td><td>1000</td><td>AU</td><td>Maximum light intensity in the middle of the plate </td><td></td> | <td>I<sub>0</sub></td><td>1000</td><td>AU</td><td>Maximum light intensity in the middle of the plate </td><td></td> | ||
Line 58: | Line 56: | ||
</div> | </div> | ||
<p> | <p> | ||
- | The simulation | + | The simulation results in static state are shown below: |
</p> | </p> | ||
<div class="floatC"> | <div class="floatC"> |
Revision as of 20:59, 24 October 2012
ODE Model
According to the previous circuit and ODE model, we listed all the differential equations and simulated this system in MATLAB with equations listed as below:
We applied the inverse square law to describe the light intensity distribution on the plate according to different radius, with a central intensity I_{0} in a region of r=1mm. Here, parameters are:
Parameter | Value | Unit | Description | Source |
a_{G} | 2 | 10^{-6}M/min | GFP production rate constant | [1] |
a_{C} | 2 | 10^{-6}M/min | CI production rate constant | [1] |
a_{L1} | 1 | 10^{-6}M/min | LacI production rate constant | [1] |
a_{L2} | 1 | 10^{-6}M/min | LacIM1 production rate constant | [1] |
b_{C} | 8.x10^{-3} | 10^{-6}M | Binding strength of CI on LacI operator | [1] |
b_{L} | 8.x10^{-1} | 10^{-6}M | Binding strength of LacI or LacIM1 on GFP operator | [1] |
b_{R} | 1.x10^{-2} | 10^{-6}M | Binding strength of Luminesensor on corresponding operator | |
r_{G} | 6.92x10^{-2} | min^{-1} | GFP dissociation rate constant | [1] |
r_{C} | 6.92x10^{-2} | min^{-1} | CI dissociation rate constant | [1] |
r_{L} | 2.31x10^{-2} | min^{-1} | LacI and LacIM1 dissociation rate constant | [1] |
I_{0} | 1000 | AU | Maximum light intensity in the middle of the plate | |
k | 500 | 10^{-6}M | Luminesensor activation rate under light | |
K | 10000 | AU | light sensitivity of Luminesensor activation |
The simulation results in static state are shown below:
Figure 1. ODE Simulation in a plate of the ring-like pattern formation.
Figure 2. ODE Simulation for the radial expression amplitude of the ring-like pattern formation.
From the Figure 1 & 2 above, we discover that, with wildtype parameters, ring-like pattern is formed based on sender-receiver communication through bio-luminescence.
Parameter Analysis
After modeling the ring-like pattern formation with wildtype parameters, we attempted to optimize it in a rational way. We have tuned the parameters both up and down, one by one, and finally discovered five parameters which predominantly influence the expression intensity, ring radius, and band width of pattern formation.
Function | Parameter | Description | Remark |
Reduce responsing time | k_{1} | Vivid lighting decay rate constant | Mainly on process from Light to Dark |
k_{3} | rate constant of monomer LexA releasing from specific binding site | ||
Enhance contrast | K_{2} | Vivid association equilibrium constant | More dimerization provides more binding opportunity |
K_{5} | dimered LexA binding equilibrium constant | More binding affinity |
Reference
- 1. Subhayu Basu et al.(2005), A synthetic multicellular system for programmed pattern formation. Nature, vol.434: 1130: 1134