Accurate predictive modelling of the growth of microbial communities requires the credible representation of the interactions of biological, chemical and mechanical processes. However, although biological and chemical processes are represented in a number of Individual-based Models (IbMs) the interaction of growth and mechanics is limited. Conversely, there are mechanically sophisticated IbMs with only elementary biology and chemistry. This study focuses on addressing these limitations by developing a flexible IbM that can robustly combine the biological, chemical and physical processes that dictate the emergent properties of a wide range of bacterial communities. This IbM is developed by creating a microbiological adaptation of the open source Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS). This innovation should provide the basis for “bottom up” prediction of the emergent behaviour of entire microbial systems. In the model presented here, bacterial growth, division, decay, mechanical contact among bacterial cells, and adhesion between the bacteria and extracellular polymeric substances are incorporated. In addition, fluid-bacteria interaction is implemented to simulate biofilm deformation and erosion. The model predicts that the surface morphology of biofilms becomes smoother with increased nutrient concentration, which agrees well with previous literature. In addition, the results show that increased shear rate results in smoother and more compact biofilms. The model can also predict shear rate dependent biofilm deformation, erosion, streamer formation and breakup.

This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International licence Newcastle University ePrints -eprint.ncl.ac.uk

The production of extracellular polymeric substance (EPS) is important for the survival of biofilms. However, EPS production is costly for bacteria and the bacterial strains that produce EPS (EPS+) grow in the same environment as non-producers (EPS−) leading to competition between these strains for nutrients and space. The outcome of this competition is likely to be dependent on factors such as initial attachment, EPS production rate, ambient nutrient levels and quorum sensing. We use an Individual-based Model (IbM) to study the competition between EPS+ and EPS− strains by varying the nature of initial colonizers which can either be in the form of single cells or multicellular aggregates. The microbes with EPS+ characteristics obtain a competitive advantage if they initially colonize the surface as smaller aggregates and are widely spread-out between the cells of EPS−, when both are deposited on the substratum. Furthermore, the results show that quorum sensing-regulated EPS production may significantly reduce the fitness of EPS producers when they initially deposit as aggregates. The results provide insights into how the distribution of bacterial aggregates during initial colonization could be a deciding factor in the competition among different strains in biofilms.

No abstract

Individual based Models (IbM) must transition from research tools to engineering tools. To make the transition we must aspire to develop large, three dimensional and physically and biologically credible models. Biological credibility can be promoted by grounding, as far as possible, the biology in thermodynamics. Thermodynamic principles are known to have predictive power in microbial ecology. However, this in turn requires a model that incorporates pH and chemical speciation. Physical credibility implies plausible mechanics and a connection with the wider environment. Here, we propose a step toward that ideal by presenting an individual based model connecting thermodynamics, pH and chemical speciation and environmental conditions to microbial growth for 5·10 5 individuals. We have showcased the model in two scenarios: a two functional group nitrification model and a three functional group anaerobic community. In the former, pH and connection to the environment had an important effect on the outcomes simulated. Whilst in the latter pH was less important but the spatial arrangements and community productivity (that is, methane production) were highly dependent on thermodynamic and reactor coupling. We conclude that if IbM are to attain their potential as tools to evaluate the emergent properties of engineered biological systems it will be necessary to combine the chemical, physical, mechanical and biological along the lines we have proposed. We have still fallen short of our ideals because we cannot (yet) calculate specific uptake rates and must develop the capacity for longer runs in larger models. However, we believe such advances are attainable. Ideally in a common, fast and modular platform. For future innovations in IbM will only be of use if they can be coupled with all the previous advances.