Team:TU-Delft/Timer Plus Sumo
From 2013.igem.org
Line 240: | Line 240: | ||
<center> | <center> | ||
<img src="http://2013.igem.org/wiki/images/5/5f/Modeling_results_timer.jpg" width="550px" height="432"> | <img src="http://2013.igem.org/wiki/images/5/5f/Modeling_results_timer.jpg" width="550px" height="432"> | ||
- | <p>Figure 2: Simulation Results</p | + | <p>Figure 2: Simulation Results</p> |
</center> | </center> | ||
<br> | <br> |
Revision as of 12:46, 4 October 2013
Timer Plus Sumo
In this section the system of Figure 1 is modeled. The structure of the timer has two repressing promoters (PcI and Ptet) and the input is the T7 promoter and the output is the protease Ulp-1. This Ulp-1 cleaves off the SUMO from the produced SUMO+peptide.
Figure 1: Circuit of the timer including sumo cleaving
By this model we are finding answers to the following questions:
- At what time point does the peptide+SUMO start to be produced?
- After how many minutes does the Ulp-1 cleave the peptide?
Differential Equations
The above circuit can be represented by the following differential equations. We assume a binary behavior of the T7 promoter. In the presence of IPTG, the T7 promoter will be active. So, we make the assumption that the T7 is binary variable with two possible states: either active 1 or inactive 0.
Parameters
For the stated parameters of the differential equations many values can be found in the literature. The degradation rates are relatively unknown, while the transcription rates vary widely. Many of these transcription rates resulted in our model in infeasible results, e.g. concentrations of millions of molecules per cell. For the PT7 a relative low transcription rate was thus selected in order to get the values in a reasonable range. The other promoters in the model we found as a fraction of the transcription rate of PT7 in the literature.
Parameter | Value | Description | Units | Reference |
c_{a} | 1020 | Translation rate per amino acid | min^{-1}#_{a}^{-1} | [7] |
c_{T7} | 4.16 | Maximum transcription rate of T7 | #_{m}/min | [2] |
c_{ptet} | 2.79 | Maximum transcription rate of Ptet | #_{m}/min | [4] |
c_{ci} | 1.79 | Maximum transcription rate of Pci | #_{m}/min | [3] |
d_{mRNA} | 0.231 | Degradation rate of mRNA | min^{-1} | [8] |
d_{TET} | 0.1386 | Degradation rate of TET | min^{-1} | [9] |
d_{CI} | 0.042 | Degradation rate of CI | min^{-1} | [9] |
d_{PEP} | 2.1*10^{-3} | Degradation rate of the peptide | min^{-1} | Assumed three times slower same as GFP |
d_{PSU} | 6.3*10^{-3} | Degradation rate of the peptide plus SUMO | min^{-1} | Assumed the same as GFP |
d_{Ulp} | 1.263*10^{-2} | Degradation rate of Ulp | min^{-1} | Assumed twice the rate of GFP |
l_{t7} | 0.002 | Leakage factor of T7 | - | Assumption |
l_{ptet} | 0.002 | Leakage factor of Ptet | - | Assumption |
l_{ci} | 0.002 | Leakage factor of Pci | - | Assumption |
k_{tet} | 6 | Dissociation constant of Ptet | #m | [10] |
k_{ci} | 20 | Dissociation constant of Pci | #m | [10] |
k_{cUlp} | 3 | Turnover rate of Ulp | min^{-1} | [6] |
n_{ci} | 3 | Hills coefficient | - | [11] |
n_{tet} | 3 | Hills coefficient | - | [11] |
s | 0 or 1 | Activation/Inactivation of T7 promoter | Binary | Assumption |
s_{ci} | 228 | Length of CI | amino acids | [12] |
s_{PSU} | 18 + 110 | Length of peptide plus SUMO | amino acids | [12] |
s_{TET} | 206 | Length of TET | amino acids | [13] |
s_{Ulp} | 233 | Length of Ulp1 | amino acids | [13] |
Variables
Variable | Description |
TET_{m} | concentration of translated TET |
PSU_{m} | concentration of translated peptide plus SUMO |
CI_{m} | concentration of translated C_{i} |
Ulp_{m} | concentration of translated Ulp |
TET | concentration of transcribed TET |
PSU | concentration of transcribed peptide plus SUMO |
CI | concentration of transcribed C_{i} |
Ulp | concentration of transcribed Ulp |
Results
Here the results are given of the simulation upon activating the T7 promoter. For starting conditions the steady state values of the concentrations are used when T7 is switched off.
Figure 2: Simulation Results
Conclusion/Discussion
In the graph of Figure 2 the behavior of the circuit can be seen and the answers to the questions can be given:
- The peptide+SUMO is being produced at time point zero, reaching the maximum number of molecules, 3400 molecules, at time point 70.
- The cleaving of the SUMO can be observed from the start probably due to the fact that the promoter is leaky. After 80 minutes more Ulp-1 is produced, leading to an increased cleaving of the SUMO.