Computation of elementary modes: a unifying framework and the new binary approach

BackgroundMetabolic pathway analysis has been recognized as a central approach to the structural analysis of metabolic networks. The concept of elementary (flux) modes provides a rigorous formalism to describe and assess pathways and has proven to be valuable for many applications. However, computing elementary modes is a hard computational task. In recent years we assisted in a multiplication of algorithms dedicated to it. We require a summarizing point of view and a continued improvement of the current methods.ResultsWe show that computing the set of elementary modes is equivalent to computing the set of extreme rays of a convex cone. This standard mathematical representation provides a unified framework that encompasses the most prominent algorithmic methods that compute elementary modes and allows a clear comparison between them. Taking lessons from this benchmark, we here introduce a new method, the binary approach, which computes the elementary modes as binary patterns of participating reactions from which the respective stoichiometric coefficients can be computed in a post-processing step. We implemented the binary approach in FluxAnalyzer 5.1, a software that is free for academics. The binary approach decreases the memory demand up to 96% without loss of speed giving the most efficient method available for computing elementary modes to date.ConclusionsThe equivalence between elementary modes and extreme ray computations offers opportunities for employing tools from polyhedral computation for metabolic pathway analysis. The new binary approach introduced herein was derived from this general theoretical framework and facilitates the computation of elementary modes in considerably larger networks.

[1]  Jason A. Papin,et al.  Extreme pathway lengths and reaction participation in genome-scale metabolic networks. , 2002, Genome research.

[2]  David Avis,et al.  A pivoting algorithm for convex hulls and vertex enumeration of arrangements and polyhedra , 1991, SCG '91.

[3]  Günter M. Ziegler,et al.  Oriented Matroids , 2017, Handbook of Discrete and Computational Geometry, 2nd Ed..

[4]  B. Palsson,et al.  Reconstructing metabolic flux vectors from extreme pathways: defining the alpha-spectrum. , 2003, Journal of theoretical biology.

[5]  Daniel A Beard,et al.  Extreme pathways and Kirchhoff's second law. , 2002, Biophysical journal.

[6]  S. Schuster,et al.  ON ELEMENTARY FLUX MODES IN BIOCHEMICAL REACTION SYSTEMS AT STEADY STATE , 1994 .

[7]  B. Palsson,et al.  Constraints-based models: regulation of gene expression reduces the steady-state solution space. , 2003, Journal of theoretical biology.

[8]  Nagiza F. Samatova,et al.  Parallel out-of-core algorithm for genome-scale enumeration of metabolic systemic pathways , 2002, Proceedings 16th International Parallel and Distributed Processing Symposium.

[9]  Jason A. Papin,et al.  Metabolic pathways in the post-genome era. , 2003, Trends in biochemical sciences.

[10]  Jason A. Papin,et al.  The genome-scale metabolic extreme pathway structure in Haemophilus influenzae shows significant network redundancy. , 2002, Journal of theoretical biology.

[11]  J. Stelling,et al.  Combinatorial Complexity of Pathway Analysis in Metabolic Networks , 2004, Molecular Biology Reports.

[12]  Steffen Klamt,et al.  FluxAnalyzer: exploring structure, pathways, and flux distributions in metabolic networks on interactive flux maps , 2003, Bioinform..

[13]  K. Ülgen,et al.  Metabolic pathway analysis of yeast strengthens the bridge between transcriptomics and metabolic networks , 2004, Biotechnology and bioengineering.

[14]  Joel Franklin,et al.  Methods of Mathematical Economics , 1980 .

[15]  Steffen Klamt,et al.  Minimal cut sets in biochemical reaction networks , 2004, Bioinform..

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

[17]  R. Heinrich,et al.  Metabolic Pathway Analysis: Basic Concepts and Scientific Applications in the Post‐genomic Era , 1999, Biotechnology progress.

[18]  Jason A. Papin,et al.  Comparison of network-based pathway analysis methods. , 2004, Trends in biotechnology.

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

[20]  H. Qian,et al.  Thermodynamic constraints for biochemical networks. , 2004, Journal of theoretical biology.

[21]  Steffen Klamt,et al.  Two approaches for metabolic pathway analysis? , 2003, Trends in biotechnology.

[22]  H. Qian,et al.  Energy balance for analysis of complex metabolic networks. , 2002, Biophysical journal.

[23]  Bernd Sturmfels,et al.  Oriented Matroids: Notation , 1999 .

[24]  D. Fell,et al.  Reaction routes in biochemical reaction systems: Algebraic properties, validated calculation procedure and example from nucleotide metabolism , 2002, Journal of mathematical biology.

[25]  R. Heinrich,et al.  The Regulation of Cellular Systems , 1996, Springer US.

[26]  Komei Fukuda,et al.  Double Description Method Revisited , 1995, Combinatorics and Computer Science.

[27]  B. Sturmfels.  Oriented Matroids , 1993 .

[28]  C. Wagner.  Nullspace Approach to Determine the Elementary Modes of Chemical Reaction Systems , 2004 .

[29]  D. Fell,et al.  Detection of elementary flux modes in biochemical networks: a promising tool for pathway analysis and metabolic engineering. , 1999, Trends in biotechnology.

[30]  Juan Carlos Nuño,et al.  METATOOL: for studying metabolic networks , 1999, Bioinform..

[31]  Stefan Schuster,et al.  Detecting strictly detailed balanced subnetworks in open chemical reaction networks , 1991 .

[32]  Jason A. Papin,et al.  Determination of redundancy and systems properties of the metabolic network of Helicobacter pylori using genome-scale extreme pathway analysis. , 2002, Genome research.

[33]  H. Raiffa,et al.  3. The Double Description Method , 1953 .

[34]  B. Palsson,et al.  Theory for the systemic definition of metabolic pathways and their use in interpreting metabolic function from a pathway-oriented perspective. , 2000, Journal of theoretical biology.

[35]  S. Schuster,et al.  Use of network analysis of metabolic systems in bioengineering , 2002 .

[36]  B. Palsson,et al.  Reconstructing metabolic flux vectors from extreme pathways: defining the α-spectrum , 2003 .

[37]  S. Schuster,et al.  Metabolic network structure determines key aspects of functionality and regulation , 2002, Nature.