Difference between revisions of "Team:ETH Zurich/Modeling/Single-cell Model"
| Line 99: | Line 99: | ||
<table> | <table> | ||
<tr> | <tr> | ||
| − | + | ||
<td> | <td> | ||
\begin{align*} | \begin{align*} | ||
| Line 110: | Line 110: | ||
&\mathop{\xrightarrow{\hspace{4em}}}_{a_\mathrm{GFP},K_{\mathrm{a,LuxRAHL}}}^{\displaystyle\mathop{\downarrow}^{\text{LuxRAHL}}} \text{GFP}\\ | &\mathop{\xrightarrow{\hspace{4em}}}_{a_\mathrm{GFP},K_{\mathrm{a,LuxRAHL}}}^{\displaystyle\mathop{\downarrow}^{\text{LuxRAHL}}} \text{GFP}\\ | ||
\end{align*} | \end{align*} | ||
| − | </td>-- | + | </td> |
| − | + | <!-- | |
<td> | <td> | ||
\begin{align*} | \begin{align*} | ||
| Line 121: | Line 121: | ||
\text{Aiia}+\text{AHL}&\mathop{\xrightarrow{\hspace{4em}}}^{K_{\mathrm{M}},v_{\mathrm{Aiia}}}\text{Aiia}\\ | \text{Aiia}+\text{AHL}&\mathop{\xrightarrow{\hspace{4em}}}^{K_{\mathrm{M}},v_{\mathrm{Aiia}}}\text{Aiia}\\ | ||
\end{align*} | \end{align*} | ||
| − | </td> | + | </td>--> |
</tr> | </tr> | ||
Revision as of 08:22, 15 September 2015
- Project
- Modeling
- Lab
- Human
Practices - Parts
- About Us
Combined Compartment Model
Overview
In this model we plan to simulate whether our system can work as an AND-GATE. Therefore we simulated the system using compartment to model the density of the E. coli as already explained in the AHL module.
First, we will simulate the model with no amplification of the lactate input to match our experimental results. In a second part, we will simulate the full model. As already done in the AHL module, we will compare three conditions:
- No degradation by AiiA, no riboregulator.
- Degradation of AHL by AiiA.
- Riboregulator controlling LuxI expression.
Results
These equations are the integration of both modules in one compartment model.
Assumptions
Here we will assume that .
Equations
Single cell model
Overview
The single cell model is provided here to simulate the combined model.
Chemical species
| Name | Description |
|---|---|
| AHL | Signaling protein, Acyl homoserine lactone (30C6-HSL) |
| LuxR | Regulator protein, that can bind to AHL to form a complex |
| LuxRAHL | Complex of LuxR and AHL, activates transcription of LuxI |
| LuxI | Autoinducer synthase |
| Aiia | AHL-lactonase, N-Acyl Homoserine Lactone Lactonase |
| Lact | Lactate |
| LacI | Lac operon repressor, DNA-binding protein, acts as a protein |
| IPTG | Isopropyl β-D-1-thiogalactopyranoside, prevents LacI from repressing the gene of interest |
| IL | Dimer formed between LacI and IPTG |
Reactions
| \begin{align*} &\mathop{\xrightarrow{\hspace{4em}}}_{a_{LacI},K_{A,appLact}}^{\displaystyle\mathop{\downarrow}^{\text{Lact}}} \text{LacI}\\ \text{IPTG} + \text{LacI} &\mathop{\mathop{\xrightarrow{\hspace{4em}}}^{\xleftarrow{\hspace{4em}}}}_{k_{\mathrm{IL}}}^{k_{\mathrm{-IL}}} \text{IL}\\ &\mathop{\xrightarrow{\hspace{4em}}}_{a_{LuxR},K_{A,appLact}}^{\displaystyle\mathop{\downarrow}^{\text{Lact}}} \text{LuxR}\\ &\mathop{\xrightarrow{\hspace{4em}}}_{a_{LuxR},K_{R,LacI}}^{\displaystyle\mathop{\bot}^{\text{LacI}}} \text{LuxR}\\ \text{AHL} + \text{LuxR} &\mathop{\mathop{\xrightarrow{\hspace{4em}}}^{\xleftarrow{\hspace{4em}}}}_{k_{\mathrm{LuxRAHL}}}^{k_{\mathrm{-LuxRAHL}}} \text{LuxRAHL}\\ &\mathop{\xrightarrow{\hspace{4em}}}_{a_\mathrm{LuxI},K_{\mathrm{a,LuxRAHL}}}^{\displaystyle\mathop{\downarrow}^{\text{LuxRAHL}}} \text{LuxI}\\ &\mathop{\xrightarrow{\hspace{4em}}}_{a_\mathrm{GFP},K_{\mathrm{a,LuxRAHL}}}^{\displaystyle\mathop{\downarrow}^{\text{LuxRAHL}}} \text{GFP}\\ \end{align*} |
Equations
Combining all of the equations from the two different modules, it yields the following system:
\begin{align*} \frac{d[LacI]}{dt}&=\frac{a_\mathrm{LacI} \cdot (\frac{[Lact]}{K_\mathrm{A,appLact}})^{n_1}}{1+(\frac{[Lact]}{K_\mathrm{A,appLact}})^{n_1}}-d_{\mathrm{LacI}}[LacI]\\ \frac{d[LuxR]}{dt}&=\frac{a_\mathrm{LuxR} \cdot (\frac{[Lact]}{K_\mathrm{A,appLact}})^{n_1}}{1+(\frac{[Lact]}{K_\mathrm{A,appLact}})^{n_1}} \cdot \frac{1}{1+(\frac{[LacI]}{K_{\mathrm{R,LacI}}\cdot (\gamma_2+1)})^{n_\mathrm{2}}}-d_{\mathrm{LuxR}}[LuxR]\\ [LuxRAHL]&= \frac{[AHL]\cdot [LuxR]}{K_{\mathrm{d,LuxRAHL}}+[AHL]}\\ \frac{d[LuxI]}{dt}&=a_{\mathrm{LuxI}}k_{\mathrm{leaky}}([LuxR]-[LuxRAHL])+\frac{a_{\mathrm{LuxI}}(\frac{[LuxRAHL]}{K_{\mathrm{A,LuxRAHL}}})^2}{1+(\frac{[LuxRAHL]}{K_{\mathrm{A,LuxRAHL}}})^2}-d_{\mathrm{LuxI}}[LuxI]\\ \frac{d[AHL]}{dt}&=a_{\mathrm{AHL}}[LuxI]-d_{\mathrm{AHL}}[AHL]-\frac{v_\mathrm{Aiia}\cdot [AHL]}{K_{\mathrm{M,AiiA}}+[AHL]}\\ \frac{d[GFP]}{dt}&=a_\mathrm{GFP}k_{\mathrm{leaky}}([LuxR]-[LuxRAHL])+\frac{a_\mathrm{GFP}(\frac{[LuxRAHL]}{K_{\mathrm{A,LuxRAHL}}})^2}{1+(\frac{[LuxRAHL]}{K_{\mathrm{A,LuxRAHL}}})^2}-d_{\mathrm{GFP}}[GFP]\\ K_\mathrm{d,LuxRAHL} &= \frac{k_\mathrm{-LuxRAHL}}{k_\mathrm{LuxRAHL}}\\ \gamma_2 &= \frac{IPTG_{tot}}{K_{IL}} \end{align*}
We would like to thank our sponsors
