Team:TU-Delft/Timer Plus Sumo
From 2013.igem.org
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
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