Team:Aberdeen Scotland/modeling/combined model

=Combined Chemotaxis, Lysis and Glue Deposition Model=

Introduction
In order to examine how the system as a whole behaves it is necessary to combine the results of modelling the gene regulatory network with a macroscopic model of interacting E.coli cells. This was done by estimating that quorum sensing becomes activated when the density of bacteria around the source of chemoattractant goes above a critical value. The radius within which quorum sensing is activated varies over time, and depends on how many bacteria there are. We combine the lysis time calculated from our stochastic simulations with this model so that this amount of time after quorum sensing is activated the bacteria lyse and deposit glue.

Results
In the simulation we assume the bacteria to be distributed throughout the pipeline so that they continually arrive at the breach by chemotaxis. Due to a lack of computing power, we performed a scaled simulation by having quorum sensing activate for a lower density of bacteria than in reality. The graph below shows the number of alive bacteria at the site of the breach:



As you can see in the above graph the number of bacteria increases initially as bacteria arrive from the surrounding area. It eventually reaches an equilibrium value as the rate of bacteria lysing equals the rate of arriving bacteria.

As the bacteria number reaches an equilibrium the effective radius within which the bacteria quorum sense also reaches an equilibrium, as illustrated below at point 3.

 

The above graph and diagram show the progression of the quorum sensing region as bacteria continue to arrive. As there is a constant rate of lysing bacteria, there is a linear increase in glue concentration at the breach, as illustrated below:



The increase in glue concentration does not stop, as we have not modelled the eventual sealing of the hole. Once the hole is sealed the chemoattractant gradient would dissipate, along with the signalling molecule IPTG. Eventually the genes for lysis and glue production would no longer be activated and the glue concentration would stop increasing.

Conclusion
This model shows, assuming that the gene regulatory network behaves as expected, that the E.coli would chemotax, activate and deposit glue. Actual experiments would need to be done with lysing cells that actually deposit glue to say for certain whether the Pico-Plumber would seal the hole, and how big the plug would be. As the glue production does not stop until the IPTG is gone, it can be assumed that the hole would actually be sealed, the size of the glue plug is however uncertain.