Team:Marburg:Safety:Modelling

From 2014.igem.org

(Difference between revisions)
 
(23 intermediate revisions not shown)
Line 5: Line 5:
{{Team:Marburg/Template:StartContinue}}
{{Team:Marburg/Template:StartContinue}}
-
In order to predict the behaviour of our SURFkiller, we created a model of our system. For this we used MATLAB software environment. The aim was to simulate cellular protein synthesis in different situations, and based on this information predict the robustness of the SURFkiller. Modelling was based on the following parameters:
+
In order to predict the behaviour of our SURFkiller, we created a model of our system. For this we used MATLAB software environment.  
 +
The aim was to simulate cellular protein synthesis in different situations, and based on this information predict the robustness of the SURFkiller.  
 +
Modelling was based on the following equations using the parameters explained in Table 1.
<html>
<html>
<table>
<table>
 +
<caption>
 +
Table 1: Parameters used in the equations modelling the SURFkiller.
 +
</caption>
<colgroup>
<colgroup>
<col width="35%" />
<col width="35%" />
Line 47: Line 52:
<td class="t1">&Kappa;<sub>x</sub></td>
<td class="t1">&Kappa;<sub>x</sub></td>
<td class="t2">Dissociation (Equilibrium) constant</td>
<td class="t2">Dissociation (Equilibrium) constant</td>
-
</tr>
 
-
<tr>
 
-
<td class="t1">n<sub>x</sub></td>
 
-
<td class="t2">Hill coefficient</td>
 
</tr>
</tr>
<tr>
<tr>
Line 64: Line 65:
'''a) Production of Antiholin'''
'''a) Production of Antiholin'''
-
<html><img src="https://static.igem.org/mediawiki/2014/8/85/Mr_killswitch_eq_1.png" /></html>
+
<html><img src="https://static.igem.org/mediawiki/2014/3/30/Mr_modelling_eq_1.png" /></html>
-
<html><img src="https://static.igem.org/mediawiki/2014/2/2a/Mr_killswitch_eq_2.png" /></html>
+
<html><img src="https://static.igem.org/mediawiki/2014/c/c8/Mr_modelling_eq_2.png" /></html>
'''b) Production of Holin'''
'''b) Production of Holin'''
-
<html><img src="https://static.igem.org/mediawiki/2014/2/2b/Mr_killswitch_eq_3.png" /></html>
+
<html><img src="https://static.igem.org/mediawiki/2014/5/55/Mr_modelling_eq_3.png" /></html>
-
<html><img src="https://static.igem.org/mediawiki/2014/0/07/Mr_killswitch_eq_4.png" /></html>
+
<html><img src="https://static.igem.org/mediawiki/2014/d/d9/Mr_modelling_eq_4.png" /></html>
'''c) Production of TetR'''  
'''c) Production of TetR'''  
-
<html><img src="https://static.igem.org/mediawiki/2014/8/88/Mr_killswitch_eq_5.png" /></html>
+
<html><img src="https://static.igem.org/mediawiki/2014/6/60/Mr_modelling_eq_5.png" /></html>
-
<html><img src="https://static.igem.org/mediawiki/2014/6/6a/Mr_killswitch_eq_6.png" /></html>
+
<html><img src="https://static.igem.org/mediawiki/2014/d/d5/Mr_modelling_eq_6.png" /></html>
'''d) Production of the ribosome-hibernation factor YvyD'''  
'''d) Production of the ribosome-hibernation factor YvyD'''  
-
<html><img src="https://static.igem.org/mediawiki/2014/3/3b/Mr_killswitch_eq_7.png" /></html>
+
<html><img src="https://static.igem.org/mediawiki/2014/9/94/Mr_modelling_eq_8.png" /></html>
-
<html><img src="https://static.igem.org/mediawiki/2014/9/9a/Mr_killswitch_eq_8.png" /></html>
+
<html><img src="https://static.igem.org/mediawiki/2014/9/9e/Mr_modelling_eq_9.png" /></html>
'''e) Production of the ribosomal protein RpL5'''
'''e) Production of the ribosomal protein RpL5'''
-
<html><img src="https://static.igem.org/mediawiki/2014/c/c7/Mr_killswitch_eq_9.png" /></html>
+
<html><img src="https://static.igem.org/mediawiki/2014/5/57/Mr_modelling_eq_10.png" /></html>
-
<html><img src="https://static.igem.org/mediawiki/2014/f/f9/Mr_killswitch_eq_10.png" /></html>
+
<html><img src="https://static.igem.org/mediawiki/2014/6/6d/Mr_modelling_eq_11.png" /></html>
-
'''Challenging different SURFkiller scenarios in silico'''
+
We assumed the transcription rate of both promoters to be ''α<sub>0</sub>''=10 nM/min representing a strong promoter. The translation rate is set to ''β''=2.3 min<sup>-1</sup> to represent the average burst size of 10 proteins/mRNA for the assumed half-life time of five minutes resulting in the degradation rate ''δ<sub>mx</sub>''=log(2)/5 min<sup>-1</sup>. The proteins except for TetR are assumed to be stable and only diluted by cell division every 30 min resulting in the degradation rate ''δ<sub>P</sub>''=log(2)/30 min<sup>-1</sup>. The additional degradation of TetR due to the lva-tag is assumed to be ''δ<sub>P</sub>''=0.018 min<sup>-1</sup> (Andersen et al., 1998).
-
To understand behaviour and robustness of SURFkiller in our model, we simulated several scenarios which would pose a challenge for the SURF killer. The efficient system would perform its function under virtually any conditions, and do it in in shortest possible time period.
+
'''Challenging different SURFkiller scenarios <i>in silico</i>'''
-
'''Scenario 1: A SURFkiller-equipped GMO leaves the lab'''
+
To understand behaviour and robustness of SURFkiller in our model, we simulated several scenarios which would pose a challenge for the SURFkiller. The efficient and secure system would kill the cells under virtually any conditions as quickly as possible.
-
The organism leaves the lab environment. Concentration of the IPTG anti-repressor quickly drops to minimal levels, as a consequence of the degradation and diffusion through the cell membrane. In this case the critical variable is the concentration of the L5 essential ribosomal protein. The results are shown in graphs below.
+
'''Scenario 1: A SURFkiller-equipped GMO leaves the laboratory'''
-
<html><div class="figure"><img src="https://static.igem.org/mediawiki/2014/f/f1/Mr_modelling_c1.jpg" width="50%" />
+
The organism leaves the laboratory environment. Concentration of the IPTG anti-repressor quickly drops to minimal levels, as a consequence of the degradation and diffusion through the cell membrane. In this case the critical variable is the concentration of the RpL5 essential ribosomal protein. The results are shown in graphs below.
-
<span class="caption">At the time point T=200 min, the mRNA concentration sharply drops, production of the L5 protein stops, and its concentration in the cells starts to decrease. If we assume that a cell needs around 2000 functional ribosomes in order to survive the cell death occurs at latest when the level of the L5 protein drops below 2000 nM. A concentration of 2000 nM in the volume of a cell, which is about 1 fL, equals about 2000 molecules. From our simulation we can see that the current number of ribosomes is below the critical number of 2000 at latest 60 minutes after the cell leaves the lab environment even if we assume that every L5 molecule leads to the formation of a functional ribosome.
+
 
 +
<html><div class="figure" style="width:90%;">
 +
<img src="https://static.igem.org/mediawiki/2014/7/73/Scence_1.png"/>
 +
<span class="caption"><b>Figure 1:</b> mRNA concentration is plotted against time <b>(left)</b>. Concentration of RpL5 plotted against time <b>(right)</b>.  
 +
The time of outbreak is t=0 (red line). The cell is killed after 63 min at latest (green line).
</span></div></html>
</span></div></html>
 +
 +
At the time point t=200 min, the mRNA concentration sharply drops, production of the RpL5 protein stops and its concentration in the cells starts to decrease. If we assume that a cell needs around 2000 functional ribosomes in order to survive the cell death occurs at latest when the level of the RpL5 protein drops below 2000 nM. A concentration of 2000 nM in the volume of a cell, which is about 1 fL, equals about 2000 molecules. From our simulation we can see that the current number of ribosomes is below the critical number of 2000 at latest 60 minutes after the cell leaves the laboratory environment even if we assume that every RpL5 molecule leads to the formation of a functional ribosome.
'''Case 2: Mutations compromising modular functions of the SURFkiller'''
'''Case 2: Mutations compromising modular functions of the SURFkiller'''
-
The efficiency of a kill-switch can be compromised with mutation that may occur on one of the promoters used in the system. SURF killer is designed to remain robust even in these situations, incorporating a secondary toxin-antitoxin system (Holin-Antiholin) that balances promoter function in our system. In this case we simulated a situation where one of the LacI promoters gets constitutive. L5 protein is always produced in this case, and doesn't cause death of the cell. However, since the toxin is also under control of the same promoter it also gets produced.
+
The efficiency of a killswitch can be compromised with mutation that may occur in one of the promoters used in the system. SURFkiller is designed to remain robust even in these situations, incorporating a secondary toxin-antitoxin system (Holin-Antiholin) that balances the promoter function in our system. In this case we simulated a situation where one of the LacI promoters gets constitutive. The essential RpL5 protein is always produced in this case, and does not lead to cell death. However, since the toxin is also under control of the same promoter it also gets produced.
-
<html><div class="figure"><img src="https://static.igem.org/mediawiki/2014/7/7a/Mr_modelling_c2.jpg" width="50%" /><span class="caption">
+
<html><div class="figure" style="width:90%;"><img src="https://static.igem.org/mediawiki/2014/2/2b/Mr_modelling_scence_2.png" />
-
As the organism leaves the lab environment, levels of antitoxin production start to sink, widening the gap between the levels of toxin and antitoxin in the organism. The critical concentration of toxin in the cell, in the case of T4 Holin, is around 1000-3000 molecules. Since number of free Holin molecules dictates death of the cell, we concentrate on the difference between toxin and antitoxin levels in the cell. Even if we take the worst case, where 3000 free T4 Holin molecules are needed, our model shows that the cell lysis will occur at latest 30 minutes after it leaves the lab environment.
+
<span class="caption">
-
</span></div></html>
+
<b>Figure 2:</b> Plotted on the left is the concentration of Holin (red) and Antiholin (blue) against time.  
 +
The right graph shows the difference in concentration between both. The cells die 32 min (green line) after the outbreak (red line).</span></div></html>
-
A variation of this scenario could be when a mutation occurs after the bacteria leave the lab but before they die due to the lack of L5.
+
As the organism leaves the laboratory environment, the level of antitoxin production start to sink, widening the gap between the level of toxin and antitoxin in the organism. The critical concentration of toxin in the cell, in the case of T4-Holin, is around 1000-3000 molecules. Since the number of free Holin molecules dictates cell death, we concentrate on the difference between toxin and antitoxin levels in the cell. Even if we take the worst case, where 3000 free T4- Holin molecules are needed, our model shows that the cell lysis will occur at latest 30 minutes after it leaves the laboratory environment.
-
<html><div class="figure"><img src="https://static.igem.org/mediawiki/2014/3/33/Mr_modelling_c3.jpg" width="50%" /><span class="caption">
+
A variation of this scenario could be when a mutation occurs after the bacteria leave the laboratory but before they die due to the lack of RpL5.
-
Figure legend:????</span></div></html>
+
-
<html><div class="figure"><img src="https://static.igem.org/mediawiki/2014/9/92/Mr_modelling_c4.jpg" width="50%" /><span class="caption">
+
<html><div class="figure" style="width:40%;"><img src="https://static.igem.org/mediawiki/2014/2/21/Mr_modelling_scence_3.png" /><span class="caption">
-
The mutation occurs at the T=250 min in the simulation depicted in graphs above, and we can see how the concentration of both L5 protein and T4 Holin rebounds after the initial drop caused by the bacteria leaving the lab. Antiholin production is not affected by this mutation, as it is located on another module in the system. Regardless of the mutation, Holin-induced cell lysis still occurs at latest 30 minutes after the mutation. Otherwise, cell dies from the lack of the L5 protein, as in the first case.
+
<b>Figure 3: Concentration of RpL5 plotted against time showing that a mutation occurs before the concentration drops below a safe level.</b></span></div></html>
 +
 
 +
<html><div class="figure" style="width:90%;"><img src="https://static.igem.org/mediawiki/2014/e/ed/Mr_modelling_scence_4.png" /><span class="caption">
 +
<b>Figure 4:</b> On the left concentration of Holin (red) and Antiholin (blue) is plotted against time. On the right graph the
 +
difference between both is shown. With the mutation occurring 50 min after the outbreak (red line) the cells are dead 31 min later (green line).
</span></div></html>
</span></div></html>
-
'''Scenario 3: Rpl5 is lost from SURFkiller'''
+
The mutation occurs at the T=250 min in the simulation depicted in graphs above, and we can see how the concentration of both the RpL5 protein and the T4-Holin rebounds after the initial drop caused by the bacteria leaving the laboratory. The Antiholin production is not affected by this mutation, as it is located on another module in the system. Regardless of the mutation, Holin-induced cell lysis still occurs at latest 30 minutes after the mutation. Otherwise, cell dies from the lack of the RpL5 protein, as in the first case.
-
For our third scenario we decided to test the unlikely case where the L5 encoding gene is removed from the operon under control of the Lac-promoter, for instance by the means of homologous recombination. It is then possible that the levels of the L5 protein will be high enough in the cell while the difference between Holin and Antiholin is kept low, even if it escapes the lab. The third security layer of our SURF killer is designed for this case, and activated if all other security measures get compromised.
+
'''Scenario 3: RpL5 is lost from SURFkiller'''
-
<html><div class="figure"><img src="https://static.igem.org/mediawiki/2014/f/fe/Mr_modelling_c5.jpg" width="50%" /><span class="caption">
+
For our third scenario we decided to test the unlikely case where the RpL5 encoding gene is removed from the operon controlled by the Lac-promoter, for instance by the means of homologous recombination. It is then possible that the levels of the RpL5 protein will be high enough in the cell while the difference between Holin and Antiholin is kept low, even if it escapes the laboratory. The third security layer of our SURFkiller is designed for this case, and activated if all other security measures get compromised.
-
Outside of the lab environment TetR protein translation stops, and its concentration starts to sink, as shown with the red line in the graph above. As its concentration sinks below the threshold level (around 20 molecules), it is unable to repress the constitutive promoter and the yvyD gene gets expressed starting synthesis of the protein, shown in blue. As YvyD levels rise, they dimerize more and more ribosomes, which in the end result in a cell death. The short period of YvyD production at the beginning is due to the fact that level of TetR protein starts below 20 molecules at the very beginning of the simulation.  
+
 
-
</span></div></html>
+
<html><div class="figure" style="width:40%;"><img src="https://static.igem.org/mediawiki/2014/a/ae/Mr_graph5.png" /><span class="caption">
 +
<b>Figure 5: Concentration of TetR (red) and YvyD (blue) plotted against time.</b></span></div></html>
 +
 
 +
Outside of the laboratory environment the TetR production stops, and its concentration starts to sink, as shown in Figure 5. As its  
 +
concentration sinks it is unable to repress the constitutive promoter and the''yvyD'' transcription slowly starts resulting in the  
 +
production of YvyD, shown in red. As the YvyD level rises, more and more ribosomes dimerize, resulting in cell death. The short period  
 +
of YvyD production at the beginning is due to the fact that in the beginning of the simulation the concentration to TetR is assumed to
 +
be zero.
'''Sensitivity of the simulation on the parameters used'''
'''Sensitivity of the simulation on the parameters used'''
-
In order to make sure that our simulation adequately describes real-world scenarios we ran a calculation on how our results depend on different choice of parameters, primarily transcription rates of mRNA and translation rates of proteins. Choice and calculation of these parameters is one of the most challenging tasks in modelling process, so this simulation sheds light on how would our SURF killer function, even with non-optimal parameters.
+
In order to make sure that our simulation adequately describes real-world scenarios we ran a calculation on how our results depend on different choice of parameters, primarily transcription rates of mRNA and translation rates of proteins. Choice and calculation of these parameters is one of the most challenging tasks in modelling process, so this simulation sheds light on how our SURFkiller would function, even with non-optimal parameters.
-
<html><div class="figure"><img src="https://static.igem.org/mediawiki/2014/9/99/Mr_modelling_c6.jpg" width="50%" /><br /><img src="https://static.igem.org/mediawiki/2014/1/1c/Mr_modelling_c7.jpg" /><span class="caption">
+
<html><div class="figure" style="width:90%;">
-
As we can see from the graphs above, if we vary the standard translation rate of 2.3 min-1 there is comparatively small change in the time needed for the cell death to occur. If the translation rate is below 0.6 the cells cannot survive even under lab conditions because the steady state number of ribosomes is below 2000. Similar is applicable to variations in transcription rate. If the assumed transcription rate of 10 nM/min doubles, the expected survival time of the cell is around 25 minutes longer.
+
<img src="https://static.igem.org/mediawiki/2014/9/93/Mr_modelling_scence_6.png" width="50%" /><span class="caption">
 +
<b>Figure 6: Time till the cells die plotted against translation or transcription rate showing the values (green line) used in the calculation of the above scenarios.</b>
</span></div></html>
</span></div></html>
 +
As we can see from the graphs above, if we vary the standard translation rate of 2.3 min <sup>-1</sup> there is comparatively small change in the time needed for the cell death to occur. If the translation rate is below 0.6 the cells cannot survive even under laboratory conditions because the steady state number of ribosomes is below 2000. Similar is applicable to variations in transcription rate. If the assumed transcription rate of 10 nM/min doubles, the expected survival time of the cell is around 25 minutes longer.
 +
 +
The modelling of our SURFkiller showed that the system indeed works and leads to a quick decrease of RpL5 upon the deprivation of IPTG. Also the further layers of security against mutations that might occur have been shown to be working. What we learned though is that it might be better to use several copies of ''yvyD'' to increase the steady state protein concentration, so that the ribosome number decreases by dimerization quickly below an essential number for survival. Our analysis of the sensitivity of the system suggest that we established a robust circuit design concerning the translation rate, which means that changes of this rate lead only to a small change in time until the cells die. However the transcription rate should be increased for example by using a stronger promoter in order to make the system even more robust.
<html><hr /></html>
<html><hr /></html>
-
Wang I., Smith DL., Young R., (2000). HOLINS: The Protein Clocks of Bacteriophage Infections. Annu. Rev. Microbiol. 2000. 54:799–825.
 
-
Young R., Bläsi U., (1995). Holins: forms and function in bacteriophage lysis. FEMS Microbiology Reviews 17 (1995) 191-205.
+
Andersen, J. B., Sternberg, C., Poulsen, L. K., Bjorn, S. P., Givskov, M., & Molin, S. (1998). New unstable variants of green fluorescent protein for studies of transient gene expression in bacteria. Applied and Environmental Microbiology, 64(6), 2240–6. Retrieved from http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=106306&tool=pmcentrez&rendertype=abstract
 +
 
 +
Lauber, M.A., Running, W.E., and Reilly, J.P. (2009) B. subtilis ribosomal proteins: structural homology and post-translational modifications. J Proteome Res 8: 4193–206 http://www.ncbi.nlm.nih.gov/pubmed/19653700. Accessed October 6, 2014.
 +
 
 +
Savva, C.G., Dewey, J.S., Deaton, J., White, R.L., Struck, D.K., Holzenburg, A., and Young, R. (2008) The holin of bacteriophage lambda forms rings with large diameter. Mol Microbiol 69: 784–793 http://www.ncbi.nlm.nih.gov/pubmed/18788120. Accessed October 14, 2014.
-
Lauber MA., Running WE., Reilly JP., (2009). B. subtilis Ribosomal Proteins: Structural Homology and Post-Translational Modifications. J. Proteome Res. 8 (2009) 4193-4206.
+
Wang, I.N., Smith, D.L., and Young, R. (2000) Holins: the protein clocks of bacteriophage infections. Annu Rev Microbiol 54: 799–825 http://www.ncbi.nlm.nih.gov/pubmed/11018145. Accessed October 13, 2014.
-
Christos G. Savva, Jill S. Dewey, John Deaton, Rebecca L. White, Douglas K. Struck, Andreas Holzenburg, Ry Young. (2008). The holin of bacteriophage lambda forms rings with large diameter. Molecular Microbiology, 69 (4), 784-793. doi: 10.1111/j.1365-2958.2008.06298.x
+
Young, R., and Bläsi, U. (1995) Holins: form and function in bacteriophage lysis. FEMS Microbiol Rev 17: 191–205 http://www.ncbi.nlm.nih.gov/pubmed/7669346. Accessed October 14, 2014.
{{Team:Marburg/Template:End}}
{{Team:Marburg/Template:End}}

Latest revision as of 02:38, 18 October 2014

Modelling cellular behaviour and robustness of the SURFkiller

In order to predict the behaviour of our SURFkiller, we created a model of our system. For this we used MATLAB software environment. The aim was to simulate cellular protein synthesis in different situations, and based on this information predict the robustness of the SURFkiller. Modelling was based on the following equations using the parameters explained in Table 1.

Table 1: Parameters used in the equations modelling the SURFkiller.
Parameter Description
[mx] mRNA concentration
α0,x Maximum transcription rate
βx Protein synthesis rate
δmx mRNA degradation rate
δp Standard dilution rate of the protein due to cell division
δR Degradation rate of TetR due to LVA tag
Κx Dissociation (Equilibrium) constant
[Χ] Protein concentration

The mathematical equations we based our model on are given as:

a) Production of Antiholin

b) Production of Holin

c) Production of TetR

d) Production of the ribosome-hibernation factor YvyD

e) Production of the ribosomal protein RpL5

We assumed the transcription rate of both promoters to be α0=10 nM/min representing a strong promoter. The translation rate is set to β=2.3 min-1 to represent the average burst size of 10 proteins/mRNA for the assumed half-life time of five minutes resulting in the degradation rate δmx=log(2)/5 min-1. The proteins except for TetR are assumed to be stable and only diluted by cell division every 30 min resulting in the degradation rate δP=log(2)/30 min-1. The additional degradation of TetR due to the lva-tag is assumed to be δP=0.018 min-1 (Andersen et al., 1998).

Challenging different SURFkiller scenarios in silico

To understand behaviour and robustness of SURFkiller in our model, we simulated several scenarios which would pose a challenge for the SURFkiller. The efficient and secure system would kill the cells under virtually any conditions as quickly as possible.

Scenario 1: A SURFkiller-equipped GMO leaves the laboratory

The organism leaves the laboratory environment. Concentration of the IPTG anti-repressor quickly drops to minimal levels, as a consequence of the degradation and diffusion through the cell membrane. In this case the critical variable is the concentration of the RpL5 essential ribosomal protein. The results are shown in graphs below.

Figure 1: mRNA concentration is plotted against time (left). Concentration of RpL5 plotted against time (right). The time of outbreak is t=0 (red line). The cell is killed after 63 min at latest (green line).

At the time point t=200 min, the mRNA concentration sharply drops, production of the RpL5 protein stops and its concentration in the cells starts to decrease. If we assume that a cell needs around 2000 functional ribosomes in order to survive the cell death occurs at latest when the level of the RpL5 protein drops below 2000 nM. A concentration of 2000 nM in the volume of a cell, which is about 1 fL, equals about 2000 molecules. From our simulation we can see that the current number of ribosomes is below the critical number of 2000 at latest 60 minutes after the cell leaves the laboratory environment even if we assume that every RpL5 molecule leads to the formation of a functional ribosome.

Case 2: Mutations compromising modular functions of the SURFkiller

The efficiency of a killswitch can be compromised with mutation that may occur in one of the promoters used in the system. SURFkiller is designed to remain robust even in these situations, incorporating a secondary toxin-antitoxin system (Holin-Antiholin) that balances the promoter function in our system. In this case we simulated a situation where one of the LacI promoters gets constitutive. The essential RpL5 protein is always produced in this case, and does not lead to cell death. However, since the toxin is also under control of the same promoter it also gets produced.

Figure 2: Plotted on the left is the concentration of Holin (red) and Antiholin (blue) against time. The right graph shows the difference in concentration between both. The cells die 32 min (green line) after the outbreak (red line).

As the organism leaves the laboratory environment, the level of antitoxin production start to sink, widening the gap between the level of toxin and antitoxin in the organism. The critical concentration of toxin in the cell, in the case of T4-Holin, is around 1000-3000 molecules. Since the number of free Holin molecules dictates cell death, we concentrate on the difference between toxin and antitoxin levels in the cell. Even if we take the worst case, where 3000 free T4- Holin molecules are needed, our model shows that the cell lysis will occur at latest 30 minutes after it leaves the laboratory environment.

A variation of this scenario could be when a mutation occurs after the bacteria leave the laboratory but before they die due to the lack of RpL5.

Figure 3: Concentration of RpL5 plotted against time showing that a mutation occurs before the concentration drops below a safe level.

Figure 4: On the left concentration of Holin (red) and Antiholin (blue) is plotted against time. On the right graph the difference between both is shown. With the mutation occurring 50 min after the outbreak (red line) the cells are dead 31 min later (green line).

The mutation occurs at the T=250 min in the simulation depicted in graphs above, and we can see how the concentration of both the RpL5 protein and the T4-Holin rebounds after the initial drop caused by the bacteria leaving the laboratory. The Antiholin production is not affected by this mutation, as it is located on another module in the system. Regardless of the mutation, Holin-induced cell lysis still occurs at latest 30 minutes after the mutation. Otherwise, cell dies from the lack of the RpL5 protein, as in the first case.

Scenario 3: RpL5 is lost from SURFkiller

For our third scenario we decided to test the unlikely case where the RpL5 encoding gene is removed from the operon controlled by the Lac-promoter, for instance by the means of homologous recombination. It is then possible that the levels of the RpL5 protein will be high enough in the cell while the difference between Holin and Antiholin is kept low, even if it escapes the laboratory. The third security layer of our SURFkiller is designed for this case, and activated if all other security measures get compromised.

Figure 5: Concentration of TetR (red) and YvyD (blue) plotted against time.

Outside of the laboratory environment the TetR production stops, and its concentration starts to sink, as shown in Figure 5. As its concentration sinks it is unable to repress the constitutive promoter and theyvyD transcription slowly starts resulting in the production of YvyD, shown in red. As the YvyD level rises, more and more ribosomes dimerize, resulting in cell death. The short period of YvyD production at the beginning is due to the fact that in the beginning of the simulation the concentration to TetR is assumed to be zero.

Sensitivity of the simulation on the parameters used

In order to make sure that our simulation adequately describes real-world scenarios we ran a calculation on how our results depend on different choice of parameters, primarily transcription rates of mRNA and translation rates of proteins. Choice and calculation of these parameters is one of the most challenging tasks in modelling process, so this simulation sheds light on how our SURFkiller would function, even with non-optimal parameters.

Figure 6: Time till the cells die plotted against translation or transcription rate showing the values (green line) used in the calculation of the above scenarios.

As we can see from the graphs above, if we vary the standard translation rate of 2.3 min -1 there is comparatively small change in the time needed for the cell death to occur. If the translation rate is below 0.6 the cells cannot survive even under laboratory conditions because the steady state number of ribosomes is below 2000. Similar is applicable to variations in transcription rate. If the assumed transcription rate of 10 nM/min doubles, the expected survival time of the cell is around 25 minutes longer.

The modelling of our SURFkiller showed that the system indeed works and leads to a quick decrease of RpL5 upon the deprivation of IPTG. Also the further layers of security against mutations that might occur have been shown to be working. What we learned though is that it might be better to use several copies of yvyD to increase the steady state protein concentration, so that the ribosome number decreases by dimerization quickly below an essential number for survival. Our analysis of the sensitivity of the system suggest that we established a robust circuit design concerning the translation rate, which means that changes of this rate lead only to a small change in time until the cells die. However the transcription rate should be increased for example by using a stronger promoter in order to make the system even more robust.


Andersen, J. B., Sternberg, C., Poulsen, L. K., Bjorn, S. P., Givskov, M., & Molin, S. (1998). New unstable variants of green fluorescent protein for studies of transient gene expression in bacteria. Applied and Environmental Microbiology, 64(6), 2240–6. Retrieved from http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=106306&tool=pmcentrez&rendertype=abstract

Lauber, M.A., Running, W.E., and Reilly, J.P. (2009) B. subtilis ribosomal proteins: structural homology and post-translational modifications. J Proteome Res 8: 4193–206 http://www.ncbi.nlm.nih.gov/pubmed/19653700. Accessed October 6, 2014.

Savva, C.G., Dewey, J.S., Deaton, J., White, R.L., Struck, D.K., Holzenburg, A., and Young, R. (2008) The holin of bacteriophage lambda forms rings with large diameter. Mol Microbiol 69: 784–793 http://www.ncbi.nlm.nih.gov/pubmed/18788120. Accessed October 14, 2014.

Wang, I.N., Smith, D.L., and Young, R. (2000) Holins: the protein clocks of bacteriophage infections. Annu Rev Microbiol 54: 799–825 http://www.ncbi.nlm.nih.gov/pubmed/11018145. Accessed October 13, 2014.

Young, R., and Bläsi, U. (1995) Holins: form and function in bacteriophage lysis. FEMS Microbiol Rev 17: 191–205 http://www.ncbi.nlm.nih.gov/pubmed/7669346. Accessed October 14, 2014.