Search

Wastewater-fed microbial fuel cells (MFCs) are a promising technology to treat low-organic carbon wastewater and recover part of the chemical energy in wastewater as electrical power. However, the interactions between electrochemically active and fermentative microorganisms cannot be easily studied in wastewater-fed MFCs because of their complex microbial communities. Defined co-culture MFCs provide a detailed understanding of such interactions. In this study, we characterize the extracellular metabolites in laboratory-scale membrane-less MFCs inoculated with Geobacter sulfurreducens and Escherichia coli co-culture and compare them with pure culture MFCs. G. sulfurreducens MFCs are sparged to maintain anaerobic conditions, while co-culture MFCs rely on E. coli for oxygen removal. G. sulfurreducens MFCs have a power output of 128 mW m(-2) , compared to 63 mW m(-2) from the co-culture MFCs. Analysis of metabolites shows that succinate production in co-culture MFCs decreases current production by G. sulfurreducens and that the removal of succinate is responsible for the increased current density in the late co-culture MFCs. Interestingly, pH adjustment is not required for co-culture MFCs but a base addition is necessary for E. coli MFCs and cultures in vials. Our results show that defined co-culture MFCs provide clear insights into metabolic interactions among bacteria while maintaining a low operational complexity., (© 2013 Wiley Periodicals, Inc.)

Uranium contamination is a serious concern at several sites motivating the development of novel treatment strategies such as the Geobacter-mediated reductive immobilization of uranium. However, this bioremediation strategy has not yet been optimized for the sustained uranium removal. While several reactive-transport models have been developed to represent Geobacter-mediated bioremediation of uranium, these models often lack the detailed quantitative description of the microbial process (e.g., biomass build-up in both groundwater and sediments, electron transport system, etc.) and the interaction between biogeochemical and hydrological process. In this study, a novel multi-scale model was developed by integrating our recent model on electron capacitance of Geobacter (Zhao et al., 2010) with a comprehensive simulator of coupled fluid flow, hydrologic transport, heat transfer, and biogeochemical reactions. This mechanistic reactive-transport model accurately reproduces the experimental data for the bioremediation of uranium with acetate amendment. We subsequently performed global sensitivity analysis with the reactive-transport model in order to identify the main sources of prediction uncertainty caused by synergistic effects of biological, geochemical, and hydrological processes. The proposed approach successfully captured significant contributing factors across time and space, thereby improving the structure and parameterization of the comprehensive reactive-transport model. The global sensitivity analysis also provides a potentially useful tool to evaluate uranium bioremediation strategy. The simulations suggest that under difficult environments (e.g., highly contaminated with U(VI) at a high migration rate of solutes), the efficiency of uranium removal can be improved by adding Geobacter species to the contaminated site (bioaugmentation) in conjunction with the addition of electron donor (biostimulation). The simulations also highlight the interactive effect of initial cell concentration and flow rate on U(VI) reduction.

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., (Published by Elsevier B.V.)

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.

There is a wide diversity of unexplored metabolism encoded in the genomes of microorganisms that have an important environmental role. Genome-scale metabolic modelling enables the individual reactions that are encoded in annotated genomes to be organized into a coherent whole, which can then be used to predict metabolic fluxes that will optimize cell function under a range of conditions. In this Review, we summarize a series of studies in which genome-scale metabolic modelling of Geobacter spp. has resulted in an in-depth understanding of their central metabolism and ecology. A similar iterative modelling and experimental approach could accelerate elucidation of the physiology and ecology of other microorganisms inhabiting a diversity of environments, and could guide optimization of the practical applications of these species.

In situ stimulation of the metabolic activity of Geobacter species through acetate amendment has been shown to be a promising bioremediation strategy to reduce and immobilize hexavalent uranium [U(VI)] as insoluble U(IV). Although Geobacter species are reducing U(VI), they primarily grow via Fe(III) reduction. Unfortunately, the biogeochemistry and the physiology of simultaneous reduction of multiple metals are still poorly understood. A detailed model is therefore required to better understand the pathways leading to U(VI) and Fe(III) reduction by Geobacter species. Based on recent experimental evidence of temporary electron capacitors in Geobacter we propose a novel kinetic model that physically distinguishes planktonic cells into electron-loaded and -unloaded states. Incorporation of an electron load-unload cycle into the model provides insight into U(VI) reduction efficiency, and elucidates the relationship between U(VI)- and Fe(III)-reducing activity and further explains the correlation of high U(VI) removal with high fractions of planktonic cells in subsurface environments. Global sensitivity analysis was used to determine the level of importance of geochemical and microbial processes controlling Geobacter growth and U(VI) reduction, suggesting that the electron load-unload cycle and the resulting repartition of the microbes between aqueous and attached phases are critical for U(VI) reduction. As compared with conventional Monod modeling approaches without inclusion of the electron capacitance, the new model attempts to incorporate a novel cellular mechanism that has a significant impact on the outcome of in situ bioremediation., (Copyright 2009 Elsevier B.V. All rights reserved.)

Background: Production or consumption of protons in growth medium during microbial metabolism plays an important role in determining the pH of the environment. Such pH changes resulting from microbial metabolism may influence the geochemical speciation of many elements in subsurface environments. Protons produced or consumed during microbial growth were measured by determining the amount of acid or base added in a 5 L batch bioreactor equipped with pH control for different species including Escherichia coli, Geobacter sulfurreducens, and Geobacter metallireducens., Results: An in silico model was used to predict the proton secretion or consumption rates and the results were compared with the data. The data was found to confirm predictions of proton consumption during aerobic growth of E. coli with acetate as the carbon source. However, in contrast to proton consumption observed during aerobic growth of E. coli with acetate, proton secretion was observed during growth of Geobacter species with acetate as the donor and Fe(III) as the extracellular electron acceptor., Conclusions: In this study, we have also shown that the final pH of the medium can be either acidic or basic depending on the choice of the electron acceptor for the same electron donor. In all cases, the in silico model could predict qualitatively the proton production/consumption rates obtained from the experimental data. Therefore, measurements of pH equivalents generated or consumed during growth can help characterize the microbial physiology further and can be valuable for optimizing practical applications such as microbial fuel cells, where growth associated pH changes can limit current generation rates.

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., (© 2009 Battelle Memorial Institute. Journal compilation © 2009 Society for Applied Microbiology and Blackwell Publishing Ltd.)

Background: Geobacter metallireducens was the first organism that can be grown in pure culture to completely oxidize organic compounds with Fe(III) oxide serving as electron acceptor. Geobacter species, including G. sulfurreducens and G. metallireducens, are used for bioremediation and electricity generation from waste organic matter and renewable biomass. The constraint-based modeling approach enables the development of genome-scale in silico models that can predict the behavior of complex biological systems and their responses to the environments. Such a modeling approach was applied to provide physiological and ecological insights on the metabolism of G. metallireducens., Results: The genome-scale metabolic model of G. metallireducens was constructed to include 747 genes and 697 reactions. Compared to the G. sulfurreducens model, the G. metallireducens metabolic model contains 118 unique reactions that reflect many of G. metallireducens' specific metabolic capabilities. Detailed examination of the G. metallireducens model suggests that its central metabolism contains several energy-inefficient reactions that are not present in the G. sulfurreducens model. Experimental biomass yield of G. metallireducens growing on pyruvate was lower than the predicted optimal biomass yield. Microarray data of G. metallireducens growing with benzoate and acetate indicated that genes encoding these energy-inefficient reactions were up-regulated by benzoate. These results suggested that the energy-inefficient reactions were likely turned off during G. metallireducens growth with acetate for optimal biomass yield, but were up-regulated during growth with complex electron donors such as benzoate for rapid energy generation. Furthermore, several computational modeling approaches were applied to accelerate G. metallireducens research. For example, growth of G. metallireducens with different electron donors and electron acceptors were studied using the genome-scale metabolic model, which provided a fast and cost-effective way to understand the metabolism of G. metallireducens., Conclusion: We have developed a genome-scale metabolic model for G. metallireducens that features both metabolic similarities and differences to the published model for its close relative, G. sulfurreducens. Together these metabolic models provide an important resource for improving strategies on bioremediation and bioenergy generation.

Geobacter species are among the most effective microorganisms known for the bioremediation of radioactive and toxic metals in contaminated subsurface environments and for converting organic compounds to electricity in microbial fuel cells. However, faster rates of electron transfer could aid in optimizing these processes. Therefore, the Optknock strain design methodology was applied in an iterative manner to the constraint-based, in silico model of Geobacter sulfurreducens to identify gene deletions predicted to increase respiration rates. The common factor in the Optknock predictions was that each resulted in a predicted increase in the cellular ATP demand, either by creating ATP-consuming futile cycles or decreasing the availability of reducing equivalents and inorganic phosphate for ATP biosynthesis. The in silico model predicted that increasing the ATP demand would result in higher fluxes of acetate through the TCA cycle and higher rates of NADPH oxidation coupled with decreases in flux in reactions that funnel acetate toward biosynthetic pathways. A strain of G. sulfurreducens was constructed in which the hydrolytic, F(1) portion of the membrane-bound F(0)F(1) (H(+))-ATP synthase complex was expressed when IPTG was added to the medium. Induction of the ATP drain decreased the ATP content of the cell by more than half. The cells with the ATP drain had higher rates of respiration, slower growth rates, and a lower cell yield. Genome-wide analysis of gene transcript levels indicated that when the higher rate of respiration was induced transcript levels were higher for genes involved in energy metabolism, especially in those encoding TCA cycle enzymes, subunits of the NADH dehydrogenase, and proteins involved in electron acceptor reduction. This was accompanied by lower transcript levels for genes encoding proteins involved in amino acid biosynthesis, cell growth, and motility. Several changes in gene expression that involve processes not included in the in silico model were also detected, including increased expression of a number of redox-active proteins, such as c-type cytochromes and a putative multicopper outer-surface protein. The results demonstrate that it is possible to genetically engineer increased respiration rates in G. sulfurreducens in accordance with predictions from in silico metabolic modeling. To our knowledge, this is the first report of metabolic engineering to increase the respiratory rate of a microorganism.

Geobacteraceae are a family of metal reducing bacteria with important applications in bioremediation and electricity generation. G. sulfurreducens is a representative of Geobacteraceae that has been extensively studied with the goal of extending the understanding of this family of organisms for optimizing their practical applications. Here, we have analyzed gene expression data from 10 experiments involving environmental and genetic perturbations and have identified putative transcription factor binding sites (TFBS) involved in regulating key aspects of metabolism. Specifically, we considered data from both a subset of 10 microarray experiments (7 of 10) and all 10 experiments. The expression data from these two sets were independently clustered, and the upstream regions of genes and operons from the clusters in both sets were used to identify TFBS using the AlignACE program. This analysis resulted in the identification of motifs upstream of several genes involved in central metabolism, sulfate assimilation, and energy metabolism, as well as genes potentially encoding acetate permease. Further, similar TFBS were identified from the analysis of both sets, suggesting that these TFBS are significant in the regulation of metabolism in G. sulfurreducens. In addition, we have utilized microarray data to derive condition specific constraints on the capacity of key enzymes in central metabolism. We have incorporated these constraints into the metabolic model of G. sulfurreducens and simulated Fe(II)-limited growth. The resulting prediction was consistent with data, suggesting that regulatory constraints are important for simulating growth phenotypes in nonoptimal environments.

Previous model-based analysis of the metabolic network of Geobacter sulfurreducens suggested the existence of several redundant pathways. Here, we identified eight sets of redundant pathways that included redundancy for the assimilation of acetate, and for the conversion of pyruvate into acetyl-CoA. These equivalent pathways and two other sub-optimal pathways were studied using 5 single-gene deletion mutants in those pathways for the evaluation of the predictive capacity of the model. The growth phenotypes of these mutants were studied under 12 different conditions of electron donor and acceptor availability. The comparison of the model predictions with the resulting experimental phenotypes indicated that pyruvate ferredoxin oxidoreductase is the only activity able to convert pyruvate into acetyl-CoA. However, the results and the modeling showed that the two acetate activation pathways present are not only active, but needed due to the additional role of the acetyl-CoA transferase in the TCA cycle, probably reflecting the adaptation of these bacteria to acetate utilization. In other cases, the data reconciliation suggested additional capacity constraints that were confirmed with biochemical assays. The results demonstrate the need to experimentally verify the activity of key enzymes when developing in silico models of microbial physiology based on sequence-based reconstruction of metabolic networks.

Geobacter sulfurreducens is a well-studied representative of the Geobacteraceae, which play a critical role in organic matter oxidation coupled to Fe(III) reduction, bioremediation of groundwater contaminated with organics or metals, and electricity production from waste organic matter. In order to investigate G. sulfurreducens central metabolism and electron transport, a metabolic model which integrated genome-based predictions with available genetic and physiological data was developed via the constraint-based modeling approach. Evaluation of the rates of proton production and consumption in the extracellular and cytoplasmic compartments revealed that energy conservation with extracellular electron acceptors, such as Fe(III), was limited relative to that associated with intracellular acceptors. This limitation was attributed to lack of cytoplasmic proton consumption during reduction of extracellular electron acceptors. Model-based analysis of the metabolic cost of producing an extracellular electron shuttle to promote electron transfer to insoluble Fe(III) oxides demonstrated why Geobacter species, which do not produce shuttles, have an energetic advantage over shuttle-producing Fe(III) reducers in subsurface environments. In silico analysis also revealed that the metabolic network of G. sulfurreducens could synthesize amino acids more efficiently than that of Escherichia coli due to the presence of a pyruvate-ferredoxin oxidoreductase, which catalyzes synthesis of pyruvate from acetate and carbon dioxide in a single step. In silico phenotypic analysis of deletion mutants demonstrated the capability of the model to explore the flexibility of G. sulfurreducens central metabolism and correctly predict mutant phenotypes. These results demonstrate that iterative modeling coupled with experimentation can accelerate the understanding of the physiology of poorly studied but environmentally relevant organisms and may help optimize their practical applications.

Highlights: [•] We model microbially-mediated iron reduction at the pore scale. [•] We simulate microbial metabolism using a genome-scale modeling approach. [•] We compare the genome-scale and Monod kinetic rate models. [•] We examine effects of pore-scale mixing on effective reaction rates. [•] Diffusion limitations on pore-scale mixing reduce effective reaction rates. [ABSTRACT FROM AUTHOR]