Team:Imperial College London/Project/Gene/Modelling

From 2011.igem.org




Module 3: Gene Guard

Containment is a serious issue concerning the release of genetically modified organisms (GMOs) into the environment. To prevent horizontal gene transfer of the genes we are expressing in our chassis, we have developed a system based on the genes encoding holin, anti-holin and endolysin. We are engineering anti-holin into the genome of our chassis, where it acts as an anti-toxin, and holin and endolysin on plasmid DNA. In the event of horizontal gene transfer with a soil bacterium, holin and endolysin will be transferred without anti-holin, rendering the recipient cell non-viable and effectively containing the Auxin Xpress and Phyto-Route genes in our chassis.




Modelling

I. INTRODUCTION

The gene guard device involves three proteins:

1) Holin is a protein forms pore complexes in the inner-membrane and this pores allow endolysin to move into periplasmic space.

2) Endolysin is an enzyme that could break down the cell wall and causes cell lysis.

3) Anti-holin binds to holin and prevent from forming pores.

Our system will have holin and endolysin genes together with axuin genes and chemoreceptor genes on the plasmids, and anti-holin genes on the genome under control of strong promoter. The production of anti-holin will prevent cell lysis of the genetic modified bacteria by deactivating holin. If the plasmid get horizontal transferred to a different cell without anti-holin on its genome(e.g. any wild type soil bacteria), the expression of holin and endolysin will induce cell lysis and hence prevent it from keeping the plasmid containing auxin gene.

II. OBJECTIVES

The key to the success of this device is the right ratio of holin and anti-holin production such that the anti-holin can deactivate all the holin produced in our genetic modified bacteria, and the holin expression level in wild type bacteria is high enough to induce cell lysis. As the gene expression is controlled by both promoter strength and ribosome binding site (RBS) strength, therefore modeling of gene guard is important for helping us to choose appropriate promoters and ribosome binding site (RBS) to inhibit and optimize cell behavior in our genetic modified cell and soil bacteria respectively.


III. DESCRIPTION

The modeling of gene guard consists two parts, one is the holin production in wild type bacteria and the other part is the deactivation of holin by anti-holin in genetic modified bacteria. Since all the genes in our system are constitutively expressed, thereby we only consider steady state.


1. Holin production in wild type soil bacteria

This part of modeling focused on deriving the relationships between promoter strength and RBS strength of holin on single plasmid. In order to model this, the following assumptions were made:

1) After plasmid transfer and duplication, there will be 50 copies of plasmids in one bacterium.

2) The literature indicates that cell lysis happens at holin concentration equals to 1000/cell[1].

Based on above assumptions, the kinetics of cell lysis in wild type bacteria was modeled using ordinary differential equations (ODEs) for the holin mRNA and holin (Equation 1).


2. Holin inhibition in Genetic Modified bacteria

In genetic modified bacteria the promoter strength and RBS strength of anti-holin on the genome should be stronger than promoter strength and RBS strength of holin on single plasmid such that all the holin produced by 300 copies of plasmids should be inhibited, where 300 is the maximum number of plasmids can be in single bacterium. In this system, the anti-holin can bind to holin to form dimer, which is defined as the deactivated holin. The mechanism in the genetic modified bacteria can be described in following ODEs in Equation 3.



IV. RESULTS

1. Holin production in wild type soil bacteria

At steady state (i.e. d[ ] /dt = 0), rearrange Equation 1 with [Holin] = 1000, we obtained an expression for P_(holin) K_(holin) as shown in Equation 2:


2. Holin inhibition in Genetic Modified bacteria

Again at steady state with d[]/dt = 0, we obtained an equation of holin concentration (Equation 4). By varying the ratio m in equation 4, we obtained that the holin concentration in genetic modified bacteria decreased to zero at ratio m ≥ 400 (see Figure 1). Then, we tested m = 400 with the anti-holin promoter we have with strength 13.1pg RNA/min/μg substrate DNA using ODEs shown in Equation 1 and 3, the result (Figure 2) shown that all the holin in the cell are inhibited as the holing concentration remains at 0 all the time, while the holin concentration in wild type bacteria reaches 1000 per cell after 5000 seconds.

Figure 1: The holin concentration against the ratio m

Figure 2: Evolution of anti-holin and hoin in genetic modified and wild type bacteria agiainst time(s)

V. MATLAB CODE



VI. REFERENCE

[1] Christos G. Savva, Jill S. Dewey, John Deaton, Rebecca L. White, Douglas K. Struck, Andreas Holzenburg and Ry Young. The holin of bacteriophage lambda forms rings with large diameter. Molecular Microbiology 69(4), 784–793. 2008.