Team:ETH Zurich/Modeling/Switch Circuit

From 2008.igem.org

Revision as of 14:27, 29 October 2008 by Nico (Talk | contribs)

Contents

Switch Circuit

Using a novel pulsing mechanism consisting of two signals – a start signal that initiates the synthesis of restriction enzymes, and a terminating signal, which switches gene expression off – we are able to delete genome fragments in vivo. To simulate this system we developed, a switch circuit which works (see Switching Curcuit) as follows: Restriction enzyme gene expression is under the control of LacI and can be induced by the addition of IPTG. In order to stop gene expression, the IPTG-sensitive LacI is replaced by IPTG-insensitive LacIS, which shuts restriction enzyme gene expression off again. The synthesis of LacIS is started by the addition of Tetracyclin to the system, which binds to the tet repressor TetR and thus de-represses the expression of the lacIS gene. This circuit is modeled by ca. 40 reactions and is simulated using ODE solver and stochastic simulations.

Circuit model.JPG

Implementation and Simulation

For the implementation, different steps having no significant effects on the aspired results have been neglected for the sake of simplicity. In the transcriptions the additional step involving the RNA-Polymerase and in the translations the step involving the ribosomes have not been taken into account. Furthermore, the effects of the dimerization of TetR as well as the impacts of dimerization and tetramerization of LacI and LacIIs have not been considered in the final implementation.

This simplified model still comprises more than 20 different species and over 30 kinetic reactions and has been implemented by using the Simbiology Toolbox in MATLAB.

diagram view of the model

To get any useful results form the model, we performed deterministic and stochastic simulations based on the Mass-Action-Kinetics. The stochastic simulations turned out to be computationally very exhaustive but generated no further significant information compared to the deterministic simulations.

Simulation Results

Sensitivity Analysis

We define the sensitivity as the change of the production of the desired fluorescence protein - which is the output of our system - depending on the change of the paramters.

Sensitivity analysis


Results and Discussion

The simulations show that our system actually should create a nice pulse-shaped expression of the fluorescent protein (FP). This expression can be started by inducing with IPTG and stopped by subsequent addition of tet into the medium. By tagging the protein it will be degraded much faster by the Clp protease, so that the overall concentration is bounded and, after activating the stop-signal, the remaining proteins disappear quickly.

Another fact, our simulations show, is that in order to get one single pulse, the tet-concentration inside the medium must not reach zero before all the IPTG has degraded too. Otherwise, there would still be IPTG in the system inhibiting the binding of LacI to the FP-promotor and leading to an unwanted expression of our protein of interest as the LacIIs degrades.
One way to overcome this problem is simply by inducing a much higher quantity of tet than IPTG, so that it takes longer for it to completely degrade or being washed away. Another way would be to wait a bit longer after the induction with IPTG, so that it has already partly vanished until stopping the FP-expression with tet.

Detailed Model

diffusion of IPTG

In order to switch on the circuit, we induce with IPTG. When IPTG is added into the medium it diffuses reversibly between the medium and the cells, where it is slowly degraded. In a chemostat extracellular IPTG is washed away.
Grafik left width=30%

diffusion of tet

The second inducer which is used in our system is tet. This one also diffuses into the cell when it is added into the medium and degrades slowly.
Grafik left width=30%

binding of TetR to LacIIs-promotor

TetR is constitutively expressed and binds to the LacIIs promotor, inhibiting its expression.
Grafik left width=30%

binding of LacI and LacIIs to GFP-promotor

LacI which is constitutively expressed and LacIIs which is under the control of a tet repressor both can bind to the GFP promotor.
Grafik left width=30%

Binding LacI to Pgfp.JPG

binding of tet to TetR

Grafik left width=30%

binding of IPTG to LacI

Grafik left width=30%

transcription and translation of LacI

Grafik left width=30%

Translation of LacI.JPG

transcription and translation of LacIIs

Grafik left width=30%

Translation of LacIIs.JPG

transcription and translation of tetR

Grafik left width=30%

Translation of tetR.JPG

transcription and translation of GFP

Grafik left width=30%

Translation of gfp.JPG

dimerization of tetR

Dimerization of tetR.JPG

dimerization and tetramerization of LacI and LacIIs

Dimerization of LacI.JPG

Dimerization of LacIIs.JPG

Tetramerization of LacI.JPG

Tetramerization of LacIIs.JPG



Parameters

In this section you can find all the parameters used in the simulation.

# Parameter name Value Units Reference/Comment
1 k_assoc(IPTG_LacI) 5.0 1/(mole*second) X
2 k_assoc(LacI) 5.0 1/(mole*second) X
3 k_assoc(LacIs) 5.0 1/(mole*second) X
4 k_assoc(tet) 5.0 1/(mole*second) X
5 k_assoc(tetR) 5.0 1/(mole*second) X
6 k_dec(IPTG) 0.0002 1/second X
7 k_dec(IPTG_ext) 0.001 1/second X
8 k_dec(LacI) 0.05 1/second X
9 k_dec(LacIs) 0.05 1/second X
10 k_dec(gfp) 0.25 1/second X
11 k_dec(mRNA_LacI) 0.05 1/second X
12 k_dec(mRNA_LacIs) 0.05 1/second X
13 k_dec(mRNA_gfp) 0.1 1/second X
14 k_dec(mRNA_tetR) 0.05 1/second X
15 k_dec(tetR) 0.05 1/second X
16 k_dec(tet) 0.0002 1/second X
17 k_dec(tet_ext) 0.001 1/second X
18 k_diff(IPTG) 0.1 1/second X
19 k_diff(tet) 0.1 1/second X
20 k_dissoc(IPTG_LacI) 1.0 1/second X
21 k_dissoc(LacI) 1.0 1/second X
22 k_dissoc(LacIs) 1.0 1/second X
23 k_dissoc(tet) 1.0 1/second X
24 k_dissoc(tetR) 1.0 1/second X
25 k_tl(LacI) 5.0 1/second X
26 k_tl(LacIs) 5.0 1/second X
27 k_tl(gfp) 5.0 1/second X
28 k_tl(tetR) 5.0 1/second X
29 k_tr(LacI) 0.1 1/second X
30 k_tr(LacIs) 0.1 1/second X
31 k_tr(gfp) 1.0 1/second X
32 k_tr(tetR) 1.0 1/second X

References

(1) "Spatiotemporal control of gene expression with pulse-generating networks", Basu et al., PNAS, 2004

(2) "Genetic circuit building blocks for cellular computation, communications, and signal processing", Weiss et al., Natural Computing, 2003

(3) "Predicting stochastic gene expression dynamics in single cells", Mettetal et al., PNAS, 2006

(4) "Engineered gene circuits", Hasty et al., Nature, 2002