Team:Imperial College/SingleCell

<!DOCTYPE html>
Single Cell Model
Click on the circuit to discover how we modeled each section

Single Cell Modelling Overview

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

Once built, the model was first used to test and compare the time taken for four different growth regulators (GP0.4, GP2, Leucine B Auxotrophy and Chloramphenicol acetyl transferase antibiotic resistance) that we shortlisted in order to determine which of them would work the fastest.

This allowed us to optimize our assembly experiments allowing us to achieve a faster route to a working prototype circuit. We decided to focus our attention on the GP2 and GP0.4 systems as they were shown to work faster.

We performed sensitivity analysis on each of our circuit designs. We did these for the parameters that we can change in lab (transcriptions rates via promoter strength, translation rates via RBS strength, copy numbers and degradation rates via the inclusion of degradation tags).

GP2
[Insert sensitivity analysis here]

GP0.4
[Insert sensitivity analysis here]

Our next process was to create a framework in which we could balance our circuit in silico. To do this we ran parameter sweeps for each of the transcription rates, translation rates and degradation rates that were indicated by the sensitivity analysis.

Transcription rates
4 pics

Translation rates
4 pics

Copy numbers
5 pics

Degradations numbers
4 pics

Communication module
We constructed four quorum systems that we considered viable choices for our system. We designed the model for the Rhl and Cin systems [Reference Bennett paper] as they have been previously shown to operate with minimal crosstalk.

Comparator module
Using STAR 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 we 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 STAR as was possible.
[RNA structure Diagram of ANTISTAR and STAR ANTISTAR]


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 [REFERENCE].

data?

Growth regulator module
We modelled 4 different growth regulator systems in silico in order to assess the speed and effectiveness of each case.

Auxotrophy (LeuB)

Antibiotic resistance (Chloramphenicol resistance)

Gene product 2

Gene Product 0.4