- Research Article
- Open Access
Hypothetical biomolecular probe based on a genetic switch with tunable symmetry and stability
- Nikolay Martyushenko^{1},
- Sigurd Hagen Johansen^{1},
- Cheol-Min Ghim^{2, 3} and
- Eivind Almaas^{1}Email author
https://doi.org/10.1186/s12918-016-0279-y
© Martyushenko et al. 2016
- Received: 23 December 2015
- Accepted: 13 April 2016
- Published: 6 June 2016
Abstract
Background
Genetic switches are ubiquitous in nature, frequently associated with the control of cellular functions and developmental programs. In the realm of synthetic biology, it is of great interest to engineer genetic circuits that can change their mode of operation from monostable to bistable, or even to multistable, based on the experimental fine-tuning of readily accessible parameters. In order to successfully design robust, bistable synthetic circuits to be used as biomolecular probes, or understand modes of operation of such naturally occurring circuits, we must identify parameters that are key in determining their characteristics.
Results
Here, we analyze the bistability properties of a general, asymmetric genetic toggle switch based on a chemical-reaction kinetic description. By making appropriate approximations, we are able to reduce the system to two coupled differential equations. Their deterministic stability analysis and stochastic numerical simulations are in excellent agreement. Drawing upon this general framework, we develop a model of an experimentally realized asymmetric bistable genetic switch based on the LacI and TetR repressors. By varying the concentrations of two synthetic inducers, doxycycline and isopropyl β-D-1-thiogalactopyranoside, we predict that it will be possible to repeatedly fine-tune the mode of operation of this genetic switch from monostable to bistable, as well as the switching rates over many orders of magnitude, in an experimental setting. Furthermore, we find that the shape and size of the bistability region is closely connected with plasmid copy number.
Conclusions
Based on our numerical calculations of the LacI-TetR asymmetric bistable switch phase diagram, we propose a generic work-flow for developing and applying biomolecular probes: Their initial state of operation should be specified by controlling inducer concentrations, and dilution due to cellular division would turn the probes into memory devices in which information could be preserved over multiple generations. Additionally, insights from our analysis of the LacI-TetR system suggest that this particular system is readily available to be employed in this kind of probe.
Keywords
- Genetic toggle switch
- Stochastic gene expression
- LacI-TetR system
- Nano-recorder
Background
Genetic switches and circuits are abundant in nature and control the regulation and genetic development programs that are at the basis of cellular function. Spanning simple feedback loops in B. subtilis competence switching [1] to complex differentiation patterns in humans [2], ultra-stable genetic switches present living organisms with a reliable mode of operation within a noisy state until a signal is received that triggers the transition to a different state of operation. The physical properties, dynamics and stability of genetic switches are important, not only because of their common occurrence in biological systems, but also for their potential use in synthetic biology in the engineering and design of more complicated cellular functions [3, 4].
An important challenge in the successful design of a genetic circuit is the fine-tuning of various circuit parameters [4]: In even simple genetic circuit topologies, it is most often necessary to fine-tune the strength and responsiveness of different circuit components in order to achieve intended functional behavior. For instance, the consequences of noise can significantly affect the function of a genetic circuit [5]. Thus, it is of great interest to develop circuits and components that are readily amenable to experimental control.
Among the early demonstrations of the engineering of synthetic circuits, Gardner et al. used molecular genetic tools to construct a bistable toggle switch consisting of two genes coding for mutually repressing proteins [6]. They also implemented control of the toggle-switch state of operation by using an inducer that could affect transcription in the circuit, and thus, affect its stability (monostable or bistable). Since then, a considerable amount of effort has been put into improving switch robustness, tunability and scalability. For instance Deans et al. combined tunability with robustness in a protein-RNA genetic switch [7], while Green et al. combined robustness with unprecedented scalability using an RNA toehold system [8]. While both protein and RNA switches are commonly used in scientific practice, the focus on robustness, being the lack of dynamics, has led to the fact that not many studies have investigated the dynamics of specific or generic genetic switches. Recently, however, Shopera et al. investigated the properties of several single and dual-positive feedback synthetic genetic switches, coupling experimental evidence to deterministic mathematical modelling, and amongst other things, generating experimental and theoretical stability diagrams [9]. Unfortunately, this study was limited by the absence of reliable kinetic information concerning the circuit system in question and by the lack of stochastic analysis, both experimental and theoretical.
Here, we discuss in detail a mathematical and computational model for a generic, asymmetric genetic toggle switch, following up on our previous reports on symmetric genetic switches [10, 11]. In these toggle switches, the active repressor species is a protein dimer. This class of switch was demonstrated to have a large dynamic range of switching frequency over a relatively small gene expression efficiency window [11].
At first we treat the general case, where the two mutually inhibiting genetic circuits are equivalent, except for a chosen set of kinetic constants making it asymmetric. For this circuit, we show that it is possible to reduce the system to two coupled non-linear differential equations that depend on only four tunable parameters (two for each gene). This is analogous to the analysis that was conducted in [11] with focus on a symmetric switch.
Next we use the framework of the general asymmetric toggle switch to develop a detailed mathematical model for an experimentally realizable switch. To make the model construct as accurate as possible, we elected to use transcription factors for which we can utilize the greatest amount of quantitative data. We decided to utilize the LacI and TetR repressors, which we will denote the LITR-switch. These repressors are especially interesting because they are amenable to fine-tuning by the addition of IPTG and doxycyline (inhibitors of lacI and TetR, and thus inducers of LacI and TetR-repressed expression respectively), allowing the degree of asymmetry caused by their differing expression strengths to be easily modified experimentally. Using values for kinetic constants derived from literature, we explored the circuit’s states of operation through both stochastic and deterministic simulations. Interestingly, some of our findings can be related to properties of the circuit previously studied in vivo by Gardner et al. [6], since one of their constructed toggle switches is quite similar in design.
Results
Processes included in the HOM2 genetic circuit description
Type of reaction | Gene 1 | Gene 2 |
---|---|---|
Repressor binding | \(\mathrm {D}^{1}_{00} + \mathrm {P}^{2}_{2} \underset {q_{12}}{\overset {k_{12}}{\rightleftharpoons }} \mathrm {D}^{1}_{10}\) | \(\mathrm {D}^{2}_{00} + \mathrm {P}^{1}_{2} \underset {q_{22}}{\overset {k_{22}}{\rightleftharpoons }} \mathrm {D}^{2}_{10}\) |
\(\mathrm {D}^{1}_{10} + \mathrm {P}^{2}_{2} \underset {q_{14}}{\overset {k_{14}}{\rightleftharpoons }} \mathrm {D}^{1}_{20}\) | \(\mathrm {D}^{2}_{10} + {P^{1}_{2}} \underset {q_{24}}{\overset {k_{24}}{\rightleftharpoons }} \mathrm {D}^{2}_{20}\) | |
RNAp binding | \(\mathrm {D}^{1}_{00} + \mathrm {R} \underset {q_{13}}{\overset {k_{13}}{\rightleftharpoons }} \mathrm {D}^{1}_{01}\) | \(\mathrm {D}^{2}_{00} + \mathrm {R} \underset {q_{23}}{\overset {k_{23}}{\rightleftharpoons }} \mathrm {D}^{2}_{01}\) |
\(\mathrm {D}^{1}_{10} + \mathrm {R} \underset {q_{15}}{\overset {k_{15}}{\rightleftharpoons }} \mathrm {D}^{1}_{11}\) | \(\mathrm {D}^{2}_{10} + \mathrm {R} \underset {q_{25}}{\overset {k_{25}}{\rightleftharpoons }} \mathrm {D}^{2}_{11}\) | |
\(\mathrm {D}^{1}_{20} + \mathrm {R} \underset {q_{17}}{\overset {k_{17}}{\rightleftharpoons }} \mathrm {D}^{1}_{21}\) | \(\mathrm {D}^{2}_{20} + \mathrm {R} \underset {q_{27}}{\overset {k_{27}}{\rightleftharpoons }} \mathrm {D}^{2}_{21}\) | |
Transcription initiation | \(\mathrm {D}^{1}_{01} \overset {\alpha _{1m}}{\rightarrow } \mathrm {E}^{1} + \mathrm {D}^{1}_{00}\) | \(\mathrm {D}^{2}_{01} \overset {\alpha _{2m}}{\rightarrow } \mathrm {E}^{2} + \mathrm {D}^{2}_{00}\) |
\(\mathrm {D}^{1}_{11} \overset {\alpha _{1m}}{\rightarrow } \mathrm {E}^{1} + \mathrm {D}^{1}_{10}\) | \(\mathrm {D}^{2}_{11} \overset {\alpha _{2m}}{\rightarrow } \mathrm {E}^{2} + \mathrm {D}^{2}_{10}\) | |
\(\mathrm {D}^{1}_{21} \overset {\alpha _{1m}}{\rightarrow } \mathrm {E}^{1} + \mathrm {D}^{1}_{20}\) | \(\mathrm {D}^{2}_{21} \overset {\alpha _{2m}}{\rightarrow } \mathrm {E}^{2} + \mathrm {D}^{2}_{20}\) | |
Elongation | \(\mathrm {E}^{1} \overset {\alpha _{1m}^{\prime }}{\rightarrow } \mathrm {M}^{1} + \mathrm {R} \) | \(\mathrm {E}^{2} \overset {\alpha _{2m}^{\prime }}{\rightarrow } \mathrm {M}^{2} + \mathrm {R} \) |
Translation | \(\mathrm {M}^{1} \overset {\alpha _{1p}}{\rightarrow } \mathrm {P}^{1} + \mathrm {M}^{1} \) | \(\mathrm {M}^{2} \overset {\alpha _{2p}}{\rightarrow } \mathrm {P}^{2} + \mathrm {M}^{2} \) |
Dimerization | \(\mathrm {P}^{1} + \mathrm {P}^{1} \underset {q_{11}}{\overset {k_{11}}{\rightleftharpoons }} \mathrm {P}^{1}_{2}\) | \(\mathrm {P}^{2} + \mathrm {P}^{2} \underset {q_{21}}{\overset {k_{21}}{\rightleftharpoons }} \mathrm {P}^{2}_{2}\) |
Degradation | \(\mathrm {M}^{1} \overset {\gamma _{1m}}{\rightarrow } \oslash \) | \(\mathrm {M}^{2} \overset {\gamma _{2m}}{\rightarrow } \oslash \) |
\(\mathrm {P}^{1} \overset {\gamma _{1p}}{\rightarrow } \oslash \) | \(\mathrm {P}^{2} \overset {\gamma _{2p}}{\rightarrow } \oslash \) | |
\(\mathrm {P}^{1}_{2} \overset {\gamma _{1p}/\sigma _{1}}{\rightarrow } \oslash \) | \(\mathrm {P}^{2}_{2} \overset {\gamma _{2p}/\sigma _{2}}{\rightarrow } \oslash \) |
where the parameters are all defined in the Methods section. The form of Eq. (1) lends itself to the following interpretation: Since the first term is always positive, we may consider it an effective synthesis rate in a birth-death process. The second term is always negative and would correspond to the decay rate in such a process. The synthesis term of a dimer is proportional to the root of its own concentration, since the dimer must be produced from the monomer and equilibrium is assumed to occur instantaneously. The first part of this term corresponds to protein synthesis at full repression, while the second represents repression by the competing agent. Furthermore, the first part of the decay term represents degradation of dimers, while the second represents that of the monomers. This result is the generalization of the symmetric HOM2 circuit in Ghim and Almaas (2009) [11] to the situation where the two interacting genes have different characteristics.
In the calculations culminating in Eq. (1), we have applied approximations (i)-(iv) described in the Methods section and the simplification that the dissociation parameters K _{ l7}=K _{ l5}. We do not believe that these significantly reduce the generality of Eq. (1). The main assumption here is the adiabatic assumptions (ii) and (iv), that the repressor binding and dimerization as well as transcriptional elongation reactions are in steady state. While making the problem analytically tractable, this assumption has been shown to have only minor effects on the results of genetic switch simulations [11]. The mentioned simplification of the two dissociation parameters equates the RNApol binding constants in a promoter singly and doubly bound by repressors. Though this is not the case in a general system, the simplification still preserves the effect of double repressor binding. Even thought the binding of the second repressor does not affect transcription initiation efficiency directly, it reduces the probability of a repressor-free promoter, thus thereby still reducing the average transcription rate. Additionally, when exploring the properties of the system, we have made the simplifying assumptions γ _{ lp }=γ _{ p } and K _{ l2}=K _{2}. The first assumption equates the protein degradation rates of the two switch-components. This is true for E.coli in exponential growth phase, as protein degradation in that case is dominated by dilution due to cell division [12]. The final assumption is that the affinities of the two repressor types for their respective promoters are the same. This assumption was implemented in order to restrict the exploration of the system to only two parameters, as seen in the next section.
Reduced system: deterministic and stochastic analysis
Parameter values used to perform a deterministic analysis of the general, reduced system based on selections made in [11]
K _{2} | K _{ l1} | r _{ l } | s _{ l } | β _{ l } | u _{ l } | σ _{ l } |
---|---|---|---|---|---|---|
20 nM | 10 nM | 25 | 0.01 | 17.5 nM^{−1} | 3.0 | 10 |
The individual gene expression efficiency from each of the promoters, described by the parameters β _{1} and β _{2}, are able to vary independently of each other. The β-variables can potentially be modified by genetic manipulation, especially by modifying the corresponding 5’UTR regions [13–15], making the β _{1} vs β _{2} stability diagram of great interest. Figure 2(b) displays the region of bistability (inside the curves) for the sequence of values s={0.002,0.01,0.03,0.1}. As promoter leakage s increases, the bistable region shows a systematic decrease.
Figure 2(c) shows the stability diagram (dotted black line) when varying β _{2} and s _{2} while keeping β _{1} and s _{1} as in Table 2. We note that, while the general shape is similar to that found for the simultaneous and identical change of β _{ l } and s _{ l } [11], the extremal part of the bistability region (large s _{2}) is much sharper in the current case.
We will later simulate the full reaction set of the LITR-system, as described by Eqs. (A.1)-(A.21) and, for the sake of brevity, do not present such simulations of the current circuit. Previously, we employed an approximate representation of the stochastic dynamics for the symmetric version of the circuit [11]: Assuming a simple birth-death process where the propensities for synthesis and decay are described by Eq. (1), this stochastic system is governed by the average rates we calculated under the adiabatic assumptions (i)-(iv), which we detail in the Methods section.
We also conduct a stochastic exploration of the stability diagram (demarcated by dotted lines) reported in Fig. 2(c) by keeping the pair of the leakage/efficiency parameters at their reference values, s _{1}=0.01β _{1}=17.5 nM^{−1}, while varying s _{2} and β _{2} as indicated by the axes in Fig. 2(c). The bistable gene-switch system has two possible starting points: either gene 1 is in an active state (with high protein copy number) or in an inactive state (with low protein copy number), and we use both conditions as starting points for two separate 10^{3} time-step and 10^{7} time-step simulations.
The results reported in Fig. 2(c) show that the bistable region determined by stochastic simulations lies within, and almost entirely fills, the one calculated using deterministic analysis. The switching rate is not of great practical relevance here, because the time is scaled relative to that of the protein decay rate. The high switching rate region that extends out of the deterministic bistable region corresponds to a situation where the two genes behave independently of one another (see next section for more detail).
LITR system: stochastic analysis
Following our analysis of the generalized reduced genetic switch, we decided to develop a realistic model for a switch that could be experimentally implemented. We simulated the circuit and looked at its dynamic properties.
For this we needed to identify two transcription factors for which enough kinetic data is available and which state of operation could be tunable with small molecule inducers. The tunability is necessary to make the circuit accessible to experimental manipulation for mapping out its stability states. The choice of small molecule inducers for the circuit manipulation (as opposed to e.g. promoter sequence alterations or genetic mutations of the repressor species) was to make the in vivo experimental system tunable and possible to re-set.
Processes included in the LacI-TetR (LITR) genetic circuit
Type of reaction | lacI (with P_{tet}) | tetR (with P_{lac}) |
---|---|---|
Repressor binding | \(\text {TetR}_{2} + \mathrm {D}_{l} \underset {q_{t3}}{\overset {k_{t3}}{\rightleftharpoons }} \mathrm {d}_{1}\text {TetR}_{2}\) | \(\text {LacI}_{4} + \mathrm {D}_{t} \underset {q_{l3}}{\overset {k_{l3}}{\rightleftharpoons }} \mathrm {d}_{1}\text {LacI}_{4}\) |
\(\mathrm {d}_{1}\text {TetR}_{2} + \mathrm {D}_{l} \underset {q_{t3}}{\overset {k_{t3}}{\rightleftharpoons }} \mathrm {d}_{2}\text {TetR}_{2}\) | \(\text {LacI}_{4i} + \mathrm {D}_{t} \underset {q_{l3}}{\overset {0.5k_{l3}}{\rightleftharpoons }} \mathrm {d}_{1}\text {LacI}_{4i}\) | |
\(\mathrm {d}_{1}\text {TetR}_{2} + \mathrm {I}_{t} {\overset {k_{t2}}{\rightarrow }} \text {TetR}_{2i} + \mathrm {D}_{l}\) | \(\mathrm {d}_{1}\text {LacI}_{4} \underset {q_{l4}}{\overset {200k_{l3}}{\rightleftharpoons }} \mathrm {d}_{2}\text {LacI}_{4}\) | |
\(\mathrm {d}_{2}\text {TetR}_{2} + \mathrm {I}_{t} {\overset {k_{t2}}{\rightarrow }} \text {TetR}_{2i} + \mathrm {d}_{1}\text {TetR}_{2}\) | \(\mathrm {d}_{1}\text {LacI}_{4i} + \mathrm {I}_{l} {\overset {k_{l5}}{\rightarrow }} \text {LacI}_{4ii} + \mathrm {D}_{t}\) | |
\(\mathrm {d}_{1}\text {LacI}_{4} + \mathrm {I}_{l}{\overset {0.5k_{l5}}{\rightarrow }} \text {LacI}_{4i} + \mathrm {D}_{t}\) | ||
\(\mathrm {d}_{1}\text {LacI}_{4} + \mathrm {I}_{l}\underset {q_{l5}}{\overset {0.5k_{l5}}{\rightleftharpoons }} \mathrm {d}_{1}\text {LacI}_{4i}\) | ||
\(\mathrm {d}_{2}\text {LacI}_{4} + \mathrm {I}_{l}\underset {q_{l5}}{\overset {k_{l5}}{\rightleftharpoons }} \mathrm {d}_{1}\text {LacI}_{4i}\) | ||
RNAp binding | \(\mathrm {D}_{l} + \mathrm {R} \underset {q_{b}}{\overset {k_{b1}}{\rightleftharpoons }} \mathrm {D}^{\mathrm {R}}_{l}\) | \(\mathrm {D}_{t} + \mathrm {R} \underset {q_{b}}{\overset {k_{b1}}{\rightleftharpoons }} \mathrm {D}^{\mathrm {R}}_{t}\) |
\(\mathrm {d}_{1}\text {LacI}_{4} + \mathrm {R} \underset {q_{b}}{\overset {k_{lb2}}{\rightleftharpoons }} \mathrm {d}_{1}\text {LacI}^{\mathrm {R}}_{4}\) | \(\mathrm {d}_{1}\text {TetR}_{2} + \mathrm {R} \underset {q_{b}}{\overset {k_{tb2}}{\rightleftharpoons }} \mathrm {d}_{2}\text {TetR}^{\mathrm {R}}_{2}\) | |
\(\mathrm {d}_{2}\text {LacI}_{4} + \mathrm {R} \underset {q_{b}}{\overset {k_{lb3}}{\rightleftharpoons }} \mathrm {d}_{2}\text {LacI}^{\mathrm {R}}_{4}\) | \(\mathrm {d}_{2}\text {TetR}_{2} + \mathrm {R} \underset {q_{b}}{\overset {k_{tb3}}{\rightleftharpoons }} \mathrm {d}_{2}\text {TetR}^{\mathrm {R}}_{2}\) | |
\(\mathrm {d}_{1}\text {LacI}_{4i} + \mathrm {R} \underset {q_{b}}{\overset {k_{lb2}}{\rightleftharpoons }} \mathrm {d}_{1}\text {LacI}^{\mathrm {R}}_{4i}\) | ||
Transcription initiation | \(\mathrm {D}^{\mathrm {R}}_{l} \overset {\alpha _{i}}{\rightarrow } \mathrm {E}^{1} + \mathrm {D}_{l}\) | \(\mathrm {D}^{\mathrm {R}}_{t} \overset {\alpha _{i}}{\rightarrow } \mathrm {E}^{2} + \mathrm {E}_{t}\) |
\(\mathrm {d}_{1}\text {LacI}^{\mathrm {R}}_{4} \overset {\alpha _{i}}{\rightarrow } \mathrm {E}^{1} + \mathrm {d}_{1}\text {LacI}_{4}\) | \(\mathrm {d}_{2}\text {TetR}^{\mathrm {R}}_{2} \overset {\alpha _{i}}{\rightarrow } \mathrm {E}^{2} + \mathrm {d}_{2}\text {TetR}_{2}\) | |
\(\mathrm {d}_{2}\text {LacI}^{\mathrm {R}}_{4} \overset {\alpha _{i}}{\rightarrow } \mathrm {E}^{1} + \mathrm {d}_{2}\text {LacI}_{4}\) | \(\mathrm {d}_{2}\text {TetR}^{\mathrm {R}}_{2} \overset {\alpha _{i}}{\rightarrow } \mathrm {E}^{2} + \mathrm {d}_{2}\text {TetR}_{2}\) | |
\(\mathrm {d}_{1}\text {LacI}^{\mathrm {R}}_{4i} \overset {\alpha _{i}}{\rightarrow } \mathrm {E}^{1} + \mathrm {d}_{1}\text {LacI}_{4i}\) | ||
Elongation | \(\mathrm {E}_{l} \overset {\alpha _{le}}{\rightarrow } \mathrm {M}_{l} + \mathrm {R} \) | \(\mathrm {E}_{t} \overset {\alpha _{te}}{\rightarrow } \mathrm {M}_{t} + \mathrm {R} \) |
Translation | \(\mathrm {M}_{l} \overset {\alpha _{lt}}{\rightarrow } \text {LacI} + \mathrm {M}_{l} \) | \(\mathrm {M}_{t} \overset {\alpha _{tt}}{\rightarrow } \text {TetR} + \mathrm {M}_{t} \) |
Dimerization | \(\text {LacI} + \text {LacI} \underset {q_{1}}{\overset {k_{1}}{\rightleftharpoons }} \text {LacI}_{2}\) | \(\text {TetR} + \text {TetR} \underset {q_{1}}{\overset {k_{1}}{\rightleftharpoons }} \text {TetR}_{2}\) |
\(\text {LacI}_{2} + \text {LacI}_{2} \underset {q_{l2}}{\overset {k_{1}}{\rightleftharpoons }} \text {LacI}_{4}\) | ||
Repressor Binding | \(\text {LacI}_{4} + \mathrm {I}_{l} \underset {q_{l5}}{\overset {k_{l5}}{\rightleftharpoons }} \text {LacI}_{4i}\) | \(\text {TetR}_{2} + \mathrm {I}_{t} \underset {q_{t2}}{\overset {k_{t2}}{\rightleftharpoons }} \text {TetR}_{2i}\) |
\(\text {LacI}_{4i} + \mathrm {I}_{l} \underset {q_{l5}}{\overset {k_{l5}}{\rightleftharpoons }} \text {LacI}_{4ii}\) | \(\text {TetR}_{2i} + \mathrm {I}_{t} \underset {q_{t2}}{\overset {k_{t2}}{\rightleftharpoons }} \text {TetR}_{2ii}\) | |
Degradation | \(\mathrm {M}_{l} \overset {\gamma _{ml}}{\rightarrow } \oslash \) | \(\mathrm {M}_{t} \overset {\gamma _{mt}}{\rightarrow } \oslash \) |
\(\text {LacI} \overset {\gamma _{p}}{\rightarrow } \oslash \) | \(\text {TetR} \overset {\gamma _{p}}{\rightarrow } \oslash \) | |
\(\text {LacI}_{2} \overset {\gamma _{p}}{\rightarrow } \oslash \) | \(\text {TetR}_{2} \overset {\gamma _{p}}{\rightarrow } \oslash \) | |
\(\text {LacI}_{4} \overset {\gamma _{p}}{\rightarrow } \oslash \) | \(\text {TetR}_{2i} \overset {\gamma _{p}}{\rightarrow } \oslash \) | |
\(\text {LacI}_{4i} \overset {\gamma _{p}}{\rightarrow } \oslash \) | \(\text {TetR}_{2ii} \overset {\gamma _{p}}{\rightarrow } \oslash \) | |
\(\text {LacI}_{4ii} \overset {\gamma _{p}}{\rightarrow } \oslash \) |
Summary of the parameter values used to perform stochastic and deterministic analysis of the realistic LacI-TetR system, based on the sources described in Methods
α _{ le } | 0.065 s^{−1} | α _{ te } | 0.11 s^{−1} | α _{ i } | 0.66 s^{−1} |
α _{ lt } | 0.01 s^{−1} | α _{ tt } | 0.017 s^{−1} | γ _{ ml } | 3.0 ×10^{−3}s^{−1} |
γ _{ mt } | 3.9 ×10^{−3}s^{−1} | γ _{ p } | 3.9 ×10^{−4}s^{−1} | k _{1} | 0.015 μm^{3} |
q _{1} | 8.8 ×10^{−5}s^{−1} | q _{ l2} | 8.8 ×10^{−8}s^{−1} | k _{ t2} | 2.4 ×10^{−4} μm^{3} |
q _{ t2} | 0.019 s^{−1} | k _{ l3} | 3.0 ×10^{−3} μm^{3} | q _{ l3} | 5.4 ×10^{−7}s^{−1} |
k _{ t3} | 4.8 ×10^{−3} μm^{3} | q _{ t3} | 2.1 ×10^{−3}s^{−1} | q _{ l4} | 3.0 ×10^{−3}s^{−1} |
k _{ l5} | 1.4 ×10^{−6} μm^{3} | q _{ l5} | 8.2 ×10^{−3}s^{−1} | k _{ b1} | 4.3 ×10^{−3} μm^{3} |
q _{ b } | 0.45 s^{−1} | k _{ l b2} | 1.2 ×10^{−5} μm^{3} | k _{ t b2} | 5.7 ×10^{−6} μm^{3} |
k _{ l b3} | 3.6 ×10^{−7} μm^{3} | k _{ t b3} | 8.6 ×10^{−6} μm^{3} |
We decided to investigate the stability of the circuit from a stochastic perspective. To this end, we simulated the full system for time-courses of 10^{7} s, each implemented with a state sampling every second. We ran two independent simulations with initial conditions being either 400 LacI tetramers or 400 TetR dimers, all other proteins and mRNA species numbers set to zero. We determined the stability by looking at the distributions of Δ=[LacI_{4}]−[TetR_{2}] over the time course of the simulation assuming three possible stability scenarios: (i) Monostability would occur if, for both initial conditions, the switch would be turned to the same side with overwhelming probability. (ii) Bistability would be determined if the initial conditions gave rise to two non-intersecting distributions of Δ. Finally, (iii) the Δ distribution could be bistable but with switching occurring between the two states (joint bistable case).
In all cases, we observe a joint bistable region with high switching rate at large inducer concentrations. In this region, the number of transcription factors not bound to inducers is very close to zero. Consequently, the lack of regulation makes the two genes virtually decoupled and constitutively expressed. Also, it appears that LacI is a more potent transcriptional repressor than TetR, because at zero [IPTG] and [DOX] the system is monostable, with lacI being the “winning” species. This agrees well with experimental observations of Gardner et al. [6]. Here, LacI was expressed under the control of a weaker promoter (P_{LtetO-1}) than the TetR (P_{trc}) in the LITR circuit. Consequently, it was necessary to further reduce its strength in order to experimentally observe bistability. In our case, both species are expressed under the same P_{R} promoter, which should give rise to strong monostability in the absence of inducers. It is worth noting that in these simulations, we use the LacI binding operator O^{sym}, which is much stronger than O^{1} used by Gardner et al. Nevertheless, when we conducted simulations of the system with appropriately reduced operator strength to explore this effect (not shown), the results showed only minor deviations. In this case we used the binding constant for the O^{1} operator taken from Garcia et al. [18]. We also observed that, if the copy number of tetR gene was set to four times that of lacI, the LITR system becomes bistable at zero inducer concentration (not shown). Thus the relative gene copy numbers could be used as a separate parameter to tune the symmetry of genetic switches.
A dilution experiment starting from any point in the diagram in Fig. 4(a) would result in an [IPTG]-[DOX] phase-space movement along a straight line towards the origin. Thus, the switch could keep both its state and bistability upon dilution if the stability diagram would be more symmetric than what is shown here. A symmetrization of the bistability phase plot could be achieved for this system by adjusting plasmid copy numbers (see above). The histogram from running the system along a titration trajectory near the diagonal in Fig. 4(a) can be seen in Fig. 4(b). Similarly, if [IPTG] and [DOX] would be set to large values within the fast-switching bistable region of Fig. 4(a), a dilution would result in a (biased) randomization of the switch values that would depend on the starting concentrations of the inducers.
LITR system: deterministic analysis
In order to investigate the agreement between stochastic modelling and deterministic analysis, we performed a deterministic stability analysis of the LITR system. The results of the analysis are given in Fig. 3. Here, the borders of the deterministically determined bistable regions are displayed as dashed black lines. It can be seen that the bistable regions of the stochastic simulation lie within the boundaries determined deterministically with the exception of the circuit with just one copy of each gene (Fig. 3(a)). Interestingly, the stochastic simulation in this case shows LacI to be a stronger transcriptional repressor than in the deterministic analysis: In panel Fig. 3(a), we see that the bistability region is shifted towards larger values of [IPTG], and it is especially evident for [DOX] ∼0. This effect is likely explained by a combination of the shallowness of the wells in the potential landscape at small inducer concentrations and our definition of bistability: The definition is based on the time that the system spends with [LacI_{4}]>[TetR_{2}] and vice versa, not on the shape of the potential landscape. This is because the potential landscape can be rugged, which increases the likelihood of detecting false wells. At larger gene copy numbers, stochastic and deterministic simulations agree very well. This is likely owing to the increased depth of the potential wells in that case. Also, the deterministic analysis does not identify the stochastically determined fast switching bistable region (red) as bistable. This is again related to our definition of bistability for the stochastic simulation. In this case, the potential well is very narrow and centered at zero (see Fig. 4b).
Discussion
In this article we have analyzed and outlined the workings of an asymmetric genetic switch. At first, we looked at a generalized example of a bistable genetic switch. While it was relatively straightforward to simplify and theorize, it was challenging to relate this specific circuit implementation to experimental realizations, because its asymmetry-controlling parameters are difficult to access directly in an experimental setting. Therefore, we conducted an investigation of a realistic circuit, the LITR system, based on a mutually repressing LacI-TetR design. While this is challenging to simplify, the computational time needed for direct stochastic and deterministic simulations of this system is not prohibitive. Since we use small molecules and the gene copy numbers to tune circuit symmetry and stability, it does not appear unrealistic to experimentally verify or disprove our findings. Additionally, by letting fluorescent reporters to be under the control of LacI and TetR, one could explore the bistability states of the [IPTG]-[DOX] plane for different gene and plasmid copy number combinations. If sufficiently automated, this could provide us with diagrams similar to those in Fig. 3.
- 1.
A chemical or physical molecular sensor is set to slightly influence the expression of either TetR or LacI.
- 2.
The concentrations of IPTG and DOX are set such that the bistable switch is entered into the fast switching regime, at a location in the phase space where a dilution curve (straight line to the origin) would be completely contained within the bistable region of the [IPTG]-[DOX] plane. Microorganisms containing the above constructs are put into the environment to be measured.
- 3.
As the inducers are diluted, the environment will influence the switch to either be biased toward the TetR- or LacI-dominated states. As the systems enter the disjoint bistable region of the IPTG-DOX plane due to dilution, the population distribution of the switch between the two states will depend on the measured parameters.
- 4.
The microorganisms may now be extracted from the environment without risk of perturbing their states, as in this region of the [IPTG]-[DOX] phase-space the bistable switch is locked into a state by the decreasing inducer concentrations.
- 5.
The states of the switches are measured.
In order to function properly, this approach requires the switch in question to posses the following properties. Firstly, the switch must be bistable in the absence of inhibitors. This is required to preserve the states of the switch throughout the extraction time. Secondly, the bistable region of the switch in inhibitor space must be roughly linear to fully contain a dilution curve. Finally, it is necessary that the bias introduced by the measurement step is small enough not to push the system out of the bistable phase and into monostability when inhibitors are diluted. In the case of the LITR switch, the first requirement is at first glance not met. However, as discussed above, it is possible to tune both the overall and the relative copy numbers of the two repressors to make the switch symmetric and bistable (with no switching) in the absence of inhibitors. The second requirement is related to the overall structure of the two-inhibitor system. In our case, a dilution curve can be plotted within the bistable region and exit into the fast switching region, but this may not necessarily be the case for other switch architectures. Finally, limiting the bias imposed by the molecular sensor on the system performance and properties could be addressed both by modifying the sensor itself, and by increasing the gene copy number in order to widen the bistable region. Overall, the above requirements put strong constraints on the choice of switch components.
Conclusions
The problem of creating a biomolecular measuring device with memory has previously been addressed by Bhomkar et al., where a recording device was made using a similar switch circuit [19]. However, its principle relied on the inhibition of bacterial division in order to prevent the reporter used for measuring from becoming diluted after the actual measurements were performed. Our method, if practical, does not require such radical precautions. With enough default switch stability, measurements could be recorded after an arbitrary amount of time or an arbitrary number of cell divisions. This stability would be ensured by the fact that each cell stores a binary value, that is relatively easy to preserve. The continuous measured value is then determined from the distribution cells between the two states of the recording device. Using stochastic and deterministic modelling with literature based reaction kinetics we demonstrated that a genetic switch based on the LacI and the TetR transcriptional repressors shows the necessary dynamics to function as such a recording device. We also derived a simplified model for a generalized asymmetric genetic switch, which showed dynamic behaviour similar to the specific case. Overall our modelling provides clues as to how a genetic switch with desirable stability and dynamic properties could be constructed in vivo, whether or not it will be used in a recording device.
Methods
Derivation of a reduced set of equations
In the following, we will simplify the equation system in Table 1 and reduce it to two coupled ordinary differential equations in the homodimers for genes 1 and 2: \([\mathrm {P}^{l}_{2}]\). In contrast to previous work [11], we will maintain that the two genes (including their transcription and translation processes) have different properties. This will allow for the study of a more realistic version of this bistable toggle switch, where the two coupled genes have different chemical kinetics.
We begin by making four assumptions: (i) There is only 1 copy of each genetic element, \(\sum _{ij}[\mathrm {D}_{ij}^{l}]=1\); (ii) All binding reactions to DNA as well as mRNA elongation are in steady state, \(\dot {[\mathrm {D}_{ij}^{l}]} =\dot {[\mathrm {E}^{l}]}=\dot {[\mathrm {M}^{l}]}=0\), where the dot over the variables stands for time derivative; (iii) The concentration of free RNA polymerase is constant, \(\dot {[\mathrm {R}]}= 0\); (iv) Dimerization reactions are in steady state, \(\dot {[\mathrm {P}^{i}_{2}]} = 0\).
with i={1,2}.
Here, we have introduced the dissociation parameter K _{ li }=q _{ li }/k _{ li } and made the approximation that (q _{ li }+α _{ lm })/k _{ li }≈q _{ li }/k _{ li }=K _{ li } for the dissociation parameters K _{ l3},K _{ l5}, and K _{ l7}.
Here, the definition of u _{ l } is based on assumption (iii) above, that \(\dot {[\mathrm {R}]}= 0\). We will interpret these gene specific parameters as: s _{ l } being a measure of promoter leakage, u _{ l } as the RNAp-promoter dissociation constant scaled by the concentration of free RNAp, [T_{ l }] as the dimensionless concentration of the active repressor molecule (homodimer), and r _{ l } as a measure of cooperativity in the repressor-DNA binding interaction.
Reduced system: derivation of simulation constants
Dimerisation of the CI-repressor — K _{1.1}
the parameter K _{1.1}=5.6 nM. Although in the same article they also showed that single-site mutations could cause a change in the dimer dissociation constant up to a value of K _{1.1}=2.7 M [20]. As this large interval has been shown experimentally it was assumed reasonable to keep using the value used in Ghim and Almaas (2009) [11], where K _{1.1}=10 nM. The dimerization constant of TetR was assumed to be equal to the one for CI for simplicity, yielding K _{2.1}=10 nM.
Binding of the CI repressor to DNA — K _{2}
Sauer reported in 1979 that K _{ d }=20 nM for CI_{2} dissociation from the O_{ R1} operator site. While this reference was not available, several others have cited the same value [21, 22]. The value was assumed to be the same for O_{ L1}, the operator site in the P_{ L }-promoter, which gives a K _{2}=20 nM.
The assumption that K _{1.2}=K _{2.2} is not really valid when comparing CI to the Tet repressor [23]. However, we chose to set them equal to one another for convenience. The parameters were much more carefully chosen for our realistic switch model. This gave K _{1.2}=K _{2.2}=20 nM.
Co-operative binding to the second operator site — r _{ l }
According to Johnson et al. [24] the co-operativity between two CI repressor dimers when binding to O_{ R1} and O_{ R2} will decrease the two dissociation constants with a factor of 2 and 12.5, respectively. As the model is made for sequential binding to the two binding sites these two factors are combined to one co-operative factor r _{2} = 25. In order to explain the operation of the Tn10 regulon, regulated by TetR, it is not necessary to include cooperativity to the binding of the repressor to the two repressor sites. However, we assume r _{2}=r _{1} for simplicity, leading to r _{ l }=25.
Concentration of free RNAp — [R]
The concentration of free RNAp was assumed to be the same as for previous modelling of the λ-phage switch [21], [R] =30 nM.
Binding of the RNAp to the P_{ L }-promoter — u _{ l }
Giladi et al. (1990) reported both the equilibrium constant for RNAp binding to the P_{ L }-promoter and the forward rate constant for the isomerization of closed to open RNAp-promoter complexes. The reported values were 8.94×10^{7} M and 4.38×10^{−3} s^{−1} respectively. These values can be used to estimate that the parameter K _{1.3}=11.2 nM [25]. The K _{1.3}-parameter was estimated using the relationship from Eq. (13). As the planned second promoter is P_{ L t e t−O } the value was assumed to be equal for the second promoter as well, which gives K _{2.3}=11.2 nM. Given our obtained RNAp concentration and Eq. (10), we get u _{ l }=3.0.
Remaining parameters — s and σ
The remaining parameters were chosen as in Ghim and Almaas (2009), giving a monomer to dimer lifetime ratio σ=10 and the leakage parameter s=0.01.
LITR system: Derivation of rate constants
Protein and mRNA decay — γ _{ p },γ _{ ml }, and γ _{ mt }
The protein degradation in an E. coli cell in exponential growth phase has been shown to be dominated by dilution due to cell division [12]. Using a division rate of once every 30 minutes, we obtain γ _{ p }=3.9×10^{−4} s^{−1}. Using the directly measured lacI mRNA half-life obtained by Bernstein et al. (2002) we determined the lacI mRNA degradation constant to be γ _{ ml }=3.0×10^{−3} s^{−1}. Combining the mRNA half-life of lacZ from [26] with the fact that tetR mRNA should be degraded three times as fast as that of lacZ [27], we obtained γ _{ mt }=3.9×10^{−3} s^{−1}.
Transcription initiation/elongation — α _{ le }, and α _{ lt }
The transcription initiation rate at both the lacI and tetR promoters was assumed to be equal to that at the λ P _{ R } promoter. We used a value measured at 37 °C [28].
Estimates of the lacI transcription elongation rate in E. coli range from 28 to 89 bp/s depending on the medium it is grown in [29]. Here we employ a very conservative estimate of 20 bp/s. Since the transcription elongation is carried out by a well-conserved molecular machinery, the same rate can be applied to both the lacI and tetR genes. Considering the lengths of the mRNA transcripts, we get the respective elongation constants of α _{ le }=6.5×10^{−2} s^{−1} and α _{ te }=1.1×10^{−1} s^{−1}.
mRNA translation to protein — α _{ lt }, and α _{ tt }
The translation rate constant of lacI was previously measured to be α _{ lt }=1.0×10^{−2} s^{−1} [30]. As no such constant was available for tetR, we chose to assume similar translation rate per amino acid. Scaling by the relative length of the tetR transcript gave us α _{ tt }=1.7×10^{−2} s^{−1}.
Protein di- and tetramerization — k _{1},q _{1}, and q _{ l2}
The association constants of the dimerization and tetramerization reactions were derived from the minimal diffusion limited estimate of the search time of a LacI tetramer for a specific site on the DNA of an E. coli cell [16]. Since no data was available for TetR, both TetR and LacI were assumed to dimerize at the same rate k _{1}=1.5×10^{−2} μm^{3}. Using the data provided in [31, 32] we derived equilibrium constants for LacI dimerization and tetramerization, from which we obtained the respective dissociation constants (using the value of k _{1} above) to be q _{1}=8.8×10^{−5} s^{−1}, and q _{ l2}=8.8×10^{−8} s^{−1}.
Repressor-inducer binding — k _{ l5},q _{ l5},k _{ t2} and q _{ t2}
The LacI-IPTG association constant was derived from the rate of dissociation of LacI from DNA after IPTG addition [16], and found to be k _{ l5}=1.4×10^{−6} μm^{3}. The dissociation constant was calculated using k _{ l5} and the respective equilibrium constant [33], giving q _{ l5}=8.2×10^{−3} s^{−1}. The TetR-doxycycline association and dissociation constants were directly measured by Kedracka et al. (2005). The values were k _{ t2}=2.4×10^{−4} μm^{3} and q _{ t2}=1.9×10^{−2} s^{−1} [34]. The constant k _{ t2} was also used as association constant for the binding of doxycycline to DNA-bound TetR. Similarly k _{ l5} was used for binding of IPTG to DNA-bound LacI (0.5k _{ l5} for reactions with two outcomes). Also q _{ l5} was used as a dissociation constant of IPTG from the DNA-LacI-IPTG complex.
DNA-LacI binding — k _{ l3},q _{ l3}, and q _{ l4}
From Levandoski et al. (1996) we know that the active species of LacI is the tetramer LacI_{4}, while the dimer is ineffective at repressing transcription. It is usual that a LacI-controlled promoter has two binding sites. One LacI_{4} is able to bind both sites. Our model goes as follows: If an IPTG molecule is bound to a doubly DNA-bound LacI_{4}, then it can displace it from one site but not from the other, but if it binds to a singly DNA-bound LacI_{4}, it would either displace it or not, depending on the dimer it binds to. LacI_{4} could exist doubly or singly bound to IPTG or to DNA as well as singly bound to both. Reactions where IPTG displaces LacI_{4} from DNA completely are assumed to be irreversible.
From the single molecule measurements of Elf et al. (2004), we calculated the association constant of LacI_{4} to DNA to be k _{ l3}=3.0×10^{−3} μm^{3}. We assumed that if the tetramer would be bound to one IPTG molecule, it would only have one dimer available for DNA binding, thus reducing this binding rate by a factor of two. We also assumed that if already bound to one binding site on the DNA, it would find the other one much faster, thus here we used 200k _{ l3}. Combining k _{ l3} with the equilibrium constant for LacI-DNA binding [32], we found the dissociation rate constant to be q _{ l3}=5.4×10^{−7} s^{−1}. It is worth mentioning that there exists a variety of different operator sites to which lacI can bind, and some studies suggest over 100 times variability in equilibrium constants between these [18]. We chose the comparatively strong binding site O ^{ s y m } for our model. It is obvious that the choice of operator would directly affect the simulation results of this specific system. Having a large degree of asymmetry would was in our view better for demonstrational purposes. However, we also tested our circuit with the O ^{1} operator site [18], yielding very similar results. Finally, we assumed the equilibrium constant to hold for binding to the second site as well, giving q _{ l4}=3.0×10^{−3} s^{−1}.
DNA-TetR binding — k _{ t3} and q _{ t3}
The active species of TetR is a dimer TetR_{2}. It can bind to a single site on the DNA, however usually TetR-controlled promoters have two TetR binding sites. We assumed that both of the binding sites have the same kinetic parameters, which were fortunately already determined [35]. They are k _{ t3}=4.8×10^{−3} μm^{3} and q _{ t3}=2.1×10^{−3} s^{−1}. Inducer doxycycline would irreversibly displace a TetR_{2} from DNA.
RNAp binding — k _{ b1},q _{ b },k _{ l b2},k _{ t b2} k _{ l b3} and k _{ t b3}
The kinetics of σ ^{70} RNAp binding to λ P _{ R } promoter has been determined. We assume that the obtained dissociation rate constant is the same whether repressors are bound to the promoter or not, while the association constant changes depending on the type and number of repressor molecules. The kinetic constants for a free promoter are k _{ b1}=4.3×10^{−3} μm^{3} and q _{ b }=4.5×10^{−1} s^{−1} [28].
We assumed that a singly DNA-bound LacI_{4} (whether bound to IPTG or not) would change the RNAp association rate as much as a bound dimer would. Using expression data for tetramerization mutants [31], we estimated that singly and doubly bound LacI_{4} would yield k _{ l b2}=1.2×10^{−5} μm^{3} and k _{ l b3}=3.6×10^{−7} μm^{3} respectively. Using the equilibrium constants for RNAp binding to a TetR-regulated promoter determined by Meier et al. (1988), we found k _{ t b2}=5.7×10^{−6} μm^{3} and k _{ t b3}=8.6×10^{−6} μm^{3} [36].
Stochastic simulation details
The stochastic simulations were performed using the software StochKit2 [37]. This provided us with a fast way of simulating large sets of stochastic equations for extended simulation periods. The stability and switching rate calculations were performed using our own custom software. The criterion for a switch was that the difference between the numbers of the active species had to cross the value of zero. The criterion for bistability was that when initiating the simulation with an excess of either protein, the histograms (of the difference in the number of the active species, e.g. [LacI_{4}]−[TetR_{2}]) would not intersect. The criterion for monostability was that for both initial conditions >95 % of the histogram values would be on the same side. Otherwise the system would be considered jointly bistable (that is, switching is possible within the given simulation time period of 10^{7}s). When generating the heat plots in Fig. 3 the simulation time was set to be 10^{7} s, and the number of samples to 10^{7}. In Fig. 2 the kinetics was much faster. Therefore the simulation time was set to 10^{3} s while the number of samples was kept the same. In Fig. 4(b-d) each line corresponds to a simulation, with the same time and sampling as in Fig. 3. The Supplementary Materials contains the initial set of equations used to derive the reduced model of the generalized asymmetric switch (Additional file 1).
Abbreviations
LITR-switch, LacI-TetR genetic switch; IPTG, isopropyl β-D-1-thiogalactopyranoside; DOX: doxycycline.
Declarations
Acknowledgements
Not applicable.
Funding
EA and NM would like to thank ERASysApp (WineSys) and Norwegian Research Council (Grant 245160), and EA and SHJ would like to thank Dept. of Biotechnology at NTNU for financial support. CMG was supported by the National Research Foundation grant funded by Korean Government (Grants 2010-0013582 and 2010-0028684).
Availability of data and materials
All necessary information is contained in the manuscript and the Supporting Material.
Authors’ contributions
CMG and EA conceived the study. EA, NM, and SHJ designed the study. NM and SHJ performed the simulations. EA and NM analyzed the data. CMG, EA and NM wrote the manuscript. All authors have read and approved the final version of the manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent to publish
Not applicable.
Ethics
Not applicable.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
- Leisner M, Kuhr JT, Radler JO, Frey E, Maier B. Kinetics of genetic switching into the state of bacterial competence. Biophys J. 2009; 96:1178–88.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang L, Walker BL, Iannaccone S, Bhatt D, Kennedy PJ, Tse WT. Bistable switches control memory and plasticity in cellular differentiation. PNAS. 2009; 106:6638–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Eldar A, Elowitz MB. Functional roles for noise in genetic circuits. Nature. 2010; 467(7312):167–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Brophy JAN, Voigt CA. Principles of genetic circuit design. Nat Methods. 2014; 11(5):508–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Nandagopal N, Elowitz M. Synthetic biology: integrated gene circuits. Science. 2008; 333:1244.View ArticleGoogle Scholar
- Gardner TS, Cantor CR, Collins JJ. Construction of a genetic toggle switch in Escherichia coli. Nature. 2000; 403:339–42.View ArticlePubMedGoogle Scholar
- Deans TL, Cantor CR, Collins JJ. A tunable genetic switch based on rnai and repressor proteins for regulating gene expression in mammalian cells. Cell. 2007; 130:363–72.View ArticlePubMedGoogle Scholar
- Green AA, Silver PA, Collins JJ, Yin P. Toehold switches: De-novo-designed regulators of gene expression. Cell. 2014; 159:925–39.View ArticlePubMedPubMed CentralGoogle Scholar
- Shopera T, Henson WR, Ng A, Lee YJ, Ng K, Moon TS. Robust, tunable genetic memory from protein sequestration combined with positive feedback. Nucleic Acids Res. 2015. [doi:10.1093/nar/gkv936].
- Ghim CM, Almaas E. Genetic noise control via protein oligomerization. BMC Syst Biol. 2008; 2(1):94.View ArticlePubMedPubMed CentralGoogle Scholar
- Ghim CM, Almaas E. Two-component genetic switch as a synthetic module with tunable stability. Phys Rev Lett. 2009; 103:028101.View ArticlePubMedGoogle Scholar
- Taniguchi Y, Choi PJ, Li G, Chen H, Babu M, Hearn J, Emili A, Xie XS. Quantifying E, coli proteome and transcriptome with single-molecule sensitivity in single cells. Science. 2011; 329:533–8.View ArticleGoogle Scholar
- Spruill LS, McDermott PJ. Role of the 5’-untranslated region in regulating translational efficiency of specific mrnas in adult cardiocytes. FASEB J. 2009; 23:2879–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Farasat I, Kushawaha M, Collens J, Easterbrook M, Guido M, Salis HM. Efficient search, mapping, and optimization of multi-protein genetic systems in diverse bacteria. Mol Syst Biol. 2014; 10:731.View ArticlePubMedPubMed CentralGoogle Scholar
- Bonde MT, Pedersen M, Klausen MS, Jensen SI, Wulff T, Harrison S, Nielsen AT, Herrgård MJ, Sommer MO. Predictable tuning of protein expression in bacteria. Nat Methods. 2016; 13:233–6.View ArticlePubMedGoogle Scholar
- Elf J, Li G, Xie XS. Probing transcription factor dynamics at the single-molecule level in a living cell. Science. 2004; 316:1191–4.View ArticleGoogle Scholar
- Clontech Laboratories Inc.tet-off and tet-on gene expression systems user manual. Mountain View, CA: Clontech Laboratories, Inc; 2012.Google Scholar
- Garcia HG, Phillips R. Quantitative dissection of the simple repression input–output function. PNAS. 2011; 108:12173–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Bhomkar P, Materi W, Wishart DS. The bacterial nanorecorder: Engineering E. coli to function as a chemical recording device. Plos ONE. 2011; 6:27559.View ArticleGoogle Scholar
- Burz DS, Beckett D, Benson N, Ackers GK. Self-assembly of bacteriophage λ cI repressor: Effects of single-site mutations on the monomer-dimer equilibrium. Nat Struct Biol. 1994; 7(6):449–55.Google Scholar
- Arkin A, Ross J, McAdams HH. Stochastic kinetic analysis of developmental pathway bifurcation in phage λ-infected Escherichia coli cells. Genetics. 1998; 149:1633–48.PubMedPubMed CentralGoogle Scholar
- Hawley DK, Johnson AD, McClure WR. Functional and physical characterization of transcription initiation complexes in the bacteriophage λ oR region. J Biol Chem. 1985; 260(14):8618–26.PubMedGoogle Scholar
- Hillen W, Berens C. Mechanisms underlying expression of tn10 encoded tetracycline resistance. Annu Rev Microbiol. 1994; 48:345–69.View ArticlePubMedGoogle Scholar
- Johnson AD, Meyer BJ, Ptashne M. Interactions between dna-bound repressors govern regulation by the λ phage repressor. Proc Natl Acad Sci U S A. 1979; 76(10):5061–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Giladi H, Gottesman M, Oppenheim AB. Integration host factor stimulates the phage λ pL promoter. J Mol Biol. 1990; 213:109–21.View ArticlePubMedGoogle Scholar
- Bernstein JA, Khodursky AB, Lin PH, Lin-Chao S, Cohen SN. Global analysis of mrna decay and abundance in Escherichia coli at single-gene resolution using two-color fluorescent dna microarrays. PNAS. 2002; 99:9697–702.View ArticlePubMedPubMed CentralGoogle Scholar
- Baumeister R, Flache P, Gabain A, Hillen W. Lack of a 5’ non-coding region in tn1721 encoded tetr mrna is associated with a low efficiency of translation and a short half-life in Escherichia coli. Nucleic Acids Res. 1991; 19:4595–600.View ArticlePubMedPubMed CentralGoogle Scholar
- Saecker RM, Tsodikov OV, McQuade KL, Schlax PE, Capp MW, Record MT. Kinetic studies and structural models of the association of E. coli σ ^{70} rna polymerase with the λ pR promoter: Large scale conformational changes in forming the kinetically significant intermediates. J Mol Biol. 2002; 319:649–71.View ArticlePubMedGoogle Scholar
- Vogel U, Jensen KF. The rna chain elongation rate in Escherichia coli depends on the growth rate. J Bacteriol. 1994; 176:2807–13.PubMedPubMed CentralGoogle Scholar
- Pedersen S. Eschericia coli ribosomes translate in vivo with variable rate. EMBO J. 1984; 3:2895–8.PubMedPubMed CentralGoogle Scholar
- Dong F, Spott S, Zimmermann O, Kisters-Woike B, MuÈller-Hill B, Barker A. Dimerisation mutants of lac repressor. i. a monomeric mutant, l251a, that binds lac operator dna as a dimer. J Mol Biol. 1999; 290:653–66.View ArticlePubMedGoogle Scholar
- Levandoski MM, Tsodikov OV, Frank DE, Melche rSE, Saecker RM, Record MT. Cooperative and anticooperative effects in binding of the first and second plasmid osym operators to a laci tetramer: Evidence for contributions of non-operator dna binding by wrapping and looping. J Mol Biol. 1996; 260:697–717.View ArticlePubMedGoogle Scholar
- Sochor MA. Kinetic studies and structural models of the association of E. coli σ ^{70} rna polymerase with the λ pR promoter: large scale conformational changes in forming the kinetically significant intermediates. PeerJ. 2014; 498:10–7717.Google Scholar
- Kedracka-Krok S, Gorecki A, Bonarek P, Wasylewski Z. Kinetic and thermodynamic studies of tet repressor-tetracycline interaction. Biochemistry. 2005; 44:1037–46.View ArticlePubMedGoogle Scholar
- Sushma A, Singh V, Singh AK, Mittal M, Datt M, Subramani B, Kumaran S. Equilibrium binding and kinetic characterization of putative tetracycline repressor family transcription regulator fad35r from mycobacterium tuberculosis. FEBS J. 2012; 279:3214–28.View ArticleGoogle Scholar
- Meier I, Wray LV, Hillen W. Differential regulation of the tn 10-encoded tetracycline resistance genes teta and tetr by the tandem tet operators o1 and o2. EMBO J. 1988; 7:567–72.PubMedPubMed CentralGoogle Scholar
- Sanft KR, Wu S, Roh M, Fu J, Lim RK, Petzold LR. Stochkit2: software for discrete stochastic simulation of biochemical systems with events. Bioinformatics. 2011; 27:2457–8.View ArticlePubMedPubMed CentralGoogle Scholar