Mobilization of broad host range plasmid from Pseudomonas putida to established biofilm of Bacillus azotoformans. II. Modeling

A strain of Pseudomonas putida that harbors plasmids RK2 and pDLB101 was exposed to a pure culture biofilm of Bacillus azotoformans grown in a rotating annular reactor. Transfer of the RK2 mobilizable pDLB101 plasmid to B. azotoformans was monitored over a 4-day period. Experimental results demonstrated that the broad host range, RSF1010 derivative pDLB101 was transferred to and expressed by B. azotoformans. In the companion article to this work, the rate of plasmid transfer was quantified as a function of the limiting nutrient, succinate, and as a function of the mechanism of transfer. A biofilm process simulation program (AQUASIM) was modified to analyze resultant experimental data. Although the AQUASIM package was not designed to simulate or predict genetic events in biofilms, modification of the rate process dynamics allowed successful modeling of plasmid transfer. For the narrow range of substrate concentrations used in these experiments, nutrient level had only a slight effect on the rate and extent of plasmid transfer in biofilms. However, further simulations using AQUASIM revealed that under nutrient poor conditions, the number of transconjugants appearing in the biofilm was limited.     

Redox-stratification controlled biofilm (ReSCoBi) for completely autotrophic nitrogen removal: the effect of co- versus counter-diffusion on reactor performance

A multi-population biofilm model for completely autotrophic nitrogen removal was developed and implemented in the simulation program AQUASIM to corroborate the concept of a redox-stratification controlled biofilm (ReSCoBi). The model considers both counter- and co-diffusion biofilm geometries. In the counter-diffusion biofilm, oxygen is supplied through a gas-permeable membrane that supports the biofilm while ammonia (NH(4)(+)) is supplied from the bulk liquid. On the contrary, in the co-diffusion biofilm, both oxygen and NH(4)(+) are supplied from the bulk liquid. Results of the model revealed a clear stratification of microbial activities in both of the biofilms, the resulting chemical profiles, and the obvious effect of the relative surface loadings of oxygen and NH(4)(+) (J(O(2))/J(NH(4)(+))) on the reactor performances. Steady-state biofilm thickness had a significant but different effect on T-N removal for co- and counter-diffusion biofilms: the removal efficiency in the counter-diffusion biofilm geometry was superior to that in the co-diffusion counterpart, within the range of 450-1,400 microm; however, the efficiency deteriorated with a further increase in biofilm thickness, probably because of diffusion limitation of NH(4)(+). Under conditions of oxygen excess (J(O(2))/J(NH(4)(+)) > 3.98), almost all NH(4)(+) was consumed by aerobic ammonia oxidation in the co-diffusion biofilm, leading to poor performance, while in the counter-diffusion biofilm, T-N removal efficiency was maintained because of the physical location of anaerobic ammonium oxidizers near the bulk liquid. These results clearly reveal that counter-diffusion biofilms have a wider application range for autotrophic T-N removal than co-diffusion biofilms.     

Modeling biofilm accumulation and mass transport in a porous medium under high substrate loading

A packed bed biofilm reactor inoculated with pure culture Pseudomonas aeruginosa was run under high substrate loading and constant flow rate conditions. The 3.1-cm-diameter cylindrical reactor was 5 cm in length and packed with 1-mm glass beads. Daily observations of biofilm thickness, influent and effluent glucose substrate concentration, and effluent dissolved and total organic carbon were made during the 13-day experiment. Biofilm thickness appeared to rech quasi-steady-state condition after 10 days. A published biofilm process simulation program (AQUASIM) was used to analyze experimental data. Comparison of observed and simulated variables revealed three distinct phases of biofilm accumulation during the experiment: an initial phase, a growth phase, and a mature biofilm phase. Different combinations of biofilm and mass transport process variables were found to be important during each phase. Biofilm detachment was highly correlated with shear at the biofilm surface during all three phases of biofilm development. (c) 1995 John Wiley & Sons, Inc.     

Storage and growth of denitrifiers in aerobic granules: part I. model development

A mathematical model, based on the Activated Sludge Model No.3 (ASM3), is developed to describe the storage and growth activities of denitrifiers in aerobic granules under anoxic conditions. In this model, mass transfer, hydrolysis, simultaneous anoxic storage and growth, anoxic maintenance, and endogenous decay are all taken into account. The model established is implemented in the well-established AQUASIM simulation software. A combination of completely mixed reactor and biofilm reactor compartments provided by AQUASIM is used to simulate the mass transport and conversion processes occurring in both bulk liquid and granules. The modeling results explicitly show that the external substrate is immediately utilized for storage and growth at feast phase. More external substrates are diverted to storage process than the primary biomass production process. The model simulation indicates that the nitrate utilization rate (NUR) of granules-based denitrification process includes four linear phases of nitrate reduction. Furthermore, the methodology for determining the most important parameter in this model, that is, anoxic reduction factor, is established.     

A Multidimensional Multispecies Continuum Model for Heterogeneous Biofilm Development

We propose a multidimensional continuum model for heterogeneous growth of biofilm systems with multiple species and multiple substrates. The new model provides a deterministic framework for the study of the interactions between several spe1cies and their effects on biofilm heterogeneity. It consists of a system of partial differential equations derived on the basis of conservation laws and reaction kinetics. The derivation and key assumptions are presented. The assumptions used are a combination of those used in the established one dimensional model, due to Wanner and Gujer, and for the viscous fluid model, of Dockery and Klapper. The work of Wanner and Gujer in particular has been extensively used through the years, and thus this new model is an extension to several spatial dimensions of an already proven working model. The model equations are solved using numerical techniques, for purposes of simulation and verification. The new model is applied to two different biofilm systems in several spatial dimensions, one of which is equivalent to a system originally studied by Wanner and Gujer. Dimensionless formulations for these two systems are given, and numerical simulation results with varying initial conditions are presented. An erratum to this article can be found at     

Linearized kinetic model of<Emphasis Type="Italic"> Listeria monocytogenes</Emphasis> biofilm growth

In this study the dynamics of biofilm formation on aluminum has been investigated. The process of cell growth has been observed using fluorescence microscopy. It has been confirmed that the process of biofilm formation can be represented as a sum of two separate processes: cell adhesion and colony proliferation. The derived set of equations describes kinetics of surface population growth and characteristic times for adsorption and combined growth processes, including characteristic time for the nutrient supply depletion. All equations contain variables based on the fundamental characteristics of bacterial population and can be easily determined from the experimental data or estimated theoretically. The developed theoretical model allows obtaining realistic values for population growth time and characteristic time for nutrient limitation occurrence during the biofilm development. Resulting equations qualitatively describe the biofilm formation process, and allow predicting microbial kinetics in the batch reactor system and determining critical values of the process parameters.     

Program Code Generator for Cardiac Electrophysiology Simulation with Automatic PDE Boundary Condition Handling

Clinical and experimental studies involving human hearts can have certain limitations. Methods such as computer simulations can be an important alternative or supplemental tool. Physiological simulation at the tissue or organ level typically involves the handling of partial differential equations (PDEs). Boundary conditions and distributed parameters, such as those used in pharmacokinetics simulation, add to the complexity of the PDE solution. These factors can tailor PDE solutions and their corresponding program code to specific problems. Boundary condition and parameter changes in the customized code are usually prone to errors and time-consuming. We propose a general approach for handling PDEs and boundary conditions in computational models using a replacement scheme for discretization. This study is an extension of a program generator that we introduced in a previous publication. The program generator can generate code for multi-cell simulations of cardiac electrophysiology. Improvements to the system allow it to handle simultaneous equations in the biological function model as well as implicit PDE numerical schemes. The replacement scheme involves substituting all partial differential terms with numerical solution equations. Once the model and boundary equations are discretized with the numerical solution scheme, instances of the equations are generated to undergo dependency analysis. The result of the dependency analysis is then used to generate the program code. The resulting program code are in Java or C programming language. To validate the automatic handling of boundary conditions in the program code generator, we generated simulation code using the FHN, Luo-Rudy 1, and Hund-Rudy cell models and run cell-to-cell coupling and action potential propagation simulations. One of the simulations is based on a published experiment and simulation results are compared with the experimental data. We conclude that the proposed program code generator can be used to generate code for physiological simulations and provides a tool for studying cardiac electrophysiology.     

Determination of pollutant diffusion coefficients in naturally formed biofilms using a single tube extractive membrane bioreactor

A novel technique has been used to determine the effective diffusion coefficients for 1,1,2-trichloroethane (TCE), a nonreacting tracer, in biofilms growing on the external surface of a silicone rubber membrane tube during degradation of 1,2-dichloroethane (DCE) by Xanthobacter autotrophicus GJ10 and monochlorobenzene (MCB) by Pseudomonas JS150. Experiments were carried out in a single tube extractive membrane bioreactor (STEMB), whose configuration makes it possible to measure the transmembrane flux of substrates. A video imaging technique (VIT) was employed for in situ biofilm thickness measurement and recording. Diffusion coefficients of TCE in the biofilms and TCE mass transfer coefficients in the liquid films adjacent to the biofilms were determined simultaneously using a resistances-in-series diffusion model. It was found that the flux and overall mass transfer coefficient of TCE decrease with increasing biofilm thickness, showing the importance of biofilm diffusion on the mass transfer process. Similar fluxes were observed for the nonreacting tracer (TCE) and the reactive substrates (MCB or DCE), suggesting that membrane-attached biofilm systems can be rate controlled primarily by substrate diffusion. The TCE diffusion coefficient in the JS150 biofilm appeared to be dependent on biofilm thickness, decreasing markedly for biofilm thicknesses of >1 mm. The values of the TCE diffusion coefficients in the JS150 biofilms <1-mm thick are approximately twice those in water and fall to around 30% of the water value for biofilms >1-mm thick. The TCE diffusion coefficients in the GJ10 biofilms were apparently constant at about the water value. The change in the diffusion coefficient for the JS150 biofilms is attributed to the influence of eddy diffusion and convective flow on transport in the thinner (<1-mm thick) biofilms.     

Towards optimum permeability reduction in porous media using biofilm growth simulations

While biological clogging of porous systems can be problematic in numerous processes (e.g., microbial enhanced oil recovery—MEOR), it is targeted during bio‐barrier formation to control sub‐surface pollution plumes in ground water. In this simulation study, constant pressure drop (CPD) and constant volumetric flow rate (CVF) operational modes for nutrient provision for biofilm growth in a porous system are considered with respect to optimum (minimum energy requirement for nutrient provision) permeability reduction for bio‐barrier applications. Biofilm growth is simulated using a Lattice‐Boltzmann (LB) simulation platform complemented with an individual‐based biofilm model (IbM). A biomass detachment technique has been included using a fast marching level set (FMLS) method that models the propagation of the biofilm–liquid interface with a speed proportional to the adjacent velocity shear field. The porous medium permeability reduction is simulated for both operational modes using a range of biofilm strengths. For stronger biofilms, less biomass deposition and energy input are required to reduce the system permeability during CPD operation, whereas CVF is more efficient at reducing the permeability of systems containing weaker biofilms. Biotechnol. Bioeng. 2009;103: 767–779. © 2009 Wiley Periodicals, Inc.     

Model parameter uncertainties in a dual-species biofilm competition model affect ecological output parameters much stronger than morphological ones

Bacterial biofilms are complex microbial depositions on immersed interfaces that form wherever the environmental conditions sustain microbial growth. Despite their name, biofilms can develop in highly irregular structures. Recently several mathematical concepts have been introduced to model these spatially structured microbial populations. Regardless of the type of model, they all have, even for microbially relatively simple systems, many parameters which generally are known at most approximately. We investigate the effect of uncertainties in model parameters on four morphological and four ecological output parameters using a nonlinear diffusion model for a biofilm in which two species compete for a shared nutrient. To this end we conduct an extensive computer simulation experiment for two different levels of data uncertainty, three different hydrodynamic conditions, and two different scenarios of bulk substrate availability. Our results indicate that input model parameter uncertainties have a much larger effect on ecological than on morphological output parameters.     

Development of an oral mucosa model to study host-microbiome interactions during wound healing

Crosstalk between the human host and its microbiota is reported to influence various diseases such as mucositis. Fundamental research in this area is however complicated by the time frame restrictions during which host-microbe interactions can be studied in vitro. The model proposed in this paper, consisting of an oral epithelium and biofilm, can be used to study microbe-host crosstalk in vitro in non-infectious conditions up to 72 h. Microbiota derived from oral swabs were cultured on an agar/mucin layer and challenged with monolayers of keratinocytes grown on plastic or collagen type I layers embedded with fibroblasts. The overall microbial biofilm composition in terms of diversity remained representative for the oral microbiome, whilst the epithelial cell morphology and viability were unaffected. Applying the model to investigate wound healing revealed a reduced healing of 30 % in the presence of microbiota, which was not caused by a reduction of the proliferation index (52.1–61.5) or a significantly increased number of apoptotic (1–1.13) or necrotic (32–30.5 %) cells. Since the model allows the separate study of the microbial and cellular exometabolome, the biofilm and epithelial characteristics after co-culturing, it is applicable for investigations within fundamental research and for the discovery and development of agents that promote wound healing.     

Modeling biocide action against biofilms

A phenomenological model of biocide action against microbial biofilms was derived. Processes incorporated in the model include bulk flow in and out of a well-mixed reactor, transport of dissolved species into the biofilm, substrate consumption by bacterial metabolism, bacterial growth, advection of cell mass within the biofilm, cell detachment from the biofilm, cell death, and biocide concentration-dependent disinfection. Simulations were performed to analyze the general behavior of the model and to perform preliminary sensitivity analysis to identify key input parameters. The model captured several general features of antimicrobial agent action against biofilms that have been observed widely by experimenters and practitioners. These included (1) rapid disinfection followed by biofilm regrowth, (2) slower detachment than disinfection, and (3) reduced susceptibility of microorganisms in biofilms. The results support the plausibility of a mechanism of biofilm resistance in which the biocide is neutralized by reaction with biofilm constituents, leading to a reduction in the bulk biocide concentration and, more significantly, biocide concentration gradients within the biofilm. Sensitivity experiments and analyses identified which input parameters influence key response variables. Each of three response variables was sensitive to each of the five input parameters, but they were most sensitive to the initial biofilm thickness and next most sensitive to the biocide disinfection rate coefficient. Statistical regression modeling produced simple equations for approximating the response variables for situations within the range of conditions covered by the sensitivity experiment. The model should be useful as a tool for studying alternative biocide control strategies. For example, the simulations suggested that a good interval between pulses of biocide is the time to minimum thickness. (c) 1996 John Wiley & Sons, Inc.     

In vitro model systems for exploring oral biofilms: From single-species populations to complex multi-species communities

Numerous in vitro biofilm model systems are available to study oral biofilms. Over the past several decades, increased understanding of oral biology and advances in technology have facilitated more accurate simulation of intraoral conditions and have allowed for the increased generalizability of in vitro oral biofilm studies. The integration of contemporary systems with confocal microscopy and 16S rRNA community profiling has enhanced the capabilities of in vitro biofilm model systems to quantify biofilm architecture and analyse microbial community composition. In this review, we describe several model systems relevant to modern in vitro oral biofilm studies: the constant depth film fermenter, Sorbarod perfusion system, drip–flow reactor, modified Robbins device, flowcells and microfluidic systems. We highlight how combining these systems with confocal microscopy and community composition analysis tools aids exploration of oral biofilm development under different conditions and in response to antimicrobial/anti-biofilm agents. The review closes with a discussion of future directions for the field of in vitro oral biofilm imaging and analysis.     

A new approach to model the spatiotemporal development of biofilm phase in porous media

Bacteria can exist within biofilms that are attached to the solid matrix of a porous medium. Under certain conditions, the biomass can fully occupy the pore space leading to reduced hydraulic conductivity and mass transport. Here, by treating biofilm as a growing, high-viscosity phase, a novel macroscopic approach to model biofilm spatial expansion and its corresponding effects on porous medium hydraulic properties is presented. The separate yet coupled flow of the water and biofilm phases is handled by using relative permeability curves that allow for biofilm movement within the porous medium and bioclogging effects. Fluid flow is governed by Darcy's law and component transport is set by the convection-diffusion equation reaction terms for each component. Here, the system of governing equations is solved by using a commercial multiphase flow reservoir simulator, which is used to validate the model against published laboratory experiments. A comparison of the model and experimental observations reveal that the model provides a reasonable means to predict biomass development in the porous medium. The results reveal that coupled flow of water and movement of biofilm, as described by relative permeability curves, is complex and has a large impact on the development of biomass and consequent bioclogging in the porous medium.     

Formation and Rapid Export of the Monochlorobimane–Glutathione Conjugate in Cultured Rat Astrocytes

Monochlorobimane (MCB) is often used to visualize glutathione (GSH) levels in cultured cells, since it is quickly converted to a fluorescent GSH conjugate (GS–MCB). To test for consequences of MCB application on the GSH metabolism of astrocytes, we have studied rat astrocyte-rich primary cultures as model system. MCB caused a concentration dependent rapid decrease in the cellular GSH content. Simultaneously, a transient accumulation of GS–MCB in the cells was observed with a maximal content 5 min after MCB application. The cellular accumulation was followed by a rapid release of GS–MCB into the medium with a maximal initial export rate of 27.9 ± 6.5 nmol h⁻¹ mg protein⁻¹. Transporters of the family of multidrug resistance proteins (Mrps) are likely to be involved in this export, since the Mrp inhibitor MK571 lowered the export rate by 60%. These data demonstrate that, due to its rapid export from astrocytes, GS–MCB is only under well-defined conditions a reliable indicator of the cellular GSH concentration and that MK571 can be used to maintain maximal GS–MCB levels in astrocytes.     

Kinetic study of the velum formation by Saccharomyces cerevisiae (beticus ssp.) during the biological aging of wines

The main objetive of this work was to evaluate and model the biofilm growth of the Saccharomyces cerevisiae (beticus ssp.) yeast during the biological aging of some types of wines. Thus, we have study how the biofilm growth, the glycerine is consumed and the acetaldehyde is produced, and how this phenomena are affected by the media ethanol concentration (0–17%, v/v), under experimental conditions similar to the industrial ones. In consequence, the growth of the S. cerevisiae (beticus ssp.) biofilm on the surface of the liquid was studied and kinetically modelled. Growth curves were fitted by using general kinetic models that include biomass and substrate inhibition factors. The alcohol content of the medium for the fastest growth rate of biofilm was found to be 4.3%, v/v. The proposed kinetic models for biomass growth, glycerine consumption and acetaldehyde formation fit well with the experimental data.The growth kinetics of S. cerevisiae beticus ssp. in biofilm phase presents a typical discontinuous microbial growth profile (with lag, exponential and stationary phases). The glycerine consumption is directly related to the substrate concentrations (ethanol and glycerine). Finally, the rate of acetaldehyde formation suggests a model associated with the rate of microbial growth, which is modified by a substrate-dependent factor. The suggested model can be used for optimization and control processes of biological aging of wines.     

Mass transfer enhancement in moving biofilm structures

Biofilms are layers of microbial cells growing on an interface and they can form highly complex structures adapted to a wide variety of environmental conditions. Biofilm streamers have a small immobile base attached to the support and a flexible tail elongated in the flow direction, which can vibrate in fast flows. Herein we report numerical results for the role of the periodical movement of biofilm streamers on the nutrient uptake and in general on the solute mass transfer enhancement due to flow-induced oscillations. We developed what to our knowledge is a novel two-dimensional fluid-structure interaction model coupled to unsteady solute mass transport and solved the model using the finite element method with a moving mesh. Results demonstrate that the oscillatory movement of the biofilm tail significantly increases the substrate uptake. The mass transfer coefficient is the highest in regions close to the streamer tip. The reason for substrate transfer enhancement is the increase in speed of tip movement relative to the surrounding liquid, thereby reducing the thickness of the mass transfer boundary layer. In addition, we show that the relative mass transfer enhancement in unsteady conditions compared with the rigid static structure is larger at higher flow velocities, and this relative increase favors a more flexible structure.     

Persistence in a Single Species CSTR Model with Suspended Flocs and Wall Attached Biofilms

We consider a mathematical model for a bacterial population in a continuously stirred tank reactor (CSTR) with wall attachment. This is a modification of the Freter model, in which we model the sessile bacteria as a microbial biofilm. Our analysis indicates that the results of the algebraically simpler original Freter model largely carry over. In a computational simulation study, we find that the vast majority of bacteria in the reactor will eventually be sessile. However, we also find that suspended biomass is relatively more efficient in removing substrate from the reactor than biofilm bacteria.