A new algorithm to find all alternate optimal flux distributions of a metabolic network

Abstract Finding all optimal solutions for a metabolic model is the challenge of metabolic modeling, but there is no practical algorithm for large scale models. A two-phase algorithm is proposed here to systematically identify all optimal solutions. In phase 1, the model is reduced using the FVA approach; in phase 2, all optimal solutions are searched by the addition of a binary variable to convert the model to an MILP problem. The proposed approach proved itself to be a more tractable method for large scale metabolic models when compared with the previously proposed algorithm. The algorithm was implemented on a metabolic model of Escherichia coli (iJR904) to find all optimal flux distributions. The model was reduced from 1076 to 80 fluxes and from 998 to 54 equations and the MILP problem was solved, resulting in 30,744 various flux distributions. For the first time, this number of optimal solutions has been reported.

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

[2]  Nagasuma R. Chandra,et al.  Flux balance analysis of biological systems: applications and challenges , 2009, Briefings Bioinform..

[3]  B. Palsson,et al.  Metabolic Flux Balancing: Basic Concepts, Scientific and Practical Use , 1994, Bio/Technology.

[4]  Hanif D. Sherali,et al.  Linear Programming and Network Flows , 1977 .

[5]  Kenneth J. Kauffman,et al.  Advances in flux balance analysis. , 2003, Current opinion in biotechnology.

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

[7]  B. Palsson,et al.  Metabolic modelling of microbes: the flux-balance approach. , 2002, Environmental microbiology.

[8]  Andrew R. Joyce,et al.  Experimental and Computational Assessment of Conditionally Essential Genes in Escherichia coli , 2006, Journal of bacteriology.

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

[10]  B. Palsson,et al.  Genome-scale in silico models of E. coli have multiple equivalent phenotypic states: assessment of correlated reaction subsets that comprise network states. , 2004, Genome research.

[11]  I. Grossmann,et al.  Recursive MILP model for finding all the alternate optima in LP models for metabolic networks , 2000 .

[12]  Motamedian Ehsan,et al.  FLUX DISTRIBUTION IN BACILLUS SUBTILIS: INSPECTION ON PLURALITY OF OPTIMAL SOLUTIONS , 2011 .

[13]  Manuel Cánovas,et al.  An insight into the role of phosphotransacetylase (pta) and the acetate/acetyl-CoA node in Escherichia coli , 2009, Microbial cell factories.

[14]  Adam M. Feist,et al.  The biomass objective function. , 2010, Current opinion in microbiology.

[15]  F. Llaneras,et al.  Stoichiometric modelling of cell metabolism. , 2008, Journal of bioscience and bioengineering.

[16]  B. Palsson,et al.  A protocol for generating a high-quality genome-scale metabolic reconstruction , 2010 .

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

[18]  Ines Thiele,et al.  Computationally efficient flux variability analysis , 2010, BMC Bioinformatics.