Team:UANL Mty-Mexico/Modelling/Biphasic Switch
From 2011.igem.org
One of the main elements of our project is the “Biphasic Switch”, as dubbed by Chris A Voigt in his review Genetic Parts to Program a Bacteria (2006).
In our model, we include the affinity change due to the cooperativity effect and the dual behavior of cI. The affinity shift towards O2 and O3 is regulated by a logic rule regarding cI concentration, so when cI concentration is below its dissociation constant towards O1, the dissociation constant towards O2 and O3 remains at basal levels. But when cI concentration rises above its O1 dissociation constant, the cooperativity effect takes place and the dissociation constants towards O2 and O3 are lowered, increasing cI affinity.
Finally, we modeled the dual behavior of cI by considering operators O2 and O3 as if they were different promoters. Thus, their effect on the change of the regulated gene’s concentration is reflected in two different elements of the differential equations, one for activation at O2 and another for repression at O3. The different dissociation constants towards O1, O2 and O3, both for basal activity and for the cooperativity effect, were calculated from the total free energies for the binding of cI to the three operators published in Ackers GK, et al., (1982) and according to the formula presented there. The effect of cI on its regulated genes is then represented as follows:
Where α is the maximum transcription rate; [cI] is the cI protein concentration; O2Kd is the dissociation constant of cI towards the O2 operator; O3Kd is the dissociation constant of cI towards operator O3; and µ is the degradation rate of the mRNA. Note that both the cI protein concentration and the dissociation constants have an exponent. This exponent is the Hill coefficient and equals two.
The figure 1 is a diagram of a Biphasic Circuit that responds to an IPTG input. The LacI gene is obviated.
In brief, after the addition of IPTG to the culture, the LacI repression begins to cease, allowing cI to be expressed. This rise of cI concentration is supposed to be dependent on the IPTG input concentration, so that at low IPTG, the induction of cI is expected to be also low and in the range of the pRM activation state. On the other hand, when the IPTG input is high enough, the high cI induction will switch pRM into its repressed state.
The pRM is controlling the expression of two genes: the reporter GFP and another transcription factor cI434, which is negatively regulating the expression of another reporter gene, YFP. This reporter gene is also regulated by cI, but in a positive manner at any concentration considered. In this way, when the pRM promoter is active, GFP and cI434 will be expressed, giving rise to the first state of our Biphasic Switch, i. e., the green state. Conversely, when pRM is inhibited by the high concentration of cI, the inhibiton by cI434 will be gradually released and the second state will be observed, i. e., the yellow state (because of YFP).
The time for the IPTG addition in all simulations is at 300 min. We considered a simulation time of 1000 min. just for the sake of clarity in the graphs, but the actual experimental times may be shorter. In all the graphs shown, the Y-axis represents protein concentration in nM units and the X-axis represents time in minutes.
The initial condition of the LacI concentration was determined in a simulation where, without any disturbance, the LacI gene was let to reach its saturation level. This saturation level is taken as an initial condition in order to avoid noisy behavior at the first minutes of simulation. Furthermore, in the wet-experiments, these simulated initial conditions may represent the incubation time previous the IPTG input, when the cells are allowed to grow, express LacI without disturbance and sustain the repression of the cI gene and of the Biphasic Switch as well. As seen, the green state is almost undetectable throughout the whole experiment, while the yellow one is activated with cI activation.
For cI to act as expected, i.e. a dual effect transcription factor, there must be some changes to the original circuitry, starting with modifications in LVA tags (all genes are considered to have one), RBS activity differences and perhaps some connectivity modifications.
The effect of RBSThe simulations in Figure 2 were performed considering that all genes have the same RBS efficiency. But in order to simulate the effect of different RBS, we assumed that protein translation is proportional to the activity of the specific RBS present in its mRNA. We decided to use the RBS Community Collection from the Registry of Standard Biological Parts (http://partsregistry.org/Ribosome_Binding_Sites/Prokaryotic/Constitutive/Community_Collection) because this set’s relative efficiencies have already been characterized.
We are aware that multiplying the maximum transcription rate times the relative RBS efficiency gives only a relative result, because the RBS characterization doesn’t report efficiencies in a nM concentration scale. The best approach would be to determine the ratio between the theoretical maximum translation rate and the translation rate for a given RBS, all in nM units. For this reason, the simulations of the RBS effects are reported in relative units, assuming that maximum transcription rate is achieved with RBS BBa_B0034, i.e. the reference RBS of the Community Collection. The RBS Relative Units are proposed as a unit measure for gene expression in our simulations, as the actual behavior of the protein concentration dynamics in nanomolar units can be also reflected in relative units, considering the basal signal as zero.
We decided to modify the RBS from cI, GFP and YFP. Protein cI was selected in order to maintain its concentration in a low range for some time, so that GFP inhibition can be delayed. Protein GFP was selected in order to increase its translation rate when using the most powerful RBS of the set. And finally, protein YFP was selected in order the noise it generates when the green state is turned on. Also, LacI LVA tag was removed, increasing its half-life. Furthermore, two LacI genes are considered in simulation D. The relative strengths of the RBS for the genes are the following:
Gene/RBS | Gene/RBS | Gene/RBS | Gene/RBS | Gene/RBS |
LacI/1.124 |
LacI/1.124 |
LacI/1.124 |
2LacI/1.124 |
2LacI/1.124 |
cI/0.6 |
cI/0.3 |
cI/0.01 |
cI/0.01 |
cI/0.01 |
GFP/1.124 |
GFP/1.124 |
GFP/1.124 |
GFP/1.124 |
GFP/1.124 |
cI434/1.124 |
cI434/1.124 |
cI434/1.124 |
cI434/1.124 |
cI434/1.124 |
YFP/0.3 |
YFP/0.01 |
YFP/0.01 |
YFP/0.01 |
YFP/0.01 |
Simulation E was performed considering an initial IPTG concentration of 285 nM, because, in a series of simulations, it was observed that GFP expression was greatly induced (8-10 nM steadily sustained) when working with IPTG in the 250-285 nM range.
The conditions of the simulation shown in figure 4 seem to fulfill for the expected behavior of the biphasic switch. That is, when IPTG concentration is null, the GFP and YFP expression remains basal and undetectable. But when it raises to a certain threshold (that was determined to be in the range of 250 – 285 nM IPTG by other, non-shown, simulations) and remains there, cI enters into its activator state and the GFP expression is turned on (to values of aprox. 10 RBS Relative Units), as well as cI434, while YFP remains basal. Then, when the IPTG concentration is incremented to a second threshold, cI starts acting as an inhibitor and turns off GFP and cI434 expression, while YFP expression is turned on to values around 10 RBS Relative Units, which are comparable to the GFP levels. Take into account that these values are calculated for a single cell, so they must be corrected for population density, as the fluorescence measured is the general signal from a sample.
In order to reset this chemically induced Biphasic Switch, the IPTG concentration must be set again to zero. This can be achieved by subculturing the IPTG induced culture into a series of fresh media cultures until no fluorescent signal is detected or by plating the culture and growing colonies in fresh media until, again, no fluorescence is detected. Another reset strategy is the Quorum Sensing connection by which a TetR gene is activated, so that it represses cI.
Introduction
The red light system is based on a chimeric photoreceptor, which is composed by a photosensible domain and an EnvZ domain. This EnvZ domain is activated through darkness-stimulated autophosphorylation and then acts as a protein kinase, which in turn, phosphorylates and activates the transcription factor OmpR. Conversely, when red light is present, EnvZ will no longer exhibit a kinase activity, instead, it will act as a phosphatase that will inactivate the phosphorylated OmpR.
When phosphorylated, OmpR will turn on transcription on the PompC promoter, which in our system controls the expression of the repressor Mnt. Mnt is in turn repressing cI, so when red-light is present, Mnt expression will be decreased because the concentration of phosphorylated OmpR declines, and cI expression will be turned-on.
The expression of the chimeric photoreceptor and OmpR is given by the gene expression equations already shown. But the phosphorylation process is expressed by a set of equations and paramaters that were adapted from Srividhya and Krishnaswamy, (2004). These equations express the binding and dissociation of OmpR to EnvZ and, in our model, it is assumed that all the OmpR proteins that bind to EnvZ are phosphorylated. The parameters have been transformed from µM/second units to nM/minute.
After red-light exposure, all the EnvZ proteins in the membrane are assumed to be de-activated. Thus, these last equations are valid only when the light input is OFF.
When the cI gene is turned on using red light, the activity of the Biphasic Switch starts. The expression of the first state (GFP) is achieved through intermittent pulses of red light and far-red light; meanwhile, the second state is achieved through continuous exposure to red light.
SimulationsA number of simulations were performed in order to determine the appropriate time proportion for the cycles of red/far-red light exposure in the intermittent pulses. The next table shows the RBS strengths and the initial conditions that were considered: