Team:Brown/Modeling/Parameters
From 2010.igem.org
Parameters
As our circuit is relatively complex, we have many factors and thus many parameters with which to describe the function of these factors mathematically.
These parameters are as follows:
- δ is a function designed to represent the presence/absence of light in a Boolean fashion. δ=0 without light and 1 with light.
- c_1 represents the rate of conformational change from LovTAP to LovTAP* when irradiated with 470nm light.
- c_2 represents the rate of SupD + T7ptag --> T7 polymerase
- α represents the maximum protein synthesis rate in nanomoles/sec.
- β represents the basal synthesis rate. We assume ∀x β_x= .01α_x
- μ represents the protein degradation rate.
- n represents a ligand’s hill coefficient.
- K_m represents the ligand’s dissociation constant
Players
The following are the players that are used in the modeling of our circuit:
Transcription factors
LovTAP
LovTAP* (After light-induced conformational change)
tetR
Mnt
AraC
LacI
CI
CI434
GAL4
Reporters
S1
S2
S3
S4
Constants
The following are the constants for which we need to find online or determine experimentally:
$\delta = \left\{ \begin{array}{ll} 0 & \mbox{without light}\\ 1 & \mbox{with light} \end{array} \right.$
c represents the rate of conformation change from LovTAP to LovTAP* under 470nm light.
Parameters
αx represents the xth synthesis rate $\left(\dfrac{\text{nanomoles}}{\text{min}}\right)$.
α = (rate of transcription) × (rate of translation)
For promoters for which reliable data is not available, we assume an average E. coli transcription speed to be 70 nucleotides/second = 4,200 nucleotides/min, and an average translation speed of 40 amino acids/second = 2400, which is then further regulated by the appropriate ribosome binding site, represented as a normalized constant. Thus, we use the following equation to calculate unknown synthesis rates[0]:
$\alpha = \left(\frac{4200}{\text{gene length}}\right)\times\left(\frac{2400\times\text{RBS Strength}}{\text{protein length}}\right) $α1 LovTAP synthesis rate constant  = 1. 06432 × 10 - 13 nanomoles. α2 tetR synthesis rate constant  = 1. 15318 × 10 - 13 nanomoles. α3 Mnt synthesis rate constant  = 6. 0562 ×  * 10 - 13 nanomoles. α4 AraC synthesis rate constant  = 5. 99989 × 10 - 14 nanomoles. α5 LacI synthesis rate constant  = 3. 94791 × 10 - 14 nanomoles. α6 CI synthesis rate constant  = 8. 93023 × 10 - 14 nanomoles. α7 CI434 synthesis rate constant  = 1. 2236 × 10 - 13 nanomoles. α8 [check]SupD synthesis rate constant  = 1. 2236 × 10 - 13 nanomoles. α9 [check]T7ptag synthesis rate constant  = 1. 2236 × 10 - 13 nanomoles. α10 [check]GAL4 synthesis rate constant  = 1. 06432 × 10 - 13 nanomoles.
α11 S1 synthesis rate constant  = 9. 68992 × 10 - 14 nanomoles. α12 S2 synthesis rate constant  = 9. 68992 × 10 - 14 nanomoles. α13 S3 synthesis rate constant  = 9. 68992 × 10 - 14 nanomoles. α14 S4 synthesis rate constant  = 9. 68992 × 10 - 14 nanomoles.
βx represents xth basal (un-induced or un-repressed) synthesis rate of a given promoter. We assume that, for all promoters, this is equal to 1% of the synthesis rate constant.
∀ x βx = 0. 1αxμx represents the degradation rate of a given protein. When the degredation rate is unknown, we assume a decay of 0.012 proteins/min (half-life of one hour) μ1 = 0. 0453 Rate of LovTAP degredation (1 / sec)μ2 = 0. 0453 Rate of LovTAP* degredation (1 / sec)μ3 = 0. 0453 tetR degredation (1 / sec)μ4 = 0. 0453 Rate of Mnt degredation (1 / sec)μ5 = 0. 0453 Rate of AraC degredation (1 / sec)μ6 = 0. 0453 Rate of LacI degredation (1 / sec)μ7 = 0. 0453 Rate of CI degredation (1 / sec)μ8 = 0. 0453 Rate of CI434 degredation (1 / sec)μ9 = 00453 Rate of SupD degredation (1 / sec)μ10 = 0. 0453 Rate of T7ptag degredation (1 / sec)μ11 = 0. 012 Rate of T7 polymerase degredation (1 / sec)μ12 = 0. 012 Rate of GAL4 degredation (1 / sec)μ13 = 0. 012 Rate of S1 degredation (1 / sec)μ14 = 0. 012 Rate of S2 degredation (1 / sec)μ15 = 0. 012 Rate of S3 degredation (1 / sec)μ16 = 0. 012 Rate of S4 degredation (1 / sec)
nx represents the hill coefficient for a given protein. n1 = 1 for LovTAP*
n2 = 3 for tetR)
n3 = 2 for CI
n4 = 2 for CI434
n5 = 2 Hill coefficient of GAL4
n6 = 2 Hill coefficient of AraC
n7 = 2 Hill coefficient of T7
n8 = 2 Hill coefficient of LacI
n9 = 1 Hill coefficient of Mnt
Kx represents, for a given protein, the ligand concentration producing half occupation (ligand concentration occupying half of the binding sites) in nanomoles. This is also the infamous "microscopic dissociation constant". Kd1 = 142 for LovTAP*
Kd2 = 0. 179 for tetR
Kd3 = 50 for CI
Kd4 = 40 for CI434
Kd5 = 0. 5 for GAL4
Kd6 = 14 for AraC
Kd7 = 2 for T7
Kd8 = 800for LacI
Kd9 = 50 for Mnt
c c1 = 20. 20, which represents the rate of light-induced LovTAP to LovTAP* conformation change (1 / sec) c2 = 20. 20 which represents the rate of the reaction SupD + T7ptag  →  T7
$\delta = \left\{ \begin{array}{ll} 0 & \mbox{without light}\\ 1 & \mbox{with light} \end{array}\right.$
For a complete list of values and sensitivities, please see our [http://brownigem.com/parameters_2010.ods spreadsheet of parameters and sensitivities].