Exploring Bacterial Carboxylate Reductases for the Reduction of Bifunctional Carboxylic Acids.

Carboxylic acid reductases (CARs) selectively reduce carboxylic acids to aldehydes using ATP and NADPH as cofactors under mild conditions. Although CARs attracts significant interest, only a few enzymes have been characterized to date, whereas the vast majority of CARs have yet to be examined. Herein the authors report that 12 bacterial CARs reduces a broad range of bifunctional carboxylic acids containing oxo-, hydroxy-, amino-, or second carboxyl groups with several enzymes showing activity toward 4-hydroxybutanoic (4-HB) and adipic acids. These CARs exhibits significant reductase activity against substrates whose second functional group is separated from the carboxylate by at least three carbons with both carboxylate groups being reduced in dicarboxylic acids. Purified CARs supplemented with cofactor regenerating systems (for ATP and NADPH), an inorganic pyrophosphatase, and an aldo-keto reductase catalyzes a high conversion (50-76%) of 4-HB to 1,4-butanediol (1,4-BDO) and adipic acid to 1,6-hexanediol (1,6-HDO). Likewise, Escherichia coli strains expressing eight different CARs efficiently reduces 4-HB to 1,4-BDO with 50-95% conversion, whereas adipic acid is reduced to a mixture of 6-hydroxyhexanoic acid (6-HHA) and 1,6-HDO. Thus, our results illustrate the broad biochemical diversity of bacterial CARs and their compatibility with other enzymes for applications in biocatalysis.

The design of long-term effective uranium bioremediation strategy using a community metabolic model.

Acetate amendment at uranium contaminated sites in Rifle, CO. leads to an initial bloom of Geobacter accompanied by the removal of U(VI) from the groundwater, followed by an increase of sulfate-reducing bacteria (SRBs) which are poor reducers of U(VI). One of the challenges associated with bioremediation is the decay in Geobacter abundance, which has been attributed to the depletion of bio-accessible Fe(III), motivating the investigation of simultaneous amendments of acetate and Fe(III) as an alternative bioremediation strategy. In order to understand the community metabolism of Geobacter and SRBs during artificial substrate amendment, we have created a genome-scale dynamic community model of Geobacter and SRBs using the previously described Dynamic Multi-species Metabolic Modeling framework. Optimization techniques are used to determine the optimal acetate and Fe(III) addition profile. Field-scale simulation of acetate addition accurately predicted the in situ data. The simulations suggest that batch amendment of Fe(III) along with continuous acetate addition is insufficient to promote long-term bioremediation, while continuous amendment of Fe(III) along with continuous acetate addition is sufficient to promote long-term bioremediation. By computationally minimizing the acetate and Fe(III) addition rates as well as the difference between the predicted and target uranium concentration, we showed that it is possible to maintain the uranium concentration below the environmental safety standard while minimizing the cost of chemical additions. These simulations show that simultaneous addition of acetate and Fe(III) has the potential to be an effective uranium bioremediation strategy. They also show that computational modeling of microbial community is an important tool to design effective strategies for practical applications in environmental biotechnology.

Genome-scale dynamic modeling of the competition between Rhodoferax and Geobacter in anoxic subsurface environments.

The advent of rapid complete genome sequencing, and the potential to capture this information in genome-scale metabolic models, provide the possibility of comprehensively modeling microbial community interactions. For example, Rhodoferax and Geobacter species are acetate-oxidizing Fe(III)-reducers that compete in anoxic subsurface environments and this competition may have an influence on the in situ bioremediation of uranium-contaminated groundwater. Therefore, genome-scale models of Geobacter sulfurreducens and Rhodoferax ferrireducens were used to evaluate how Geobacter and Rhodoferax species might compete under diverse conditions found in a uranium-contaminated aquifer in Rifle, CO. The model predicted that at the low rates of acetate flux expected under natural conditions at the site, Rhodoferax will outcompete Geobacter as long as sufficient ammonium is available. The model also predicted that when high concentrations of acetate are added during in situ bioremediation, Geobacter species would predominate, consistent with field-scale observations. This can be attributed to the higher expected growth yields of Rhodoferax and the ability of Geobacter to fix nitrogen. The modeling predicted relative proportions of Geobacter and Rhodoferax in geochemically distinct zones of the Rifle site that were comparable to those that were previously documented with molecular techniques. The model also predicted that under nitrogen fixation, higher carbon and electron fluxes would be diverted toward respiration rather than biomass formation in Geobacter, providing a potential explanation for enhanced in situ U(VI) reduction in low-ammonium zones. These results show that genome-scale modeling can be a useful tool for predicting microbial interactions in subsurface environments and shows promise for designing bioremediation strategies.

Direct coupling of a genome-scale microbial in silico model and a groundwater reactive transport model.

The activity of microorganisms often plays an important role in dynamic natural attenuation or engineered bioremediation of subsurface contaminants, such as chlorinated solvents, metals, and radionuclides. To evaluate and/or design bioremediated systems, quantitative reactive transport models are needed. State-of-the-art reactive transport models often ignore the microbial effects or simulate the microbial effects with static growth yield and constant reaction rate parameters over simulated conditions, while in reality microorganisms can dynamically modify their functionality (such as utilization of alternative respiratory pathways) in response to spatial and temporal variations in environmental conditions. Constraint-based genome-scale microbial in silico models, using genomic data and multiple-pathway reaction networks, have been shown to be able to simulate transient metabolism of some well studied microorganisms and identify growth rate, substrate uptake rates, and byproduct rates under different growth conditions. These rates can be identified and used to replace specific microbially-mediated reaction rates in a reactive transport model using local geochemical conditions as constraints. We previously demonstrated the potential utility of integrating a constraint-based microbial metabolism model with a reactive transport simulator as applied to bioremediation of uranium in groundwater. However, that work relied on an indirect coupling approach that was effective for initial demonstration but may not be extensible to more complex problems that are of significant interest (e.g., communities of microbial species and multiple constraining variables). Here, we extend that work by presenting and demonstrating a method of directly integrating a reactive transport model (FORTRAN code) with constraint-based in silico models solved with IBM ILOG CPLEX linear optimizer base system (C library). The models were integrated with BABEL, a language interoperability tool. The modeling system is designed in such a way that constraint-based models targeting different microorganisms or competing organism communities can be easily plugged into the system. Constraint-based modeling is very costly given the size of a genome-scale reaction network. To save computation time, a binary tree is traversed to examine the concentration and solution pool generated during the simulation in order to decide whether the constraint-based model should be called. We also show preliminary results from the integrated model including a comparison of the direct and indirect coupling approaches and evaluated the ability of the approach to simulate field experiment.

Coupling a genome-scale metabolic model with a reactive transport model to describe in situ uranium bioremediation.

The increasing availability of the genome sequences of microorganisms involved in important bioremediation processes makes it feasible to consider developing genome-scale models that can aid in predicting the likely outcome of potential subsurface bioremediation strategies. Previous studies of the in situ bioremediation of uranium-contaminated groundwater have demonstrated that Geobacter species are often the dominant members of the groundwater community during active bioremediation and the primary organisms catalysing U(VI) reduction. Therefore, a genome-scale, constraint-based model of the metabolism of Geobacter sulfurreducens was coupled with the reactive transport model HYDROGEOCHEM in an attempt to model in situ uranium bioremediation. In order to simplify the modelling, the influence of only three growth factors was considered: acetate, the electron donor added to stimulate U(VI) reduction; Fe(III), the electron acceptor primarily supporting growth of Geobacter; and ammonium, a key nutrient. The constraint-based model predicted that growth yields of Geobacter varied significantly based on the availability of these three growth factors and that there are minimum thresholds of acetate and Fe(III) below which growth and activity are not possible. This contrasts with typical, empirical microbial models that assume fixed growth yields and the possibility for complete metabolism of the substrates. The coupled genome-scale and reactive transport model predicted acetate concentrations and U(VI) reduction rates in a field trial of in situ uranium bioremediation that were comparable to the predictions of a calibrated conventional model, but without the need for empirical calibration, other than specifying the initial biomass of Geobacter. These results suggest that coupling genome-scale metabolic models with reactive transport models may be a good approach to developing models that can be truly predictive, without empirical calibration, for evaluating the probable response of subsurface microorganisms to possible bioremediation approaches prior to implementation.