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 :

Millifluidic compartments used for the diffusion model

Species

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

According to Fick's law, the flow of external AHL diffusing into the cell is $$ \sigma \mathcal{A} ([AHL_{ext}-[AHL_{int}]]) $$ where $\sigma$ is the membrane permeability and \mathcal{A} the area of the membrane. Thus the rate of external AHL diffusing into all cells of a bead is $$Diff(AHL_{int})=\frac{\sigma \mathcal{A}}{V_{int}} ([AHL_{ext}-[AHL_{int}]])$$ and the rate of internal AHL diffusing out of all cells of a bead into the bead is $$Diff(AHL_{ext})=\frac{\sigma \mathcal{A}}{V_{ext}} ([AHL_{ext}-[AHL_{int}]])$$

$$ \begin{align*} \frac{d[LuxAHLext]}{dt} &= N \ V_{E. coli}\ D_m (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