<!DOCTYPE html>
The first stage of our modelling process was to construct a single cell in silico model of our circuit. Our model was built using mass action kinetics in Simbiology (Matlab toolbox) and built up reaction by reaction.
We first separated the models into 3 modules: Quorum Communication, STAR-antiSTAR Comparator and Growth Regulation.
We constructed the four quorum systems that we considered viable choices for our system (cin, rhl, lux and las) to allow us to directly compare the expected behaviour and plan our growth module experiments accordingly. We designed the overall model for the Rhl and Cin systems (Chen et al., 2015) as they have been previously shown to operate with minimal crosstalk.
We used numbers obtained from Chen et al for C4 and C14 production. This is a high level production term that ignores parts of the central dogma. We made this assumption due to the limited data on the enzymatic kinetics of the autoinducer synthases.
The Rhl system consists of an autoinducer synthase RhlI (that produces C4 AHL) and a transcriptional regulator RhlR that dimerises to activate a pRhl promoter which activates STAR transcription.
The Cin system consists of an autoinducer synthase CinI (that produces C14 AHL) and a transcriptional regulator CinR that dimerizes to activate a pCin promoter which induces Anti-STAR production.
Communication module:
The results indicate that the species within this module reach steady states at time points unique to that species.
Figure 1: Production of C4 AHLs and C14 AHLs against time
Figure 2: Production of C4R and C14R regulatory proteins against time
Figure 3: C4:C4R and C14:C14R complex formation against time
Using STAR (Short Transcription Activating RNA) technology, we were able to develop a novel method of comparing the sizes of two populations from their quorum signal concentrations.
We used RNAstruct developed by Matthews Lab to help aid the development of the ANTISTAR. This software allowed us to determine the secondary structure and free energy to optimize the way in which our ANTISTAR sequence was designed. This was done so that our ANTISTAR sequence would have as high an affinity to the STAR sequence as was possible.
Figure 4: Secondary structure of our ANTISTAR sequence
To calculate the kinetics of the RNA interactions that occur within this module we adapted a method developed by Eric Winfree known as DNA strand displacement kinetics (Zhang and Winfree, 2009).
Comparator module:
Induced production of STAR and ANTISTAR was shown to be much higher than their basal production.
Figure 5: Basal vs Induced STAR production
Figure 6:Basal vs induced AntiSTAR production
STAR:ANTI STAR Complex formation occurs continuously whereas STAR: STAR Target complex formation is dependant on the copy number of the STAR target (in this case it was 30 represented by a psB3k3 plasmid backbone.
Figure 7:STAR:AntiSTAR Complex formation
Figure 8:STAR:STAR Target Complex formation
We constructed 4 different models for each of our growth regulators.
LeuB codes for an enzyme in the biosynthetic pathway of leucine in E. coli. The gene has been used before as a control method in co-culture. When using LeuB, we would need to use a strain that is auxotrophic for leucine and a growth medium lacking leucine. This system would require an inverter for it to operate so we incorporated a Tet inverter into the model.
LeuB module:
The figure shows successful repression of LeuB by the Tet system.
Figure 9: Effect of tet repressor on the production of LeuB
Antibiotic resistance (Chloramphenicol resistance)
Cat, encoding Chloramphenicol acetyltransferase, which is an enzyme that confers resistance to the antibiotic chloramphenicol. We decided to use chloramphenicol because it is a bacteriostatic rather than a bactericidal antibiotic. When using Cat, chloramphenicol would be added to the growth medium. This system would require an inverter for it to operate so we incorporated a Tet inverter into the model.
Chloramphenicol resistance module:
Figure 10: Effect of tet repressor on the production of CAT
Gp2 is a gene from the E. coli bacteriophage T7 phage, which slows down cell growth by binding reversibly to the E. coli RNA polymerase complex, thus inhibiting transcription. T7 phage infection is characterised by the hindrance of bacterial growth, and Gp2 has been suggested as a potential antimicrobial agent.
Reactions
GP 2 module:
Simulation Results
Figure 11: Basal vs induced production of Gp2
Gene Product 0.4
Gp0.4 is another T7 phage gene which inhibits growth by binding to the FtsZ ring during mitosis, preventing cytokinesis (the final stage of cell division where the two daughter cells separate).
Reactions
GP 0.4 module:
Simulation Results
Figure 12: Basal vs Induced production of GP 0.4
TEXT GOES HERE
Works Cited
Chen, Y., Kim, J., Hirning, A., Josi, K. and Bennett, M. (2015). Emergent genetic oscillations in a synthetic microbial consortium. Science, 349(6251), pp.986-989.
Zhang, D. and Winfree, E. (2009). Control of DNA Strand Displacement Kinetics Using Toehold Exchange. J. Am. Chem. Soc., 131(47), pp.17303-17314.