Co-evolution of strain design methods based on flux balance and elementary mode analysis

More than a decade ago, the first genome-scale metabolic models for two of the most relevant microbes for biotechnology applications, Escherichia coli and Saccaromyces cerevisiae, were published. Shortly after followed the publication of OptKnock, the first strain design method using bilevel optimization to couple cellular growth with the production of a target product. This initiated the development of a family of strain design methods based on the concept of flux balance analysis. Another family of strain design methods, based on the concept of elementary mode analysis, has also been growing. Although the computation of elementary modes is hindered by computational complexity, recent breakthroughs have allowed applying elementary mode analysis at the genome scale. Here we review and compare strain design methods and look back at the last ten years of in silico strain design with constraint-based models. We highlight some features of the different approaches and discuss the utilization of these methods in successful in vivo metabolic engineering applications.Novo Nordisk UK Research Foundation(NORTE-07-0124-FEDER-000028

Genome-scale constraint-based modeling of Geobacter metallireducens

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\u27 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

Constraint-based modeling analysis of the metabolism of two Pelobacter species

BACKGROUND: Pelobacter species are commonly found in a number of subsurface environments, and are unique members of the Geobacteraceae family. They are phylogenetically intertwined with both Geobacter and Desulfuromonas species. Pelobacter species likely play important roles in the fermentative degradation of unusual organic matters and syntrophic metabolism in the natural environments, and are of interest for applications in bioremediation and microbial fuel cells. RESULTS: In order to better understand the physiology of Pelobacter species, genome-scale metabolic models for Pelobacter carbinolicus and Pelobacter propionicus were developed. Model development was greatly aided by the availability of models of the closely related Geobacter sulfurreducens and G. metallireducens. The reconstructed P. carbinolicus model contains 741 genes and 708 reactions, whereas the reconstructed P. propionicus model contains 661 genes and 650 reactions. A total of 470 reactions are shared among the two Pelobacter models and the two Geobacter models. The different reactions between the Pelobacter and Geobacter models reflect some unique metabolic capabilities such as fermentative growth for both Pelobacter species. The reconstructed Pelobacter models were validated by simulating published growth conditions including fermentations, hydrogen production in syntrophic co-culture conditions, hydrogen utilization, and Fe(III) reduction. Simulation results matched well with experimental data and indicated the accuracy of the models. CONCLUSIONS: We have developed genome-scale metabolic models of P. carbinolicus and P. propionicus. These models of Pelobacter metabolism can now be incorporated into the growing repertoire of genome scale models of the Geobacteraceae family to aid in describing the growth and activity of these organisms in anoxic environments and in the study of their roles and interactions in the subsurface microbial community

Role of Geobacter sulfurreducens Outer Surface c-Type Cytochromes in Reduction of Soil Humic Acid and Anthraquinone-2,6-Disulfonate▿

Deleting individual genes for outer surface c-type cytochromes in Geobacter sulfurreducens partially inhibited the reduction of humic substances and anthraquinone-2,6,-disulfonate. Complete inhibition was obtained only when five of these genes were simultaneously deleted, suggesting that diverse outer surface cytochromes can contribute to the reduction of humic substances and other extracellular quinones

Novel strategy for three-dimensional real-time imaging of microbial fuel cell communities: monitoring the inhibitory effects of proton accumulation within the anode biofilm

Harvesting electricity from the environment, organic wastes, or renewable biomass with microbial fuel cells (MFCs) is an appealing strategy, but the destructive sampling required to investigate the anode-associated biofilms has hampered research designed to better understand and optimize microbe–anode interactions. Therefore, a MFC that permits real-time imaging of the anode biofilm with confocal scanning laser microscopy was developed. In this new MFC Geobacter sulfurreducens, an organism closely related to those often found on MFC anodes and capable of high current densities, produced current comparable to that previously reported with other MFC designs. G. sulfurreducens engineered to produce the fluorescent protein mcherry to facilitate real-time imaging produced current comparable to wild-type cells. Introducing C-SNARF-4, a pH-sensitive fluoroprobe, into the anode chamber revealed strong pH gradients within the anode biofilms. The pH decreased with increased proximity to the anode surface and from the exterior to the interior of biofilm pillars. Near the anode surface pH levels were as low as 6.1 compared to ca. 7 in the external medium. Various controls demonstrated that the proton accumulation was associated with current production. Dropping the pH of culture medium from 7 to 6 severely limited the growth of G. sulfurreducens. These results demonstrate that it is feasible to non-destructively monitor the activity of anode biofilms in real time and suggest that the accumulation of protons that are released from organic matter oxidation within anode biofilms can limit current production

Genome-Scale Dynamics 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

Purification and Characterization of OmcZ, an Outer-Surface, Octaheme c-Type Cytochrome Essential for Optimal Current Production by Geobacter sulfurreducens

Previous studies have demonstrated that Geobacter sulfurreducens requires the c-type cytochrome OmcZ, which is present in large (OmcZL; 50-kDa) and small (OmcZS; 30-kDa) forms, for optimal current production in microbial fuel cells. This protein was further characterized to aid in understanding its role in current production. Subcellular-localization studies suggested that OmcZS was the predominant extracellular form of OmcZ. N- and C-terminal amino acid sequence analysis of purified OmcZS and molecular weight measurements indicated that OmcZS is a cleaved product of OmcZL retaining all 8 hemes, including 1 heme with the unusual c-type heme-binding motif CX14CH. The purified OmcZS was remarkably thermally stable (thermal-denaturing temperature, 94.2°C). Redox titration analysis revealed that the midpoint reduction potential of OmcZS is approximately −220 mV (versus the standard hydrogen electrode [SHE]) with nonequivalent heme groups that cover a large reduction potential range (−420 to −60 mV). OmcZS transferred electrons in vitro to a diversity of potential extracellular electron acceptors, such as Fe(III) citrate, U(VI), Cr(VI), Au(III), Mn(IV) oxide, and the humic substance analogue anthraquinone-2,6-disulfonate, but not Fe(III) oxide. The biochemical properties and extracellular localization of OmcZ suggest that it is well suited for promoting electron transfer in current-producing biofilms of G. sulfurreducens