Yield prediction and stoichiometry of multi-step biodegradation reactions involving oxygenation

Competitive and substrate limited environments drive metabolic heterogeneity for comammox Nitrospira

Nitrospira has been revealed as a high versatile genus. Although previously considered only responsible for the conversion of nitrite to nitrate, now we know that Nitrospira can perform complete ammonia oxidation to nitrate too (comammox). Comammox activity was firstly reported as dominant in extremely limited oxygen environments, where anaerobic ammonia oxidation was also occurring (anammox). To explain the comammox selection, we developed an Individual-based Model able to describe Nitrospira and anammox growth in suspended flocs assembled in a dynamic nitrogen and oxygen-limiting environment. All known and hypothesized nitrogen transformations of Nitrospira were considered: ammonia and nitrite oxidation, comammox, nitrate-reducing ammonia oxidation, and anaerobic nitrite-reducing ammonia oxidation. Through bioenergetics analysis, the growth yield associated to each activity was estimated. The other kinetic parameters necessary to describe growth were calibrated according to the reported literature values. Our modeling results suggest that even extremely low oxygen concentrations (~1.0 µM) allow for a proportional growth of anammox versus Nitrospira similar to the one experimentally observed. The strong oxygen limitation was followed by a limitation of ammonia and nitrite, because anammox, without strong competitors, were able to grow faster than Nitrospira depleting the environment in nitrogen. These substrate limitations created an extremely competitive environment that proved to be decisive in the community assembly of Nitrospira and anammox. Additionally, a diversity of metabolic activities for Nitrospira was observed in all tested conditions, which in turn, explained the transient nitrite accumulation observed in aerobic environments with higher ammonia availability.

Stoichiometric ratios for biotics and xenobiotics capture effective metabolic coupling to re(de)fine biodegradation

Preserving human and environmental health requires anthropogenic pollutants to be biologically degradable. Depending on concentration, both nutrients and pollutants induce and activate metabolic capacity in the endemic bacterial consortium, which in turn aids their degradation. Knowledge on such 'acclimation' is rarely implemented in risk assessment cost-effectively. As a result, an accurate description of the mechanisms and kinetics of biodegradation remains problematic. In this study, we defined a yield 'effectivity', comprising the effectiveness at which a pollutant (substrate) enhances its own degradation by inducing (biomass) cofactors involved therein. Our architecture for calculation represents the interplay between concentration and metabolism via both stoichiometric and thermodynamic concepts. The calculus for yield 'effectivity' is biochemically intuitive, implicitly embeds co-metabolism and distinguishes 'endogenic' from 'exogenic' substances' reflecting various phenomena in biodegradation and bio-transformation studies. We combined data on half-lives of pollutants/nutrients in wastewater and surface water with transition-state rate theory to obtain also experimental values for effective yields. These quantify the state of acclimation: the portion of biodegradation kinetics attributable to (contributed by) 'natural metabolism', in view of similarity to natural substances. Calculated and experimental values showed statistically significant correspondence. Particularly, carbohydrate metabolism and nucleic acid metabolism appeared relevant for acclimation (R² = 0.11-0.42), affecting rates up to 10^4.9(±0.7) times: under steady-state acclimation, a compound stoichiometrically identical to carbohydrates or nucleic acids, is 10^3.2 to 10^4.9 times faster aerobically degraded than a compound marginally similar. Our new method, simulating (contribution by) the state of acclimation, supplements existing structure-biodegradation and kinetic models for predicting biodegradation in wastewater and surface water. The accuracy of prediction may increase when characterizing nutrients/co-metabolites in terms of, e.g., elemental analysis. We discuss strengths and limitations of our approach by comparison to empirical and mechanism-based methods.

Growth of nitrite-oxidizing Nitrospira and ammonia-oxidizing Nitrosomonas in marine recirculating trickling biofilter reactors

Aerobic ammonia and nitrite oxidation reactions are fundamental biogeochemical reactions contributing to the global nitrogen cycle. Although aerobic nitrite oxidation yields 4.8-folds less Gibbs free energy (∆G_r ) than aerobic ammonia oxidation in the NH₄ ⁺ -feeding marine recirculating trickling biofilter reactors operated in the present study, nitrite-oxidizing and not ammonia-oxidizing Nitrospira (sublineage IV) outnumbered ammonia-oxidizing Nitrosomonas (relative abundance; 53.8% and 7.59% respectively). CO₂ assimilation efficiencies during ammonia or nitrite oxidation were 0.077 μmol-¹⁴ CO₂ /μmol-NH₃ and 0.053-0.054 μmol-¹⁴ CO₂ /μmol-NO₂ ^- respectively, and the difference between ammonia and nitrite oxidation was much smaller than the difference of ∆G_r . Free-energy efficiency of nitrite oxidation was higher than ammonia oxidation (31%-32% and 13% respectively), and high CO₂ assimilation and free-energy efficiencies were a determinant for the dominance of Nitrospira over Nitrosomonas. Washout of Nitrospira and Nitrosomonas from the trickling biofilter reactors was also examined by quantitative PCR assay. Normalized copy numbers of Nitrosomonas amoA were 1.5- to 1.7-folds greater than Nitrospira nxrB and 16S rRNA gene in the reactor effluents. Nitrosomonas was more susceptible for washout than Nitrospira in the trickling biofilter reactors, which was another determinant for the dominance of Nitrospira in the trickling biofilter reactors.

Bioenergetics analysis of ammonia-oxidizing bacteria and the estimation of their maximum growth yield

The currently accepted biochemistry and bioenergetics of ammonia-oxidizing bacteria (AOB) show an inefficient metabolism: only 53.8% of the energy released when a mole of ammonia is oxidised and less than two of the electrons liberated can be directed to the autotrophic anabolism. However, paradoxically, AOB seem to thrive in challenging conditions: growing readily in virtually most aerobic environment, yet limited AOB exist in pure culture. In this study, a comprehensive model of the biochemistry of the metabolism of AOB is presented. Using bioenergetics calculations and selecting the minimum estimation for the energy dissipated in each of the metabolic steps, the model predicts the highest possible true yield of 0.16 gBio/gN and a yield of 0.13 gBio/gN when cellular maintenance is considered. Observed yields should always be lower than these values but the range of experimental values in literature vary between 0.04 and 0.45 gBio/gN. In this work, we discuss if this variance of observed values for AOB growth yield could be understood if other non-considered alternative energy sources are present in the biochemistry of AOB. We analyse how the predicted maximum growth yield of AOB changes considering co-metabolism, the use of hydroxylamine as a substrate, the abiotic oxidation of NO, energy harvesting in the monooxygenase enzyme or the use of organic carbon sources.

Engineering the Dark Food Chain

Meeting global food needs in the face of climate change and resource limitation requires innovative approaches to food production. Here, we explore incorporation of new dark food chains into human food systems, drawing inspiration from natural ecosystems, the history of single cell protein, and opportunities for new food production through wastewater treatment, microbial protein production, and aquaculture. The envisioned dark food chains rely upon chemoautotrophy in lieu of photosynthesis, with primary production based upon assimilation of CH₄ and CO₂ by methane- and hydrogen-oxidizing bacteria. The stoichiometry, kinetics, and thermodynamics of these bacteria are evaluated, and opportunities for recycling of carbon, nitrogen, and water are explored. Because these processes do not require light delivery, high volumetric productivities are possible; because they are exothermic, heat is available for downstream protein processing; because the feedstock gases are cheap, existing pipeline infrastructure could facilitate low-cost energy-efficient delivery in urban environments. Potential life-cycle benefits include: a protein alternative to fishmeal; partial decoupling of animal feed from human food; climate change mitigation due to decreased land use for agriculture; efficient local cycling of carbon and nutrients that offsets the need for energy-intensive fertilizers; and production of high value products, such as the prebiotic polyhydroxybutyrate.

Microbial growth yield estimates from thermodynamics and its importance for degradation of pesticides and formation of biogenic non-extractable residues

In biodegradation studies with isotope-labelled pesticides, fractions of non-extractable residues (NER) remain, but their nature and composition is rarely known, leading to uncertainty about their risk. Microbial growth leads to incorporation of carbon into the microbial mass, resulting in biogenic NER. Formation of microbial mass can be estimated from the microbial growth yield, but experimental data is rare. Instead, we suggest using prediction methods for the theoretical yield based on thermodynamics. Recently, we presented the Microbial Turnover to Biomass (MTB) method that needs a minimum of input data. We have estimated the growth yield of 40 organic chemicals (31 pesticides) using the MTB and two existing methods. The results were compared to experimental values, and the sensitivity of the methods was assessed. The MTB method performed best for pesticides. Having the theoretical yield and using the released CO₂ as a measure for microbial activity, we predicted a range for the formation of biogenic NER. For the majority of the pesticides, a considerable fraction of the NER was estimated to be biogenic. This novel approach provides a theoretical foundation applicable to the evaluation and prediction of biogenic NER formation during pesticide degradation experiments, and may also be employed for the interpretation of NER data from regulatory studies.

Thermodynamic analysis of biodegradation pathways

Microorganisms provide a wealth of biodegradative potential in the reduction and elimination of xenobiotic compounds in the environment. One useful metric to evaluate potential biodegradation pathways is thermodynamic feasibility. However, experimental data for the thermodynamic properties of xenobiotics is scarce. The present work uses a group contribution method to study the thermodynamic properties of the University of Minnesota Biocatalysis/Biodegradation Database. The Gibbs free energies of formation and reaction are estimated for 914 compounds (81%) and 902 reactions (75%), respectively, in the database. The reactions are classified based on the minimum and maximum Gibbs free energy values, which accounts for uncertainty in the free energy estimates and a feasible concentration range relevant to biodegradation. Using the free energy estimates, the cumulative free energy change of 89 biodegradation pathways (51%) in the database could be estimated. A comparison of the likelihood of the biotransformation rules in the Pathway Prediction System and their thermodynamic feasibility was then carried out. This analysis revealed that when evaluating the feasibility of biodegradation pathways, it is important to consider the thermodynamic topology of the reactions in the context of the complete pathway. Group contribution is shown to be a viable tool for estimating, a priori, the thermodynamic feasibility and the relative likelihood of alternative biodegradation reactions. This work offers a useful tool to a broad range of researchers interested in estimating the feasibility of the reactions in existing or novel biodegradation pathways.

Thermodynamic electron equivalents model for bacterial yield prediction: modifications and comparative evaluations

Modifications are made to an earlier thermodynamic model (TEEM1) for prediction of maximum microbial yields from aerobic and anaerobic as well as heterotrophic and autotrophic growth. The revised model (TEEM2) corrects for lower yields found with aerobic oxidations of organic compounds where an oxygenase is involved and with growth on single-carbon (C1) compounds. TEEM1 and TEEM2 are based on energy release and consumption as determined from the reduction potential or Gibbs free energy of (1/2)-reaction reduction equations together with losses of energy during energy transfer. Energy transfer efficiency is a key parameter needed to make predictions with TEEM2, and was determined through evaluations with extensive data sets on aerobic heterotrophic yield available in the literature. For compounds following normal catabolic pathways, the best-fit value for energy transfer efficiency was 0.37, which permitted accurate predictions of growth with a precision of 15%-20% as determined by standard deviation. Using the same energy transfer efficiency, a similar precision, but somewhat less accuracy was found for organic compounds where oxidation involves an oxygenase (estimates 8% too high) and for C1 compounds (estimates 17% too high). In spite of the somewhat lower accuracy, the TEEM2 modifications resulted in improved predictions over TEEM1 and the comparison models.

Predicting the accumulation of harmful metabolic byproducts during the treatment of VOC emissions in suspended growth bioreactors

A predicting model is proposed to evaluate metabolic byproducts accumulation and process performance in suspended growth reactors treating air emissions contaminated with volatile organic compounds (VOCs). The model presented integrates a multistep kinetic model and a general mechanistic model describing bioreactor operation. This integrated model is based on general equations modeling, both mass transport and the mechanisms underlying pollutant biotransformation and byproducts accumulation, and can be applied to a wide range of operating conditions (VOC substrate, O2, and nutrients limitation). The model was tested for predicting benzyl alcohol (BA) accumulation in a chemostat reactor treating toluene. BA accumulates in Pseudomonas putida F1 cultures degrading toluene as a result of methyl monooxygenation reaction parallel to the main TOD degradation pathway. The operational conditions leading to BA accumulation are evaluated through simulations assays. Simulation results indicate that BA accumulation occurs when other substrates rather than toluene are limiting. Therefore, operation under toluene limitation is highly recommended to ensure not only the detoxification goals but also to avoid potential mutagenic effects of BA over the microbial culture.

Expanded thermodynamic true yield prediction model: adjustments and limitations

Bacterial yield prediction is critical for bioprocess optimization and modeling of natural biological systems. In previous work, an expanded thermodynamic true yield prediction model was developed through incorporating carbon balance and nitrogen balance along with electron balance and energy balance. In the present work, the application of the expanded model is demonstrated in multiple growth situations (aerobic heterotrophs, anoxic, anaerobic heterotrophs, and autolithotrophs). Two adjustments are presented that enable improved prediction when additional information regarding the environmental conditions (pH) or degradation pathway (requirement for oxygenase- or oxidase-catalyzed reactions) is known. A large data set of reported yields is presented and considered for suitability in model validation. Significant uncertainties of literature-reported yield values are described. Evaluation of the model with experimental yield values shows good predictive ability. However, the wide range in reported yields and the variability introduced into the prediction by uncertainty in model parameters, limits comprehensive validation. Our results suggest that the uncertainty of the experimental data used for validation limits further improvement of thermodynamic prediction models.

Bacterial growth yields on EDTA, NTA, and their biodegradation intermediates

Ethylenediaminetetraacetic acid (EDTA) and nitrilotriacetic acid (NTA) are widely used anthropogenic chelating agents for control of metal speciation and are ubiquitous in natural waters and wastewaters. This is the first report of systematic measurement of the growth yields of a mixed culture (BNC1-BNC2) on EDTA and its biodegradation intermediates, and of Aminobacter aminovorans (aka Chelatobacter heintzii) ATCC 29600 on NTA and its biodegradation intermediates. The yields measured for BNC1-BNC2 co-culture were 75.0 g of cell dry weight (CDW) (mole of EDTA)(-1), 68.6 g of CDW (mole of ED3 A)(-1), 51.2 g of CDW (mole of N,N'-EDDA)(-1), 34.5 g of CDW (mole of ED)(-1), 26.3 g of CDW (mole of IDA)(-1), 12.2 g of CDW (mole of glycine)(-1), and 9.7 g of CDW (mole of glyoxylate)(-1). The yields measured for A. aminovorans were 44.3 g of CDW (mole of NTA)(-1), 37.9 g of CDW (mole of IDA)(-1), 15.2 g of CDW (mole of glycine)(-1), and 10.4 g of CDW (mole of glyoxylate)(-1). The biodegradation pathways of EDTA, NTA, and several of their metabolic intermediates include reactions catalyzed by oxygenase enzymes, which may reduce energy available for cell synthesis. Comparison of measured yields with predicted yields indicates that the effect of oxygenase reaction on cell yield can be quantified experimentally as well as modeled based on thermodynamics.

Direct estimation of the oxygen requirements of Achromobacter xylosoxidans for aerobic degradation of monoaromatic hydrocarbons (BTEX) in a bioscrubber

The O2 requirements for biomass production and supplying maintenance energy demands during the degradation of both benzene and ethylbenzene by Achromobacter xylosoxidans Y234 were measured using a newly proposed technique involving a bioscrubber. Using this approach, relevant microbial parameter estimates were directly and simultaneously obtained via linear regression of pseudo steady-state data. For benzene and ethylbenzene, the biomass yield on O2, Y(X/O2), was estimated on a cell dry weight (CDW) basis as 1.96 +/- 0.25 mg CDW mgO2(-1) and 0.98 +/- 0.17 mg CDW mgO2(-1), while the specific rate of O2 consumption for maintenance, m(O2), was estimated as 0.041 +/- 0.008 mgO(2) mg CDW(-1) h(-1) and 0.053 +/- 0.022 mgO(2) mg CDW(-1) h(-1), respectively.

Thermodynamics of microbial growth and metabolism: An analysis of the current situation

This paper attempts to review in how far thermodynamic analysis can be used to understand and predict the performance of microorganisms with respect to growth and bio-product synthesis. In the first part, a simple thermodynamic model of microbial growth is developed which explains the relationship between the driving force for growth in terms of Gibbs energy dissipation and biomass yield. From the currently available literature, it appears that the Gibbs energy dissipation per C-mol of biomass grown, which represents the driving force for chemotrophic growth, may have been adapted by evolutionary processes to strike a reasonable compromise between metabolic rate and growth efficiency. Based on empirical correlations of the C-molar Gibbs energy dissipation, the wide variety of biomass yields observed in nature can be explained and roughly predicted. This type of analysis may be highly useful in environmental applications, where such wide variations occur. It is however not able to predict biomass yields in very complex systems such as mammalian cells nor is it able to predict or to assess bio-product or recombinant protein yields. For this purpose, a much more sophisticated treatment that accounts for individual metabolic pathways separately is required. Based on glycolysis as a test example, it is shown in the last part that simple thermodynamic analysis leads to erroneous conclusions even in well-known, simple cases. Potential sources for errors have been analyzed and can be used to identify the most important needs for future research.

Expanded thermodynamic model for microbial true yield prediction

Thermodynamic methods to predict true yield and stoichiometry of bacterial reactions have been widely used in biotechnology and environmental engineering. However, yield predictions are often inaccurate for certain simple organic compounds. This work evaluates an existing method and identifies the cause of prediction errors for compounds with low degree of reductance of carbon. For these compounds, carbon, not energy or reducing equivalents, constrains growth. Existing thermodynamically-based models do not account for the potential of carbon-limited growth. The improved method described here consists of four balances: carbon balance, nitrogen balance, electron balance, and energy balance. Two efficiency terms, K1 and K2 are defined and estimated from a priori analysis. The results show that K1 and K2 are nearly the same in value so that only one coefficient, K = 0.41 is used in the modified model. Comparisons with observed yields show that use of the new model and parameters results in significantly improved yield estimation based on inclusion of the carbon balance. The average estimation error is less than 6% for the data set presented.

Phenol and catechol biodegradation by the haloalkaliphile Halomonas campisalis: influence of pH and salinity

Removal of aromatic compounds from alkaline and/or saline industrial wastewater is an environmental concern for industry. In addition, aromatics may be accumulating in soda lakes, unique natural systems, where the fate and toxicity of these contaminants is unknown. To determine the feasibility of aromatic compound biodegradation in saline and alkaline conditions, the effect of pH and salinity on the biodegradation of phenol as a model aromatic waste compound by the haloalkaliphilic bacterium Halomonas campisalis was examined. Phenol was degraded as a source of carbon and energy at pH 8-11 and 0-150 g/L NaCl. Metabolic intermediates catechol, cis,cis-muconate, and (+)-muconolactone were identified, thus indicating that phenol was degraded via the beta-ketoadipate metabolic pathway. Although phenol and catechol were completely degraded in all cases, small amounts of cis,cis-muconate accumulated proportionally to increases in pH. There was no noticeable influence of salinity on cis,cis-muconate accumulation except at 0 g/L NaCl where it was completely degraded. These results indicate that it may be feasible to use haloalkaliphiles forthe treatment of aromatics present in saline and/or alkaline systems. This is the first report of phenol and catechol biodegradation under combined saline and alkaline conditions.