Team:ETH Zurich/modeling/whole

From 2014.igem.org

(Difference between revisions)
(With Leakiness)
 
(179 intermediate revisions not shown)
Line 1: Line 1:
{{:Team:ETH Zurich/tpl/head|Whole cell model}}
{{:Team:ETH Zurich/tpl/head|Whole cell model}}
 +
{{:Team:ETH Zurich/tpl/read-more}}
 +
{{:Team:ETH Zurich/tpl/fortables}}
 +
<center>
<center>
-
{{:Team:ETH Zurich/tpl/scrollbutton|Model|green}}
+
{{:Team:ETH Zurich/tpl/scrollbutton|Model|blue}}
-
{{:Team:ETH Zurich/tpl/scrollbutton2|Species|blue}}
+
{{:Team:ETH Zurich/tpl/scrollbutton2|Species|green}}
{{:Team:ETH Zurich/tpl/scrollbutton3|Reactions|red}}
{{:Team:ETH Zurich/tpl/scrollbutton3|Reactions|red}}
-
{{:Team:ETH Zurich/tpl/scrollbutton3|Equations|green}}
+
{{:Team:ETH Zurich/tpl/scrollbutton|Equations|blue}}
-
{{:Team:ETH Zurich/tpl/scrollbutton2|Leakiness|red}}
+
{{:Team:ETH Zurich/tpl/scrollbutton2|Results|green}}
-
{{:Team:ETH Zurich/tpl/scrollbutton|Cross-talk|blue}}
+
{{:Team:ETH Zurich/tpl/scrollbuttontworows|Alternate|Design|red}}
</center>
</center>
Line 13: Line 16:
== Model ==
== Model ==
-
The whole cell model is the combination of the Quorum sensing, Integrase and XOR modules. The model shows the behaviour of the a single cell in response to incoming signals. The model enables us to understand the effect of leakiness, cross-talk and their combinations on the whole system.
+
The whole cell model is the combination of the Quorum sensing, Integrase and XOR modules. The model shows the behaviour of a single cell in response to incoming signals. The model enables us to understand the effect of leakiness, cross-talk and their combinations on the whole system. All the simulations here are described for a red cell .i.e. a cell producing LasI and GFP when it is ON.
-
<br/>
+
 
-
<br/>
+
 
 +
[[File:ETHZ_WholeCellRedCell.png|800px|center|thumb|'''Figure 1''' The four cases with different combinations of input for a red cell.]]
<html></article></html>
<html></article></html>
   
   
-
<html><article id="Species"></html>
+
<html><article id="Species" class='carousel'></html>
 +
 
== Chemical Species ==
== Chemical Species ==
{| class="wikitable"
{| class="wikitable"
|-
|-
-
! Name
+
! '''Name'''
! Description
! Description
|-
|-
-
|Lux-AHL
+
|'''LuxAHL'''
|30C6-HSL is an acyl homoserine lactone which mainly binds to LuxR.
|30C6-HSL is an acyl homoserine lactone which mainly binds to LuxR.
|-
|-
-
|LuxR
+
|'''LuxR'''
-
|Constitutively expressed regulator protein that can bind Lux-AHL and stimulate transcription of Bxb1.
+
|Constitutively expressed regulator protein that can bind LuxAHL and stimulate transcription of Bxb1.
|-
|-
-
|RLux
+
|'''RLux'''
-
|LuxR and Lux-AHL complex which can dimerize.
+
|LuxR and LuxAHL complex which can dimerize.
|-
|-
-
|DRLux  
+
|'''DRLux '''
|Dimerized form of RLux.
|Dimerized form of RLux.
|-
|-
-
|mRNA<sub>Bxb1</sub>
+
|'''mRNA<sub>Bxb1</sub>'''
|mRNA of the Bxb1 integrase being transcribed by the Lux promoter.
|mRNA of the Bxb1 integrase being transcribed by the Lux promoter.
|-
|-
-
|Bxb1
+
|'''Bxb1'''
|Serine integrase that can fold into two conformations - Bxb1a and Bxb1b. We chose to use a common connotation for both conformations - Bxb1.
|Serine integrase that can fold into two conformations - Bxb1a and Bxb1b. We chose to use a common connotation for both conformations - Bxb1.
|-
|-
-
|Las-AHL
+
|'''DBxb1'''
 +
|Dimerized form of Bxb1. We chose to use a common connotation for both homodimers, DBxb1a and DBxb1b.
 +
|-
 +
|'''LasAHL'''
|30C12-HSL is an acyl homoserine lactone which mainly binds to LasR.  
|30C12-HSL is an acyl homoserine lactone which mainly binds to LasR.  
|-
|-
-
|LasR
+
|'''LasR'''
-
|Constitutively expressed regulator protein that can bind Las-AHL and stimulate transcription of ΦC31.
+
|Constitutively expressed regulator protein that can bind LasAHL and stimulate transcription of ΦC31.
|-
|-
-
|RLas
+
|'''RLas'''
-
|LasR and Las-AHL complex which can dimerize.
+
|LasR and LasAHL complex which can dimerize.
|-
|-
-
|DRLas
+
|'''DRLas'''
|Dimerized form of RLas.
|Dimerized form of RLas.
|-
|-
-
|mRNA<sub>ΦC31</sub>
+
|'''mRNA<sub>ΦC31</sub>'''
|mRNA of the ΦC31 integrase being transcribed by the Lux promoter.
|mRNA of the ΦC31 integrase being transcribed by the Lux promoter.
|-
|-
-
|ΦC31
+
|'''ΦC31'''
|Serine integrase that can fold into two conformations - ΦC31a and ΦC31b. We chose to use a common connotation for both conformations - ΦC31.
|Serine integrase that can fold into two conformations - ΦC31a and ΦC31b. We chose to use a common connotation for both conformations - ΦC31.
|-
|-
-
|SI<sub>Bxb1</sub>
+
|'''DΦC31'''
 +
|Dimerized form of ΦC31. We chose to use a common connotation for both homodimers, DΦC31a and DΦC31b.
 +
|-
 +
|'''SI<sub>Bxb1</sub>'''
|Inactive DNA binding site for Bxb1. No dimer is bound to this site.
|Inactive DNA binding site for Bxb1. No dimer is bound to this site.
|-
|-
-
|SA<sub>Bxb1</sub>
+
|'''SA<sub>Bxb1</sub>'''
|Active DNA binding site for Bxb1. A dimer is bound to this site.
|Active DNA binding site for Bxb1. A dimer is bound to this site.
|-
|-
-
|SF<sub>Bxb1</sub>
+
|'''SF<sub>Bxb1</sub>'''
|Flipped DNA binding site for Bxb1. The site generated after recombination by integrase.
|Flipped DNA binding site for Bxb1. The site generated after recombination by integrase.
|-
|-
-
|SI<sub>ΦC31</sub>
+
|'''SI<sub>ΦC31</sub>'''
|Inactive DNA binding site for ΦC31. No dimer is bound to this site.
|Inactive DNA binding site for ΦC31. No dimer is bound to this site.
|-
|-
-
|SA<sub>ΦC31</sub>
+
|'''SA<sub>ΦC31</sub>'''
|Active DNA binding site for ΦC31. A dimer is bound to this site.
|Active DNA binding site for ΦC31. A dimer is bound to this site.
|-
|-
-
|SF<sub>ΦC31</sub>
+
|'''SF<sub>ΦC31</sub>'''
|Flipped DNA binding site for ΦC31. The site generated after recombination by integrase.
|Flipped DNA binding site for ΦC31. The site generated after recombination by integrase.
|-
|-
-
|T<sub>on,i<sub>
+
|'''T<sub>on,i<sub>'''
|The number of terminators which are blocking the transcription of GFP and LuxI/LasI initially.
|The number of terminators which are blocking the transcription of GFP and LuxI/LasI initially.
|-
|-
-
|T<sub>offBxb1</sub>
+
|'''T<sub>offBxb1</sub>'''
|The number of terminators turned off by recombination due to Bxb1. Favours the transcription of GFP and LuxI.
|The number of terminators turned off by recombination due to Bxb1. Favours the transcription of GFP and LuxI.
|-
|-
-
|T<sub>offΦC31</sub>
+
|'''T<sub>offΦC31</sub>'''
|The number of terminators turned off by recombination due to ΦC31. Favours the transcription of GFP and LuxI.
|The number of terminators turned off by recombination due to ΦC31. Favours the transcription of GFP and LuxI.
|-
|-
-
|T<sub>on,f</sub>
+
|'''T<sub>on,f</sub>'''
|The number of terminators blocking the transcription of GFP and LuxI/LasI after recombination by  Bxb1 and ΦC31. No transcription.
|The number of terminators blocking the transcription of GFP and LuxI/LasI after recombination by  Bxb1 and ΦC31. No transcription.
|-
|-
-
|mRNA<sub>GFP</sub>
+
|'''mRNA<sub>GFP</sub>'''
|mRNA for Green fluorescent protein which is produced when the cells are ON.
|mRNA for Green fluorescent protein which is produced when the cells are ON.
|-
|-
-
|GFP
+
|'''GFP'''
|Green fluorescent protein which is produced when the cells are ON.
|Green fluorescent protein which is produced when the cells are ON.
|-
|-
-
|mRNA<sub>LuxI</sub>
+
|'''mRNA<sub>LuxI</sub>'''
|mRNA for LuxI which is produced when the cells are ON.
|mRNA for LuxI which is produced when the cells are ON.
|-
|-
-
|LuxI
+
|'''LuxI'''
-
|Enzyme catalysing the production of Lux-AHL from SAM and ACP.
+
|Enzyme catalysing the production of LuxAHL from SAM and ACP.
|-
|-
-
|mRNA<sub>LasI</sub>
+
|'''mRNA<sub>LasI</sub>'''
|mRNA for LasI which is produced when the cell are ON.
|mRNA for LasI which is produced when the cell are ON.
|-
|-
-
|LasI
+
|'''LasI'''
-
|Enzyme catalysing  the production of Las-AHL from SAM and ACP.
+
|Enzyme catalysing  the production of LasAHL from SAM and ACP.
|}
|}
 +
 +
{{:Team:ETH_Zurich/tpl/topbutton|blue}}
 +
{{:Team:ETH_Zurich/tpl/rmbutton|blue|Species}}
<html></article></html>
<html></article></html>
-
<html><article id="Reactions"></html>
+
<html><article id="Reactions" class='carousel'></html>
 +
 
==Reactions==
==Reactions==
-
;For the Lux system:
 
$$ \begin{align}
$$ \begin{align}
-
&\rightarrow LuxR \\
+
\emptyset&\rightarrow LuxR \\
-
Lux-AHL+LuxR & \leftrightarrow RLux\\
+
LuxAHL+LuxR & \leftrightarrow RLux\\
RLux+RLux &\leftrightarrow DRLux\\
RLux+RLux &\leftrightarrow DRLux\\
DRLux+P_{luxOFF} & \leftrightarrow P_{luxON}\\
DRLux+P_{luxOFF} & \leftrightarrow P_{luxON}\\
P_{luxON}&\rightarrow P_{luxON}+mRNA_{Bxb1}\\
P_{luxON}&\rightarrow P_{luxON}+mRNA_{Bxb1}\\
mRNA_{Bxb1}&\rightarrow Bxb1\\
mRNA_{Bxb1}&\rightarrow Bxb1\\
-
AHL &\rightarrow \\
+
Bxb1 + Bxb1 &\leftrightarrow DBxb1 \\
-
LuxR &\rightarrow  \\
+
DBxb1 + SI_{Bxb1} & \leftrightarrow SA_{Bxb1}\\
-
RLux &\rightarrow\\
+
LuxAHL &\rightarrow \emptyset\\
-
DRLux &\rightarrow\\
+
LuxR &\rightarrow  \emptyset\\
-
mRNA_{Bxb1} &\rightarrow\\
+
RLux &\rightarrow \emptyset\\
-
Bxb1 &\rightarrow
+
DRLux &\rightarrow \emptyset\\
-
\end{align}$$
+
mRNA_{Bxb1} &\rightarrow \emptyset\\
-
 
+
Bxb1 &\rightarrow \emptyset\\
-
;For the Las system:
+
DBxb1 &\rightarrow \emptyset\\
-
\begin{align}
+
\\
-
&\rightarrow LasR \\
+
\emptyset&\rightarrow LasR \\
-
Las-AHL+LasR & \leftrightarrow RLas \\
+
LasAHL+LasR & \leftrightarrow RLas \\
RLas+RLas & \leftrightarrow DRLas\\
RLas+RLas & \leftrightarrow DRLas\\
DRLas+P_{LasOFF} & \leftrightarrow P_{LasON}\\
DRLas+P_{LasOFF} & \leftrightarrow P_{LasON}\\
P_{LasON}&\rightarrow P_{LasON}+mRNA_{\phi C31}\\
P_{LasON}&\rightarrow P_{LasON}+mRNA_{\phi C31}\\
-
Las-AHL &\rightarrow \\
+
mRNA_{\phi C31}&\rightarrow \phi C31\\
-
LasR &\rightarrow \\
+
\phi C31 + \phi C31 &\leftrightarrow D\phi C 31 \\
-
RLas &\rightarrow\\
+
D\phi C 31 + SI_{\phi C31} & \leftrightarrow SA_{\phi C31}\\
-
DRLas &\rightarrow\\
+
LasAHL &\rightarrow  \emptyset\\
-
mRNA_{\phi C31} &\rightarrow \\
+
LasR &\rightarrow \emptyset\\
-
\phi C31 &\rightarrow \\
+
RLas &\rightarrow \emptyset\\
-
\end{align}
+
DRLas &\rightarrow \emptyset\\
 +
mRNA_{\phi C31} &\rightarrow \emptyset\\
 +
\phi C31 &\rightarrow \emptyset\\
 +
D\phi C31 &\rightarrow \emptyset\\
 +
\\
 +
SA_{Bxb1}+SA_{Bxb1}+T_{on,i}& \rightarrow T_{offBxb1}+ SF_{Bxb1}+SF_{Bxb1}\\
 +
T_{offBxb1} &\rightarrow T_{offBxb1} + mRNA_{GFP} + mRNA_{LasI} \\
 +
SA_{\phi C31}+SA_{\phi C31}+T_{on,i}& \rightarrow T_{off\phi C31}+SF_{\phi C31}+SF_{\phi C31}\\
 +
T_{off\phi C31} &\rightarrow T_{off\phi C31} + mRNA_{GFP} + mRNA_{LasI} \\
 +
SA_{Bxb1}+SA_{Bxb1}+T_{off,\phi C31}& \rightarrow T_{on,f}+SF_{Bxb1}+SF_{Bxb1}\\
 +
SA_{\phi C31}+SA_{\phi C31}+T_{offBxb1}& \rightarrow T_{on,f}+SF_{\phi C31}+SF_{\phi C31}\\
 +
mRNA_{GFP} &\rightarrow GFP\\
 +
mRNA_{LasI} &\rightarrow LasI\\
 +
mRNA_{GFP} &\rightarrow \emptyset\\
 +
mRNA_{LasI}&\rightarrow \emptyset\\
 +
GFP &\rightarrow \emptyset\\
 +
LasI&\rightarrow \emptyset\\
 +
\\
 +
S_{LasI}+LasI & \rightarrow LasAHL\\
 +
\end{align}$$
 +
 
 +
 
 +
{{:Team:ETH_Zurich/tpl/topbutton|red}}
 +
{{:Team:ETH_Zurich/tpl/rmbutton|red|Reactions}}
<html></article></html>
<html></article></html>
-
<html><article id=""></html>
+
<html><article id="Equations" class='carousel'></html>
 +
 
 +
==Equations==
 +
The following equations were used for the whole cell model.
 +
The production of Bxb1 and ΦC31 is reduced to a hill function of RLux an RLas respectively. For details, check [https://2014.igem.org/Team:ETH_Zurich/modeling/qs Quorum sensing module].
 +
 
 +
The equation is for a system which produces LasI and GFP as outputs. The LasI catalyses the production of LasAHL from SAM and ACP. Since, we do not know the amount of ACP or SAM present in the cell, we assume the rate of production of LasAHL to depend only on LasI concentration.
 +
 
-
Applying mass action kinetic laws, we obtain the following set of differential equations.
 
$$\begin{align*}
$$\begin{align*}
-
\frac{d[Lux-AHL]}{dt} &= k_{-RLux}[R_{Lux}]-k_{RLux}[Lux-AHL][LuxR]-d_{Lux-AHL}[Lux-AHL]\\
+
\frac{d[LuxAHL]}{dt} &= -d_{LuxAHL}[LuxAHL]\\
-
\frac{d[LuxR]}{dt} &= \alpha_{LuxR} -k_{RLux}[Lux-AHL][LuxR] + k_{-RLux}[RLux] - d_{LuxR}[LuxR] \\
+
\frac{d[LuxR]}{dt} &= \alpha_{LuxR} -k_{RLux}[LuxAHL][LuxR] + k_{-RLux}[RLux] - d_{LuxR}[LuxR] \\
-
\frac{d[RLux]}{dt} &= k_{RLux}[Lux-AHL][LuxR] - k_{-RLux}[RLux]  - 2 k_{DRLux} [RLux]^2 + 2 k_{-DRLux} [DRLux] - d_{RLux} [RLux] \\  
+
\frac{d[RLux]}{dt} &=  k_{RLux}[LuxAHL][LuxR] - k_{-RLux}[RLux] - d_{RLux} [RLux] \\
-
\frac{d[DRLux]}{dt} &= k_{DRLux} [RLux]^2 - k_{-DRLux} [DRLux] - d_{DRLux} [DRLux] \\
+
\frac{d[mRNA_{Bxb1}]}{dt} &= L_{P_{Lux}} + \frac{K_{mRNA_{Bxb1}}[RLux]^2}{K_{mLux}^2 + [RLux]^2 }- d_{mRNA_{Bxb1}} [mRNA_{Bxb1}]\\
-
\frac{d[P_{LuxON}]}{dt} &= k_{P_{LuxON}} [P_{LuxOFF}][DRLux] - k_{-P_{LuxON}} [P_{LuxON}]\\
+
\frac{d[Bxb1]}{dt} &=  k_{Bxb1} [mRNA_{Bxb1}] -2 k_{DBxb1}[Bxb1]^2+ 2 k_{-DBxb1}[DBxb1] - d_{Bxb1}[Bxb1]\\
-
\frac{d[mRNA_{Bxb1}]}{dt} &= L_{P_{Lux}} + k_{mRNA_{Bxb1}} [P_{LuxON}] - d_{mRNA_{Bxb1}} [mRNA_{Bxb1}]\\
+
\frac{d[DBxb1]}{dt}&= k_{DBxb1}[Bxb1]^2-k_{-DBxb1}[DBxb1]-k_{SABxb1}[DBxb1][SI_{Bxb1}]+k_{-SABxb1}[SA_{Bxb1}]-d_{DBxb1}[DBxb1]\\
-
\frac{d[Bxb1]}{dt} &=  k_{mRNA_{Bxb1}} [mRNA_{Bxb1}] - d_{Bxb1}[Bxb1]\\
+
\frac{d[SI_{Bxb1}]}{dt}&=-k_{SABxb1}[DBxb1][SI_{Bxb1}]+k_{-SABxb1}[SA_{Bxb1}]\\
 +
\frac{d[SA_{Bxb1}]}{dt}&=k_{SABxb1}[DBxb1][SI_{Bxb1}]-k_{-SABxb1}[SA_{Bxb1}] - 2 k_{ToffBxb1} [SA_{Bxb1}]^2 [T_{on,i}] - 2 k_{-Toff\phi C31} [SA_{Bxb1}]^2 [T_{off\phi C31}]\\
 +
\frac{d[SF_{Bxb1}]}{dt}&=  2 k_{ToffBxb1} [SA_{Bxb1}]^2 [T_{on,i}] + 2 k_{-Toff\phi C31} [SA_{Bxb1}]^2 [T_{off\phi C31}]\\
 +
\\
 +
\frac{d[LasAHL]}{dt} &= -d_{LasAHL}[LasAHL] + DF * k_{LasAHL} [LasI]\\
 +
\frac{d[LasR]}{dt} &= \alpha_{LasR} -k_{RLas}[LasAHL][LasR] + k_{-RLas}[RLas] - d_{LasR}[LasR] \\
 +
\frac{d[RLas]}{dt} &=  k_{RLas}[LasAHL][LasR] - k_{-RLas}[RLas]  - d_{RLas} [RLas] \\
 +
\frac{d[mRNA_{\phi c31}]}{dt} &= L_{P_{Las}} + \frac{K_{mRNA_{\phi C31}}[RLas]^2}{K_{mLas}^2 + [RLas]^2 } - d_{mRNA_{\phi c31}} [mRNA_{\phi c31}]\\
 +
\frac{d[\phi c31]}{dt} &=  k_{\phi c31} [mRNA_{\phi c31}] -2 k_{D\phi c31}[\phi c31]^2+ 2 k_{-D\phi c31}[D\phi c31] - d_{\phi c31}[\phi c31]\\
 +
\frac{d[D\phi c31]}{dt}&= k_{D\phi c31}[\phi c31]^2-k_{-D\phi c31}[D\phi c31]-k_{SA\phi c31}[D\phi c31][SI_{\phi c31}]+k_{-SA\phi c31}[SA_{\phi c31}]-d_{D\phi c31}[D\phi c31]\\
 +
\frac{d[SI_{\phi c31}]}{dt}&=-k_{SA\phi c31}[D\phi c31][SI_{\phi c31}]+k_{-SA\phi c31}[SA_{\phi c31}]\\
 +
\frac{d[SA_{\phi C31}]}{dt}&=k_{SA\phi C31}[D\phi C31][SI_{\phi C31}]-k_{-SA\phi C31}[SA_{\phi C31}] - 2 k_{Toff\phi C31} [SA_{\phi C31}]^2 [T_{on,i}] - 2 k_{-ToffBxb1} [SA_{\phi C31}]^2 [T_{offBxb1}]\\
 +
\frac{d[SF_{\phi C31}]}{dt}&=  2 k_{Toff\phi C31} [SA_{\phi C31}]^2 [T_{on,i}] + 2 k_{-ToffBxb1} [SA_{\phi C31}]^2 [T_{offBxb1}]\\
 +
\\
 +
\frac{d[T_{on,i}]}{dt}&= -k_{ToffBxb1}[T_{on,i}][SA_{Bxb1}]^2 -k_{Toff\phi C31} [T_{on,i}][SA_{\phi C31}]^2 \\
 +
\frac{d[T_{offBxb1}]}{dt}&=k_{ToffBxb1} [T_{on,i}][SA_{Bxb1}]^2 -  k_{-ToffBxb1} [T_{offBxb1}] [SA_{\phi C31}]^2\\
 +
\frac{d[T_{off\phi C31}]}{dt}&=k_{Toff\phi C31} [T_{on,i}][SA_{\phi C31}]^2 -  k_{-Toff\phi C31} [T_{offBxb1}] [SA_{Bxb1}]^2\\
 +
\frac{d[T_{on,f}]}{dt}&=  k_{-ToffBxb1} [T_{offBxb1}] [SA_{\phi C31}]^2 + k_{-Toff\phi C31} [T_{offBxb1}] [SA_{Bxb1}]^2\\
 +
\\
 +
\frac{d[mRNA_{LasI}]}{dt}&=  k_{mRNA_{LasI}} (\frac{[T_{offBxb1}]}{\theta + [T_{offBxb1}]} +  \frac{[T_{off\phi C31}]}{\theta + [T_{off\phi C31}]}) - d_{mRNA_{LasI}}[mRNA_{LasI}]\\
 +
\frac{d[LasI]}{dt}&=  k_{LasI} [mRNA_{LasI}] - d_{LasI}[LasI]\\
 +
\frac{d[mRNA_{GFP}]}{dt}&=  k_{mRNA_{GFP}} (\frac{[T_{offBxb1}]}{\theta + [T_{offBxb1}]} +  \frac{[T_{off\phi C31}]}{\theta + [T_{off\phi C31}]}) - d_{mRNA_{GFP}}[mRNA_{GFP}]\\
 +
\frac{d[GFP]}{dt}&=  k_{GFP} [mRNA_{GFP}] - d_{GFP}[GFP]\\
\end{align*}$$
\end{align*}$$
-
The same holds true for the Las system.
+
{{:Team:ETH_Zurich/tpl/topbutton|blue}}
 +
{{:Team:ETH_Zurich/tpl/rmbutton|blue|Equations}}
<html></article></html>
<html></article></html>
-
<html><article id="Parameters"></html>
+
<html><article id="Results"></html>
-
== Characterization: K<sub>mLux</sub> and K<sub>mLas</sub> ==
+
== Dynamics ==
-
=== Data ===
+
=== Ideal case ===
-
For the Quorum sensing module we used established experimentally determined parameters for the rate of formation of RLux (reference). Since, in the literature the other parameters were estimated or fitted to their data, we decided to determine the parameters specific to our system. Hence, we used our data for the remaining parameters. Our data was mainly a transfer function of normalized GFP concentration as a function of input Lux-AHL concentrations. (link to data)
+
[[File:ETHZ_IdealCase.png|center|800px|thumb|'''Figure 2''' The four cases for the whole cell under ideal conditions. On [https://2014.igem.org/Team:ETH_Zurich/project/overview#Implementation_in_E._coli our project overview], you can see  see in action how the circuit responds to different inputs.]]
-
=== Assumptions ===
 
 +
Under ideal conditions, there is no leakiness or cross-talk. The Figure 2 shows the production of GFP as a function of four different combinations of inputs in an ideal whole cell. No GFP is produced when there is no LasAHL and LuxAHL or when both are present. GFP is produced only when one of the input AHLs are present, thus emulating an XOR system. In the case with 0.5 nM of LuxAHL and no LasAHL as input, the rate of production of GFP reduces after six hours. This could be attributed to the delayed production of LasAHL from LasI. In the case with only LasAHL as input, there is a positive feedback of LasAHL.
-
From the original set of reactions, we reduce the rate of production of mRNA<sub>Bxb1</sub> as a Hill function of RLux instead of Mass action kinetics in terms of P<sub>LuxON</sub>  and P<sub>LuxOFF</sub>.
 
 +
[[File:ETHZ_00TerminatorwithLeakiness.png|center|600px|thumb|'''Figure 3''' The dynamics of flipping of the terminators when the cell receives 0.5 nM LuxAHL and 0 nM of LasAHL and produces GFP and LasI]]
-
Assumption A
 
-
We assumed that the dimerization of RLux to DRLux is quick. Quasi steady state approximation (QSSA)  as follows
+
In the case above, the cell receives 0.5 nM of LuxAHL and 10 nM LasAHL. This results in production of Bxb1 and ΦC31 respectively. Each terminator is flanked by binding sites for each of these integrases. These integrases flip the terminators by binding to their respective sites. Since, both integrases are present, the terminator is flipped twice thus eventually turning the system OFF. This is summarized in the figure above.Then, they flip the terminators again thus turning the system OFF.  
-
$$\frac{d[DRLux]}{dt} = k_{DRLux} [RLux]^2 - k_{-DRLux} [DRLux] - d_{DRLux} [DRLux] \approx 0\\$$
 
-
Assumption B
+
At first the systems are ON either due to flipping by ΦC31 or Bxb1. This results in the production of LasI which catalyses the production of LasAHL. The LasAHL further increases the production of ΦC31 which flips the already flipped terminators by Bxb1. The remaining Bxb1 flips the already flipped terminators by ΦC31. Thus, the system eventually turns OFF.Even though there is more Bxb1 flipped terminators than those flipped by ΦC31. The dynamics of T<sub>oni</sub>, T<sub>offBxb1</sub>, T<sub>offΦC31</sub> and T<sub>onf</sub> is summarised in the figure above.
-
Further, from literature, we found that DRLux is specific to DNA and the dissociation constant is low (k<sub>m</sub> = 0.1nM) {Reference}.  Therefore, we using QSSA again,
 
-
$$\frac{d[P_{LuxON}]}{dt} = k_{P_{LuxON}} [P_{LuxOFF}][DRLux] - k_{-P_{LuxON}} [P_{LuxON}] \approx 0\\$$
 
-
Solving, we get the rate of production of mRNA<sub>Bxb1</sub> as
+
[[File:ETHZ_XORWholeCell.png|center|500px|thumb|'''Figure 4''' Predicted XOR behaviour for the whole cell model without leakiness]]
-
$$\frac{d[mRNA_{Bxb1}]}{dt} = L_{P_{Lux}} + \frac{K_{mRNA_{Bxb1}}RLux^2}{K_{mLux} + RLux^2 }- d_{mRNA_{Bxb1}} [mRNA_{Bxb1}]\\$$
+
-
where
 
-
$$K_{mLux} = \frac {k_{-P_{LuxON}}}{k_{P_{LuxON}}}.\frac {k_{-DRLux} + d_{DRLux}}{k_{DRLux}}$$
+
The figure above, shows the predicted XOR gate performance by a single cell model with no basal leakiness for different concentrations of LasAHL and LuxAHL.
-
is a lumped parameter which we fitted to our data.
+
=== With Leakiness ===
-
Similarly, lumped parameter K<sub>mLas</sub> was derived for the las system and fitted to a transfer function of normalized GFP concentration as a function of input Las-AHL.
 
-
<!--Further,
+
[[File:ETHZ_Leakiness.png|center|800px|thumb|'''Figure 5''' The four cases for the whole cell with basal leakiness]]
-
$$K_{mRNA_{Bxb1}} = k_{mRNA_{Bxb1}} P_{LuxTOT}$$ which we assumed. -->
+
-
=== Parameter fitting ===
 
-
We used MEIGO Toolbox to fit the parameters to the transfer functions.
+
From our experimental data, we observed some basal leakiness for P<sub>Lux</sub> and P<sub>Las</sub> even after using riboregulators. From the model we see that, this small basal leakiness is amplified downstream. The basal leakiness results in production of integrases which further act on the XOR module and cause the switching of the terminator. Thus, there is GFP and LasI produced, and the LasI produced further catalyses the production of more LasAHL. Thus, we observe some GFP even without inputs.
-
=== Range of validity of the assumptions ===
+
<!--[[File:ETHZ_00TerminatorwithLeakiness.png|center|500px|thumb|'''Figure 6''' No inputs and only basal leakiness]]-->
-
These assumptions hold true for all input Lux-AHL and Las-AHL concentrations.  
+
-
<html></article></html>
+
<!--The figure above summarizes the predicted effect of basal leakiness on the flipping of the terminator. The basal leakiness results in production of Bxb1 and ΦC31 which result in flipping of the terminator. In this case, since the cell produces LasI there is increased production of LasAHL. The LasAHL produced induces the production of ΦC31 which further causes flipping of all terminators flanked by the ΦC31 sites. Thus, by 200 minutes almost all ΦC31 sites are inactive and the cell will stay ON. -->
 +
However, if we measure the fluorescence at around 400 mins, the effect of leakiness becomes negligible: in the case with null inputs, the level of the output concentration, GFP, is negligible compared to the level of the GFP concentration produced when the cell receives only one type of AHL. At this point, our model behaves like an XOR gate. Therefore, one of the solutions we propose is to kill or freeze the cells in each row after 6.5 hours.
-
<html><article id="Leakiness"></html>
+
 
 +
<!-- Alternatively, we propose a modified construct, where production of LasI is regulated by a weaker promoter and is induced by a protein, whose production is coupled with production of GFP. -->
 +
 
 +
Initially from our model we observed that the feedback was rapid and hence, the amplification was much higher. However, from literature <sup>[[Team:ETH_Zurich/project/references|[9]]]</sup> we see that the XOR module is relatively slow. We were able to correct this by modelling transcription and translation steps. The delay introduced seems more realistic although we do not have our own experimental data to validate the same. Further, we use a dilution factor (DF) which represents the density of the cells in the bead. By choosing and appropriate DF we can get a more delayed response.
 +
 
 +
$$DF = \frac{No. of cells * V_{E.coli}}{V_{Bead}}$$
 +
 
 +
 
 +
[[File:ETHZ_XORWholeCellwithLeakiness.png|center|800px|thumb|'''Figure 6''' Predicted XOR behaviour for the whole cell model with leakiness]].
 +
 
 +
=== With Leakiness and Crosstalk ===
 +
 
 +
We modelled cross-talk by fitting Hill-functions to experimental data.  For each quorum sensing module, there are two levels of [https://2014.igem.org/Team:ETH_Zurich/expresults cross-talk]. At the first level, we could have any of the two AHLs binding to a given regulator thus activating the promoter while at the second level there is cross talk between regulators and their respective native promoters. Thus, for activation of a given promoter, we can have four Hill functions corresponding to the combinations of inducer and regulator interactions. Since we are using two quorum sensing modules activating two different promoters, we have eight Hill functions.
 +
 
 +
As all the measurements were in terms of fluorescence, the V<sub>max</sub> for each Hill function was also in terms of fluorescence. In order to observe the effect of cross-talk we normalized the V<sub>max</sub> of non native interactions (in Lux system, P<sub>Lux</sub> being activated by LasAHL-LuxR, LasAHL-LasR or LuxAHL-LasR complexes) with the V<sub>max</sub> of the native interaction (LuxR binding to LuxAHL and activating P<sub>Lux</sub>). These ratios acted as weights for the effect of a non-native interaction on the promoter.
 +
 
 +
 
 +
[[File:ETHZ_DynamicsWithCrossTalk.png|center|800px|thumb|'''Figure 7''' The four cases for the whole cell with basal leakiness and cross-talk]]
 +
 
 +
 
 +
For an ideal whole cell model, it was necessary  to have no or minimum leakiness and cross-talk. However, in reality this was not the case. Although we were able to reduce leakiness significantly using riboregulators, we still had the issue of cross-talk. However, we observe that the system acts as an XOR gate till about 300 mins even with cross talk.
 +
In Figure 7 we observe the four cases with both cross-talk and leakiness.
 +
 
 +
 
 +
[[File:ETHZ_WikiXORCTDiffLeakiness.png|center|500px|thumb|'''Figure 8''' The four cases for the whole cell with basal leakiness and cross-talk]]
 +
 
 +
Figure 8, shows the predicted XOR gate performance by a single cell model with basal leakiness and cross talk for different concentrations of LasAHL and LuxAHL at the end of 270 min. The system turns OFF when the input concentration of either AHLs is greater than 100 nM due to the effect of cross-talk.
<html></article></html>
<html></article></html>
 +
<html><article id="Alternate"></html>
-
<html><article id="Cross-talk"></html>
+
== Alternate Design ==
 +
We collaborated with the [https://2014.igem.org/Team:Edinburgh Edinburgh iGEM team 2014]. Their project is focused on communication between bacteria using metabolites as communicating molecules. Combining one metabolite molecule with another quorum sensing molecule could be a solution to avoid cross-talk. We derived a model on the molecular level of their problem (see the [https://2014.igem.org/Team:ETH_Zurich/modeling/qs#Alternate_Design Alternate Design section on the quorum sensing module page]). We inserted the equations in our system, putting B (metabolite) instead of LasAHL as communicating molecule.
 +
<br/>
 +
<br/>
 +
Here, we present the simulation results. The four possibilities of the XOR gate are presented in the following figure. The system was the one with a feedback corresponding to the metabolite pathway. Thus, when the only input is Lux-AHL, the cell will produce the metabolite B and create a feedback loop, which will switch off the gate. This behvior can indeed be observed on the following figure.
 +
<br/>
 +
<br/>
 +
[[File:ETH_Zurich_Collaboration_Four_Cases.jpg|800px|center|thumb|'''Figure 9''' Dynamic simulations of the deterministic single-cell model. One input is a metabolite concentration. Cross-talk is then supposed to be assumed.]]
 +
<br/>
 +
<br/>
 +
It exhibits an XOR behavior. It could be possible to combine quorum sensing and metabolite wiring to make the pattern appear.
<html></article></html>
<html></article></html>
{{:Team:ETH Zurich/tpl/foot}}
{{:Team:ETH Zurich/tpl/foot}}

Latest revision as of 02:39, 18 October 2014

iGEM ETH Zurich 2014