Identifying quantitative operation principles in metabolic pathways: a systematic method for searching feasible enzyme activity patterns leading to cellular adaptive responses

BackgroundOptimization methods allow designing changes in a system so that specific goals are attained. These techniques are fundamental for metabolic engineering. However, they are not directly applicable for investigating the evolution of metabolic adaptation to environmental changes. Although biological systems have evolved by natural selection and result in well-adapted systems, we can hardly expect that actual metabolic processes are at the theoretical optimum that could result from an optimization analysis. More likely, natural systems are to be found in a feasible region compatible with global physiological requirements.ResultsWe first present a new method for globally optimizing nonlinear models of metabolic pathways that are based on the Generalized Mass Action (GMA) representation. The optimization task is posed as a nonconvex nonlinear programming (NLP) problem that is solved by an outer-approximation algorithm. This method relies on solving iteratively reduced NLP slave subproblems and mixed-integer linear programming (MILP) master problems that provide valid upper and lower bounds, respectively, on the global solution to the original NLP. The capabilities of this method are illustrated through its application to the anaerobic fermentation pathway in Saccharomyces cerevisiae. We next introduce a method to identify the feasibility parametric regions that allow a system to meet a set of physiological constraints that can be represented in mathematical terms through algebraic equations. This technique is based on applying the outer-approximation based algorithm iteratively over a reduced search space in order to identify regions that contain feasible solutions to the problem and discard others in which no feasible solution exists. As an example, we characterize the feasible enzyme activity changes that are compatible with an appropriate adaptive response of yeast Saccharomyces cerevisiae to heat shockConclusionOur results show the utility of the suggested approach for investigating the evolution of adaptive responses to environmental changes. The proposed method can be used in other important applications such as the evaluation of parameter changes that are compatible with health and disease states.

[1]  Pio A. Aguirre,et al.  Global Optimal Structures of Heat Exchanger Networks by Piecewise Relaxation , 2007 .

[2]  I. Chou,et al.  Recent developments in parameter estimation and structure identification of biochemical and genomic systems. , 2009, Mathematical biosciences.

[3]  Madhukar S. Dasika,et al.  OptCircuit: An optimization based method for computational design of genetic circuits , 2008, BMC Systems Biology.

[4]  Eberhard O Voit,et al.  Theoretical Biology and Medical Modelling Identification of Metabolic System Parameters Using Global Optimization Methods , 2022 .

[5]  Enrique Herrero,et al.  PROTEINS: Structure, Function, and Bioinformatics 57:481–492 (2004) Predictive Reconstruction of the Mitochondrial Iron-Sulfur Cluster Assembly Metabolism. II. Role of Glutaredoxin Grx5 , 2022 .

[6]  Michael A Savageau,et al.  Hysteretic and graded responses in bacterial two‐component signal transduction , 2008, Molecular microbiology.

[7]  A Sorribas,et al.  Analysis and prediction of the effect of uncertain boundary values in modeling a metabolic pathway. , 2000, Biotechnology and bioengineering.

[8]  J. Scott Turner,et al.  The Tinkerer's Accomplice: How Design Emerges from Life Itself , 2007 .

[9]  Eberhard O Voit,et al.  A systems-theoretical framework for health and disease: inflammation and preconditioning from an abstract modeling point of view. , 2009, Mathematical biosciences.

[10]  Lawrence M. Seiford,et al.  Recent developments in dea : the mathematical programming approach to frontier analysis , 1990 .

[11]  K. F. Tipton,et al.  Biochemical systems analysis: A study of function and design in molecular biology , 1978 .

[12]  M. Savageau.  Comparison of classical and autogenous systems of regulation in inducible operons , 1974, Nature.

[13]  Marta Cascante,et al.  Multicriteria optimization of biochemical systems by linear programming: application to production of ethanol by Saccharomyces cerevisiae. , 2003, Biotechnology and bioengineering.

[14]  Alberto Marín-Sanguino,et al.  Optimization of biochemical systems by linear programming and general mass action model representations. , 2003, Mathematical biosciences.

[15]  F. Chaplen,et al.  Developments in metabolic engineering. , 1997, Current opinion in biotechnology.

[16]  Eberhard O Voit,et al.  Theoretical Biology and Medical Modelling , 2022 .

[17]  Michael A. Savageau,et al.  Optimal design of feedback control by inhibition , 1974, Journal of Molecular Evolution.

[18]  Enrique Herrero,et al.  Comprehensive Transcriptional Analysis of the Oxidative Response in Yeast* ♦ , 2008, Journal of Biological Chemistry.

[19]  M A Savageau,et al.  Irreversibility in unbranched pathways: preferred positions based on regulatory considerations. , 2001, Biophysical journal.

[20]  E O Voit,et al.  Assessment of effects of experimental imprecision on optimized biochemical systems. , 2001, Biotechnology and bioengineering.

[21]  Manuel Cánovas,et al.  Modeling, optimization and experimental assessment of continuous L-(-)-carnitine production by Escherichia coli cultures. , 2002, Biotechnology and bioengineering.

[22]  E. Klipp,et al.  Integrative model of the response of yeast to osmotic shock , 2005, Nature Biotechnology.

[23]  M. Savageau.  Biochemical systems analysis. II. The steady-state solutions for an n-pool system using a power-law approximation. , 1969, Journal of theoretical biology.

[24]  Ignacio E. Grossmann,et al.  Logic-based outer approximation for globally optimal synthesis of process networks , 2005, Comput. Chem. Eng..

[25]  Eberhard O. Voit,et al.  Simulation and validation of modelled sphingolipid metabolism in Saccharomyces cerevisiae , 2005, Nature.

[26]  Ignacio E. Grossmann,et al.  An improved piecewise outer-approximation algorithm for the global optimization of MINLP models involving concave and bilinear terms , 2008, Comput. Chem. Eng..

[27]  A. Bennett.  The Origin of Species by means of Natural Selection; or the Preservation of Favoured Races in the Struggle for Life , 1872, Nature.

[28]  Enrique Herrero,et al.  Predictive reconstruction of the mitochondrial iron–sulfur cluster assembly metabolism: I. The role of the protein pair ferredoxin–ferredoxin reductase (Yah1–Arh1) , 2004, Proteins.

[29]  D Kell,et al.  Making cells work--metabolic engineering for everyone. , 1997, Trends in biotechnology.

[30]  E. Voit.  Design principles and operating principles: the yin and yang of optimal functioning. , 2003, Mathematical biosciences.

[31]  A. Sorribas,et al.  Mathematical formalisms based on approximated kinetic representations for modeling genetic and metabolic pathways , 2008, Biotechnology & genetic engineering reviews.

[32]  Michael A Savageau,et al.  Signalling network with a bistable hysteretic switch controls developmental activation of the σF transcription factor in Bacillus subtilis , 2006, Molecular microbiology.

[33]  Eberhard O. Voit,et al.  Computational Analysis of Biochemical Systems: A Practical Guide for Biochemists and Molecular Biologists , 2000 .

[34]  D. Botstein,et al.  Genomic expression programs in the response of yeast cells to environmental changes. , 2000, Molecular biology of the cell.

[35]  R. Raman,et al.  Modelling and computational techniques for logic based integer programming , 1994 .

[36]  Michael A Savageau,et al.  Evidence of selection for low cognate amino acid bias in amino acid biosynthetic enzymes , 2005, Molecular microbiology.

[37]  M. Cascante,et al.  Comparative characterization of the fermentation pathway of Saccharomyces cerevisiae using biochemical systems theory and metabolic control analysis: model definition and nomenclature. , 1995, Mathematical biosciences.

[38]  Christodoulos A. Floudas,et al.  Deterministic global optimization - theory, methods and applications , 2010, Nonconvex optimization and its applications.

[39]  Carolyn R. Bertozzi,et al.  Methods and Applications , 2009 .

[40]  E. O. Voit,et al.  Biochemical systems analysis of genome-wide expression data , 2000, Bioinform..

[41]  C. Darwin,et al.  On the Tendency of Species to form Varieties; and on the Perpetuation of Varieties and Species by Natural Means of Selection , 1858 .

[42]  Ignacio E. Grossmann,et al.  Retrospective on optimization , 2004, Comput. Chem. Eng..

[43]  E. Lander,et al.  Remodeling of yeast genome expression in response to environmental changes. , 2001, Molecular biology of the cell.

[44]  Rui Alves,et al.  Comparative analysis of prototype two‐component systems with either bifunctional or monofunctional sensors: differences in molecular structure and physiological function , 2003, Molecular microbiology.

[45]  Julio R. Banga,et al.  Optimization in computational systems biology , 2008, BMC Systems Biology.

[46]  D. Botstein,et al.  Cluster analysis and display of genome-wide expression patterns. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[47]  Michael A Savageau,et al.  Evolution of enzymes in a series is driven by dissimilar functional demands , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[48]  Eberhard O. Voit,et al.  System estimation from metabolic time-series data , 2008, Bioinform..

[49]  E. Koonin.  Darwinian evolution in the light of genomics , 2008, Nucleic acids research.

[50]  M. Savageau.  Biochemical systems analysis. II. The steady-state solutions for an n-pool system using a power-law approximation. , 1969, Journal of theoretical biology.

[51]  J E Bailey,et al.  Strategies and Challenges in Metabolic Engineering a , 1990, Annals of the New York Academy of Sciences.

[52]  W. S. Hlavacek,et al.  Rules for coupled expression of regulator and effector genes in inducible circuits. , 1996, Journal of molecular biology.

[53]  Eberhard O Voit,et al.  Yield optimization of regulated metabolic systems using deterministic branch‐and‐reduce methods , 2008, Biotechnology and bioengineering.

[54]  Michael A. Savageau,et al.  Significance of autogenously regulated and constitutive synthesis of regulatory proteins in repressible biosynthetic systems , 1975, Nature.

[55]  Michael A Savageau,et al.  Quantitative evolutionary design of glucose 6-phosphate dehydrogenase expression in human erythrocytes , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[56]  M. Savageau.  Biochemical systems analysis. III. Dynamic solutions using a power-law approximation , 1970 .

[57]  E O Voit,et al.  Optimization in integrated biochemical systems , 1992, Biotechnology and bioengineering.

[58]  A Sorribas,et al.  Mathematical models of purine metabolism in man. , 1998, Mathematical biosciences.

[59]  Rui Alves,et al.  Integrating Bioinformatics and Computational Biology: Perspectives and Possibilities for In Silico Network Reconstruction in Molecular Systems Biology , 2008 .

[60]  Rui Alves,et al.  Use of physiological constraints to identify quantitative design principles for gene expression in yeast adaptation to heat shock , 2006, BMC Bioinformatics.

[61]  Julio Vera,et al.  Detection of potential enzyme targets by metabolic modelling and optimization: Application to a simple enzymopathy , 2007, Bioinform..

[62]  W. S. Hlavacek,et al.  Completely uncoupled and perfectly coupled gene expression in repressible systems. , 1997, Journal of molecular biology.