Team:TUDelft/Modeling Cascade
From 2009.igem.org
(Difference between revisions)
(→Modeling Negative Cascade) |
|||
Line 1: | Line 1: | ||
{{Template:TUDelftiGEM2009}} | {{Template:TUDelftiGEM2009}} | ||
- | =Modeling | + | =Modeling the Transcriptional Cascade= |
[[Image:Negative_Feedforward.jpg|center| <center> Negative cascade assembly and overview </center> | thumb | 550px]]<br> | [[Image:Negative_Feedforward.jpg|center| <center> Negative cascade assembly and overview </center> | thumb | 550px]]<br> | ||
+ | |||
+ | A full description of the Transcriptional Cascade can be found [https://2009.igem.org/Team:TUDelft/Synthetic_Transcriptional_Cascade here]. | ||
+ | |||
+ | ==ODEs== | ||
+ | |||
+ | The kinetic equations were written out in a Matlab script for both transcription and translation. | ||
+ | |||
[[Image:TUD_eq_cas.png]] | [[Image:TUD_eq_cas.png]] | ||
- | [[Image: | + | |
+ | {| border="1" align="center" | ||
+ | | <b>Symbol</b> || <b>Definition</b> | ||
+ | |- align="left" | ||
+ | | k<sub>IPTGin</sub>, k<sub>IPTGout</sub> || rate constants | ||
+ | |- align="left" | ||
+ | | k<sub>50IPTG</sub>, k<sub>50LacI</sub>, k<sub>50TetR</sub>, k<sub>50CI</sub> || dissociation constants | ||
+ | |- align="left" | ||
+ | | d<sub>mRNA</sub> || mRNA degradation rate | ||
+ | |- align="left" | ||
+ | | d<sub>TetR</sub>, d<sub>CI</sub>, d<sub>RFP</sub>, d<sub>GFP</sub> || protein degradation rates | ||
+ | |- align="left" | ||
+ | | a<sub>pLac</sub>, a<sub>pTet</sub>, a<sub>λp</sub> || transcription leakage (%) | ||
+ | |- align="left" | ||
+ | | c<sub>pLac</sub>, c<sub>pTet</sub>, c<sub>λp</sub> || maximum transcription rates | ||
+ | |- align="left" | ||
+ | | α<sub>1</sub>, α<sub>2</sub>, α<sub>3</sub>, α<sub>4</sub> || translation rates | ||
+ | |- align="left" | ||
+ | | n<sub>IPTG</sub>, n<sub>LacI</sub>, n<sub>TetR</sub>, n<sub>CI</sub> || Hill coefficients | ||
+ | |- align="left" | ||
+ | | [X]<sub>mRNA</sub> || concentration of X mRNA | ||
+ | |- align="left" | ||
+ | |} | ||
+ | |||
+ | |||
+ | ==Sensitivity== | ||
+ | |||
+ | {| border="1" align="center" | ||
+ | | <b>Parameter</b> || <b>Normalized Sensitivity</b> | ||
+ | |- align="left" | ||
+ | | k<sub>IPTGin</sub>, k<sub>IPTGout</sub> || | ||
+ | |- align="left" | ||
+ | | k<sub>50IPTG</sub>, k<sub>50LacI</sub>, k<sub>50TetR</sub>, k<sub>50CI</sub> || | ||
+ | |- align="left" | ||
+ | | d<sub>mRNA</sub> || | ||
+ | |- align="left" | ||
+ | | d<sub>TetR</sub>, d<sub>CI</sub>, d<sub>RFP</sub>, d<sub>GFP</sub> || | ||
+ | |- align="left" | ||
+ | | a<sub>pLac</sub>, a<sub>pTet</sub>, a<sub>λp</sub> || | ||
+ | |- align="left" | ||
+ | | c<sub>pLac</sub>, c<sub>pTet</sub>, c<sub>λp</sub> || | ||
+ | |- align="left" | ||
+ | | α<sub>1</sub>, α<sub>2</sub>, α<sub>3</sub>, α<sub>4</sub> || | ||
+ | |- align="left" | ||
+ | | n<sub>IPTG</sub>, n<sub>LacI</sub>, n<sub>TetR</sub>, n<sub>CI</sub> || | ||
+ | |- align="left" | ||
+ | | [X]<sub>mRNA</sub> || | ||
+ | |- align="left" | ||
+ | |} | ||
+ | |||
+ | ==Parameter Sweeps== | ||
+ | |||
+ | {|border = "0" | ||
+ | |- | ||
+ | |rowspan="2"| | ||
+ | <gallery> | ||
+ | Image:test.jpg|'''1111 :''' parameter vs parameter | ||
+ | |||
+ | Image:test.jpg|'''1111 :''' parameter vs parameter | ||
+ | |||
+ | Image:test.jpg|'''1111 : ''' parameter vs parameter | ||
+ | |||
+ | </gallery><gallery> | ||
+ | |||
+ | Image:test.jpg|'''1111 : ''' parameter vs parameter | ||
+ | |||
+ | Image:test.jpg|'''111 : ''' parameter vs parameter | ||
+ | |||
+ | Image:test.jpg|'''1111 : ''' parameter vs parameter | ||
+ | |||
+ | |||
+ | </gallery> | ||
+ | |} | ||
+ | |||
+ | ==Stability== | ||
+ | |||
+ | Jacobian | ||
+ | |||
+ | {|border = "0" | ||
+ | |- | ||
+ | |rowspan="2"| | ||
+ | <gallery> | ||
+ | Image:test.jpg|'''1111 :''' parameter vs parameter | ||
+ | |||
+ | Image:test.jpg|'''1111 :''' parameter vs parameter | ||
+ | |||
+ | Image:test.jpg|'''1111 : ''' parameter vs parameter | ||
+ | |||
+ | </gallery><gallery> | ||
+ | |||
+ | Image:test.jpg|'''1111 : ''' parameter vs parameter | ||
+ | |||
+ | Image:test.jpg|'''111 : ''' parameter vs parameter | ||
+ | |||
+ | Image:test.jpg|'''1111 : ''' parameter vs parameter | ||
+ | |||
+ | |||
+ | </gallery> | ||
+ | |} | ||
+ | |||
+ | ==Design Recommendations== | ||
+ | |||
+ | Based on the results of the simulations, a series of recommendations were given to the delay team to aid them in choosing parts which would maximize the delay time. | ||
+ | |||
+ | # Significant transcription leakages greatly shorten the delay time. Attempt to minimize leakages. Leakage of λp is a far bigger problem than pTet leakage. | ||
+ | # Use a weak promoter and a weak RBS on the last stage (λp) of the cascade. | ||
+ | # A weak pLac promoters is favorable. | ||
+ | # A strong pTet promoter is favorable. | ||
+ | # A strong RBS on CI gene is favorable. | ||
+ | # A weak RBS on TetR gene is favorable. | ||
+ | # A weak RBS on the endonuclease is favorable although a strong RBS can be used for the GFP gene. | ||
+ | # When choosing RBS and promoter strengths avoid the red areas on the stability plots. | ||
+ | |||
{{Template:TUDelftiGEM2009_end}} | {{Template:TUDelftiGEM2009_end}} |
Revision as of 16:43, 29 September 2009
Modeling the Transcriptional Cascade
A full description of the Transcriptional Cascade can be found here.
ODEs
The kinetic equations were written out in a Matlab script for both transcription and translation.
Symbol | Definition |
kIPTGin, kIPTGout | rate constants |
k50IPTG, k50LacI, k50TetR, k50CI | dissociation constants |
dmRNA | mRNA degradation rate |
dTetR, dCI, dRFP, dGFP | protein degradation rates |
apLac, apTet, aλp | transcription leakage (%) |
cpLac, cpTet, cλp | maximum transcription rates |
α1, α2, α3, α4 | translation rates |
nIPTG, nLacI, nTetR, nCI | Hill coefficients |
[X]mRNA | concentration of X mRNA |
Sensitivity
Parameter | Normalized Sensitivity |
kIPTGin, kIPTGout | |
k50IPTG, k50LacI, k50TetR, k50CI | |
dmRNA | |
dTetR, dCI, dRFP, dGFP | |
apLac, apTet, aλp | |
cpLac, cpTet, cλp | |
α1, α2, α3, α4 | |
nIPTG, nLacI, nTetR, nCI | |
[X]mRNA |
Parameter Sweeps
|
Stability
Jacobian
|
Design Recommendations
Based on the results of the simulations, a series of recommendations were given to the delay team to aid them in choosing parts which would maximize the delay time.
- Significant transcription leakages greatly shorten the delay time. Attempt to minimize leakages. Leakage of λp is a far bigger problem than pTet leakage.
- Use a weak promoter and a weak RBS on the last stage (λp) of the cascade.
- A weak pLac promoters is favorable.
- A strong pTet promoter is favorable.
- A strong RBS on CI gene is favorable.
- A weak RBS on TetR gene is favorable.
- A weak RBS on the endonuclease is favorable although a strong RBS can be used for the GFP gene.
- When choosing RBS and promoter strengths avoid the red areas on the stability plots.