SynDustry Fuse. Produce. Use.
Quantitative evolutionary stability analysis and design guidelines for kill switches
1. Concentration dynamics of kill switches
This chapter summarizes the first working phase of the iGEM Marburg 2016 modeling work: The continuous ODE model presented previously is used with the biologically relevant parameters to compute the concentrations' dynamics \(\vec{X}(t)\) for all times. This shows that the designs are comparable in terms of the toxin expression levels. The latter three designs are comparable in terms of the number of model parameters.
At first, the lab environment is simulated in order to evaluate the toxin expression during the kill switch's OFF state. After that, the organism is released into wild life conditions – no inducers present – and the toxin expression of the activated kill switch is compared to the OFF state. From there, a toxin threshold is found.
This sections serves as verification that the non-mutated kill switch designs work as intended.
(a) BNU China 2014
The BNU China 2014 kill switch (find design here) is simulated and the results are shown in figure 1. One observes a significant increase of the toxin concentration once the organism is released into wild life conditions which lacks an inducer. The lethal toxin threshold is defined as based on the simulation results. This is the mean of the kill switch's OFF and ON state. It is consistently drawn as dashed red line in the following.
Interestingly, one can observe a noticeable delay between the inducer turned off and the toxin turning on. This is due to the signal traveling through the network at finite speed.
The next kill switch designs represent different generic kill switch topologies and abstract from specific biological substances to substances that behave just like their biological counterparts. This enables to focus on the bare kill switch design (see algorithm description for further details here).
(b) Parallel (OR)
This design utilizes two redundant toxin regulating branches (find design here). Being connected by an OR gate, this provides a backup mechanism in case one branch is disabled through mutations. The result can be found in figure 2.
(c) Serial
The Serial kill switch design is simulated (find design here). Figure 3 shows the kill switch typical behavior.
Since this design introduces a longer serial cascade than the BNU China 2014 design, the delay between inducer deactivation and toxin expression is longer.
(d) Parallel (AND)
A parallel design that connects both branches with an AND gate is studied (find design here). Figure 4 shows that both branches need to be functional to trigger the kill switch's ON state.
As evident from figure 4, the deactivation of one of the two branches leads to the kill switch being in the OFF state permanently - thus, the branches are indeed linked by an AND gate.
2. Evolutionary stability
For different incubation times and mutation rates, the escape probability of the BNU China 2014 kill switch design is studied. Figure 5 shows the results from which two intuitive results are confirmed:
- The higher the mutation rate, the larger the escape probability: As the likelihood for a mutation to occur increases, the probability of the linear cascade to be destroyed through one of these mutations increases as well.
- Longer incubation times lead to higher escape probabilities: This is due to more possible attempts for mutations to destroy the kills switch.
(b) Time scaling by mutation rate
Scaling the incubation time with the mutation rate lets the escape probability curves collapse onto one master curve. This can be explained: Hitting a certain threshold takes more incubation time when a low mutation rate is applied but it takes less incubation generations when working with a high mutation rate. The master curve for the BNU China 2014 kill switch can be found in figure 6.
Therefore, the mutation rate is arbitrary for quantifying the evolutionary stability of a kill switch design. Computations have been performed for \(\eta\in\{0.4081875 \cdot 10^{-3}, 0.816375 \cdot 10^{-3}, 1.63275 \cdot 10^{-3}\}\) but only \(\eta=1.63275 \cdot 10^{-3}\) is shown in the following. Other mutation rates are equivalent and the results show the same behavior.
(c) Design guidelines for a kill switch
We are providing design guidelines for a kill switch that is robust against evolutionary stability and show quantitative evidence for our suggestion. A kill switch design can be understood as a network with the toxin being a node that is controlled through regulation pathways.
The BNU China 2014 kill switch design uses a single regulation pathway to control its toxin MazF. From intuition it is clear that a kill switch with multiple parallel regulation pathways for its toxin – with only one required to control the toxin expression – should be less prone to being destroyed by mutations. Let a mutation destroy one of the parallel regulation pathways: There are still other pathways that back up the destroyed one and ensure a working kill switch.
We utilize three comparable kill switches that differ in their network topology. Figure 7 shows the escape probabilities for the Parallel (OR), Serial and Parallel (AND) designs. Parallel (OR) – implementing a topology that allows destroyed branches to be backed up – shows by far the smallest escape probability. Where a single mutation is capable of destroying the Parallel (AND) and Serial design, a destroyed pathway can be compensated by the functional regulation pathway for the Parallel (OR) design.
Concluding, we suggest to implement a redundant toxin regulation for kill switches. This significantly increases the evolutionary stability and decreases the escape probability of genetically modified organisms into the environment. Biotechnological applications ranging from small to industrial scale can be made safer with having this in mind.
However, introducing additional topology such as multiple redundant toxin regulations increases a kill switches complexity. There are more points of contact for mutations to destroy the topology. Translating this in the model we use, more complex kill switches have more ODE parameters that can be mutated. Therefore, the escape probability is expected to somehow scale with a kill switch's complexity.
This remaining question of a trade-off between kill switch complexity and desired network topology to find a robust kill switch is studied in the last chapter.
(d) Comparison between BNU China 2014 and Parallel (OR) kill switch
We compare the escape probabilities of the BNU China 2014 design and the hypothetical Parallel (OR) kill switch. The latter introduces on the one hand redundant toxin regulation but on the other hand an increased complexity. That manifests in almost double the number of ODE parameters in comparison to the BNU China 2014 design.
Figure 8 shows that the stability-optimized Parallel (OR) kill switch is more robust against evolution and mutations than the BNU China 2014 topology. The gain in evolutionary stability through the additional regulation units overcompensates the introduced complexity. Finally, this leads to a more robust kill switch when using a redundant toxin regulation.
This result has various applications in real life. Redundant toxin regulation not only ensures projects in Synthetic Biology to be safer but also enables longer incubation times until a certain threshold in escape probability is reached. Organisms can therefore be longer incubated without becoming evolutionary instable. This simplifies industrial applications and increases overall efficiency in production processes.
The code framework is available upon request (Martin Lellep, mailto:martin.lellep@physik.uni-marburg.de). We provide possible outlooks on the following website.