Identifiability, Sensitivity, and Genetic Algorithms in Bacterial Biofilm Selection Models

Bacteria often develop distinct phenotypes to adapt to environmental stress. In particular, they can produce biofilms, dense communities of bacteria that live in a complex extracellular matrix. While previous studies have investigated how bacterial biofilms are regulated under laboratory conditions, they have not considered (1) the data requirements necessary to estimate model parameters and (2) how bacteria respond to recurring stressors in their natural habitats. To address (1), we adapted a mechanistic population model to explore the dynamics of biofilm formation in the presence of predator stress, using synthetic data. We used a Maximum Likelihood Estimation framework to measure crucial parameters underpinning the biofilm formation dynamics. We used genetic algorithms to propose an optimal data collection schedule that minimised parameter identifiability confidence interval widths. Our sensitivity analysis revealed that, within the explored regimes, we could simplify the binding dynamics and eliminate biofilm detachment. To address (2), we proposed a structured version of our model to capture the long-term behaviour and evolutionary selection. In our extended model, the subpopulations feature different maximal rates of biofilm formation. We compared the selection under different predator types and amounts and identified key parameters that affected the speed of selection via sensitivity analysis.