Line 293: | Line 293: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{1}\]</th> | <th scope="row">\[๐_{1}\]</th> | ||
โ | <td> | + | <td>The transcription rate of\[๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>1.955</td> | <td>1.955</td> | ||
Line 299: | Line 299: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{d1}\]</th> | <th scope="row">\[๐_{d1}\]</th> | ||
โ | <td> | + | <td>The degradation rate of\[๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>0.002</td> | <td>0.002</td> | ||
Line 305: | Line 305: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{2}\]</th> | <th scope="row">\[๐_{2}\]</th> | ||
โ | <td> | + | <td>The transcription rate of\[๐๐
๐๐ด_{Csy4}\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>1.116</td> | <td>1.116</td> | ||
Line 311: | Line 311: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{d2}\]</th> | <th scope="row">\[๐_{d2}\]</th> | ||
โ | <td> | + | <td>The degradation rate of\[๐๐
๐๐ด_{Csy4}\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>0.241</td> | <td>0.241</td> | ||
Line 317: | Line 317: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{p2}\]</th> | <th scope="row">\[๐_{p2}\]</th> | ||
โ | <td> | + | <td>The translation rate of\[Protein_{๐ถ๐ ๐ฆ4}\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>1.134</td> | <td>1.134</td> | ||
Line 323: | Line 323: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{dp2}\]</th> | <th scope="row">\[๐_{dp2}\]</th> | ||
โ | <td> | + | <td>The degradation rate of\[Protein_{๐ถ๐ ๐ฆ4}\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>6.547</td> | <td>6.547</td> | ||
Line 329: | Line 329: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{on}\]</th> | <th scope="row">\[๐_{on}\]</th> | ||
โ | <td> | + | <td>The binding constant of \[Protein_{๐ถ๐ ๐ฆ4}\] and \[๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}\]</td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>23995.469</td> | <td>23995.469</td> | ||
Line 335: | Line 335: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{off}\]</th> | <th scope="row">\[๐_{off}\]</th> | ||
โ | <td> | + | <td>The dissociation constant of \[Protein_{๐ถ๐ ๐ฆ4}\] and \[๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}\]</td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>2.703</td> | <td>2.703</td> | ||
Line 341: | Line 341: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{dc1}\]</th> | <th scope="row">\[๐_{dc1}\]</th> | ||
โ | <td> | + | <td>The degradation rate of\[Protein_{๐ถ๐ ๐ฆ4}-crRNA-RBS-mRNA_{๐๐๐}\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>0.024</td> | <td>0.024</td> | ||
Line 347: | Line 347: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{obs}\]</th> | <th scope="row">\[๐_{obs}\]</th> | ||
โ | <td> | + | <td>The cleavage rate of \[Protein_{๐ถ๐ ๐ฆ4}\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>0.327</td> | <td>0.327</td> | ||
Line 353: | Line 353: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{d3}\]</th> | <th scope="row">\[๐_{d3}\]</th> | ||
โ | <td> | + | <td>The degradation rate of \[๐๐
๐๐ด_{gfp}\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>0.472</td> | <td>0.472</td> | ||
Line 359: | Line 359: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{dc2}\]</th> | <th scope="row">\[๐_{dc2}\]</th> | ||
โ | <td> | + | <td>The degradation rate of \[Protein_{๐ถ๐ ๐ฆ4}-crRNA\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>0.024</td> | <td>0.024</td> | ||
Line 365: | Line 365: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{p1}\]</th> | <th scope="row">\[๐_{p1}\]</th> | ||
โ | <td> | + | <td>The translation rate of \[Protein_{gfp}\]</td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>1.711</td> | <td>1.711</td> | ||
Line 371: | Line 371: | ||
<tr> | <tr> | ||
<th scope="row">\[๐_{dp1}\]</th> | <th scope="row">\[๐_{dp1}\]</th> | ||
โ | <td> | + | <td>The degradation rate of \[Protein_{gfp}\] </td> |
<td>\[h^{-1}\]</td> | <td>\[h^{-1}\]</td> | ||
<td>0.479</td> | <td>0.479</td> |
Revision as of 22:51, 11 October 2018
Design
Reactions
We can describe our miniToe system to be followings:
\[\rightarrow [๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}]\]
\[\rightarrow [๐๐
๐๐ด_{Csy4}]\]
\[[๐๐
๐๐ด_{๐ถ๐ ๐ฆ4}]\rightarrow [๐๐
๐๐ด_{๐ถ๐ ๐ฆ4}] + [๐๐๐๐ก๐๐๐_{๐ถ๐ ๐ฆ4}]\]
\[[Protein_{๐ถ๐ ๐ฆ4}]+[crRNA-RBS-mRNA_{gfp}]\leftrightarrow [๐๐๐๐ก๐๐๐_{๐ถ๐ ๐ฆ4}-crRNA-RBS-๐๐
๐๐ด_{๐๐๐}]\]
\[[Protein_{๐ถ๐ ๐ฆ4}-crRNA-RBS-mRNA_{gfp}]\rightarrow [๐๐
๐๐ด_{gfp}] + [๐๐๐๐ก๐๐๐_{๐ถ๐ ๐ฆ4}-๐๐๐
๐๐ด]\]
\[[๐๐
๐๐ด_{gfp}]\rightarrow [๐๐
๐๐ด_{gfp}] + [๐๐๐๐ก๐๐๐_{gfp}]\]
\[[crRNA-RBS-mRNA_{gfp}]\rightarrow \emptyset\]
\[[๐๐
๐๐ด_{๐ถ๐ ๐ฆ4}]\rightarrow \emptyset\]
\[[๐๐๐๐ก๐๐๐_{๐ถ๐ ๐ฆ4}]\rightarrow \emptyset\]
\[[Protein_{๐ถ๐ ๐ฆ4}-crRNA-RBS-mRNA_{gfp}]\rightarrow \emptyset\]
\[[๐๐
๐๐ด_{gfp}]\rightarrow \emptyset\]
\[[๐๐๐๐ก๐๐๐_{๐ถ๐ ๐ฆ4}-๐๐๐
๐๐ด]\rightarrow \emptyset\]
\[[๐๐๐๐ก๐๐๐_{gfp}]\rightarrow \emptyset\]
Two equations, describing the functional binding and cleavage of Csy4 protein in biology, and three parameters kon koff kobs describing the same things in mathematics, are the core of our model.
ODEs
To simulate the dynamics of GFP, we use ordinary differential equations to model the reactions above. And ODEs are given as follows:
\[\frac{d[๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}]}{dt}=๐_{1}-๐_{d1}[๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}]\]
\[-๐_{on}[Protein_{๐ถ๐ ๐ฆ4}][๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}]\]
\[+๐_{off}[Protein_{๐ถ๐ ๐ฆ4}-๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}]\]
\[\frac{d[๐๐
๐๐ด_{๐ถ๐ ๐ฆ4}]}{dt}=๐_{2}-๐_{d2}[๐๐
๐๐ด_{๐ถ๐ ๐ฆ4}]\]
\[\frac{d[Protein_{๐ถ๐ ๐ฆ4}]}{dt}=๐_{p2}[๐๐
๐๐ด_{๐ถ๐ ๐ฆ4}]-๐_{dp2}[Protein_{๐ถ๐ ๐ฆ4}]\]
\[-๐_{on}[Protein_{๐ถ๐ ๐ฆ4}][๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}]\]
\[+๐_{off}[Protein_{๐ถ๐ ๐ฆ4}-๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}]\]
\[\frac{d[๐๐๐๐ก๐๐๐_{๐ถ๐ ๐ฆ4}-crRNA-RBS-๐๐
๐๐ด_{๐๐๐}]}{dt}=๐_{on}[Protein_{๐ถ๐ ๐ฆ4}][๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}]\]
\[-๐_{of}[Protein_{๐ถ๐ ๐ฆ4}][Protein_{๐ถ๐ ๐ฆ4}-crRNA-RBS-mRNA_{gfp}]\]
\[-๐_{d1}[Protein_{๐ถ๐ ๐ฆ4}][Protein_{๐ถ๐ ๐ฆ4}-crRNA-RBS-mRNA_{gfp}]\]
\[-๐_{obs}[Protein_{๐ถ๐ ๐ฆ4}][Protein_{๐ถ๐ ๐ฆ4}-crRNA-RBS-mRNA_{gfp}]\]
\[\frac{d[๐๐
๐๐ด_{๐๐๐}]}{dt}=๐_{obs}[๐๐๐๐ก๐๐๐_{๐ถ๐ ๐ฆ4}-crRNA-RBS-๐๐
๐๐ด_{๐๐๐}]-๐_{d3}[๐๐
๐๐ด_{gfp}]\]
\[\frac{d[๐๐๐๐ก๐๐๐_{๐ถ๐ ๐ฆ4}-crRNA]}{dt}=๐_{obs}[๐๐๐๐ก๐๐๐_{๐ถ๐ ๐ฆ4}-crRNA-RBS-๐๐
๐๐ด_{๐๐๐}]-๐_{dc2}[๐๐๐๐ก๐๐๐_{๐ถ๐ ๐ฆ4}-๐๐๐
๐๐ด]\]
\[-๐_{on}[Protein_{๐ถ๐ ๐ฆ4}][๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}]\]
\[+๐_{off}[Protein_{๐ถ๐ ๐ฆ4}-๐๐๐
๐๐ด โ ๐
๐ต๐ โ ๐๐
๐๐ด_{๐๐๐}]\]
\[\frac{d[Protein_{gfp}]}{dt}=๐_{p1}[๐๐
๐๐ด_{๐๐๐}]-๐_{dp1}[Protein_{gfp}]\]
For the readability, the complex symbol is simplified as:
\[\frac{d[A]}{dt}=๐_{1}-๐_{d1}[A]-๐_{on}[C][A]+๐_{off}[D]\]
\[\frac{d[B]}{dt}=๐_{2}-๐_{d2}[B]\]
\[\frac{d[C]}{dt}=๐_{p2}[B]-๐_{dp2}[C]-๐_{on}[C][A]+๐_{off}[D]\]
\[\frac{d[D]}{dt}=๐_{on}[C][A]-๐_{dp2}[C]-๐_{off}[D]-๐_{dc1}[D]-๐_{obs}[D]\]
\[\frac{d[E]}{dt}=๐_{obs}[D]-๐_{d3}[E]\]
\[\frac{d[F]}{dt}=๐_{obs}[D]-๐_{d3}[F]\]
\[\frac{d[G]}{dt}=๐_{p1}[E]-๐_{dp1}[G]\]
Data Processing
The leak in the experiment is an big problem in estimating parameters in our ODEs model, so we processing the data by following formula๏ผ \[Data(without leak)=Data(+IPTG)-Data(-IPTG)\]
By doing this, we can reduce some factor which may be influence estimation not just the leak, but also some background noise. So we can get more precise parameters of the Csy4.
Species, symbols and parameters
Species | Symbol | Initial value | Units |
---|---|---|---|
\[๐๐๐ ๐๐ด โ ๐ ๐ต๐ โ ๐๐ ๐๐ด_{๐๐๐}\] | A | 15 | \[mol/L\] |
\[๐๐ ๐๐ด_{Csy4}\] | B | 0 | \[mol/L\] |
\[Protein_{๐ถ๐ ๐ฆ4}\] | C | 0 | \[mol/L\] |
\[Protein_{๐ถ๐ ๐ฆ4}-crRNA-RBS-mRNA_{๐๐๐}\] | D | 0 | \[mol/L\] |
\[๐๐ ๐๐ด_{gfp}\] | E | 0 | \[mol/L\] |
\[Protein_{๐ถ๐ ๐ฆ4}-๐๐๐ ๐๐ด\] | F | 0 | \[mol/L\] |
\[Protein_{gfp}\] | G | 0 | \[mol/L\] |
Because the [๐๐๐ ๐๐ด โ ๐ ๐ต๐ โ ๐๐ ๐๐ด๐๐๐] is under controlled by a constitutive promoter, so we set the initial concentration to 15mol/L .
The other parament we used in the ODEs is listed in the following table:
Parameters | Definition | Units | Value |
---|---|---|---|
\[๐_{1}\] | The transcription rate of\[๐๐๐ ๐๐ด โ ๐ ๐ต๐ โ ๐๐ ๐๐ด_{๐๐๐}\] | \[h^{-1}\] | 1.955 |
\[๐_{d1}\] | The degradation rate of\[๐๐๐ ๐๐ด โ ๐ ๐ต๐ โ ๐๐ ๐๐ด_{๐๐๐}\] | \[h^{-1}\] | 0.002 |
\[๐_{2}\] | The transcription rate of\[๐๐ ๐๐ด_{Csy4}\] | \[h^{-1}\] | 1.116 |
\[๐_{d2}\] | The degradation rate of\[๐๐ ๐๐ด_{Csy4}\] | \[h^{-1}\] | 0.241 |
\[๐_{p2}\] | The translation rate of\[Protein_{๐ถ๐ ๐ฆ4}\] | \[h^{-1}\] | 1.134 |
\[๐_{dp2}\] | The degradation rate of\[Protein_{๐ถ๐ ๐ฆ4}\] | \[h^{-1}\] | 6.547 |
\[๐_{on}\] | The binding constant of \[Protein_{๐ถ๐ ๐ฆ4}\] and \[๐๐๐ ๐๐ด โ ๐ ๐ต๐ โ ๐๐ ๐๐ด_{๐๐๐}\] | \[h^{-1}\] | 23995.469 |
\[๐_{off}\] | The dissociation constant of \[Protein_{๐ถ๐ ๐ฆ4}\] and \[๐๐๐ ๐๐ด โ ๐ ๐ต๐ โ ๐๐ ๐๐ด_{๐๐๐}\] | \[h^{-1}\] | 2.703 |
\[๐_{dc1}\] | The degradation rate of\[Protein_{๐ถ๐ ๐ฆ4}-crRNA-RBS-mRNA_{๐๐๐}\] | \[h^{-1}\] | 0.024 |
\[๐_{obs}\] | The cleavage rate of \[Protein_{๐ถ๐ ๐ฆ4}\] | \[h^{-1}\] | 0.327 |
\[๐_{d3}\] | The degradation rate of \[๐๐ ๐๐ด_{gfp}\] | \[h^{-1}\] | 0.472 |
\[๐_{dc2}\] | The degradation rate of \[Protein_{๐ถ๐ ๐ฆ4}-crRNA\] | \[h^{-1}\] | 0.024 |
\[๐_{p1}\] | The translation rate of \[Protein_{gfp}\] | \[h^{-1}\] | 1.711 |
\[๐_{dp1}\] | The degradation rate of \[Protein_{gfp}\] | \[h^{-1}\] | 0.479 |
Simulation
With the help of Simbiology toolbox in Matlab๏ผwe simulate our miniToe system for 30h, the result can be seen in the Fig.1.
Fig.1 The dynamics of GFP by model prediction
Fig.2 The comparison between experimental data and simulation data
Discussion
Combining the biology and math, we now discuss why the dynamics of GFP is like the curve in the Fig.1. In order to explain in detail, we plot the dynamics of all species in the miniToe system in Fig.3.
Fig.3 The dynamics of all species in the miniToe system
Sensitivities Analysis
After building the ODE model, we try to do something more deeply to our miniToe system by analyze the sensitivity of parameters. Fig.4 shows the numerical integration of sensitivities of parameters in 30 hours.
Fig.4 The numerical integration of sensitivities of parameters in 30h
Enlighten by the sensitivity analysis, we give a prediction curve that shows that what will happen in the GFP expression curve if we change the Csy4, and it can be seen in the Fig.5.
Fig.5 The curve of GFP with the changing cleavage rate
Fig.6 The relationship between the stable expression level of GFP and the claevage rate
Contact Us : oucigem@163.com | ©2018 OUC IGEM.All Rights Reserved. | โฆโฆโฆโฆ