Team:ETH Zurich/modeling/diffmodel
From 2014.igem.org
(→Equations) |
(→Equations) |
||
Line 69: | Line 69: | ||
\frac{d[LuxR]}{dt} &= \alpha_{LuxR} -k_{RLux}[LuxAHL_{int}][LuxR] + k_{-RLux}[RLux] - d_{LuxR}[LuxR] \\ | \frac{d[LuxR]}{dt} &= \alpha_{LuxR} -k_{RLux}[LuxAHL_{int}][LuxR] + k_{-RLux}[RLux] - d_{LuxR}[LuxR] \\ | ||
\frac{d[RLux]}{dt} &= k_{RLux}[LuxAHL_{int}][LuxR] - k_{-RLux}[RLux] - 2 k_{DRLux} [RLux]^2 + 2 k_{-DRLux} [DRLux] - d_{RLux} [RLux] \\ | \frac{d[RLux]}{dt} &= k_{RLux}[LuxAHL_{int}][LuxR] - k_{-RLux}[RLux] - 2 k_{DRLux} [RLux]^2 + 2 k_{-DRLux} [DRLux] - d_{RLux} [RLux] \\ | ||
- | \frac{d[mRNA_{GFP}]}{dt} = \frac{k_{mRNA_{GFP}}[RLux]^2}{K_{mLux}^2 + [RLux]^2}- d_{mRNA_{GFP}} [mRNA_{GFP}]\\ | + | \frac{d[mRNA_{GFP}]}{dt} = \frac{k_{mRNA_{GFP}}[RLux]^2}{K_{mLux}^2 + [RLux]^2}- d_{mRNA_{GFP}} [mRNA_{GFP}]\\ |
\frac{d[GFP]}{dt} &= k_{GFP} [mRNA_{GFP}] - d_{GFP}[GFP]\\ | \frac{d[GFP]}{dt} &= k_{GFP} [mRNA_{GFP}] - d_{GFP}[GFP]\\ | ||
\end{align*} | \end{align*} |
Revision as of 09:09, 15 October 2014
Diffusion model
Model
The main experiment for investigating diffusion is a propagation of the whole pattern through the chip via the quorum sensing module only. We have used beads containing cells which are able to sense luxAHL, produce GFP when they sense it, and amplify the signal for the next row. The combination of the quorum sensing module with diffusion enables to check that cells will amplify the signal enough from one row to the next one, and to check what would be the time scale of the pattern formation.
We used a reaction-diffusion model to combine quorum sensing reactions and diffusion. Here are the different compartments considered :
Species
Name | Description |
---|---|
LuxAHLint | LuxAHL inside cells |
LuxAHLext | LuxAHL outside cells |
LuxR | Constitutively expressed regulator protein that can bind LuxAHL and stimulate transcription of GFP. |
RLux | LuxR and LuxAHL complex which can dimerize. |
mRNAGFP | mRNA for Green fluorescent protein which is produced when the cells are ON. |
GFP | Green fluorescent protein which is produced when the cells are ON. |
mRNALuxI | mRNA for LuxI which is produced when the cells are ON. |
LuxI | Enzyme catalysing the production of LuxAHL from SAM and ACP. |
The beads are the compartments containing the cells. The number of cells per bead is called N. The fraction of volume occupied by cells to total bead volume is $$\frac{V_{\text{internal to cells}}}{V_{\text{bead}}}=N \ \ V_{E. coli}$$
Reactions
$$ \begin{align} \emptyset&\rightarrow LuxR \\ LuxAHLext & \leftrightarrow LuxAHLint\\ LuxAHLint+LuxR & \leftrightarrow RLux\\ RLux+RLux &\leftrightarrow DRLux\\ DRLux+P_{luxOFF} & \leftrightarrow P_{luxON}\\ P_{luxON}&\rightarrow P_{luxON}+mRNA_{GFP}\\ mRNA_{GFP}&\rightarrow GFP \end{align}$$
Equations
$$ \begin{align*} \frac{d[LuxAHLext]}{dt} &= N V_{E. coli} Dm (LuxAHL_{int}-LuxAHL_{ext}) -d_{LuxAHLext}[LuxAHL_{ext}]\\ frac{d[LuxAHLint]}{dt} &= Dm (LuxAHL_{ext}-LuxAHL_{int}) + k_{-RLux}[R_{Lux}]-k_{RLux}[LuxAHL_{int}][LuxR]-d_{LuxAHLint}[LuxAHL_{int}]\\ \frac{d[LuxR]}{dt} &= \alpha_{LuxR} -k_{RLux}[LuxAHL_{int}][LuxR] + k_{-RLux}[RLux] - d_{LuxR}[LuxR] \\ \frac{d[RLux]}{dt} &= k_{RLux}[LuxAHL_{int}][LuxR] - k_{-RLux}[RLux] - 2 k_{DRLux} [RLux]^2 + 2 k_{-DRLux} [DRLux] - d_{RLux} [RLux] \\ \frac{d[mRNA_{GFP}]}{dt} = \frac{k_{mRNA_{GFP}}[RLux]^2}{K_{mLux}^2 + [RLux]^2}- d_{mRNA_{GFP}} [mRNA_{GFP}]\\ \frac{d[GFP]}{dt} &= k_{GFP} [mRNA_{GFP}] - d_{GFP}[GFP]\\ \end{align*} $$
Results