# American Institute of Mathematical Sciences

June  2017, 14(3): 625-653. doi: 10.3934/mbe.2017036

## Mathematical analysis of a quorum sensing induced biofilm dispersal model and numerical simulation of hollowing effects

 1 Biomedical Physics, Dept. Physics, Ryerson University, 350 Victoria Street Toronto, ON, M5B 2K3, Canada 2 Institute for Mathematics and Scientific Computing, University of Graz, Heinrichstr. 36,8010 Graz, Austria 3 Dept. Mathematics and Statistics, University of Guelph, 50 Stone Road East, ON, N1G 2W1, Canada

* Corresponding author: Blessing O. Emerenini

Received  August 20, 2015 Accepted  October 26, 2016 Published  December 2016

We analyze a mathematical model of quorum sensing induced biofilm dispersal. It is formulated as a system of non-linear, density-dependent, diffusion-reaction equations. The governing equation for the sessile biomass comprises two non-linear diffusion effects, a degeneracy as in the porous medium equation and fast diffusion. This equation is coupled with three semi-linear diffusion-reaction equations for the concentrations of growth limiting nutrients, autoinducers, and dispersed cells. We prove the existence and uniqueness of bounded non-negative solutions of this system and study the behavior of the model in numerical simulations, where we focus on hollowing effects in established biofilms.

Citation: Blessing O. Emerenini, Stefanie Sonner, Hermann J. Eberl. Mathematical analysis of a quorum sensing induced biofilm dispersal model and numerical simulation of hollowing effects. Mathematical Biosciences & Engineering, 2017, 14 (3) : 625-653. doi: 10.3934/mbe.2017036
##### References:

show all references

##### References:
Schematic of the biofilm system cf [18]: The aqueous phase is the domain $\Omega_1(t) = \{x\in \Omega: M(t;x) =0 \}$, the biofilm phase $\Omega_2(t) = \{x\in \Omega: M(t;x) >0 \}$. These regions change over time as the biofilm grows. Biofilm colonies form attached to the substratum, which is a part of the boundary of the domain
2-D structural representation of the microbial floc growth for autoinducer production rate $\alpha=30.7$ and maximum dispersal rate $\eta_1=0.6$ for selected time instances t. Color coded is the biomass density $M$, iso-lines of the autoinducer concentration $A$ are plotted in grayscale
1-D Spatial representation of the development and dispersal of bacterial cells from the microbial floc: The snapshots are taken at different computational time $t$, with an autoinducer production rate $\alpha=30.7$ and a dispersal rate of $\eta_1=0.6$
Temporal plots of simulations computed for a non-quorum sensing producing microfloc (Non-QS) and a quorum sensing producing microfloc using seven different constitutive autoinducer production rate $\alpha = \{92.0,46.0,30.7,23.0,18.4,15.3,13.1\}$ and fixed maximum dispersal rate $\eta_1=0.6$. Shown are (a) the total sessile biomass fraction $M_{tot}$ in the floc, (b) the floc size $\omega$ (c) dispersed cells $N_{tot}$, (d) relative variation $R$, and (e) signal concentration $A_{ave}$
2-D structural representation of the microbial biofilm growth for autoinducer production rate $\alpha=30.7$ and maximum dispersal rate $\eta_1=0.6$ for selected time instances $t$. Color coded is the biomass density $M$, iso-lines of the autoinducer concentration $A$ are plotted in grayscale
2-D structural representation of the microbial biofilm growth for autoinducer production rate $\alpha=92.0$ and maximum dispersal rate $\eta_1=0.6$ for selected time instances $t$. Color coded is the biomass density $M$, iso-lines of the autoinducer concentration $A$ are plotted in grayscale
Temporal plots of simulations computed for a non-quorum sensing producing biofilm (Non-QS) and a quorum sensing producing biofilm using seven different constitutive autoinducer production rate $\alpha = \{92.0,46.0,30.7,23.0,18.4,15.3,13.1\}$ and fixed maximum dispersal rate $\eta_1=0.6$. Shown are (a) the total sessile biomass fraction $M_{tot}$ in the floc, (b) the floc size $\omega$ (c) dispersed cells $N_{tot}$, (d) relative variation $R$, and (e) signal concentration $A_{ave}$
Comparison of the sessile biomass $M_{tot}$ and the dispersed cells $N_{tot}$ under different boundary conditions for the signal molecule $A$: Homogenous Dirichlet conditions and Neumann conditions. The left panel is for a microbial floc while the right panel is for a biofilm
Parameters used in the numerical simulations
 Parameter Description Value Source $k_1$ half saturation concentration (growth) $0.4$ [44] $k_2$ lysis rate $0.067$ assumed $\sigma$ nutrient consumption rate $793.65$ [19] $\eta_1$ maximum dispersal rate varied [18] $\lambda$ quorum sensing abiotic decay rate $0.02218$ [39] $\alpha$ constitutive autoinducer production rate varied - $\beta$ induced autoinducer production rate $10 \times \alpha$ [19] $m$ degree of polymerization $2.5$ [19] $d_1$ constant diffusion coefficients for $N$ $4.1667$ assumed $d_2$ constant diffusion coefficients for $C$ $4.1667$ [15] $d_3$ constant diffusion coefficients for $A$ $3.234$ [15] $d$ biomass motility coefficient $4.2 \times 10^{-8}$ [13] $a$ biofilm diffusion exponent $4.0$ [13] $b$ biofilm diffusion exponent $4.0$ [13] $L$ system length $1.0$ [15] $H$ system height $1.0$ assumed
 Parameter Description Value Source $k_1$ half saturation concentration (growth) $0.4$ [44] $k_2$ lysis rate $0.067$ assumed $\sigma$ nutrient consumption rate $793.65$ [19] $\eta_1$ maximum dispersal rate varied [18] $\lambda$ quorum sensing abiotic decay rate $0.02218$ [39] $\alpha$ constitutive autoinducer production rate varied - $\beta$ induced autoinducer production rate $10 \times \alpha$ [19] $m$ degree of polymerization $2.5$ [19] $d_1$ constant diffusion coefficients for $N$ $4.1667$ assumed $d_2$ constant diffusion coefficients for $C$ $4.1667$ [15] $d_3$ constant diffusion coefficients for $A$ $3.234$ [15] $d$ biomass motility coefficient $4.2 \times 10^{-8}$ [13] $a$ biofilm diffusion exponent $4.0$ [13] $b$ biofilm diffusion exponent $4.0$ [13] $L$ system length $1.0$ [15] $H$ system height $1.0$ assumed
 [1] Omid Nikan, Seyedeh Mahboubeh Molavi-Arabshai, Hossein Jafari. Numerical simulation of the nonlinear fractional regularized long-wave model arising in ion acoustic plasma waves. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020466 [2] H. M. Srivastava, H. I. Abdel-Gawad, Khaled Mohammed Saad. Oscillatory states and patterns formation in a two-cell cubic autocatalytic reaction-diffusion model subjected to the Dirichlet conditions. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020433 [3] Christian Beck, Lukas Gonon, Martin Hutzenthaler, Arnulf Jentzen. On existence and uniqueness properties for solutions of stochastic fixed point equations. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020320 [4] Yichen Zhang, Meiqiang Feng. A coupled $p$-Laplacian elliptic system: Existence, uniqueness and asymptotic behavior. Electronic Research Archive, 2020, 28 (4) : 1419-1438. doi: 10.3934/era.2020075 [5] Leilei Wei, Yinnian He. A fully discrete local discontinuous Galerkin method with the generalized numerical flux to solve the tempered fractional reaction-diffusion equation. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020319 [6] Shao-Xia Qiao, Li-Jun Du. Propagation dynamics of nonlocal dispersal equations with inhomogeneous bistable nonlinearity. Electronic Research Archive, , () : -. doi: 10.3934/era.2020116 [7] Parikshit Upadhyaya, Elias Jarlebring, Emanuel H. Rubensson. A density matrix approach to the convergence of the self-consistent field iteration. Numerical Algebra, Control & Optimization, 2021, 11 (1) : 99-115. doi: 10.3934/naco.2020018 [8] Stefano Bianchini, Paolo Bonicatto. Forward untangling and applications to the uniqueness problem for the continuity equation. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020384 [9] Jiaquan Liu, Xiangqing Liu, Zhi-Qiang Wang. Sign-changing solutions for a parameter-dependent quasilinear equation. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020454 [10] Bahaaeldin Abdalla, Thabet Abdeljawad. Oscillation criteria for kernel function dependent fractional dynamic equations. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020443 [11] Wenbin Li, Jianliang Qian. Simultaneously recovering both domain and varying density in inverse gravimetry by efficient level-set methods. Inverse Problems & Imaging, , () : -. doi: 10.3934/ipi.2020073 [12] Mark F. Demers. Uniqueness and exponential mixing for the measure of maximal entropy for piecewise hyperbolic maps. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 217-256. doi: 10.3934/dcds.2020217 [13] Soniya Singh, Sumit Arora, Manil T. Mohan, Jaydev Dabas. Approximate controllability of second order impulsive systems with state-dependent delay in Banach spaces. Evolution Equations & Control Theory, 2020  doi: 10.3934/eect.2020103 [14] Zexuan Liu, Zhiyuan Sun, Jerry Zhijian Yang. A numerical study of superconvergence of the discontinuous Galerkin method by patch reconstruction. Electronic Research Archive, 2020, 28 (4) : 1487-1501. doi: 10.3934/era.2020078 [15] Vieri Benci, Marco Cococcioni. The algorithmic numbers in non-archimedean numerical computing environments. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020449 [16] Weiwei Liu, Jinliang Wang, Yuming Chen. Threshold dynamics of a delayed nonlocal reaction-diffusion cholera model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020316 [17] Abdelghafour Atlas, Mostafa Bendahmane, Fahd Karami, Driss Meskine, Omar Oubbih. A nonlinear fractional reaction-diffusion system applied to image denoising and decomposition. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020321 [18] Pierre-Etienne Druet. A theory of generalised solutions for ideal gas mixtures with Maxwell-Stefan diffusion. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020458 [19] Yongxiu Shi, Haitao Wan. Refined asymptotic behavior and uniqueness of large solutions to a quasilinear elliptic equation in a borderline case. Electronic Research Archive, , () : -. doi: 10.3934/era.2020119 [20] Marion Darbas, Jérémy Heleine, Stephanie Lohrengel. Numerical resolution by the quasi-reversibility method of a data completion problem for Maxwell's equations. Inverse Problems & Imaging, 2020, 14 (6) : 1107-1133. doi: 10.3934/ipi.2020056

2018 Impact Factor: 1.313