MapMaker and PathTracer for tracking carbon in genome-scale metabolic models.

Constraint-based reconstruction and analysis (COBRA) modeling results can be difficult to interpret given the large numbers of reactions in genome-scale models. While paths in metabolic networks can be found, existing methods are not easily combined with constraint-based approaches. To address this limitation, two tools (MapMaker and PathTracer) were developed to find paths (including cycles) between metabolites, where each step transfers carbon from reactant to product. MapMaker predicts carbon transfer maps (CTMs) between metabolites using only information on molecular formulae and reaction stoichiometry, effectively determining which reactants and products share carbon atoms. MapMaker correctly assigned CTMs for over 97% of the 2,251 reactions in an Escherichia coli metabolic model (iJO1366). Using CTMs as inputs, PathTracer finds paths between two metabolites. PathTracer was applied to iJO1366 to investigate the importance of using CTMs and COBRA constraints when enumerating paths, to find active and high flux paths in flux balance analysis (FBA) solutions, to identify paths for putrescine utilization, and to elucidate a potential CO2 fixation pathway in E. coli. These results illustrate how MapMaker and PathTracer can be used in combination with constraint-based models to identify feasible, active, and high flux paths between metabolites.

[1]  R. Mahadevan,et al.  The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. , 2003, Metabolic engineering.

[2]  B. Palsson,et al.  An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR) , 2003, Genome Biology.

[3]  Oliver Kohlbacher,et al.  MetaRoute: fast search for relevant metabolic routes for interactive network navigation and visualization , 2008, Bioinform..

[4]  J. van Helden,et al.  Metabolic pathfinding using RPAIR annotation. , 2009, Journal of molecular biology.

[5]  B. Palsson,et al.  Systems approach to refining genome annotation , 2006, Proceedings of the National Academy of Sciences.

[6]  Masanori Arita.  The metabolic world of Escherichia coli is not small. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Joshua J. Hamilton,et al.  Quantitative assessment of thermodynamic constraints on the solution space of genome-scale metabolic models. , 2013, Biophysical journal.

[8]  Jennifer L Reed,et al.  FOCAL: an experimental design tool for systematizing metabolic discoveries and model development , 2012, Genome Biology.

[9]  B. Palsson,et al.  Constraining the metabolic genotype–phenotype relationship using a phylogeny of in silico methods , 2012, Nature Reviews Microbiology.

[10]  Yves Deville,et al.  NeAT: a toolbox for the analysis of biological networks, clusters, classes and pathways , 2008, Nucleic Acids Res..

[11]  Francisco J. Planes,et al.  Path finding methods accounting for stoichiometry in metabolic networks , 2011, Genome Biology.

[12]  Peter D. Karp,et al.  Pathway Tools version 13.0: integrated software for pathway/genome informatics and systems biology , 2015, Briefings Bioinform..

[13]  Tom M. Conrad,et al.  Omic data from evolved E. coli are consistent with computed optimal growth from genome-scale models , 2010, Molecular systems biology.

[14]  D. Fell,et al.  The small world inside large metabolic networks , 2000, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[15]  Ali R. Zomorrodi,et al.  Mathematical optimization applications in metabolic networks. , 2012, Metabolic engineering.

[16]  Vinay Satish Kumar,et al.  GrowMatch: An Automated Method for Reconciling In Silico/In Vivo Growth Predictions , 2009, PLoS Comput. Biol..

[17]  Emmanuel Barillot,et al.  BiNoM 2.0, a Cytoscape plugin for accessing and analyzing pathways using standard systems biology formats , 2013, BMC Systems Biology.

[18]  Costas D Maranas,et al.  OptStrain: a computational framework for redesign of microbial production systems. , 2004, Genome research.

[19]  Francisco J. Planes,et al.  Refining carbon flux paths using atomic trace data , 2014, Bioinform..

[20]  C. Schilling,et al.  Flux coupling analysis of genome-scale metabolic network reconstructions. , 2004, Genome research.

[21]  Jeffrey D Orth,et al.  What is flux balance analysis? , 2010, Nature Biotechnology.

[22]  S. Wodak,et al.  Inferring meaningful pathways in weighted metabolic networks. , 2006, Journal of molecular biology.

[23]  Janet M Thornton,et al.  Analysis of metabolic networks using a pathway distance metric through linear programming. , 2003, Metabolic engineering.

[24]  Robert Giegerich,et al.  PathFinder: reconstruction and dynamic visualization of metabolic pathways , 2002, Bioinform..

[25]  Daniel Machado,et al.  Systematic Evaluation of Methods for Integration of Transcriptomic Data into Constraint-Based Models of Metabolism , 2014, PLoS Comput. Biol..

[26]  Michel Dumontier,et al.  CO: A chemical ontology for identification of functional groups and semantic comparison of small molecules , 2005, FEBS letters.

[27]  Costas D Maranas,et al.  Microbial 1-butanol production: Identification of non-native production routes and in silico engineering interventions. , 2010, Biotechnology journal.

[28]  Adam M. Feist,et al.  A comprehensive genome-scale reconstruction of Escherichia coli metabolism—2011 , 2011, Molecular systems biology.

[29]  Dipak Barua,et al.  An Automated Phenotype-Driven Approach (GeneForce) for Refining Metabolic and Regulatory Models , 2010, PLoS Comput. Biol..

[30]  R. Albert,et al.  The large-scale organization of metabolic networks , 2000, Nature.

[31]  G. Fuchs.  Alternative pathways of carbon dioxide fixation: insights into the early evolution of life? , 2011, Annual review of microbiology.

[32]  J. Keasling,et al.  Catabolite Repression of the Propionate Catabolic Genes in Escherichia coli and Salmonella enterica: Evidence for Involvement of the Cyclic AMP Receptor Protein , 2005, Journal of bacteriology.