Team:NTNU Trondheim/Model

Introduction

Bacteria produce insignificant amount of biofilm when they are alone. When bacterias 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 3 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 contains 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 3 dimensional grid, K is a half-saturation constant, Z_QSM 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 V_max and Y_max 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 and Discussion