Team:USTC-China/Project/modeling
From 2011.igem.org
(→Results) |
(→Results) |
||
(14 intermediate revisions not shown) | |||
Line 85: | Line 85: | ||
[[File:Z2.png|320px|margin-left:4px]] | [[File:Z2.png|320px|margin-left:4px]] | ||
[[File:S1.png|320px]] | [[File:S1.png|320px]] | ||
+ | <p> '''Fig. 1'''</p> | ||
<p>Following pictures show the switch's performances under different max transcription rate of ci mRNA (t = 43200). 1 means 0.1*k1, 2 means 0.2*k1, and 3 means 0.3*k1</p> | <p>Following pictures show the switch's performances under different max transcription rate of ci mRNA (t = 43200). 1 means 0.1*k1, 2 means 0.2*k1, and 3 means 0.3*k1</p> | ||
Line 90: | Line 91: | ||
[[File:Delta2.png|250px]] | [[File:Delta2.png|250px]] | ||
[[File:Delta3.png|250px]] | [[File:Delta3.png|250px]] | ||
- | + | <p> '''Fig. 2'''</p> | |
<p> | <p> | ||
Following pictures show how cheZ concentrations change on time under conditions of different AHL and theophylline concentrations (t = 43200).</p> | Following pictures show how cheZ concentrations change on time under conditions of different AHL and theophylline concentrations (t = 43200).</p> | ||
Line 96: | Line 97: | ||
[[File:Z a3 t1000.png|320px]] | [[File:Z a3 t1000.png|320px]] | ||
[[File:Z a10 t1000.png|320px]] | [[File:Z a10 t1000.png|320px]] | ||
+ | <p> '''Fig. 3'''</p> | ||
<p>Following pictures show how a single | <p>Following pictures show how a single | ||
Line 107: | Line 109: | ||
[[File:A3 t1000.png|320px]] | [[File:A3 t1000.png|320px]] | ||
[[File:A10 t1000.png|320px]] | [[File:A10 t1000.png|320px]] | ||
+ | <p> '''Fig. 4'''</p> | ||
<p>Following pictures show how a single | <p>Following pictures show how a single | ||
Line 119: | Line 122: | ||
[[File:Z t306.png|320px]] | [[File:Z t306.png|320px]] | ||
[[File:A10.png|320px]] | [[File:A10.png|320px]] | ||
+ | <p> '''Fig. 5'''</p> | ||
+ | |||
<p>AHL = 3 nM</p> | <p>AHL = 3 nM</p> | ||
| | ||
[[File:Z t226.png|320px]] | [[File:Z t226.png|320px]] | ||
[[File:A3nm.png|320px]] | [[File:A3nm.png|320px]] | ||
+ | <p> '''Fig. 6'''</p> | ||
<p> | <p> | ||
Line 133: | Line 139: | ||
[[File:Try 24.png|320px]] | [[File:Try 24.png|320px]] | ||
[[File:Output2.gif|320px]] | [[File:Output2.gif|320px]] | ||
- | + | <p> '''Fig. 7'''</p> | |
<p> | <p> | ||
Line 164: | Line 170: | ||
[[File:Try 21.png|320px]] | [[File:Try 21.png|320px]] | ||
[[File:Output~8.gif|320px]] | [[File:Output~8.gif|320px]] | ||
+ | <p> '''Fig. 8'''</p> | ||
<p>For theophylline gradient (25,500)</p> | <p>For theophylline gradient (25,500)</p> | ||
Line 169: | Line 176: | ||
[[File:Try 19.png|320px]] | [[File:Try 19.png|320px]] | ||
[[File:Output3.gif|320px]] | [[File:Output3.gif|320px]] | ||
+ | <p> '''Fig. 9'''</p> | ||
+ | <p>When adding intrinsic noises (fluctuations of promoters) into our system, it seems that the result is improved a little.</p> | ||
+ | | ||
+ | [[File:Pos 13.png|320px]] | ||
+ | <p> '''Fig. 10'''</p> | ||
<br/> | <br/> | ||
- | <p>' | + | <p>After several days' work, it seems that we have improved the performance of our system in some aspects.</p> |
- | + | <p>First of all, the negative feedback of regulatory ci is moved and a positive feedback of regulatory ci is added. This design is described as this, aptamer-cheZ is added downstream of ci434-gfp, and lasI is added downstream of ci-rfp and an AHL positively regulated promoter is in use. We can see from the following picture that the AHL switch performance is enhanced.</p> | |
- | <p> | + | <p>Following two pictures shows how steady-state cheZ concentration changes under different values of AHL concentration and delta. The first picture showing below shows the AHL switch performance, delta is 0.3 here. For the second picture, delta means max transcription rate of regulatory ci is set as delta*k1, AHL concentration is 10 nm here. Compared to previous negative feedback (Fig. 1), the switch performance is improved. Simulation time is 43200.</p> |
- | <p>The | + | |
- | <p> | + | [[File:Steady3.png|320px]] |
- | <p>In | + | [[File:Delta1000.png|320px]] |
- | <p> | + | <p> '''Fig. 11'''</p> |
- | < | + | <p>In the situation of following picture, cell number growth rate is constant (0.97), and initial cell number is 20. Noises are included.</p> |
+ | | ||
+ | [[File:Pos 10.png|320px]] | ||
+ | <p> '''Fig. 12'''</p> | ||
+ | |||
+ | <br/> | ||
+ | |||
<html ><a href= "#results "><u>Results</u></a> </html> | <html ><a href= "#results "><u>Results</u></a> </html> | ||
<html><a href= "#odes "> <u>Odes</u> </a> </html> | <html><a href= "#odes "> <u>Odes</u> </a> </html> |
Latest revision as of 10:04, 28 October 2011
![](https://static.igem.org/mediawiki/2011/3/31/%E6%A0%87%E9%A2%98%E6%A0%8F%E6%9C%80%E7%BB%88%E6%8E%A8%E8%8D%90%E6%9C%89%E5%AD%97%E7%89%88.png)
Contents |
Results
We have been working hard to model our system.
Before presenting our results, we would like to say that lasR protein is assumed to be always over production in our system, although it is not in our original design, it really helps a lot to simplify our modeling work.
Suggestions and questions are welcomed, for contact, you can email to PJH_panda(jhpanda@mail.ustc.edu.cn) or any other member of our team.
All codes are written in Python26(http://www.python.org/getit/releases/2.6/), you can download our codes here File:IGEM USTC 2011.zip, you will need Scipy(http://www.scipy.org/Download), Numpy(http://www.scipy.org/Download), and Matplotlib(http://matplotlib.sourceforge.net/) to run our codes.
Now see what we have done
Firstly, we would like to present some basic character of our system. For this part, the genetic design is a Toggle switch with a regulatory ci. We use lasI/lasR repressible promoter here.
In a single cell, the system acts like a fast switch upon AHL concentration, on the other hand, cheZ ouput is a slow switch of theophylline concentration.
We can see that cheZ concentration rises rapidly when AHL concentration rises from 3nM to 4nM. For the switch is so fast that we can use quorum system in our design. CheZ output rises relatively rapid when theophylline concentration is around 210nM.
For the system without quorum part, If AHL concentration is low, regulate ci protein accumulates, the toggle switch is in the ci state, cheZ production is very low, so the cell can rarely move. And when AHL concentration is high, not enough regulate ci protein is produced, the toggle switch is in the ci434 state, cheZ production is high, so the cell can spend less time tumbling, and will move to the direction where cheZ concentration rises.
Following pictures show how steady-state cheZ concentrations change under conditions of different AHL and theophylline concentrations (t = 64800).
Fig. 1
Following pictures show the switch's performances under different max transcription rate of ci mRNA (t = 43200). 1 means 0.1*k1, 2 means 0.2*k1, and 3 means 0.3*k1
Fig. 2
Following pictures show how cheZ concentrations change on time under conditions of different AHL and theophylline concentrations (t = 43200).
Fig. 3
Following pictures show how a single cell moves on time under conditions of different AHL and theophylline concentrations in one single round of simulation (t = 72000).
Fig. 4
Following pictures show how a single cell moves and how cheZ concentration changes by time under conditions of different AHL and gradient of theophylline concentration (30,600) in one single round of simulation (t = 72000).
AHL = 10 nM
Fig. 5
AHL = 3 nM
Fig. 6
Secondly, we want to see the movement of 100 cells (or more) without quorum part under 10nM AHL and theophylline gradient (25,500)
Assume that 100 cells are around the center for 0.25mm distance at most. Following pictures shows the position of the 100 cells without quorum part under 10nM AHL and gradient of theophylline concentration (25,500) after t = 86400. And a little movie is shown
Fig. 7
Thirdly, after the previous two modelling part, we believe that if we add quorum part into our system, the results will be acceptable.
For the whole system with quorum part, in the beginning, we assume that the concentrations of all related macromolecules are zero. Cells begin to divide, cell number is growing, and AHL is accumulating, when AHL concentration is high enough, toggle switch in some of the cells turns to the ci434 state, and these cells is moving out. After some time, AHL concentration in these cells will drop, and the toggle switch will most likely turn to the ci state, cells stop moving.
For this part we assume that initial cell number is 20, and numbers of the cells grow under the rate of k=0.97*exp(-t/5), the unit of time is hour.
In the beginning, max transcription rate of regulatory ci mRNA is set as k1. However, we found cheZ concentration keeps dropping during simulation time. So we set it to be 0.3*k1, this will not change the basic character of bistable switch. From the first part of our modeling, AHL acts as a fast switch between 3nM and 4nM. If we set max transcription rate of regulatory ci mRNA as 0.3*k1, the switch is between 1nM and 2nM.
For no gradient of theophylline
Fig. 8
For theophylline gradient (25,500)
Fig. 9
When adding intrinsic noises (fluctuations of promoters) into our system, it seems that the result is improved a little.
Fig. 10
After several days' work, it seems that we have improved the performance of our system in some aspects.
First of all, the negative feedback of regulatory ci is moved and a positive feedback of regulatory ci is added. This design is described as this, aptamer-cheZ is added downstream of ci434-gfp, and lasI is added downstream of ci-rfp and an AHL positively regulated promoter is in use. We can see from the following picture that the AHL switch performance is enhanced.
Following two pictures shows how steady-state cheZ concentration changes under different values of AHL concentration and delta. The first picture showing below shows the AHL switch performance, delta is 0.3 here. For the second picture, delta means max transcription rate of regulatory ci is set as delta*k1, AHL concentration is 10 nm here. Compared to previous negative feedback (Fig. 1), the switch performance is improved. Simulation time is 43200.
Fig. 11
In the situation of following picture, cell number growth rate is constant (0.97), and initial cell number is 20. Noises are included.
Fig. 12
Results Odes Parameters References
Odes
Odes and equations
C1, regulate ci mRNA concentration
C2, regulate ci protein concentration
C3, bistable ci mRNA concentration
C4, total ci protein concentration
C5, bistable ci434 mRNA concentration
C6, bistable ci434 protein concentration
C7, cheZ mRNA concentration
C8, cheZ protein concentration
C9, I mRNA concentration
C10, I protein concentration
C11, AHL produced
A, AHL concentration
For regulatory ci mRNA and protein
For Bistable Part
For cheZ mRNA and protein
For Quorum Sensing Part
The real AHL concentration is not the AHL produced. Considering about fast diffusivity of AHL, we assume that external AHL concentration is zero, although it is not quite exact, but it can describe the character of AHL well. That is, for cell i, the compact of cell j is,
Tumbling frequency of a cell is determined by its cheYp concentration
CheYp concentration is related to cheZ concentration, and for the reaction is so fast, we assume that current cheZ concentration determines current cheYp concentration
We assume that cells move at the speed of 0.025um/s on our semisolid plate. For the movement of a single cell, if the cell tumbles, it will stay of the current second, and change its direction later, if not, the cell keeps on moving in one direction.
For the gradient of theophylline,
For cell number growth, here, we assume that the growth rate drops by time, the unit of time is hour
Results
Odes
Parameters
References
Parameters
Parameters
Name | Value | Ref. |
k1, max transcription rate of regulatory ci mRNA | 0.0933 nM/s | 1 |
k2, translation rate of regulatory ci protein | 0.0072 /s | 1 |
k3, max transcription rate of bistable ci mRNA | 0.0933 nM/s | 1 |
k4, translation rate of regulatory ci protein | 0.0072 /s | 1 |
k'4, translation rate of bistable ci protein | 0.048 /s | 1 |
k5, max transcription rate of ci434 mRNA | 0.0987 nM/s | 1 |
k6, translation rate of ci434 protein | 0.0845 /s | 1 |
k7, max transcription rate of cheZ mRNA | 0.0834 nM/s | 1, 2 |
k8, translation rate of cheZ protein | 0.1869 /s | 1, 2, 4, 7 |
k9, max transcription rate of lasI mRNA | 0.014 nM/s | 3 |
k10, translation rate of lasI protein | 0.016 /s | 3 |
k11, AHL synthesis rate | 0.06 /s | 3 |
K1, Kd between AHL and Plux promoter | 1.6 nM | 10 |
K3, Kd between ci protein and bistable ci promoter | 40 nM | 1 |
K'3, Kd between ci434 protein and bistable ci promoter | 50 nM | 1 |
K5, Kd between ci protein and bistable ci434 promoter | 40 nM | 1 |
K7, Kd between theophylline and RNA aptamer | 210 nM | 4 |
n1, Hill co-effiency of AHL and Plux promter | 1.6 | 10 |
n3, Hill co-effiency of ci protein and bistable ci promoter | 4 | 1 |
n'3, Hill co-effiency of ci434 protein and bistable ci promoter | 2 | 1 |
n5, Hill co-effiency of ci protein and ci434 promoter | 4 | 1 |
n7, Hill co-effiency of theophylline and RNA aptamer | 3 | 7 |
γ1, Degradation rate of regulatory ci mRNA | 0.00434 /s | 1 |
γ2, Degradation rate of regulatory ci protein | 0.000935 /s | 1 |
γ3, Degradation rate of bistbale ci mRNA | 0.00434 /s | 1 |
γ4, Degradation rate of total ci protein | 0.000935 /s | 1 |
γ5, Degradation rate of ci434 mRNA | 0.00434 /s | 1 |
γ6, Degradation rate of ci434 protein | 0.000935 /s | 1 |
γ7, Degradation rate of cheZ mRNA | 0.00434 /s | 1 |
γ8, Degradation rate of cheZ protein | 0.00434 /s | 1 |
γ9, Degradation rate of lasI mRNA | 0.006 /s | 3 |
γ10, Degradation rate of lasI protein | 0.0001 /s | 3 |
γ11, Degradation rate of AHL | 0.00038 /s | 3 |
diffu_a, diffusion constant of AHL | 10^6 /mm^2 | assumption |
diffu_t, diffusion constant of theophylline | 500, 600 /mm^2 | assumption |
da, diffusivity constant of AHL | 0.23 /s | 3 |
[SetYp], wild type ecoli cheYp concentration | 4.4 uM | 8 |
ny, Hill co-effiency of cheYp and cell tumble rate | 5.5 | 9 |
ky, cheYp phosphorylation rate | 3*10^7 /(M*s) | 8 |
k-y, cheYp dephosphorylation rate | 5*10^5 /(M*s) | 8 |
[Tp], concentration of wild type chemotaxis phosphorylated receptor Tar | 0.038 uM | 8 |
v, cell moving speed | 0.025 mm/s | 2, assumption |
k, cell number growth rate | 0.97 /h | 6 |
Results
Odes
Parameters
References
References
1. the wiki of iGEM09 Peking U team, https://2009.igem.org/Team:PKU_Beijing/Modeling;
2. the E. coli Statistics, http://gchelpdesk.ualberta.ca/CCDB/cgi-bin/STAT_NEW.cgi;
3. Systems analysis of a quorum sensing network: Design constraints imposed by the functional requirements, network topology and kinetic constants, Systems Biology Group, Bioinformatics Institute, 30 Biopolis Str., Singapore 138671, Singapore, doi:10.1016/j.biosystems.2005.04.006;
4. Molecular interactions and metal binding in the theophylline-binding core of an RNA aptamer, G R Zimmermann, C L Wick, T P Shields, R D Jenison, and A Pardi, http://rnajournal.cshlp.org/content/6/5/659.short
5. A Minimal Model of Metabolism-Based Chemotaxis, Matthew D. Egbert, Xabier E. Barandiaran, Ezequiel A. Di Paolo, http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1001004
6. Programmed population control by cell–cell communication and regulated killing, Lingchong You, Robert Sidney Cox, III, Ron Weiss Frances, H. Arnold, http://www.nature.com/nature/journal/v428/n6985/abs/nature02491.html
7. Guiding bacteria with small molecules and RNA, Shana Topp and Justin P. Gallivan. http://pubs.acs.org/doi/abs/10.1021/ja0692480.
8. A model of excitation and adaptation in bacterial chemotaxis, Peter A. Spiro, John S. Parkinson, and Hans G. Othmer. http://www.pnas.org/content/94/14/7263.full.
9. Origins of Individual Swimming Behavior in Bacteria, Matthew D. Levin, http://www.cell.com/biophysj/abstract/S0006-3495(98)77777-X.
10. Conversion of the Vibrio fischeri Transcriptional activator, LuxR, to a repressor: Egland, Kristi A. and Greenberg, E. P. Journal of Bacteriology. Feb 2000, p. 805-811: http://jb.asm.org/cgi/content/abstract/182/3/805.