Search

Genome-scale metabolic models comprehensively describe an organisms metabolism and can be tailored using omics data to model condition-specific physiology. The quality of context-specific models is impacted by (i) choice of algorithm and parameters and (ii) alternate context-specific models that equally explain the -omics data. Here we quantify the influence of alternate optima on microbial and mammalian model extraction using GIMME, iMAT, MBA, and mCADRE. We find that metabolic tasks defining an organisms phenotype must be explicitly and quantitatively protected. The scope of alternate models is strongly influenced by algorithm choice and the topological properties of the parent genome-scale model with fatty acid metabolism and intracellular metabolite transport contributing much to alternate solutions in all models. mCADRE extracted the most reproducible context-specific models and models generated using MBA had the most alternate solutions. There were fewer qualitatively different solutions generated by GIMME in E. coli, but these increased substantially in the mammalian models. Screening ensembles using a receiver operating characteristic plot identified the best-performing models. A comprehensive evaluation of models extracted using combinations of extraction methods and expression thresholds revealed that GIMME generated the best-performing models in E. coli, whereas mCADRE is better suited for complex mammalian models. These findings suggest guidelines for benchmarking -omics integration algorithms and motivate the development of a systematic workflow to enumerate alternate models and extract biologically relevant context-specific models.

Kinetic models of metabolic networks offer the promise of quantitative phenotype prediction. The mechanistic characterization of enzyme catalyzed reactions allows for tracing the effect of perturbations in metabolite concentrations and reaction fluxes in response to genetic and environmental perturbation that are beyond the scope of stoichiometric models. In this study, we develop a two-step computational pipeline for the rapid parameterization of kinetic models of metabolic networks using a curated metabolic model and available 13C-labeling distributions under multiple genetic and environmental perturbations. The first step involves the elucidation of all intracellular fluxes in a core model of E. coli containing 74 reactions and 61 metabolites using 13C-Metabolic Flux Analysis (13C-MFA). Here, fluxes corresponding to the mid-exponential growth phase are elucidated for seven single gene deletion mutants from upper glycolysis, pentose phosphate pathway and the Entner-Doudoroff pathway. The computed flux ranges are then used to parameterize the same (i.e., k-ecoli74) core kinetic model for E. coli with 55 substrate-level regulations using the newly developed K-FIT parameterization algorithm. The K-FIT algorithm employs a combination of equation decomposition and iterative solution techniques to evaluate steady-state fluxes in response to genetic perturbations. k-ecoli74 predicted 86% of flux values for strains used during fitting within a single standard deviation of 13C-MFA estimated values. By performing both tasks using the same network, errors associated with lack of congruity between the two networks are avoided, allowing for seamless integration of data with model building. Product yield predictions and comparison with previously developed kinetic models indicate shifts in flux ranges and the presence or absence of mutant strains delivering flux towards pathways of interest from training data significantly impact predictive capabilities. Using this workflow, the impact of completeness of fluxomic datasets and the importance of specific genetic perturbations on uncertainties in kinetic parameter estimation are evaluated. [ABSTRACT FROM AUTHOR]

Synechococcus elongatus UTEX 2973 (Synechococcus 2973) has the shortest reported doubling time (2.1 h) among cyanobacteria, making it a promising platform for the solar-based production of biochemicals. In this meta-analysis, its intracellular flux distribution was recomputed using genome-scale isotopic nonstationary 13C-metabolic flux analysis given the labeling dynamics of 13 metabolites reported in an earlier study. To achieve this, a genome-scale mapping model, namely imSyu593, was constructed using the imSyn617 mapping model for Synechocystis sp. PCC 6803 (Synechocystis 6803) as the starting point encompassing 593 reactions. The flux elucidation revealed nearly complete conversion (greater than 96%) of the assimilated carbon into biomass in Synechococcus 2973. In contrast, Synechocystis 6803 achieves complete conversion of only 86% of the assimilated carbon. This high biomass yield was enabled by the reincorporation of the fixed carbons lost in anabolic and photorespiratory pathways in conjunction with flux rerouting through a nondecarboxylating reaction such as phosphoketolase. This reincorporation of lost CO2 sustains a higher flux through the photorespiratory C2 cycle that fully meets the glycine and serine demands for growth. In accordance with the high carbon efficiency drive, acetyl-coenzyme A was entirely produced using the carbon-efficient phosphoketolase pathway. Comparison of the Synechococcus 2973 flux map with that of Synechocystis 6803 revealed differences in the use of Calvin cycle and photorespiratory pathway reactions. The two species used different reactions for the synthesis of metabolites such as fructose-6-phosphate, glycine, sedoheptulose-7-phosphate, and Ser. These findings allude to a highly carbon-efficient metabolism alongside the fast carbon uptake rate in Synechococcus 2973, which explains its faster growth rate. [ABSTRACT FROM AUTHOR]