Team:Hamburg/Modeling/Networks
Genetic toggle switch for PezT expression
Our GroEL-PezT Biobrick shall enable heat-inducible cell lysis. Conducting heat induction would release our produced miRNA2911, which was triggered by light before. One major problem here is an activation of the GroEL promotor and thus PezT caused cell lysis at normal temperature (some kind of leakiness), therefore not reaching sufficient miRNA product. To stabilise our cell lysis event, we can employ a gene regulatory network called the genetic toggle switch [1]. It consists of a double negative feedback loop leading to a bistability network with memory function. This means the following:
- the network can only attain two possible steady states (concentrations of all components): PezT-expression-state versus no-PezT-expression state
- when one of the two possible steady states is once reached, this network state is naturally stable (memory function)
- each state is induced by certain, different sets of start conditions
- it is possible to switch the network state, when using inducer molecules (IPTG, cIts)
The construction of the network can be seen in the figure below. The repressors lacI and clts conduct the double negative feedback, and can themselves be inactivated through IPTG and heat, respectively. The bistability arises from the fact that either lacI or clts, but never both at the same time, can have a high concentration, and each of them promotes its own production. The inducers (6h 2mM IPTG and 35min temperature increase from 32°C to 42°C) are only needed for setting the desired state once, then the state is maintained in a cell without external influence.
To analyse this network mathematically, a representation of it through ordinary differential equations (ODEs) can be obtained [1]. It models the repressor concentrations of lacI ($u$) and cIts ($v$). The kinetic parameters were measured in experiments.
$\dot{u} = \frac{\alpha_1}{1+v^{\beta}} - u$
$\dot{v} = \frac{\alpha_2}{1+\tilde{u}^{\gamma}} - v$
Here $\tilde{u}$ means
$\tilde{u} = \frac{u}{\big(1+\frac{IPTG}{K}\big)^{\eta}}$
representing the active fraction of lacI (IPTG unbound fraction) which can be deduced by a conservation law for lacI forms. Further, the parameters are
$\alpha_1 = 156.25 \qquad \alpha_2 = 15.6\\ \beta = 2.5 \qquad \gamma = 1\\ \eta = 2.0015 \qquad K = 2.9618 \cdot 10^{-5}$
DLE analysis
It is now interesting to see which start concentrations of lacI and cIts (or initial conditions for the ODE) would yield a certain network state. Direct finite-time Lyapunov Exponents (DLEs) are a calculation method based on ODEs for the prediction of those processes. They measure the rate of exponential separation between network trajectories [2], thus revealing regions of initial conditions leading to one certain network state and regions separating two different states. We implemented a code for DLE analysis in Matlab (R2015a). In the figure below, the general bistability of this network and the initial conditions leading to a certain network state (besides the threshold with higher DLEs) can be discerned.
Following our practical workflow of the induction of the no-cell-lysis-state with IPTG to ensure sufficient miRNA production, we need to know an IPTG value where no initial conditions exist yielding a cell lysis state. In other words: Examining the IPTG value where the network gets monostable. Then a stable miRNA production is possible and the cell lysis can be triggered by heat induction, when we want it. The figure below for an IPTG value of 5e-5 M shows the required monostability, the threshold moved to the dimension axis for lacI.
- Gardner, T. S., Cantor, C. R. & Collins, J. J. Construction of a genetic toggle switch in Escherichia coli. Nature 403, 339–342 (2000).
- Aldridge, B., Haller, G., Sorger, P. & Lauffenburger, D. Direct Lyapunov exponent analysis enables parametric study of transient signalling governing cell behaviour. IEE Proc - Syst Biol 153, 425 (2006).