This Article 
 Bibliographic References 
 Add to: 
Manipulating the Steady State of Metabolic Pathways
May/June 2011 (vol. 8 no. 3)
pp. 732-747
Bin Song, University of Florida-Gainesville, FL
İ. Esra Büyüktahtakin, University of Arizona-Tucson, AZ
Sanjay Ranka, University of Florida-Gainesville, FL
Tamer Kahveci, University of Florida-Gainesville, FL
Metabolic pathways show the complex interactions among enzymes that transform chemical compounds. The state of a metabolic pathway can be expressed as a vector, which denotes the yield of the compounds or the flux in that pathway at a given time. The steady state is a state that remains unchanged over time. Altering the state of the metabolism is very important for many applications such as biomedicine, biofuels, food industry, and cosmetics. The goal of the enzymatic target identification problem is to identify the set of enzymes whose knockouts lead the metabolism to a state that is close to a given goal state. Given that the size of the search space is exponential in the number of enzymes, the target identification problem is very computationally intensive. We develop efficient algorithms to solve the enzymatic target identification problem in this paper. Unlike existing algorithms, our method works for a broad set of metabolic network models. We measure the effect of the knockouts of a set of enzymes as a function of the deviation of the steady state of the pathway after their knockouts from the goal state. We develop two algorithms to find the enzyme set with minimal deviation from the goal state. The first one is a traversal approach that explores possible solutions in a systematic way using a branch and bound method. The second one uses genetic algorithms to derive good solutions from a set of alternative solutions iteratively. Unlike the former one, this one can run for very large pathways. Our experiments show that our algorithms' results follow those obtained in vitro in the literature from a number of applications. They also show that the traversal method is a good approximation of the exhaustive search algorithm and it is up to 11 times faster than the exhaustive one. This algorithm runs efficiently for pathways with up to 30 enzymes. For large pathways, our genetic algorithm can find good solutions in less than 10 minutes.

[1] K. Backman, M.J. O'Connor, A. Maruya, E. Rudd, D. McKay, R. Balakrishnan, M. Radjai, V. DiPasquantonio, D. Shoda, and R. Hatch, "Genetic Engineering of Metabolic Pathways Applied to the Production of Phenylalanine," Annals of the New York Academy of Sciences, vol. 589, pp. 16-24, 1990.
[2] H.P.J. Bonarius, G. Schmid, and J. Tramper, "Flux Analysis of Underdetermined Metabolic Networks: The Quest for the Missing Constraints," Trends in Biotechnology, vol. 15, pp. 308-314, 1997.
[3] A.P. Burgard, P. Pharkya, and C.D. Maranas, "Optknock: A Bilevel Programming Framework for Identifying Gene Knockout Strategies for Microbial Strain Optimization," Biotechnology and Bioeng., vol. 84, pp. 647-657, 2003.
[4] G.Q. Chen and Q. Wu, "The Application of Polyhydroxyalkanoates as Tissue Engineering Materials," Biomaterials, vol. 26, pp. 6565-6578, 2005.
[5] H.H. Chen, "Heart Failure: A State of Brain Natriuretic Peptide Deficiency or Resistance or Both," J. Am. College of Cardiology, vol. 49, no. 10, pp. 1089-1091, 2007.
[6] T.C. Chou, "Theoretical Basis, Experimental Design, and Computerized Simulation of Synergism and Antagonism in Drug Combination Studies," Pharmacological Rev., vol. 58, no. 3, pp. 621-681, 2006.
[7] R.K. De, M. Das, and S. Mukhopadhyay, "Incorporation of Enzyme Concentrations into FBA and Identification of Optimal Metabolic Pathways," BMC Systems Biology, vol. 2, no. 65, 2008.
[8] J.S. Edwards and B.O. Palsson, "The Escherichia coli MG1655 in Silico Metabolic Genotype: Its Definition, Characteristics, and Capabilities," Proc. Nat'l Academy of Sciences USA, vol. 97, pp. 5528-5533, 2000.
[9] J. Fernandes, J.M. Saudubray, G.v.d. Berghe, and J.H. Walter, Inborn Metabolic Diseases: Diagnosis and Treatment. Springer, 2006.
[10] R.A. Fisher, "Frequency Distribution of the Values of the Correlation Coefficient in Samples from an Indefinitely Large Population," Biometrika, vol. 10, pp. 507-521, 1915.
[11] J. Forster, I. Famili, P. Fu, B.O. Palsson, and J. Nielsen, "Genome-Scale Reconstruction of the Saccharomyces Cerevisiae Metabolic Network," Genome Research, vol. 13, pp. 244-253, 2003.
[12] S.R. Goldsmith, "Type 5 Phosphodiesterase Inhibition in Heart Failure: The Next Step," J. Am. College of Cardiology, vol. 50, pp. 2145-2147, 2007.
[13] E. Horowitz, S. Sahni, and D. Mehta, "Fundamentals of Data Structures in C++," Silicon Press, 2007.
[14] J.C. BurnettJr., "Modulating cGMP in Heart Failure," BMC Pharmacology, vol. 7, no. S14, 2007.
[15] T. Kahveci, "Np Completeness for Optimal Enzyme Combination Identification," technical report, CISE Dept., Univ. of Florida, Jan. 2008.
[16] M. Kanehisa, "A Database for Post-Genome Analysis," Trends in Genetics, vol. 13, no. 9, pp. 375-376, 1997.
[17] M. Kanehisa and S. Goto, "KEGG: Kyoto Encyclopedia of Genes and Genomes," Nucleic Acids Research, vol. 28, no. 1, pp. 27-30, Jan. 2000.
[18] K.J. Kauffman, P. Prakash, and J.S. Edwards, "Advances in Flux Balance Analysis," Current Opinion in Biotechnology, vol. 14, no. 5, pp. 491-496, 2003.
[19] S. Klamt and E.D. Gilles, "Minimal Cut Sets in Biochemical Reaction Networks," Bioinformatics, vol. 20, no. 2, pp. 226-234, 2004.
[20] S. Loewe, "The Problem of Synergism and Antagonism of Combined Drugs," Arzneimittelforschung, vol. 3, pp. 285-290, 1953.
[21] P. Pacher, A. Nivorozhkin, and C. Szabo, "Therapeutic Effects of Xanthine Oxidase Inhibitors: Renaissance Half a Century after the Discovery of Allopurinol," Pharmacological Rev., vol. 58, no. 1, pp. 87-114, 2006.
[22] B.O. Palsson, Systems Biology: Properties of Reconstructed Networks. Cambridge Univ. Press, 2006.
[23] K.R. Patil, I. Rocha, J. Forster, and J. Nielsen, "Evolutionary Programming as a Platform for in Silico Metabolic Engineering," BMC Bioinformatics, vol. 6, no. 308, 2005.
[24] M. Peschel and W. Mende, The Predator-Prey Model: Do We Live in a Volterra World? Akademie-Verlag, 1986.
[25] S. Philipp, J. Monti, I. Pagel, T. Langenickel, T. Notter, F. Ruschitzka, T. Luscher, R. Dietz, and R. Willenbrock, "Treatment with Darusentan over 21 Days Improved cGMP Generation in Patients with Chronic Heart Failure," Clinical Science, vol. 103, pp. 249-253, 2002.
[26] N.D. Price, J.L. Reed, J.A. Papin, S.L. Wiback, and B.O. Palsson, "Network-Based Analysis of Metabolic Regulation in the Human Red Blood Cell," J. Theoretical Biology, vol. 225, pp. 185-194, 2008.
[27] R. Ramakrishna, J.S. Edwards, A. McCulloch, and B.O. Palsson, "Fluxbalance Analysis of Mitochondrial Energy Metabolism: Consequences of Systemic Stoichiometric Constraints," Am. J. Physiology-Regulatory, Integrative and Comparative Physiology, vol. 280, pp. R695-R704, 2007.
[28] D. Rodriguez-Amaya, "Food Carotenoids: Analysis, Composition and Alterations during Storage and Processing of Foods," Forum of Nutrition, vol. 56, pp. 35-37, 2003.
[29] M.A. Savageau and E.O. Voit, "Recasting Nonlinear Differential Equations as S-Systems: A Canonical Nonlinear Form," Math. Biosciences, vol. 87, pp. 83-115, 1987.
[30] C. Scriver, A.L. Beaudet, D. Valle, W.S. Sly, B. Vagelstein, B. Childs, and K.W. Kinzler, The Online Metabolic and Molecular Bases of Inherited Disease. McGraw-Hill, 2007.
[31] B. Song, P. Sridhar, T. Kahveci, and S. Ranka, "Double Iterative Optimization for Metabolic Network-Based Drug Target Identification," Int'l J. Data Mining and Bioinformatics, vol. 3, pp. 124-144, 2007.
[32] P. Sridhar, T. Kahveci, and S. Ranka, "An Iterative Algorithm for Metabolic Network-Based Drug Target Identification," Proc. Pacific Symp. Biocomputing, 2007.
[33] P. Sridhar, B. Song, T. Kahveci, and S. Ranka, "OPMET: A Metabolic Network-Based Algorithm for Optimal Drug Target Identification," Proc. Pacific Symp. Biocomputing, 2008.
[34] G.J. Vazque, M.J. Pettinari, and B.S. Mendez, "Evidence of an Association between Poly(3-Hydroxybutyrate) Accumulation and Phosphotransbutyrylase Expression in Bacillus Megaterium," Int'l Microbiology, vol. 6, no. 2, pp. 127-129, 2003.
[35] A.I. Vogel, A.R. Tatchell, B.S. Furnis, A.J. Hannaford, and P.W.G. Smith, Vogel's Textbook of Practical Organic Chemistry, fifth ed., Prentice Hall, 1996.
[36] E.O. Voit, Computational Analysis of Biochemical Systems: A Practical Guide for Biochemists and Molecular Biologists. Cambridge Univ. Press, 2000.
[37] E.O. Voit, "Metabolic Modeling: A Tool of Drug Discovery in the Post-Genomic Era," Drug Discovery Today, vol. 7, no. 11, pp. 621-628, May 2002.
[38] Y.T. Yang, G.N. Bennett, and K.Y. San, "Genetic and Metabolic Engineering," Electronic J. Biotechnology, vol. 1, no. 3, pp. 134-141, 1998.

Index Terms:
Metabolic pathway, steady state, traversal approach, genetic algorithm.
Bin Song, İ. Esra Büyüktahtakin, Sanjay Ranka, Tamer Kahveci, "Manipulating the Steady State of Metabolic Pathways," IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol. 8, no. 3, pp. 732-747, May-June 2011, doi:10.1109/TCBB.2010.41
Usage of this product signifies your acceptance of the Terms of Use.