Team:ZJU-China/Modeling
From 2014.igem.org
(One intermediate revision not shown) | |||
Line 82: | Line 82: | ||
<p><a href="https://2014.igem.org/Team:ZJU-China/SSR">background link</a></p> | <p><a href="https://2014.igem.org/Team:ZJU-China/SSR">background link</a></p> | ||
<p><b>description: </b></p> | <p><b>description: </b></p> | ||
- | <p>In this part, what we want to do is to find out the probability of the recombination of gene of interest through simple molecular dynamics simulation. Although this simulation is quite simple, it | + | <p>In this part, what we want to do is to find out the probability of the recombination of gene of interest through simple molecular dynamics simulation. Although this simulation is quite simple, it certainly can tell us something right in some aspects within a certain accuracy.</p> |
The most important things for simulation are initial conditions and boundary conditions. Next, I will describe the initial conditions and boundary conditions in detail.</p> | The most important things for simulation are initial conditions and boundary conditions. Next, I will describe the initial conditions and boundary conditions in detail.</p> | ||
<p><b>Initial conditions:</b> </p> | <p><b>Initial conditions:</b> </p> | ||
Line 189: | Line 189: | ||
<tr> | <tr> | ||
<td><img src="https://static.igem.org/mediawiki/2014/4/47/ZJU_left_arow.png"> </img></td><td> <a href="https://2014.igem.org/Team:ZJU-China/Solution">Previous: Solution</a></td> | <td><img src="https://static.igem.org/mediawiki/2014/4/47/ZJU_left_arow.png"> </img></td><td> <a href="https://2014.igem.org/Team:ZJU-China/Solution">Previous: Solution</a></td> | ||
- | <td width= | + | <td width=700px></td> |
<td><a href="https://2014.igem.org/Team:ZJU-China/Results">Next: Results</a> </td><td><img src="https://static.igem.org/mediawiki/2014/1/19/ZJU_right_arow.png" > </img> </td> | <td><a href="https://2014.igem.org/Team:ZJU-China/Results">Next: Results</a> </td><td><img src="https://static.igem.org/mediawiki/2014/1/19/ZJU_right_arow.png" > </img> </td> | ||
</tr> | </tr> |
Latest revision as of 03:33, 18 October 2014
1. The Whole Genetic Pathways
ODE equations:
Before recombination:
Formular.1 |
After combination, if combination succeeds.
Formular.2 |
Formular.3 |
After putting in Ara
Formular.4 |
Formulary:
Take GFP for example
Name | description |
---|---|
mgfp | The number of GFP mRNA |
pgfp | The number of GFP protein |
Npla | The number of plasmid |
αgfp | The maximal transcription rate of GFP |
α0gfp | The leak of the promoter |
αmgfp | The degradation rate of mRNA |
βmgfp | The translate rate of mRNA |
βpgfp | The degradation rate of GFP protein |
2. Recombination
description:
In this part, what we want to do is to find out the probability of the recombination of gene of interest through simple molecular dynamics simulation. Although this simulation is quite simple, it certainly can tell us something right in some aspects within a certain accuracy.
The most important things for simulation are initial conditions and boundary conditions. Next, I will describe the initial conditions and boundary conditions in detail.Initial conditions:
What is initial condition? Simply, initial condition is the condition when your simulation starts. More simply, initial condition is that you know every molecular coordinate as well as velocity if needs.Boundary conditions:
What is boundary condition? E coli has a boundary, when the molecule runs out of its boundary, we should adjust it back in the E coli. In this simulation, periodic boundary condition is used.
Some basic biology facts and simulation parameter choice:
figure.1 E.coli cell |
- As shown above, the shape of E coli is similar to a cylinder. So in our simulation, we regard E coli as a cylinder whose radius is 0.5 micrometer, height is 2 micrometer.
- 2.By looking up some online information, we find the average velocity of protein in cells is about 10 , we estimate the average velocity of gene of interest fragment is the same order of magnitude of the protein for their mass is the same order of magnitude.
- 3.E coli replicate its chromosome in 40 minutes, the proceed rate of replication fork is about 10^5 bp/min. A fragment about 1kb needs 0.6s.
results |
3. Bistable Switch
ODE equtions:
Formular.5 |
Formular.6 |
Formulary:
Name | description |
---|---|
[ ] | [ ] stands for the concentration |
kc | Inversion rate constant |
kdi | dissociation equilibrium constant of int dimer-recombination site complex |
ki | dissociation equilibrium constant of int-int dimer |
kdix | dissociation equilibrium constant of int-xis dimer complex on a recombination site |
αset | The transcription rate of input set |
αreset | The transcription rate of input reset |
αI | The maximal transcription rate of int |
αX | The maximal transcription rate of xis |
γI | The degradation rate of int |
γX | The degradation rate of xis |
kd | The dissociation equilibrium constant |
Parameter non-dimensionalization
We nondimensionalize all concentration and time units,in terms of Ki and Kc-1.Kdi=Kdix=Ki.γi=γx=KiKc.
figure.2The response to set input |
Source code Download:ZJU_Modeling.zip
Reference
[1]Mosberg, J. A., M. J. Lajoie and G. M. Church (2010). "Lambda red recombineering in Escherichia coli occurs through a fully single-stranded intermediate." Genetics 186(3): 791-799.
[2]Bonnet, J., P. Subsoontorn and D. Endy (2012). "Rewritable digital data storage in live cells via engineered control of recombination directionality." Proceedings of the National Academy of Sciences of the United States of America 109(23): 8884-8889.
Previous: Solution | Next: Results |