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.

We present NUFEB (Newcastle University Frontiers in Engineering Biology), a flexible, efficient, and open source software for simulating the 3D dynamics of microbial communities. The tool is based on the Individual-based Modelling (IbM) approach, where microbes are represented as discrete units and their behaviour changes over time due to a variety of processes. This approach allows us to study population behaviours that emerge from the interaction between individuals and their environment. NUFEB is built on top of the classical molecular dynamics simulator LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator), which we extended with IbM features. A wide range of biological, physical and chemical processes are implemented to explicitly model microbial systems, with particular emphasis on biofilms. NUFEB is fully parallelised and allows for the simulation of large numbers of microbes (107 individuals and beyond). The parallelisation is based on a domain decomposition scheme that divides the domain into multiple sub-domains which are distributed to different processors. NUFEB also offers a collection of post-processing routines for the visualisation and analysis of simulation output. In this article, we give an overview of NUFEB’s functionalities and implementation details. We provide examples that illustrate the type of microbial systems NUFEB can be used to model and simulate.

No abstract

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.

Biofilms occur in a broad range of environments under heterogeneous physicochemical conditions, such as in bioremediation plants, on surfaces of biomedical implants, and in the lungs of cystic fibrosis patients. In these scenarios, biofilms are subjected to shear forces, but the mechanical integrity of these aggregates often prevents their disruption or dispersal. Biofilms' physical robustness is the result of the multiple biopolymers secreted by constituent microbial cells which are also responsible for numerous biological functions. A better understanding of the role of these biopolymers and their response to dynamic forces is therefore crucial for understanding the interplay between biofilm structure and function. In this paper, we review experimental techniques in rheology, which help quantify the viscoelasticity of biofilms, and modeling approaches from soft matter physics that can assist our understanding of the rheological properties. We describe how these methods could be combined with synthetic biology approaches to control and investigate the effects of secreted polymers on the physical properties of biofilms. We argue that without an integrated approach of the three disciplines, the links between genetics, composition, and interaction of matrix biopolymers and the viscoelastic properties of biofilms will be much harder to uncover.