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.

Reactions

$$ \begin{align} &\rightarrow LuxR \\ Lux-AHL+LuxR & \leftrightarrow RLux\\ RLux+RLux &\leftrightarrow DRLux\\ DRLux+P_{luxOFF} & \leftrightarrow P_{luxON}\\ P_{luxON}&\rightarrow P_{luxON}+mRNA_{Bxb1}\\ mRNA_{Bxb1}&\rightarrow Bxb1\\ Bxb1 + Bxb1 &\leftrightarrow DBxb1 \\ DBxb1 + SI_{Bxb1} & \leftrightarrow SA_{Bxb1}\\ Lux-AHL &\rightarrow \\ LuxR &\rightarrow \\ RLux &\rightarrow\\ DRLux &\rightarrow\\ mRNA_{Bxb1} &\rightarrow\\ Bxb1 &\rightarrow\\ DBxb1 &\rightarrow\\ \\ &\rightarrow LasR \\ Las-AHL+LasR & \leftrightarrow RLas \\ RLas+RLas & \leftrightarrow DRLas\\ DRLas+P_{LasOFF} & \leftrightarrow P_{LasON}\\ P_{LasON}&\rightarrow P_{LasON}+mRNA_{\phi C31}\\ mRNA_{\phi C31}&\rightarrow \phi C31\\ \phi C31 + \phi C31 &\leftrightarrow D\phi C 31 \\ D\phi C 31 + SI_{\phi C31} & \leftrightarrow SA_{\phi C31}\\ Las-AHL &\rightarrow \\ LasR &\rightarrow \\ RLas &\rightarrow\\ DRLas &\rightarrow\\ mRNA_{\phi C31} &\rightarrow \\ \phi C31 &\rightarrow \\ D\phi C31 &\rightarrow \\ \\ 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 \\ mRNA_{LasI}&\rightarrow\\ GFP &\rightarrow \\ LasI&\rightarrow\\ \\ S_{LasI}+LasI & \rightarrow Las-AHL\\ \end{align}$$

Equations

The following equations were used for the whole cell model.

$$\begin{align*} \frac{d[Lux-AHL]}{dt} &= k_{-RLux}[R_{Lux}]-k_{RLux}[Lux-AHL][LuxR]-d_{Lux-AHL}[Lux-AHL]\\ \frac{d[LuxR]}{dt} &= \alpha_{LuxR} -k_{RLux}[Lux-AHL][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[DRLux]}{dt} &= k_{DRLux} [RLux]^2 - k_{-DRLux} [DRLux] - d_{DRLux} [DRLux] \\ \frac{d[P_{LuxON}]}{dt} &= k_{P_{LuxON}} [P_{LuxOFF}][DRLux] - k_{-P_{LuxON}} [P_{LuxON}]\\ \frac{d[mRNA_{Bxb1}]}{dt} &= L_{P_{Lux}} + k_{mRNA_{Bxb1}} [P_{LuxON}] - d_{mRNA_{Bxb1}} [mRNA_{Bxb1}]\\ \frac{d[Bxb1]}{dt} &= k_{mRNA_{Bxb1}} [mRNA_{Bxb1}] -2 k_{DBxb1}[Bxb1]^2+ 2 k_{-DBxb1}[DBxb1] - d_{Bxb1}[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[SA_{Bxb1}]}{dt}&=k_{SABxb1}[DBxb1][SI_{Bxb1}]-k_{-SABxb1}[SA_{Bxb1}]\\ \\ \frac{d[Las-AHL]}{dt} &= k_{-RLas}[R_{Las}]-k_{RLas}[Las-AHL][LasR]-d_{Las-AHL}[Las-AHL]\\ \frac{d[LasR]}{dt} &= \alpha_{LasR} -k_{RLas}[Las-AHL][LasR] + k_{-RLas}[RLas] - d_{LasR}[LasR] \\ \frac{d[RLas]}{dt} &= k_{RLas}[Las-AHL][LasR] - k_{-RLas}[RLas] - 2 k_{DRLas} [RLas]^2 + 2 k_{-DRLas} [DRLas] - d_{RLas} [RLas] \\ \frac{d[DRLas]}{dt} &= k_{DRLas} [RLas]^2 - k_{-DRLas} [DRLas] - d_{DRLas} [DRLas] \\ \frac{d[P_{LasON}]}{dt} &= k_{P_{LasON}} [P_{LasOFF}][DRLas] - k_{-P_{LasON}} [P_{LasON}]\\ \frac{d[mRNA_{\phi c31}]}{dt} &= L_{P_{Las}} + k_{mRNA_{\phi c31}} [P_{LasON}] - d_{mRNA_{\phi c31}} [mRNA_{\phi c31}]\\ \frac{d[\phi c31]}{dt} &= k_{mRNA_{\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[SA_{\phi c31}]}{dt}&=k_{SA\phi c31}[D\phi c31][SI_{\phi c31}]-k_{-SA\phi c31}[SA_{\phi c31}]\\ \end{align*}$$

The same holds true for the Las system.

Dynamics

Ideal case

With Leakiness

With Leakiness and Crosstalk