bAccurately predicting the interactions between microbial metabolism and the physical subsurface environment is necessary to enhance subsurface energy development, soil and groundwater cleanup, and carbon management. This study was an initial attempt to confirm the metabolic functional roles within an in silico model using environmental proteomic data collected during field experiments. Shotgun global proteomics data collected during a subsurface biostimulation experiment were used to validate a genome-scale metabolic model of Geobacter metallireducens-specifically, the ability of the metabolic model to predict metal reduction, biomass yield, and growth rate under dynamic field conditions. The constraint-based in silico model of G. metallireducens relates an annotated genome sequence to the physiological functions with 697 reactions controlled by 747 enzymecoding genes. Proteomic analysis showed that 180 of the 637 G. metallireducens proteins detected during the 2008 experiment were associated with specific metabolic reactions in the in silico model. When the field-calibrated Fe(III) terminal electron acceptor process reaction in a reactive transport model for the field experiments was replaced with the genome-scale model, the model predicted that the largest metabolic fluxes through the in silico model reactions generally correspond to the highest abundances of proteins that catalyze those reactions. Central metabolism predicted by the model agrees well with protein abundance profiles inferred from proteomic analysis. Model discrepancies with the proteomic data, such as the relatively low abundances of proteins associated with amino acid transport and metabolism, revealed pathways or flux constraints in the in silico model that could be updated to more accurately predict metabolic processes that occur in the subsurface environment.
Microbial environments have been engineered for enhanced oil recovery (7,23,29,33,42) and remediation of soil and groundwater contaminated by chlorinated hydrocarbons, metals, and radionuclides (4,24,30,40,47,51,63) in the past 30 years and have recently been studied to manage carbon dioxide (46). Mathematical models of microbial processes have been used for experimental interpretation, design, and prediction in recent decades (3,10,20,22,39,48,49,53,56). These processes were very complex to model because there was no way to measure the detailed metabolisms inside the system. They were usually represented by specific terminal electron acceptor process (TEAP) reactions with fixed stoichiometry throughout the simulation and were traditionally modeled by Monod kinetics. These models are sufficient under nonlimiting conditions of nutrients. However, the kinetic parameters of these models, which were calibrated but not directly measurable, cannot reflect the sophisticated mechanisms developed by microorganisms to adapt to the changing environment through the regulation of metabolic pathways, such as their ability to respond to nutrient gradients and environmental stress (5, 6).With the ad...