Circuit containing hydrolases

A seven-species model was implemented to simulate the behaviour of our multicellular sender–receiver system. The model consist of one partial differential equation for AHL and 6 ordinary differential equations with Hill functions that captured the activation of protein synthesis as a function of the concentration of the signalling molecule.

Video: GusA expression levels over 24 hours, concentration in mol/m³. Distance between colonies: 1.5 cm.

For the agar plate and mine cells modules, we use the system of equations and parameters set of the previous simulation.

Mine Cells

The PDEs for the states involved in the sender module are given below:

\begin{align} \frac{d[LuxI]}{dt} =\alpha_{LuxI} - d_{LuxI} \cdot [LuxI]\\ \end{align} \begin{align} \frac{d[NagZ]}{dt} =\alpha_{NagZ} - d_{NagZ} \cdot [NagZ]\\ \end{align} \begin{align} \frac{d[AHL]}{dt}= DF \cdot \left(\alpha_{AHL} \cdot LuxI-d_{AHL,i} \cdot [AHL]\right) + C_{agar} \cdot D_{AHL} \nabla^{2} AHL\\ \end{align}

System of differential equations for mine cells.

Agar Plate

The PDE for OHHL in the agar plate is given below:

Equation 2: Processes taking place on the agar plate: diffusion and decay of OHHL.

Receiver Cells

Receivers are engineered to respond differently to two OHHL concentration levels. Basically, cells should be capable of produce a visible response and in a reasonable amount of time, because the player will need to discriminate between the presence of 0, 1 or 2 adjacent mines. To achieve this goal, we incorporate two enzymatic reporters acting as filters (GusA and AES). The expression of these enzymes is under the control of pLux promoters with different sensitivities, the wild type and a promoter mutant (G1 mutant); The affinity constants for the promoters were obtained from our experimental data. An importan feature of these enzymes is that they can catalyze the hydrolysis of various chromogenic compounds to give rise to a relatively quick coloured response.

The intracellular species of interest in the receiver cells module include: LuxR, OHHL, LuxR/OHHL complex (denoted as R) and the hydrolases (GusA and AES).

Equation System 3 : System of Differential equations for the receiver cells.

In addition to new proteins incorporated to the circuit, it is important to emphasize that the grid was changed to a three neighbours setup.

Results

In figure 1 and 3 are shown the results regarding GusA expression levels. From the time course plot (Fig 1.), it can be pointed out that the activation of the gene is consistent with the OHHL concentration that is predicted to reach a receiver colony. Additionally, one can say that the incubations of the plates should not be longer than 17 hours, otherwise receiver cells that are more that one cell away from a mine can start synthesizing the enzyme.

Figure 1: Averaged GusA concentration for a colony with 0, 1 or 2 adjacent mines cells, respectively. Distance between colonies: 1.5 cm.

Figure 2: Averaged OHHL concentration sensed by a colony with 0, 1 or 2 adjacent mines cells, respectively. The horizontal line indicates the activation coefficient by the dimer LuxR/OHHL (K_R1). Distance between colonies: 1.5 cm.

Figure 3: GusA expression profile after 16 hours. Display: GusA concentration in mM. Distance between colonies: 1.5 cm.

Figure 4: 1D averaged OHHL concentration gradient after 11 hours, radial direction

Regarding the expression of AES under the control of pLuxR variant (Fig 4.), there is not activation before 15 hours only basal expression. Should be noted that the predicted concentration of enzyme is very low if it is compared to GusA.