- Poster presentation
- Open Access
Elementary flux modes – state-of-the-art implementation and scope of application
BMC Systems Biologyvolume 1, Article number: P2 (2007)
Elementary flux modes discribe all possible flux distributions in a metabolic network under steady state conditions, acting as a minimal constructive basis of the solution space. The algorithms to compute EFMs originate from computational geometry, where the problem is known as enumeration of extreme rays of polyhedral cones. Due to combinatorial complexity, the algorithms scale poorly and various improvements have been proposed. We present a selection of the most important contributions including a new dynamic adjacent ray enumeration approach. The new approach is based on candidate narrowing, which has recently been proposed by the authors. However, the present variant is new and optimizes the method significantly, resulting to the best of our knowledge in the most efficient algorithm for elementary flux mode computation known today.
As a second aspect, we focus on the application of elementary modes, placing it into the context of constraint-based modeling approaches. Therefore, we compare pathway analysis (aggregating EFMs and related methods) with flux balance analysis (FBA) and minimization of metabolic adjustment (MoMA). A critical evaluation of the methods assesses properties such as computational costs, network size limitations, applicability and practical benefit.
Computing the flux cone
Thermodynamic constraints (irreversible reactions) and the steady state assumption (constant concentration of internal metabolites) define the flux cone, the solution space for possible flux values of the metabolic network under the given constraints. Elementary modes are the extreme rays (edges) of the flux cone. Every valid flux vector is a positive linear combination of elementary modes.
Pairs of adjacent elementary modes are used to generate new modes. Finding adjacent modes is one of the most crucial aspects of the algorithm.
With pattern trees, adjacent modes can be enumerated efficiently, demanding memory that is linear in the number of intermediary modes. A multithreaded version of the candidate narrowing procedure is presented, which allows exploiting modern Dual-Core processors.
See Figure 2.
Marco Terzer, Jörg Stelling: Accelerating the computation of elementary modes using pattern trees. WABI, volume 4175 of Lecture Notes in Computer Science. 2006, 333-343. Springer
Klamt S, Gagneur J, von Kamp A: Algorithmic approaches for computing elementary modes in large biochemical reaction networks. IEE Proc Systems Biol. 152: 249-55. 10.1049/ip-syb:20050035.