Predicting biofilm deformation with a viscoelastic phase-field model: Modeling and experimental studies

Biofilms commonly develop in flowing aqueous environments, where the flow causes the biofilm to deform. Because biofilm deformation affects the flow regime, and because biofilms behave as complex heterogeneous viscoelastic materials, few models are able to predict biofilm deformation. In this study, a phase-field (PF) continuum model coupled with the Oldroyd-B constitutive equation was developed and used to simulate biofilm deformation. The accuracy of the model was evaluated using two types of biofilms: a synthetic biofilm, made from alginate mixed with bacterial cells, and a Pseudomonas aeruginosa biofilm. Shear rheometry was used to experimentally determine the mechanical parameters for each biofilm, used as inputs for the model. Biofilm deformation under fluid flow was monitored experimentally using optical coherence tomography. The comparison between the experimental and modeling geometries, for selected horizontal cross sections, after fluid-driven deformation was good. The relative errors ranged from 3.2 to 21.1% for the synthetic biofilm and from 9.1 to 11.1% for the P. aeruginosa biofilm. This is the first demonstration of the effectiveness of a viscoelastic PF biofilm model. This model provides an important tool for predicting biofilm viscoelastic deformation. It also can benefit the design and control of biofilms in engineering systems.     

Chemical and antimicrobial treatments change the viscoelastic properties of bacterial biofilms

Changes in the viscoelastic material properties of bacterial biofilms resulting from chemical and antimicrobial treatments were measured by rheometry. Colony biofilms of Staphylococcus epidermidis or a mucoid Pseudomonas aeruginosa were subjected to a classical creep test performed using a parallel plate rheometer. Data were fit to the 4-parameter Burger model to quantify the material properties. Biofilms were exposed to the chloride salts of several common mono-, di-, and tri- valent cations, and to urea, industrial biocides, and antibiotics. Many of these treatments resulted in statistically significant alterations in the material properties of the biofilm. Multivalent cations stiffened the P. aeruginosa biofilm, while ciprofloxacin and glutaraldehyde weakened it. Urea, rifampin, and a quaternary ammonium biocide weakened the S. epidermidis biofilm. In general, there was no correspondence between the responses of the two different types of biofilms to a particular treatment. These results underscore the distinction between the killing power of an antimicrobial agent and its ability to alter biofilm mechanical properties and thereby influence biofilm removal. Understanding biofilm rheology and how it is affected by chemical treatment could lead to improvements in biofilm control.     

Chemical and antimicrobial treatments change the viscoelastic properties of bacterial biofilms

Changes in the viscoelastic material properties of bacterial biofilms resulting from chemical and antimicrobial treatments were measured by rheometry. Colony biofilms of Staphylococcus epidermidis or a mucoid Pseudomonas aeruginosa were subjected to a classical creep test performed using a parallel plate rheometer. Data were fit to the 4-parameter Burger model to quantify the material properties. Biofilms were exposed to the chloride salts of several common mono-, di-, and tri- valent cations, and to urea, industrial biocides, and antibiotics. Many of these treatments resulted in statistically significant alterations in the material properties of the biofilm. Multivalent cations stiffened the P. aeruginosa biofilm, while ciprofloxacin and glutaraldehyde weakened it. Urea, rifampin, and a quaternary ammonium biocide weakened the S. epidermidis biofilm. In general, there was no correspondence between the responses of the two different types of biofilms to a particular treatment. These results underscore the distinction between the killing power of an antimicrobial agent and its ability to alter biofilm mechanical properties and thereby influence biofilm removal. Understanding biofilm rheology and how it is affected by chemical treatment could lead to improvements in biofilm control.     

Simulation of growth and detachment in biofilm systems under defined hydrodynamic conditions

Detachment from biofilms was evaluated using a mixed culture biofilm grown on primary wastewater in a tube reactor. The growth of biofilms and the detachment of biomass from biofilms are strongly influenced by hydrodynamic conditions. In a long-term study, three biofilms were cultivated in a biofilm tube reactor. The conducted experiments of biofilm growth and detachment can be divided into three phases: 1) an exponential phase with a rapid increase of the biofilm thickness, 2) a quasi-steady-state with spontaneous fluctuation of the biofilm thickness between 500 and 1,200 microm in the investigated biofilm systems, and 3) a washout experiment with increased shear stress in three to four steps after several weeks of quasi-steady-state. Whereas the biofilm thickness during the homogeneous growth phase can be regarded constant throughout the reactor, it was found to be very heterogeneous during the quasi-steady-state and the washout experiments. Growth and detachment during all three phases was simulated with the same one-dimensional biofilm model. For each of the three phases, a different detachment rate model was used. During the homogeneous growth phase, detachment was modeled proportional to the biofilm growth rate. During the quasi-steady-state phase, detachment was described by random detachment events assuming a base biofilm thickness. Finally, the washout experiment was simulated with detachment being a function of the biofilm thickness before the increase of the shear stress.     

Optical method for long-term and large-scale monitoring of spatial biofilm development

A method was developed that allows biofilm monitoring on the square centimeter scale over extended periods of time. The method is based on image acquisition using a desktop scanner and subsequent image analysis. It was shown that results from grey level analysis are highly correlated with physical properties of the biofilm like average biomass and biofilm thickness. The scanner method was applied to monitor overall biofilm growth, detachment, and surface roughness during two 3 and 4 week long experiments. Two significantly different growth dynamics during the biofilm development could be identified, depending on the biofilm history. Surface roughness on transects in flow direction was always higher than on transects perpendicular to the flow, reflecting the anisotropic characteristics of biofilms growing in a flow field.     

Real time monitoring of biofilm development under flow conditions in porous media

Biofilm growth can impact the effectiveness of industrial processes that involve porous media. To better understand and characterize how biofilms develop and affect hydraulic properties in porous media, both spatial and temporal development of biofilms under flow conditions was investigated in a translucent porous medium by using Pseudomonas fluorescens HK44, a bacterial strain genetically engineered to luminesce in the presence of an induction agent. Real-time visualization of luminescent biofilm growth patterns under constant pressure conditions was captured using a CCD camera. Images obtained over 8 days revealed that variations in bioluminescence intensity could be correlated to biofilm cell density and hydraulic conductivity. These results were used to develop a real-time imaging method to study the dynamic behavior of biofilm evolution in a porous medium, thereby providing a new tool to investigate the impact of biological fouling in porous media under flow conditions.     

Viscoelastic fluid description of bacterial biofilm material properties

A mathematical model describing the constitutive properties of biofilms is required for predicting biofilm deformation, failure, and detachment in response to mechanical forces. Laboratory observations indicate that biofilms are viscoelastic materials. Likewise, current knowledge of biofilm internal structure suggests modeling biofilms as associated polymer viscoelastic systems. Supporting experimental results and a system of viscoelastic fluid equations with a linear Jeffreys viscoelastic stress-strain law are presented here. This system of equations is based on elements of associated polymer physics and is also consistent with presented and previous experimental results. A number of predictions can be made. One particularly interesting result is the prediction of an elastic relaxation time on the order of a few minutes-biofilm disturbances on shorter time scales produce an elastic response, biofilm disturbances on longer time scales result in viscous flow, i.e., nonreversible biofilm deformation. Although not previously recognized, evidence of this phenomenon is in fact present in recent experimental results.     

The sociobiology of biofilms

Biofilms are densely packed communities of microbial cells that grow on surfaces and surround themselves with secreted polymers. Many bacterial species form biofilms, and their study has revealed them to be complex and diverse. The structural and physiological complexity of biofilms has led to the idea that they are coordinated and cooperative groups, analogous to multicellular organisms. We evaluate this idea by addressing the findings of microbiologists from the perspective of sociobiology, including theories of collective behavior (self-organization) and social evolution. This yields two main conclusions. First, the appearance of organization in biofilms can emerge without active coordination. That is, biofilm properties such as phenotypic differentiation, species stratification and channel formation do not necessarily require that cells communicate with one another using specialized signaling molecules. Second, while local cooperation among bacteria may often occur, the evolution of cooperation among all cells is unlikely for most biofilms. Strong conflict can arise among multiple species and strains in a biofilm, and spontaneous mutation can generate conflict even within biofilms initiated by genetically identical cells. Biofilms will typically result from a balance between competition and cooperation, and we argue that understanding this balance is central to building a complete and predictive model of biofilm formation.     

Mapping of Bacterial Biofilm Local Mechanics by Magnetic Microparticle Actuation

Most bacteria live in the form of adherent communities forming three-dimensional material anchored to artificial or biological surfaces, with profound impact on many human activities. Biofilms are recognized as complex systems but their physical properties have been mainly studied from a macroscopic perspective. To determine biofilm local mechanical properties, reveal their potential heterogeneity, and investigate their relation to molecular traits, we have developed a seemingly new microrheology approach based on magnetic particle infiltration in growing biofilms. Using magnetic tweezers, we achieved what was, to our knowledge, the first three-dimensional mapping of the viscoelastic parameters on biofilms formed by the bacterium Escherichia coli. We demonstrate that its mechanical profile may exhibit elastic compliance values spread over three orders of magnitude in a given biofilm. We also prove that heterogeneity strongly depends on external conditions such as growth shear stress. Using strains genetically engineered to produce well-characterized cell surface adhesins, we show that the mechanical profile of biofilm is exquisitely sensitive to the expression of different surface appendages such as F pilus or curli. These results provide a quantitative view of local mechanical properties within intact biofilms and open up an additional avenue for elucidating the emergence and fate of the different microenvironments within these living materials.     

Heterogeneity of biofilms in rotating annular reactors: Occurrence, structure, and consequences

A rotating annular reactor (Roto Torque) was used for qualitative and quantitative studied on biofilm heterogeneity. In contrast to the classic image of biofilms as smooth, homogeneous layers of biomass on a substratum, studies using various pure and mixed cultures consistently revealed more-dimensional structures that resembled dunes and ridges, among others. These heterogeneities were categorized and their underlying causes analyzed. Contrary to expectations, motility of the microorganisms not a decisive factor in determining biofilm homogeneity. Small Variations in substratum geometry homogeneity. Small variations in substratum geometry and flow patterns were clearly reflected in the biofilm pattern. Nonhomogeneous flow and shear patterns in the reactor, together with inadequate mixing resulted in significant, position-dependent differences in surface growth. It was therefore not possible to take representative samples of the attached biomass. Like many other types of reactors, the Roto Torque reactor is valuable for qualitative and morphological biofilm experiments but less suitable for quantitative physiological and kinetics studies using attached microorganisms. (c) 1994 John Wiley & Sons, Inc.     

A model of fluid-biofilm interaction using a Burger material law

A two-dimensional finite element model of the biofilm response to flow was developed. The numerical code sequentially coupled the fluid dynamics of turbulent, incompressible flow with the mechanical response of a single hemispherical biofilm cluster (approximately 100 microm) attached to the flow boundary. A non-linear Burger material law was used to represent the viscoelastic response of a representative microbial biofilm. This constitutive law was incorporated into the numerical model as a Prony series representation of the biofilm's relaxation modulus. Model simulations illuminated interesting details of this fluid-structure interaction. Simulations revealed that softer biofilms (characterized by lower elastic moduli) were highly susceptible to lift forces and consequently were subject to even greater drag forces found higher in the velocity field. A bimodal deformation path due to the two Burger relaxation times was also observed in several simulations. This suggested that interfacial biofilm may be most susceptible to hydrodynamically induced detachment during the initial relaxation time. This result may prove useful in developing removal strategies. Additionally, plots of lift versus drag suggested that the deformation paths taken by viscoelastic biofilms are largely insensitive to specific material coefficients. Softer biofilms merely seem to follow the same path (as a stiffer biofilm) at a faster rate. These relationships may be useful in estimating the hydrodynamic forces acting on an attached biofilm based on changes in scale and cataloged material properties.     

Applying the digital image correlation method to estimate the mechanical properties of bacterial biofilms subjected to a wall shear stress

A digital image correlation (DIC) method was applied to characterize the mechanical behavior of Pseudomonas aeruginosa biofilms in response to wall shear stress using digital video micrographs taken from biofilm flow cells. The appearance of the biofilm in the transmitted light photomicrographs presented a natural texture which was highly conducive to random encoding for DIC. The displacement fields were calculated for two biofilm specimens. The DIC method concurred with previous analysis showing that biofilms exhibit viscoelastic behavior, but had the advantage over simple length measurements of longitudinal strain that it could precisely measure local strains in length (x) and width (y) within biofilm clusters with a 2 μm resolution as a function of time and wall shear stress. It was concluded that DIC was more accurate at measuring elastic moduli than simple length measurements, but that time-lapse 3D images would enable even more accurate estimates to be performed.     

Population dynamics and in situ kinetics of nitrifying bacteria in autotrophic nitrifying biofilms as determined by real-time quantitative PCR

Population dynamics of ammonia-oxidizing bacteria (AOB) and uncultured Nitrospira-like nitrite-oxidizing bacteria (NOB) dominated in autotrophic nitrifying biofilms were determined by using real-time quantitative polymerase chain reaction (RTQ-PCR) and fluorescence in situ hybridization (FISH). Although two quantitative techniques gave the comparable results, the RTQ-PCR assay was easier and faster than the FISH technique for quantification of both nitrifying bacteria in dense microcolony-forming nitrifying biofilms. Using this RTQ-PCR assay, we could successfully determine the maximum specific growth rate (mu = 0.021/h) of uncultured Nitrospira-like NOB in the suspended enrichment culture. The population dynamics of nitrifying bacteria in the biofilm revealed that once they formed the biofilm, the both nitrifying bacteria grew slower than in planktonic cultures. We also calculated the spatial distributions of average specific growth rates of both nitrifying bacteria in the biofilm based on the concentration profiles of NH4+, NO2-, and O2, which were determined by microelectrodes, and the double-Monod model. This simple model estimation could explain the stratified spatial distribution of AOB and Nitrospira-like NOB in the biofilm. The combination of culture-independent molecular techniques and microelectrode measurements is a very powerful approach to analyze the in situ kinetics and ecophysiology of nitrifying bacteria including uncultured Nitrospira-like NOB in complex biofilm communities.     

Hydrodynamic deformation and removal of Staphylococcus epidermidis biofilms treated with urea, chlorhexidine, iron chloride, or DispersinB

The force-deflection and removal characteristics of bacterial biofilm were measured by two different techniques before and after chemical, or enzymatic, treatment. The first technique involved time lapse imaging of a biofilm grown in a capillary flow cell and subjected to a brief shear stress challenge imparted through increased fluid flow. Biofilm removal was determined by calculating the reduction in biofilm area from quantitative analysis of transmission images. The second technique was based on micro-indentation using an atomic force microscope. In both cases, biofilms formed by Staphylococcus epidermidis were exposed to buffer (untreated control), urea, chlorhexidine, iron chloride, or DispersinB. In control experiments, the biofilm exhibited force-deflection responses that were similar before and after the same treatment. The biofilm structure was stable during the post-treatment shear challenge (1% loss). Biofilms treated with chlorhexidine became less deformable after treatment and no increase in biomass removal was seen during the post-treatment shear challenge (2% loss). In contrast, biofilms treated with urea or DispersinB became more deformable and exhibited significant biofilm loss during the post-treatment flow challenge (71% and 40%, respectively). During the treatment soak phase, biofilms exposed to urea swelled. Biofilms exposed to iron chloride showed little difference from the control other than slight contraction during the treatment soak. These observations suggest the following interpretations: (1) chemical or enzymatic treatments, including those that are not frankly antimicrobial, can alter the cohesion of bacterial biofilm; (2) biocidal treatments (e.g., chlorhexidine) do not necessarily weaken the biofilm; and (3) biofilm removal following treatment with agents that make the biofilm more deformable (e.g., urea, DispersinB) depend on interaction between the moving fluid and the biofilm structure. Measurements such as those reported here open the door to development of new technologies for controlling detrimental biofilms by targeting biofilm cohesion rather than killing microorganisms.     

3D finite element model of biofilm detachment using real biofilm structures from CLSM data

In this work, a three‐dimensional model of fluid–structure interactions (FSI) in biofilm systems is developed in order to simulate biofilm detachment as a result of mechanical processes. Therein, fluid flow past the biofilm surface results in a mechanical load on the structure which in turn causes internal stresses in the biofilm matrix. When the strength of the matrix is exceeded parts of the structure are detached. The model is used to investigate the influence of several parameters related to the mechanical strength of the biofilm matrix, Young's modulus, Reynolds number, and biofilm structure on biofilm detachment. Variations in biofilm strength and flow conditions significantly influence the simulation outcome. With respect to structural properties the model is widely independent from a change of Young's modulus. A further result of this work indicates that the change of biofilm structure due to growth or other processes will significantly change the stress distribution in the biofilm and thereby the detachment rate. An increase of the mechanical load by increasing fluid flow results in a flat surface of the remaining biofilm structure. It is concluded that the change of structure during biofilm development is the key determinant in terms of the detachment behavior. Biotechnol. Bioeng. 2009;103: 177–186. © 2008 Wiley Periodicals, Inc.     

Biofilm deformation in response to fluid flow in capillaries

Biofilms are complex mixtures of microorganisms and extracellular matrix that exist on many wetted surfaces. Recently, magnetic resonance microscopy has been used to measure fluid velocities near biofilms that are attached to the walls of capillary channels. These velocity measurements showed unexpectedly high secondary velocities (i.e., high velocity magnitudes perpendicular to the direction of bulk flow and perpendicular to the surface that the biofilm is attached), and the presence of high secondary velocities near a biofilm could increase the delivery of substrates to the biofilm. A mathematical model, based on the immersed boundary method, is used here to examine the physical interaction between a biofilm and a moving fluid in a capillary and to analyze possible factors that may contribute to the elevated secondary velocities observed experimentally. The simulation predicts the formation of a recirculation downstream of a biofilm, and this recirculation deforms and lifts the biofilm upward from the surface to which the biofilm is attached. Changing the mechanical properties (i.e., stiffness) of the biofilm impacts both the lifting of the biofilm and the magnitude of the secondary velocities. The maximum lifting of the biofilm occurs when the biofilm properties are similar to previous experimental measurements, which indicates that the mechanical properties of the biofilm may be tuned for the generation of maximum secondary velocity magnitude and transport of substrates to the biofilm.     

The Role of Exopolymers Produced by Sphingomonas paucimobilis in Biofilm Formation and Composition

Exopolymers have been associated with the initial adhesion of bacteria, which is the primary step for biofilm formation. Moreover, the polymeric matrix of biofilms has a considerable influence on some of the most important physical and physiological properties of biofilms. The role of extracellular polymers in biofilm formation was studied using three mutants of Sphingomonas paucimobilis with increasing capabilities for exopolymer production. The physical, biochemical and physiological properties of three different layers of each biofilm were determined. The layers were detached by submitting the biofilm to increasing shear stress. The results revealed that the presence of exopolymers in the growth medium was essential for biofilm formation. The mutant producing the highest amount of exopolymer formed very thick biofilms, while the biofilms formed by the medium exopolymer producer were on average 8 times thinner. The lowest exopolymer producer did not form biofilm. In both types of biofilms, exopolymer density increased with depth, although this tendency was more significant in thinner biofilms. Cell distribution was also more heterogeneous in thinner biofilms, exhibiting a greater accumulation of cells in the inner layers. The thicker biofilms had very low activity in the inner layer. This was related to a high accumulation of proteins and DNA in this layer due to cell lysis and hydrolytic activity. Activity in the thin biofilm was constant throughout its depth, suggesting that there was no nutrient limitation. The production of exopolymers by each cell was constant throughout the depth of the biofilms, although it was greater in the case of the higher producer.     

Biofilm formation and interactions of bacterial strains found in wastewater treatment systems

Biofilm formation and adherence properties of 13 bacterial strains commonly found in wastewater treatment systems were studied in pure and mixed cultures using a crystal violet microtiter plate assay. Four different culture media were used, wastewater, acetate medium, glucose medium and diluted nutrient broth. The medium composition strongly affected biofilm formation. All strains were able to form pure culture biofilms within 24 h in at least one of the tested culture media and three strains were able to form biofilm in all four culture media, namely Acinetobacter calcoaceticus ATCC 23055, Comamonas denitrificans 123 and Pseudomonas aeruginosa MBL 0199. The adherence properties assessed were initial adherence, cell surface hydrophobicity, and production of amyloid fibers and extracellular polymeric substances. The growth of dual-strain biofilms showed that five organisms formed biofilm with all 13 strains while seven formed no or only weak biofilm when cocultured. In dual-strain cultures, strains with different properties were able to complement each other, giving synergistic effects. Strongest biofilm formation was observed when a mixture of all 13 bacteria were grown together. These results on attachment and biofilm formation can serve as a tool for the design of tailored systems for the degradation of municipal and industrial wastewater.     

A novel in vitro flat‐bed perfusion biofilm model for determining the potential antimicrobial efficacy of topical wound treatments

Aims: To develop an in vitro flat‐bed perfusion biofilm model that could be used to determine the antimicrobial efficacy of topically applied treatments. Methods and Results: Pseudomonas aeruginosa and Staphylococcus aureus biofilms were grown within continuously perfused cellulose matrices. Enumeration of the biofilm density and eluate was performed at various sampling times, enabling determination of the biofilm growth rate. Two antimicrobial wound dressings were applied to the surface of mature biofilms and periodically sampled. To enable real‐time imaging of biofilm growth and potential antimicrobial kinetics, a bioluminescent Ps. aeruginosa biofilm was monitored using low‐light photometry. Target species produced reproducible steady‐state biofilms at a density of c. 10⁷ per biofilm support matrix, after 24‐h perfusion. Test dressings elicited significant antimicrobial effects, producing differing kill kinetic profiles. There was a good correlation between photon and viable count data. Conclusions: The model enables determination of the antimicrobial profile of topically applied treatments against target species biofilms, accurately differentiating bactericidal from bacteriostatic effects. Moreover, these effects could be monitored in real time using bioluminescence. Significance and Impact of the Study: This is the first in vitro biofilm model which can assess the antimicrobial potential of topical therapies in a dynamic growth environment.