Introduction
Bacteria produce insignificant amount of biofilm when they are alone. When bacteria are grouped up however, they produce a biofilm at a much higher rate per bacteria than they did alone. This is attributed to the quorom sensing molecule (QSM). When the bacteria senses a high concentration of QSM it starts producing biofilm at a much higher rate than when it is low.
We have produced a computational model inspired by Fozard et al. 2012 article, "Inhibition of quorom sensing in a computational biofilm simulation". We use a three dimensional grid of boxes, which contain information about bacteria and concentrations of molecules. By setting up a number of differential equations, one can from an initial state of bacteria and concentrations learn how they evolve over time. In our project we are interested in finding how many bacteria is activated by the quorom sensing molecules and how the biofilm is produced. The code is available on github.
Theory
The boxes used in the simulation contain information of particles and concentration of QSM and substrate. Each particle can either be a collection of bacteria, or an exopolysaccharide (EPS) particle. The concentration was updated according to
where i is the index to the box and q is the type of molecules, substrate and QSM. D is the diffusion constant, l is the side-length of the boxes in the three dimensional grid, K is a half-saturation constant, ZQSM is the production rate of QSM in up- and down-regulated states, M is the mass of the bacteria particle, and n is the number of up- and down regulated bacteria. The mass M of each bacteria is modeled to be proportional to the difference in substrate uptake and maintenance rate,
where Vmax and Ymax is the maximum substrate uptake and maximum yield respectively and m is the maintenance rate. The production of exopolysaccharide (EPS), or biofilm, is depdendent on the number of up- and down-regulated cells,
At each time-step, a pressure value p was calculated for all the boxes,
The number of particles to be moved out of a box at index i to a neighbouring box is dependent on the pressure difference,
where N is the total amount of particles in the box. The relative probability of moving to a specific box was chosen to be linearly dependent on the pressure difference,
At each timestep, each cell had a probability Q± per minute to convert from down-regulated to up-regulated,
Results
Figure 1 shows how the bacteria count increased over time, and how many of them were up-regulated to produce more EPS and QSM. Three simulations were run, one without refilling the system with substrate, one refilling after 2 hours and one after 4 hours. Figure 2 shows how the amount of EPS increased over time for the three refilling scenarios run.
From Figure 2 we see that the amount of EPS produced is only slightly changed when refilling with new substrate. Therefore we chose to simplify the experimental procedure by not refilling the substrate.