A three-dimensional, stochastic simulation of biofilm growth and transport-related factors that affect structure

Modeling Microbial Community Networks: Methods and Tools for Studying Microbial Interactions

Microbial interactions function as a fundamental unit in complex ecosystems. By characterizing the type of interaction (positive, negative, neutral) occurring in these dynamic systems, one can begin to unravel the role played by the microbial species. Towards this, various methods have been developed to decipher the function of the microbial communities. The current review focuses on the various qualitative and quantitative methods that currently exist to study microbial interactions. Qualitative methods such as co-culturing experiments are visualized using microscopy-based techniques and are combined with data obtained from multi-omics technologies (metagenomics, metabolomics, metatranscriptomics). Quantitative methods include the construction of networks and network inference, computational models, and development of synthetic microbial consortia. These methods provide a valuable clue on various roles played by interacting partners, as well as possible solutions to overcome pathogenic microbes that can cause life-threatening infections in susceptible hosts. Studying the microbial interactions will further our understanding of complex less-studied ecosystems and enable design of effective frameworks for treatment of infectious diseases.

Graphene-Based Coating to Mitigate Biofilm Development in Marine Environments

Due to its several economic and ecological consequences, biofouling is a widely recognized concern in the marine sector. The search for non-biocide-release antifouling coatings has been on the rise, with carbon-nanocoated surfaces showing promising activity. This work aimed to study the impact of pristine graphene nanoplatelets (GNP) on biofilm development through the representative marine bacteria Cobetia marina and to investigate the antibacterial mechanisms of action of this material. For this purpose, a flow cytometric analysis was performed and a GNP/polydimethylsiloxane (PDMS) surface containing 5 wt% GNP (G5/PDMS) was produced, characterized, and assessed regarding its biofilm mitigation potential over 42 days in controlled hydrodynamic conditions that mimic marine environments. Flow cytometry revealed membrane damage, greater metabolic activity, and endogenous reactive oxygen species (ROS) production by C. marina when exposed to GNP 5% (w/v) for 24 h. In addition, C. marina biofilms formed on G5/PDMS showed consistently lower cell count and thickness (up to 43% reductions) than PDMS. Biofilm architecture analysis indicated that mature biofilms developed on the graphene-based surface had fewer empty spaces (34% reduction) and reduced biovolume (25% reduction) compared to PDMS. Overall, the GNP-based surface inhibited C. marina biofilm development, showing promising potential as a marine antifouling coating.

Assessment of the Antibiofilm Performance of Chitosan-Based Surfaces in Marine Environments

Marine biofouling is a natural process often associated with biofilm formation on submerged surfaces, creating a massive economic and ecological burden. Although several antifouling paints have been used to prevent biofouling, growing ecological concerns emphasize the need to develop new and environmentally friendly antifouling approaches such as bio-based coatings. Chitosan (CS) is a natural polymer that has been widely used due to its outstanding biological properties, including non-toxicity and antimicrobial activity. This work aims to produce and characterize poly (lactic acid) (PLA)-CS surfaces with CS of different molecular weight (Mw) at different concentrations for application in marine paints. Loligo opalescens pens, a waste from the fishery industry, were used as a CS source. The antimicrobial activity of the CS and CS-functionalized surfaces was assessed against Cobetia marina, a model proteobacterium for marine biofouling. Results demonstrate that CS targets the bacterial cell membrane, and PLA-CS surfaces were able to reduce the number of culturable cells up to 68% compared to control, with this activity dependent on CS Mw. The antifouling performance was corroborated by Optical Coherence Tomography since PLA-CS surfaces reduced the biofilm thickness by up to 36%, as well as the percentage and size of biofilm empty spaces. Overall, CS coatings showed to be a promising approach to reducing biofouling in marine environments mimicked in this work, contributing to the valorization of fishing waste and encouraging further research on this topic.

Challenges of biofilm control and utilization: lessons from mathematical modelling

This article reviews modern applications of mathematical descriptions of biofilm formation. The focus is on theoretically obtained results which have implications for areas including the medical sector, food industry and wastewater treatment. Examples are given as to how models have contributed to the overall knowledge on biofilms and how they are used to predict biofilm behaviour. We conclude that the use of mathematical models of biofilms has demonstrated over the years the ability to significantly contribute to the vast field of biofilm research. Among other things, they have been used to test various hypotheses on the nature of interspecies interactions, viability of biofilm treatment methods or forces behind observed biofilm pattern formations. Mathematical models can also play a key role in future biofilm research. Many models nowadays are analysed through computer simulations and continue to improve along with computational capabilities. We predict that models will keep on providing answers to important challenges involving biofilm formation. However, further strengthening of the ties between various disciplines is necessary to fully use the tools of collective knowledge in tackling the biofilm phenomenon.

Cellular automata model for human articular chondrocytes migration, proliferation and cell death: An in vitro validation

Articular cartilage is characterized by low cell density of only one cell type, chondrocytes, and has limited self-healing properties. When articular cartilage is affected by traumatic injuries, a therapeutic strategy such as autologous chondrocyte implantation is usually proposed for its treatment. This approach requires in vitro chondrocyte expansion to yield high cell number for cell transplantation. To improve the efficiency of this procedure, it is necessary to assess cell dynamics such as migration, proliferation and cell death during culture. Computational models such as cellular automata can be used to simulate cell dynamics in order to enhance the result of cell culture procedures. This methodology has been implemented for several cell types; however, an experimental validation is required for each one. For this reason, in this research a cellular automata model, based on random-walk theory, was devised in order to predict articular chondrocyte behavior in monolayer culture during cell expansion. Results demonstrated that the cellular automata model corresponded to cell dynamics and computed-accurate quantitative results. Moreover, it was possible to observe that cell dynamics depend on weighted probabilities derived from experimental data and cell behavior varies according to the cell culture period. Thus, depending on whether cells were just seeded or proliferated exponentially, culture time probabilities differed in percentages in the CA model. Furthermore, in the experimental assessment a decreased chondrocyte proliferation was observed along with increased passage number. This approach is expected to having other uses as in enhancing articular cartilage therapies based on tissue engineering and regenerative medicine.

Copper Corrosion and Biocorrosion Events in Premise Plumbing

Corrosion of copper pipes may release high amounts of copper into the water, exceeding the maximum concentration of copper for drinking water standards. Typically, the events with the highest release of copper into drinking water are related to the presence of biofilms. This article reviews this phenomenon, focusing on copper ingestion and its health impacts, the physicochemical mechanisms and the microbial involvement on copper release, the techniques used to describe and understand this phenomenon, and the hydrodynamic effects. A conceptual model is proposed and the mathematical models are reviewed.

Continuum and discrete approach in modeling biofilm development and structure: a review

The scientific community has recognized that almost 99% of the microbial life on earth is represented by biofilms. Considering the impacts of their sessile lifestyle on both natural and human activities, extensive experimental activity has been carried out to understand how biofilms grow and interact with the environment. Many mathematical models have also been developed to simulate and elucidate the main processes characterizing the biofilm growth. Two main mathematical approaches for biomass representation can be distinguished: continuum and discrete. This review is aimed at exploring the main characteristics of each approach. Continuum models can simulate the biofilm processes in a quantitative and deterministic way. However, they require a multidimensional formulation to take into account the biofilm spatial heterogeneity, which makes the models quite complicated, requiring significant computational effort. Discrete models are more recent and can represent the typical multidimensional structural heterogeneity of biofilm reflecting the experimental expectations, but they generate computational results including elements of randomness and introduce stochastic effects into the solutions.

Influence of Nutrient Availability and Quorum Sensing on the Formation of Metabolically Inactive Microcolonies Within Structurally Heterogeneous Bacterial Biofilms: An Individual-Based 3D Cellular Automata Model

The resistance of bacterial biofilms to antibiotic treatment has been attributed to the emergence of structurally heterogeneous microenvironments containing metabolically inactive cell populations. In this study, we use a three-dimensional individual-based cellular automata model to investigate the influence of nutrient availability and quorum sensing on microbial heterogeneity in growing biofilms. Mature biofilms exhibited at least three structurally distinct strata: a high-volume, homogeneous region sandwiched between two compact sections of high heterogeneity. Cell death occurred preferentially in layers in close proximity to the substratum, resulting in increased heterogeneity in this section of the biofilm; the thickness and heterogeneity of this lowermost layer increased with time, ultimately leading to sloughing. The model predicted the formation of metabolically dormant cellular microniches embedded within faster-growing cell clusters. Biofilms utilizing quorum sensing were more heterogeneous compared to their non-quorum sensing counterparts, and resisted sloughing, featuring a cell-devoid layer of EPS atop the substratum upon which the remainder of the biofilm developed. Overall, our study provides a computational framework to analyze metabolic diversity and heterogeneity of biofilm-associated microorganisms and may pave the way toward gaining further insights into the biophysical mechanisms of antibiotic resistance.

Understanding the mechanisms of solid-water reactions through analysis of surface topography

The topography of a reactive surface contains information about the reactions that form or modify the surface and, therefore, it should be possible to characterize reactivity using topography parameters such as surface area, roughness, or fractal dimension. As a test of this idea, we consider a two-dimensional (2D) lattice model for crystal dissolution and examine a suite of topography parameters to determine which may be useful for predicting rates and mechanisms of dissolution. The model is based on the assumption that the reactivity of a surface site decreases with the number of nearest neighbors. We show that the steady-state surface topography in our model system is a function of, at most, two variables: the ratio of the rate of loss of sites with two neighbors versus three neighbors (d(2)/d(3)) and the ratio of the rate of loss of sites with one neighbor versus three neighbors (d(1)/d(3)). This means that relative rates can be determined from two parameters characterizing the topography of a surface provided that the two parameters are independent of one another. It also means that absolute rates cannot be determined from measurements of surface topography alone. To identify independent sets of topography parameters, we simulated surfaces from a broad range of d(1)/d(3) and d(2)/d(3) and computed a suite of common topography parameters for each surface. Our results indicate that the fractal dimension D and the average spacing between steps, E[s], can serve to uniquely determine d(1)/d(3) and d(2)/d(3) provided that sufficiently strong correlations exist between the steps. Sufficiently strong correlations exist in our model system when D>1.5 (which corresponds to D>2.5 for real 3D reactive surfaces). When steps are uncorrelated, surface topography becomes independent of step retreat rate and D is equal to 1.5. Under these conditions, measures of surface topography are not independent and any single topography parameter contains all of the available mechanistic information about the surface. Our results also indicate that root-mean-square roughness cannot be used to reliably characterize the surface topography of fractal surfaces because it is an inherently noisy parameter for such surfaces with the scale of the noise being independent of length scale.

Cellular automata-based modelling and simulation of biofilm structure on multi-core computers

The article presents a mathematical model of biofilm growth for aerobic biodegradation of a toxic carbonaceous substrate. Modelling of biofilm growth has fundamental significance in numerous processes of biotechnology and mathematical modelling of bioreactors. The process following double-substrate kinetics with substrate inhibition proceeding in a biofilm has not been modelled so far by means of cellular automata. Each process in the model proposed, i.e. diffusion of substrates, uptake of substrates, growth and decay of microorganisms and biofilm detachment, is simulated in a discrete manner. It was shown that for flat biofilm of constant thickness, the results of the presented model agree with those of a continuous model. The primary outcome of the study was to propose a mathematical model of biofilm growth; however a considerable amount of focus was also placed on the development of efficient algorithms for its solution. Two parallel algorithms were created, differing in the way computations are distributed. Computer programs were created using OpenMP Application Programming Interface for C++ programming language. Simulations of biofilm growth were performed on three high-performance computers. Speed-up coefficients of computer programs were compared. Both algorithms enabled a significant reduction of computation time. It is important, inter alia, in modelling and simulation of bioreactor dynamics.

Modeling MIC copper release from drinking water pipes

Copper is used for household drinking water distribution systems given its physical and chemical properties that make it resistant to corrosion. However, there is evidence that, under certain conditions, it can corrode and release unsafe concentrations of copper to the water. Research on drinking water copper pipes has developed conceptual models that include several physical-chemical mechanisms. Nevertheless, there is still a necessity for the development of mathematical models of this phenomenon, which consider the interaction among physical-chemical processes at different spatial scales. We developed a conceptual and a mathematical model that reproduces the main processes in copper release from copper pipes subject to stagnation and flow cycles, and corrosion is associated with biofilm growth on the surface of the pipes. We discuss the influence of the reactive surface and the copper release curves observed. The modeling and experimental observations indicated that after 10h stagnation, the main concentration of copper is located close to the surface of the pipe. This copper is associated with the reactive surface, which acts as a reservoir of labile copper. Thus, for pipes with the presence of biofilm the complexation of copper with the biomass and the hydrodynamics are the main mechanisms for copper release.

The evolution of antibiotic susceptibility and resistance during the formation of Escherichia coli biofilms in the absence of antibiotics

Background

Explanations for bacterial biofilm persistence during antibiotic treatment typically depend on non-genetic mechanisms, and rarely consider the contribution of evolutionary processes.

Results

Using Escherichia coli biofilms, we demonstrate that heritable variation for broad-spectrum antibiotic resistance can arise and accumulate rapidly during biofilm development, even in the absence of antibiotic selection.

Conclusions

Our results demonstrate the rapid de novo evolution of heritable variation in antibiotic sensitivity and resistance during E. coli biofilm development. We suggest that evolutionary processes, whether genetic drift or natural selection, should be considered as a factor to explain the elevated tolerance to antibiotics typically observed in bacterial biofilms. This could be an under-appreciated mechanism that accounts why biofilm populations are, in general, highly resistant to antibiotic treatment.

DIFFUSION IN BIOFILMS RESPIRING ON ELECTRODES

The goal of this study was to measure spatially and temporally resolved effective diffusion coefficients (D(e)) in biofilms respiring on electrodes. Two model electrochemically active biofilms, Geobacter sulfurreducens PCA and Shewanella oneidensis MR-1, were investigated. A novel nuclear magnetic resonance microimaging perfusion probe capable of simultaneous electrochemical and pulsed-field gradient nuclear magnetic resonance (PFG-NMR) techniques was used. PFG-NMR allowed noninvasive, nondestructive, high spatial resolution in situ D(e) measurements in living biofilms respiring on electrodes. The electrodes were polarized so that they would act as the sole terminal electron acceptor for microbial metabolism. We present our results as both two-dimensional D(e) heat maps and surface-averaged relative effective diffusion coefficient (D(rs)) depth profiles. We found that 1) D(rs) decreases with depth in G. sulfurreducens biofilms, following a sigmoid shape; 2) D(rs) at a given location decreases with G. sulfurreducens biofilm age; 3) average D(e) and D(rs) profiles in G. sulfurreducens biofilms are lower than those in S. oneidensis biofilms-the G. sulfurreducens biofilms studied here were on average 10 times denser than the S. oneidensis biofilms; and 4) halting the respiration of a G. sulfurreducens biofilm decreases the D(e) values. Density, reflected by D(e), plays a major role in the extracellular electron transfer strategies of electrochemically active biofilms.

Inhibition of quorum sensing in a computational biofilm simulation

Bacteria communicate through small diffusible molecules in a process known as quorum sensing. Quorum-sensing inhibitors are compounds which interfere with this, providing a potential treatment for infections associated with bacterial biofilms. We present an individual-based computational model for a developing biofilm. Cells are aggregated into particles for computational efficiency, but the quorum-sensing mechanism is modelled as a stochastic process on the level of individual cells. Simulations are used to investigate different treatment regimens. The response to the addition of inhibitor is found to depend significantly on the form of the positive feedback in the quorum-sensing model; in cases where the model exhibits bistability, the time at which treatment is initiated proves to be critical for the effective prevention of quorum sensing and hence potentially of virulence.

Anisotropic nutrient transport in three-dimensional single species bacterial biofilms

The ability for a biofilm to grow and function is critically dependent on the nutrient availability, and this in turn is dependent on the structure of the biofilm. This relationship is therefore an important factor influencing biofilm maturation. Nutrient transport in bacterial biofilms is complex; however, mathematical models that describe the transport of particles within biofilms have made three simplifying assumptions: the effective diffusion coefficient (EDC) is constant, the EDC is that of water, and/or the EDC is isotropic. Using a Monte Carlo simulation, we determined the EDC, both parallel to and perpendicular to the substratum, within 131 real, single species, three-dimensional biofilms that were constructed from confocal laser scanning microscopy images. Our study showed that diffusion within bacterial biofilms was anisotropic and depth dependent. The heterogeneous distribution of bacteria varied between and within species, reducing the rate of diffusion of particles via steric hindrance. In biofilms with low porosity, the EDCs for nutrient transport perpendicular to the substratum were significantly lower than the EDCs for nutrient transport parallel to the substratum. Here, we propose a reaction-diffusion model to describe the nutrient concentration within a bacterial biofilm that accounts for the depth dependence of the EDC.

A novel planar flow cell for studies of biofilm heterogeneity and flow-biofilm interactions

Biofilms are microbial communities growing on surfaces, and are ubiquitous in nature, in bioreactors, and in human infection. Coupling between physical, chemical, and biological processes is known to regulate the development of biofilms; however, current experimental systems do not provide sufficient control of environmental conditions to enable detailed investigations of these complex interactions. We developed a novel planar flow cell that supports biofilm growth under complex two-dimensional fluid flow conditions. This device provides precise control of flow conditions and can be used to create well-defined physical and chemical gradients that significantly affect biofilm heterogeneity. Moreover, the top and bottom of the flow chamber are transparent, so biofilm growth and flow conditions are fully observable using non-invasive confocal microscopy and high-resolution video imaging. To demonstrate the capability of the device, we observed the growth of Pseudomonas aeruginosa biofilms under imposed flow gradients. We found a positive relationship between patterns of fluid velocity and biofilm biomass due to faster microbial growth under conditions of greater local nutrient influx, but this relationship eventually reversed because high hydrodynamic shear leads to the detachment of cells from the surface. These results reveal that flow gradients play a critical role in the development of biofilm communities. By providing new capability for observing biofilm growth, solute and particle transport, and net chemical transformations under user-specified environmental gradients, this new planar flow cell system has broad utility for studies of environmental biotechnology and basic biofilm microbiology, as well as applications in bioreactor design, environmental engineering, biogeochemistry, geomicrobiology, and biomedical research.

The onion model, a simple neutral model for the evolution of diversity in bacterial biofilms

Bacterial biofilms are particularly resistant to a wide variety of antimicrobial compounds. Their persistence in the face of antibiotic therapies causes significant problems in the treatment of infectious diseases. Seldom have evolutionary processes like genetic drift and mutation been invoked to explain how resistance to antibiotics emerges in biofilms, and we lack a simple and tractable model for the genetic and phenotypic diversification that occurs in bacterial biofilms. Here, we introduce the 'onion model', a simple neutral evolutionary model for phenotypic diversification in biofilms. We explore its properties and show that the model produces patterns of diversity that are qualitatively similar to observed patterns of phenotypic diversity in biofilms. We suggest that models like our onion model, which explicitly invoke evolutionary process, are key to understanding biofilm resistance to bactericidal and bacteriostatic agents. Elevated phenotypic variance provides an insurance effect that increases the likelihood that some proportion of the population will be resistant to imposed selective agents and may thus enhance persistence of the biofilm. Accounting for evolutionary change in biofilms will improve our ability to understand and counter diseases that are caused by biofilm persistence.

In situ effective diffusion coefficient profiles in live biofilms using pulsed-field gradient nuclear magnetic resonance

Diffusive mass transfer in biofilms is characterized by the effective diffusion coefficient. It is well documented that the effective diffusion coefficient can vary by location in a biofilm. The current literature is dominated by effective diffusion coefficient measurements for distinct cell clusters and stratified biofilms showing this spatial variation. Regardless of whether distinct cell clusters or surface-averaging methods are used, position-dependent measurements of the effective diffusion coefficient are currently: (1) invasive to the biofilm, (2) performed under unnatural conditions, (3) lethal to cells, and/or (4) spatially restricted to only certain regions of the biofilm. Invasive measurements can lead to inaccurate results and prohibit further (time-dependent) measurements which are important for the mathematical modeling of biofilms. In this study our goals were to: (1) measure the effective diffusion coefficient for water in live biofilms, (2) monitor how the effective diffusion coefficient changes over time under growth conditions, and (3) correlate the effective diffusion coefficient with depth in the biofilm. We measured in situ two-dimensional effective diffusion coefficient maps within Shewanella oneidensis MR-1 biofilms using pulsed-field gradient nuclear magnetic resonance methods, and used them to calculate surface-averaged relative effective diffusion coefficient (D(rs)) profiles. We found that (1) D(rs) decreased from the top of the biofilm to the bottom, (2) D(rs) profiles differed for biofilms of different ages, (3) D(rs) profiles changed over time and generally decreased with time, (4) all the biofilms showed very similar D(rs) profiles near the top of the biofilm, and (5) the D(rs) profile near the bottom of the biofilm was different for each biofilm. Practically, our results demonstrate that advanced biofilm models should use a variable effective diffusivity which changes with time and location in the biofilm.

Contributions of biofilm and suspended sludge to nitrogen transformation and nitrous oxide emission in hybrid sequencing batch system

Hybrid system combines the nature of suspended growth and attached growth has been widely applied to wastewater treatment. In this research, the contributions to N transformation and N2O emission by biofilm and suspended sludge in the hybrid sequencing-batch reactor for a simultaneous nitrification, denitrification and phosphorus removal process were investigated. For the hybrid system, nitrification occurred mostly in the suspended sludge, while the biofilm played the major role in denitrification. The interaction of the biofilm and the suspended sludge in the same reactor resulted in a better overall nitrogen removal performance with simultaneous nitrification and denitrification. However, N2O emission was the main end product of nitrogen removal for the hybrid system; while it was N2 for the biofilm. The relative low N2O emissions from the pure biofilm and the pure suspended sludge corresponded to the relatively high nitrate at the end of the aeration period compared with the hybrid system.

Biofilms and their relevance to veterinary medicine

Bacteria are renowned for their ability to tolerate and adapt to a wide range of adverse environmental conditions. The primary mechanism that facilitates these adaptations is thought to be the capacity to form and maintain biofilms. Within a biofilm, bacteria become attached to a surface where they exist in complex communities which are able to interact with each other through intracellular communication and thus rapidly adapt to changing environments. The organisms within biofilms are notorious for their resistance towards the host immune response and antibacterial agents compared to their free-living planktonic counterparts. Consequently, biofilms are of significant importance to both clinical and veterinary science. However, although bacterial infections are widely reported in animals their association with biofilms is rarely discussed. The aim of this review is to look at the characteristics of biofilm infections in humans and to relate this knowledge to veterinary science in order to assess their relevance in this area.

Motility influences biofilm architecture in Escherichia coli

Eight Escherichia coli strains were studied in minimal medium with a continuous flow system using confocal microscopy. K12 wild-type strains ATCC 25404 and MG1655 formed the best biofilms ( approximately 43 microm thick, 21 to 34% surface coverage). JM109, DH5alpha, and MG1655 motA formed intermediate biofilms ( approximately 13 microm thick, 41 to 58% surface coverage). BW25113, MG1655 qseB, and MG1655 fliA had poor biofilms (surface coverage less than 5%). The best biofilm-formers, ATCC 25404 and MG1655, displayed the highest motility, whereas the worst biofilm former, BW25113, was motility-impaired. The differences in motility were due to differences in expression of the motility loci qseB, flhD, fliA, fliC, and motA (e.g., qseB expression in MG1655 was 139-fold higher than BW25113 and 209-fold higher than JM109). Motility affected the biofilm architecture as those strains which had poor motility (E. coli JM109, E. coli MG1655 motA, and DH5alpha) formed flatter microcolonies compared with MG1655 and ATCC 25404, which had more dramatic vertical structures as a result of their enhanced motility. The presence of flagella was also found to be important as qseB and fliA mutants (which lack flagella) had less biofilm than the isogenic paralyzed motA strain (threefold less thickness and 15-fold less surface coverage).

Spatial periodicity of Escherichia coli K-12 biofilm microstructure initiates during a reversible, polar attachment phase of development and requires the polysaccharide adhesin PGA

Using fast Fourier transform (FFT) analysis, we previously observed that cells within Escherichia coli biofilm are organized in nonrandom or periodic spatial patterns. Here, we developed a gravity displacement assay for examining cell adherence and used it to quantitatively monitor the formation of two distinct forms of cell attachment, temporary and permanent, during early biofilm development. Temporarily attached cells were mainly surface associated by a cell pole; permanent attachments were via the lateral cell surface. While temporary attachment precedes permanent attachment, both forms can coexist in a population. Exposure of attached cells to gravity liberated an unattached population capable of rapidly reassembling a new monolayer, composed of temporarily attached cells, and possessing periodicity. A csrA mutant, which forms biofilm more vigorously than its wild-type parent, exhibited an increased proportion of permanently attached cells and a form of attachment that was not apparent in the parent strain, permanent polar attachment. Nevertheless, it formed periodic attachment patterns. In contrast, biofilm mutants with altered lipopolysaccharide synthesis (waaG) exhibited increased cell-cell interactions, bypassed the polar attachment step, and produced FFT spectra characteristic of aperiodic cell distribution. Mutants lacking the polysaccharide adhesin beta-1,6-N-acetyl-d-glucosamine (DeltapgaC) also exhibited aperiodic cell distribution, but without apparent cell-cell interactions, and were defective in forming permanent attachments. Thus, spatial periodicity of biofilm microstructure is genetically determined and evident during the formation of temporary cell surface attachments.

Cell population dynamics modulate the rates of tissue growth processes

The development and testing of a discrete model describing the dynamic process of tissue growth in three-dimensional scaffolds is presented. The model considers populations of cells that execute persistent random walks on the computational grid, collide, and proliferate until they reach confluence. To isolate the effect of population dynamics on tissue growth, the model assumes that nutrient and growth factor concentrations remain constant in space and time. Simulations start either by distributing the seed cells uniformly and randomly throughout the scaffold, or from an initial condition designed to simulate the migration and cell proliferation phase of wound healing. Simulations with uniform seeding show that cell migration enhances tissue growth by counterbalancing the adverse effects of contact inhibition. This beneficial effect, however, diminishes and disappears completely for large migration speeds. By contrast, simulations with the "wound" seeding mode show a continual enhancement of tissue regeneration rates with increasing cell migration speeds. We conclude that cell locomotory parameters and the spatial distribution of seed cells can have profound effects on the dynamics of the process and, consequently, on the pattern and rates of tissue growth. These results can guide the design of experiments for testing the effectiveness of biomimetic modifications for stimulating tissue growth.

A general description of detachment for multidimensional modelling of biofilms

A general method for describing biomass detachment in multidimensional biofilm modelling is introduced. Biomass losses from processes acting on the entire surface of the biofilm, such as erosion, are modelled using a continuous detachment speed function F(det). Discrete detachment events, i.e. sloughing, are implicitly derived from simulations. The method is flexible to allow F(det) to take several forms, including expressions dependent on any state variables such as the local biofilm density. This methodology for biomass detachment was integrated with multidimensional (2D and 3D) particle-based multispecies biofilm models by using a novel application of the level set method. Application of the method is illustrated by trends in the dynamics of biofilms structure and activity derived from simulations performed on a simple model considering uniform biomass (case study I) and a model discriminating biomass composition in heterotrophic active mass, extracellular polymeric substances (EPS) and inert mass (case study II). Results from case study I demonstrate the effect of applied detachment forces as a fundamental factor influencing steady-state biofilm activity and structure. Trends from experimental observations reported in literature were correctly described. For example, simulation results indicated that biomass sloughing is reduced when erosion forces are increased. Case study II illustrates the application of the detachment methodology to systems with non-uniform biomass composition. Simulations carried out at different bulk concentrations of substrate show changes in biofilm structure (in terms of shape, density and spatial distribution of biomass components) and activity (in terms of oxygen and substrate consumption) as a consequence of either oxygen-limited or substrate-limited growth.

A framework for multidimensional modelling of activity and structure of multispecies biofilms

Concepts from previous biofilm models were integrated to create a framework for the implementation of multidimensional (2D and 3D) multispecies biofilm models. The framework is here described at three levels: (i) mathematical representation of the processes involved in biofilm formation, (ii) numerical implementation into a computer program (freely available from our website http://www.biofilms.bt.tudelft.nl/frameworkMaterial) and (iii) using the program for the creation of biofilm models with multiple bacterial and solute species. An improved version of the individual-based modelling (IbM) that allows structured biomass was used. In this approach biomass composition may be discriminated into any number of particulate species, including extracellular polymeric substances (EPS) for which specific functionality was included. Detachment is also included, described as occurring at the biofilm surface with variable local rates derived from functions of state variables. The application of this modelling framework to a multispecies system with structured biomass is illustrated in a case study where the competition between an organism capable of accumulating polyhydroxybutyrate (PHB, an internal storage compound) and an EPS-producing organism in a two-species biofilm is analysed. Results illustrate that biofilms enriched in PHB-producing organisms may be obtained by supplying substrate intermittently in feast/famine cycles.