Direct calculation of elementary flux modes satisfying several biological constraints in genome-scale metabolic networks

MOTIVATION The concept of Elementary Flux Mode (EFM) has been widely used for the past 20 years. However, its application to genome-scale metabolic networks (GSMNs) is still under development because of methodological limitations. Therefore, novel approaches are demanded to extend the application of EFMs. A novel family of methods based on optimization is emerging that provides us with a subset of EFMs. Because the calculation of the whole set of EFMs goes beyond our capacity, performing a selective search is a proper strategy. RESULTS Here, we present a novel mathematical approach calculating EFMs fulfilling additional linear constraints. We validated our approach based on two metabolic networks in which all the EFMs can be obtained. Finally, we analyzed the performance of our methodology in the GSMN of the yeast Saccharomyces cerevisiae by calculating EFMs producing ethanol with a given minimum carbon yield. Overall, this new approach opens new avenues for the calculation of EFMs in GSMNs. AVAILABILITY AND IMPLEMENTATION Matlab code is provided in the supplementary online materials CONTACT fplanes@ceit.es. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

[1]  Francisco J. Planes,et al.  A critical examination of stoichiometric and path-finding approaches to metabolic pathways , 2008, Briefings Bioinform..

[2]  Susumu Goto,et al.  KEGG for integration and interpretation of large-scale molecular data sets , 2011, Nucleic Acids Res..

[3]  W. M. Ingledew,et al.  Effects of particulate materials and osmoprotectants on very-high-gravity ethanolic fermentation by Saccharomyces cerevisiae , 1994, Applied and environmental microbiology.

[4]  B. Palsson,et al.  Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. , 2003, Genome research.

[5]  Ping Ji,et al.  Decomposing flux distributions into elementary flux modes in genome-scale metabolic networks , 2011, Bioinform..

[6]  Angel Rubio,et al.  Computing the shortest elementary flux modes in genome-scale metabolic networks , 2009, Bioinform..

[7]  Peter D. Karp,et al.  EcoCyc: a comprehensive database of Escherichia coli biology , 2010, Nucleic Acids Res..

[8]  Christoph Kaleta,et al.  Computing Elementary Flux Modes in Genome-scale Metabolic Networks , 2009, GCB.

[9]  Jens Nielsen,et al.  Economic and environmental impacts of microbial biodiesel , 2013, Nature Biotechnology.

[10]  Bernhard O. Palsson,et al.  BiGG: a Biochemical Genetic and Genomic knowledgebase of large scale metabolic reconstructions , 2010, BMC Bioinformatics.

[11]  Leen Stougie,et al.  A note on the complexity of finding and enumerating elementary modes , 2010, Biosyst..

[12]  Adam M. Feist,et al.  A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information , 2007, Molecular systems biology.

[13]  Christoph Kaleta,et al.  In Silico Evidence for Gluconeogenesis from Fatty Acids in Humans , 2011, PLoS Comput. Biol..

[14]  Francisco J. Planes,et al.  Exploring metabolic pathways in genome-scale networks via generating flux modes , 2011, Bioinform..

[15]  Jörg Stelling,et al.  Large-scale computation of elementary flux modes with bit pattern trees , 2008, Bioinform..

[16]  Pedro Mendes,et al.  Yeast 5 – an expanded reconstruction of the Saccharomyces cerevisiae metabolic network , 2012, BMC Systems Biology.

[17]  A. Lehninger.  Principles of Biochemistry , 1984 .

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

[19]  Eugénio C. Ferreira,et al.  Random sampling of elementary flux modes in large-scale metabolic networks , 2012, Bioinform..

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

[21]  H. Kitano.  Systems Biology: A Brief Overview , 2002, Science.

[22]  Stefan Schuster,et al.  Systems biology Metatool 5.0: fast and flexible elementary modes analysis , 2006 .

[23]  Alexander Bockmayr,et al.  A new constraint-based description of the steady-state flux cone of metabolic networks , 2009, Discret. Appl. Math..

[24]  lexander,et al.  THE GENERALIZED SIMPLEX METHOD FOR MINIMIZING A LINEAR FORM UNDER LINEAR INEQUALITY RESTRAINTS , 2012 .

[25]  Angel Rubio,et al.  Advances in network-based metabolic pathway analysis and gene expression data integration , 2015, Briefings Bioinform..

[26]  L. Vogler.  Lehninger Principles Of Biochemistry 5th Edition , 2016 .

[27]  Angel Rubio,et al.  Selection of human tissue-specific elementary flux modes using gene expression data , 2013, Bioinform..

[28]  Desmond S. Lun,et al.  Analysis of complex metabolic behavior through pathway decomposition , 2011, BMC Systems Biology.

[29]  D. Fell,et al.  A general definition of metabolic pathways useful for systematic organization and analysis of complex metabolic networks , 2000, Nature Biotechnology.