Evaluating alternate biokinetic models for trace pollutant cometabolism

On the reliable estimation of sequential Monod kinetic parameters

While dozens of studies have attempted to estimate the Monod kinetic parameters of microbial reductive dechlorination, published values in the literature vary by 2-6 orders of magnitude. This lack of consensus can be attributed in part to limitations of both experimental design and parameter estimation techniques. To address these issues, Hamiltonian Monte Carlo was used to produce more than one million sets of realistic simulated microcosm data under a variety of experimental conditions. These data were then employed in model fitting experiments using a number of parameter estimation algorithms for determining Monod kinetic parameters. Analysis of data from conventional triplicate microcosms yielded parameter estimates characterized by high collinearity, resulting in poor estimation accuracy and precision. Additionally, confidence intervals computed by commonly used classical regression analysis techniques contained true parameter values much less frequently than their nominal confidence levels. Use of an alternative experimental design, requiring the same number of analyses as conventional experiments but comprised of microcosms with varying initial chlorinated ethene concentrations, is shown to result in order-of-magnitude decreases in parameter uncertainty. A Metropolis algorithm which can be run on a typical personal computer is demonstrated to return more reliable parameter interval estimates.

Bioremediation of chlorinated ethenes contaminated groundwater and the reactive transport modeling - A review

Improper disposal of chlorinated ethenes (CEs), a class of widely used solvents in chemical manufacturing and cleaning industries, often leads to severe groundwater contamination. In situ bioremediation of CE-contaminated groundwater has received continuous attention in recent years. The reactive transport simulation is a valuable tool for planning and designing in situ bioremediation systems. This paper presents a detailed and comprehensive review on the main biotransformation pathways of CEs in aquifers, the mathematical modeling of bioremediation processes, and the available software developed for the simulation of reactive transport of CEs over past three decades. The aim of this research is to provide guidance on the selection of appropriate models and software suitable for systems of varying scales, and to discern prevailing research trends while identifying areas worthy of further study. This paper provides a detailed summary of the equations, parameters, and applications of existing biotransformation models from literature studies, highlighting the operation, benefits, and limitations of software available for CEs reactive transport simulations. Lastly, the support of reactive transport simulation programs for the design of full-scale in situ bioremediation systems was elucidated. Further research is needed for incorporating the effects of key subsurface environmental factors on biodegradation processes into models and balancing model complexity with computer data processing power to better support the development and application of reactive transport modeling software.

Unraveling biogeochemical complexity through better integration of experiments and modeling

The evolution of groundwater quality in natural and contaminated aquifers is affected by complex interactions between physical transport and biogeochemical reactions. Identifying and quantifying the processes that control the overall system behavior is the key driver for experimentation and monitoring. However, we argue that, in contrast to other disciplines in earth sciences, process-based computer models are currently vastly underutilized in the quest for understanding subsurface biogeochemistry. Such models provide an essential avenue for quantitatively testing hypothetical combinations of interacting, complex physical and chemical processes. If a particular conceptual model, and its numerical counterpart, cannot adequately reproduce observed experimental data, its underlying hypothesis must be rejected. This quantitative process of hypothesis testing and falsification is central to scientific discovery. We provide a perspective on how closer interactions between experimentalists and numerical modelers would enhance this scientific process, and discuss the potential limitations that are currently holding us back. We also propose a data-model nexus involving a greater use of numerical process-based models for a more rigorous analysis of experimental observations while also generating the basis for a systematic improvement in the design of future experiments.

Cometabolism of 17α-ethynylestradiol by nitrifying bacteria depends on reducing power availability and leads to elevated nitric oxide formation

17α-ethynylestradiol (EE2) is a priority emerging contaminant (EC) in diverse environments that can be cometabolized by ammonia oxidizing bacteria (AOB). However, its transformation kinetics and the underlying molecular mechanism are unclear. In this study, kinetic parameters, including maximum specific EE2 transformation rate, EE2 half-saturation coefficient, and EE2transformation capacity of AOBwere obtained by using the model AOB strain, Nitrosomonas europaea 19718. The relationship between EE2 cometabolism and ammonia oxidation was divided into three phases according to reducing power availability, namely "activation", "coupling", and "saturation". Specifically, there was a universal lag of EE2 transformation after ammonia oxidation was initiated, suggesting that sufficient reducing power (approximately 0.95 ± 0.06 mol NADH/L) was required to activate EE2 cometabolism. Interestingly, nitric oxide emission increased by 12 ± 2% during EE2 cometabolism, along with significantly upregulated nirK cluster genes. The findings are of importance to understanding the cometabolic behavior and mechanism of EE2 in natural and engineered environments. Maintaining relatively high and stable reducing power supply from ammonia oxidation can potentially improve the cometabolic removal of EE2 and other ECs during wastewater nitrification processes.

Aerobic cometabolism of 1,4-dioxane by isobutane-utilizing microorganisms including Rhodococcus rhodochrous strain 21198 in aquifer microcosms: Experimental and modeling study

Aerobic cometabolism of the emerging contaminant 1,4-dioxane (1,4-D) by isobutane-utilizing microorganisms was assessed in pure culture and aquifer microcosm studies. The bacterium Rhodococcus rhodochrous strain ATCC 21198 transformed low, environmentally-relevant concentrations of 1,4-D when grown on isobutane. Microcosms were constructed with aquifer solids from Fort Carson, Colorado, a site contaminated with 1,4-D and trichloroethene (TCE). Multiple additions of isobutane and 1,4-D over 300 days were transformed in microcosms biostimulated with isobutane and microcosms bioaugmented with strain 21198. Results showed that, over time and with sufficient inorganic nutrients, biostimulation of native microorganisms with isobutane was just as effective as bioaugmentation with strain 21198 to achieve 1,4-D transformation in the microcosms. The presence of TCE at 0.2 mg/L did not inhibit 1,4-D transformation, though TCE itself was not readily transformed. An iterative process was used to determine kinetic parameter values to fit Michaelis-Menten/Monod models to experimental data for simultaneous isobutane utilization, biomass growth, and cometabolic transformation of 1,4-D. Parameter optimization resulted in good model fit to the data over multiple transformations of isobutane and 1,4-D in both short- and long-term experiments. Results suggest low concentrations of 1,4-D studied in the microcosms were cometabolically transformed according to a pseudo first-order rate of 0.37 L/mg TSS/day of 21198. Isobutane consumption was modeled with a maximum rate of 2.58 mg/mg TSS/day and a half saturation constant of 0.09 mg/L. 1,4-D transformation was competitively inhibited by the presence of isobutane and transformation rates were significantly reduced when inorganic nutrients were limiting. Simulations of the repeated additions found a first-order microbial endogenous decay coefficient of 0.03 day^-1 fit the alternating periods of active transformation and stagnation between isobutane and 1,4-D additions over approximately one year. The model fitting process highlighted the importance of determining kinetic parameters from data representing low concentrations typically found in the environment.

Roles of ammonia-oxidizing bacteria in improving metabolism and cometabolism of trace organic chemicals in biological wastewater treatment processes: A review

While there has been a significant recent improvement in the removal of pollutants in natural and engineered systems, trace organic chemicals (TrOCs) are posing a major threat to aquatic environments and human health. There is a critical need for developing potential strategies that aim at enhancing metabolism and/or cometabolism of these compounds. Recently, knowledge regarding biodegradation of TrOCs by ammonia-oxidizing bacteria (AOB) has been widely developed. This review aims to delineate an up-to-date version of the ecophysiology of AOB and outline current knowledge related to biodegradation efficiencies of the frequently reported TrOCs by AOB. The paper also provides an insight into biodegradation pathways by AOB and transformation products of these compounds and makes recommendations for future research of AOB. In brief, nitrifying WWTFs (wastewater treatment facilities) were superior in degrading most TrOCs than non-nitrifying WWTFs due to cometabolic biodegradation by the AOB. To fully understand and/or enhance the cometabolic biodegradation of TrOCs by AOB, recent molecular research has focused on numerous crucial factors including availability of the compounds to AOB, presence of growth substrate (NH₄-N), redox potentials, microorganism diversity (AOB and heterotrophs), physicochemical properties and operational parameters of the WWTFs, molecular structure of target TrOCs and membrane-based technologies, may all significantly impact the cometabolic biodegradation of TrOCs. Still, further exploration is required to elucidate the mechanisms involved in biodegradation of TrOCs by AOB and the toxicity levels of formed products.

Transcriptional response of Mycobacterium sp. strain A1-PYR to multiple polycyclic aromatic hydrocarbon contaminations

Cometabolism mechanisms of organic pollutants in environmental microbes have not been fully understood. In this study, a global analysis of Mycobacterium sp. strain A1-PYR transcriptomes on different PAH substrates (single or binary of pyrene (PYR) and phenanthrene (PHE)) was conducted. Comparative results demonstrated that expression levels of 23 PAH degradation enzymes were significantly higher in the binary substrate than in the PYR-only one. These enzymes constituted an integrated enzymatic system to actualize all transformation steps of PYR, and most of their encoded genes formed a novel gene cascade in the genome of strain A1-PYR. The roles of different genotypes of enzymes in PYR cometabolism were also discriminated even though all of their gene sequences were presented in the genome of this strain. NidAB and PdoA2B2 instead of NidA3B3 served the initial oxidization of PAHs, and PcaL replaced PcaCD to catalyze the formation of 3-oxoadipate. Novel genes associated with PYR cometabolism was also predicted by the relationships between their transcription profiles and PYR removal. The results showed that ABC-type transporters probably played important roles in the transport of PAHs and their metabolites through cell membrane, and [4Fe-4S] ferredoxin might be essential for dioxygenases (NidAB and PdoA2B2) to achieve oxidative activities. This study provided molecular insight in that microbial degrader subtly cometabolized recalcitrant PAHs with relatively more degradable ones.

Why are organic micropollutants not fully biotransformed? A mechanistic modelling approach to anaerobic systems

Biotransformation of most organic micropollutants (OMPs) during wastewater treatment is not complete and an unexplained steady decrease of the biotransformation rate with time is reported for many OMPs in different biological processes. To minimize and accurately predict the emission of OMPs into the environment, the mechanisms and limitations behind their biotransformations should be clarified. Aiming to achieve this objective, the present study follows a mechanistic modelling approach, based on the formulation of four models according to different biotransformation hypotheses: Michaelis-Menten kinetics, chemical equilibrium between the parent compound and the transformation product (TP), enzymatic inhibition by the TP, and a limited compound bioavailability due to its sequestration in the solid phase. These models were calibrated and validated with kinetic experiments performed in two different anaerobic systems: continuous reactors enriched with methanogenic biomass and batch assays with anaerobic sludge. Model selection was conducted according to model suitability criteria (goodness of fitting the experimental data, confidence of the estimated parameters, and model parsimony) but also considering mechanistic evidences. The findings suggest that reversibility of the biological reactions and/or sequestration of compounds are likely the causes preventing the complete biotransformation of OMPs, and biotransformation is probably limited by thermodynamics rather than by kinetics. Taking into account its simplicity and broader applicability spectrum, the reversible biotransformation is the proposed model to explain the incomplete biotransformation of OMPs.

Modeling of Pharmaceutical Biotransformation by Enriched Nitrifying Culture under Different Metabolic Conditions

Pharmaceutical removal could be significantly enhanced through cometabolism during nitrification processes. To date, pharmaceutical biotransformation models have not considered the formation of transformation products associated with the metabolic type of microorganisms. Here we report a comprehensive model to describe and evaluate the biodegradation of pharmaceuticals and the formation of their biotransformation products by enriched nitrifying cultures. The biotransformation of parent compounds was linked to the microbial processes via cometabolism induced by ammonium-oxidizing bacteria (AOB) growth, metabolism by AOB, cometabolism by heterotrophs (HET) growth, and metabolism by HET in the model framework. The model was calibrated and validated using experimental data from pharmaceutical biodegradation experiments at realistic levels, taking two pharmaceuticals as examples, i.e., atenolol and acyclovir. Results demonstrated the good predictive performance of the established biotransformation model under different metabolic conditions, as well as the reliability of the established model in predicting different pharmaceutical biotransformations. The linear positive correlation between ammonia oxidation rate and pharmaceutical degradation rate confirmed the major role of cometabolism induced by AOB in the pharmaceutical removal. Dissolved oxygen was also revealed to be capable of regulating the pharmaceutical biotransformation cometabolically, and the substrate competition between ammonium and pharmaceuticals existed especially at high ammonium concentrations.

Identification of Biokinetic Models Using the Concept of Extents

The development of a wide array of process technologies to enable the shift from conventional biological wastewater treatment processes to resource recovery systems is matched by an increasing demand for predictive capabilities. Mathematical models are excellent tools to meet this demand. However, obtaining reliable and fit-for-purpose models remains a cumbersome task due to the inherent complexity of biological wastewater treatment processes. In this work, we present a first study in the context of environmental biotechnology that adopts and explores the use of extents as a way to simplify and streamline the dynamic process modeling task. In addition, the extent-based modeling strategy is enhanced by optimal accounting for nonlinear algebraic equilibria and nonlinear measurement equations. Finally, a thorough discussion of our results explains the benefits of extent-based modeling and its potential to turn environmental process modeling into a highly automated task.

Impact of external carbon dose on the removal of micropollutants using methanol and ethanol in post-denitrifying Moving Bed Biofilm Reactors

Addition of external carbon sources to post-denitrification systems is frequently used in wastewater treatment plants to enhance nitrate removal. However, little is known about the fate of micropollutants in post-denitrification systems and the influence of external carbon dosing on their removal. In this study, we assessed the effects of two different types and availability of commonly used carbon sources -methanol and ethanol- on the removal of micropollutants in biofilm systems. Two laboratory-scale moving bed biofilm reactors (MBBRs), containing AnoxKaldnes K1 carriers with acclimated biofilm from full-scale systems, were operated in continuous-flow using wastewater dosed with methanol and ethanol, respectively. Batch experiments with 22 spiked pharmaceuticals were performed to assess removal kinetics. Acetyl-sulfadiazine, atenolol, citalopram, propranolol and trimethoprim were easily biotransformed in both MBBRs (biotransformations rate constants k_bio between 1.2 and 12.9 L g_biomass^-1 d^-1), 13 compounds were moderately biotransformed (rate constants between 0.2 and 2 L g_biomass^-1 d^-1) and 4 compounds were recalcitrant. The methanol-dosed MBBR showed higher k_bio (e.g., 1.5-2.5-fold) than in the ethanol-dosed MBBR for 9 out of the 22 studied compounds, equal k_bio for 10 compounds, while 3 compounds (i.e., targeted sulfonamides) were biotransformed faster in the ethanol-dosed MBBR. While biotransformation of most of the targeted compounds followed first-order kinetics, removal of venlafaxine, carbamazepine, sulfamethoxazole and sulfamethizole could be described with a cometabolic model. Analyses of the microbial composition in the biofilms using 16S rRNA amplicon sequencing revealed that the methanol-dosed MBBR contained higher microbial richness than the one dosed with ethanol, suggesting that improved biotransformation of targeted compounds could be associated with higher microbial richness. During continuous-flow operation, at conditions representative of full-scale denitrification systems (hydraulic residence time = 2 h), the removal efficiencies of micropollutants were below 35% in both MBBRs, with the exception of atenolol and trimethoprim (>80%). Overall, this study demonstrated that MBBRs used for post-denitrification could be optimized to enhance the biotransformation of a number of micropollutants by accounting for optimal carbon sources and extended residence time.

Removal of Antibiotics in Biological Wastewater Treatment Systems-A Critical Assessment Using the Activated Sludge Modeling Framework for Xenobiotics (ASM-X)

Many scientific studies present removal efficiencies for pharmaceuticals in laboratory-, pilot-, and full-scale wastewater treatment plants, based on observations that may be impacted by theoretical and methodological approaches used. In this Critical Review, we evaluated factors influencing observed removal efficiencies of three antibiotics (sulfamethoxazole, ciprofloxacin, tetracycline) in pilot- and full-scale biological treatment systems. Factors assessed include (i) retransformation to parent pharmaceuticals from e.g., conjugated metabolites and analogues, (ii) solid retention time (SRT), (iii) fractions sorbed onto solids, and (iv) dynamics in influent and effluent loading. A recently developed methodology was used, relying on the comparison of removal efficiency predictions (obtained with the Activated Sludge Model for Xenobiotics (ASM-X)) with representative measured data from literature. By applying this methodology, we demonstrated that (a) the elimination of sulfamethoxazole may be significantly underestimated when not considering retransformation from conjugated metabolites, depending on the type (urban or hospital) and size of upstream catchments; (b) operation at extended SRT may enhance antibiotic removal, as shown for sulfamethoxazole; (c) not accounting for fractions sorbed in influent and effluent solids may cause slight underestimation of ciprofloxacin removal efficiency. Using tetracycline as example substance, we ultimately evaluated implications of effluent dynamics and retransformation on environmental exposure and risk prediction.

Biotransformation of pharmaceuticals by ammonia oxidizing bacteria in wastewater treatment processes

Pharmaceutical residues could potentially pose detrimental effects on aquatic ecosystems and human health, with wastewater treatment being one of the major pathways for pharmaceuticals to enter into the environment. Enhanced removal of pharmaceuticals by ammonia oxidizing bacteria (AOB) has been widely observed in wastewater treatment processes. This article reviews the current knowledge on the biotransformation of pharmaceuticals by AOB. The relationship between the pharmaceuticals removal and nitrification process was revealed. The important role of AOB-induced cometabolism on the biotransformation of pharmaceuticals as well as their transformation products and pathways was elucidated. Kinetics and mathematical models describing the biotransformation of pharmaceuticals by AOB were also reviewed. The results highlighted the high degradation capabilities of AOB toward some refractory pharmaceuticals, with their degradations being clearly related to the nitrification rate and their transformation products being identified, which may exhibit similar or higher ecotoxicological impacts compared to the parent compound.

Kinetics of aerobic cometabolic biodegradation of chlorinated and brominated aliphatic hydrocarbons: A review

This review analyses kinetic studies of aerobic cometabolism (AC) of halogenated aliphatic hydrocarbons (HAHs) from 2001-2015 in order to (i) compare the different kinetic models proposed, (ii) analyse the estimated model parameters with a focus on novel HAHs and the identification of general trends, and (iii) identify further research needs. The results of this analysis show that aerobic cometabolism can degrade a wide range of HAHs, including HAHs that were not previously tested such as chlorinated propanes, highly chlorinated ethanes and brominated methanes and ethanes. The degree of chlorine mineralization was very high for the chlorinated HAHs. Bromine mineralization was not determined for studies with brominated aliphatics. The examined research period led to the identification of novel growth substrates of potentially high interest. Decreasing performance of aerobic cometabolism were found with increasing chlorination, indicating the high potential of aerobic cometabolism in the presence of medium- and low-halogenated HAHs. Further research is needed for the AC of brominated aliphatic hydrocarbons, the potential for biofilm aerobic cometabolism processes, HAH-HAH mutual inhibition and the identification of the enzymes responsible for each aerobic cometabolism process. Lastly, some indications for a possible standardization of future kinetic studies of HAH aerobic cometabolism are provided.