Systemic properties of ensembles of metabolic networks: application of graphical and statistical methods to simple unbranched pathways

MOTIVATION Mathematical models are the only realistic method for representing the integrated dynamic behavior of complex biochemical networks. However, it is difficult to obtain a consistent set of values for the parameters that characterize such a model. Even when a set of parameter values exists, the accuracy of the individual values is questionable. Therefore, we were motivated to explore statistical techniques for analyzing the properties of a given model when knowledge of the actual parameter values is lacking. RESULTS The graphical and statistical methods presented in the previous paper are applied here to simple unbranched biosynthetic pathways subject to control by feedback inhibition. We represent these pathways within a canonical nonlinear formalism that provides a regular structure that is convenient for randomly sampling the parameter space. After constructing a large ensemble of randomly generated sets of parameter values, the structural and behavioral properties of the model with these parameter sets are examined statistically and classified. The results of our analysis demonstrate that certain properties of these systems are strongly correlated, thereby revealing aspects of organization that are highly probable independent of selection. Finally, we show how specification of a given behavior affects the distribution of acceptable parameter values.

[1]  R. Somogyi,et al.  The gene expression matrix: towards the extraction of genetic network architectures , 1997 .

[2]  Rodney L. Levine,et al.  Modulation by molecular interactions , 1985 .

[3]  A Sorribas,et al.  Metabolic pathway characterization from transient response data obtained in situ: parameter estimation in S-system models. , 1993, Journal of theoretical biology.

[4]  L. Glass.  Classification of biological networks by their qualitative dynamics. , 1975, Journal of theoretical biology.

[5]  A. J,et al.  Cell Behaviour as a Dynamic Attractor in the Intracellular Signalling System , 1999 .

[6]  M A Savageau,et al.  Concepts relating the behavior of biochemical systems to their underlying molecular properties. , 1971, Archives of biochemistry and biophysics.

[7]  P. Brown,et al.  Exploring the metabolic and genetic control of gene expression on a genomic scale. , 1997, Science.

[8]  Paul Bratley,et al.  Algorithm 659: Implementing Sobol's quasirandom sequence generator , 1988, TOMS.

[9]  P. Törönen,et al.  Analysis of gene expression data using self‐organizing maps , 1999, FEBS letters.

[10]  J. Clegg,et al.  Properties and metabolism of the aqueous cytoplasm and its boundaries. , 1984, The American journal of physiology.

[11]  K. Abe,et al.  Design of practically stable n-dimensional feedback systems : a state-space approach , 1996 .

[12]  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.

[13]  O. J. Dunn,et al.  Applied statistics: analysis of variance and regression , 1975 .

[14]  M A Savageau,et al.  Accuracy of alternative representations for integrated biochemical systems. , 1987, Biochemistry.

[15]  D. Botstein,et al.  The transcriptional program of sporulation in budding yeast. , 1998, Science.

[16]  P. Srere,et al.  Binding of citrate synthase and malate dehydrogenase to mitochondrial inner membranes: tissue distribution and metabolite effects. , 1984, Biochemical and biophysical research communications.

[17]  A. Huitson,et al.  Applied Statistics: Analysis of Variance and Regression , 1976 .

[18]  P. Srere,et al.  Complexes of sequential metabolic enzymes. , 1987, Annual review of biochemistry.

[19]  S. Kauffman.  Homeostasis and Differentiation in Random Genetic Control Networks , 1969, Nature.

[20]  Stuart A. Kauffman,et al.  The origins of order , 1993 .

[21]  M A Savageau,et al.  Subunit structure of regulator proteins influences the design of gene circuitry: analysis of perfectly coupled and completely uncoupled circuits. , 1995, Journal of molecular biology.

[22]  C. Yanofsky,et al.  Role of regulatory features of the trp operon of Escherichia coli in mediating a response to a nutritional shift , 1994, Journal of bacteriology.

[23]  H. S. Marinho,et al.  Lipid peroxidation in mitochondrial inner membranes. I. An integrative kinetic model. , 1996, Free radical biology & medicine.

[24]  M. Savageau.  Biochemical Systems Analysis: A Study of Function and Design in Molecular Biology , 1976 .

[25]  Ian R. Petersen,et al.  Robust observability for a class of time-varying discrete-time uncertain systems , 1996 .

[26]  D Feng,et al.  Cut set analysis of compartmental models with applications to experiment design. , 1991, The American journal of physiology.

[27]  K R Godfrey,et al.  Structural identifiability of the parameters of a nonlinear batch reactor model. , 1992, Mathematical biosciences.

[28]  J. Diamond,et al.  How do biological systems discriminate among physically similar ions? , 1975, The Journal of experimental zoology.

[29]  R. D. Bliss,et al.  Role of feedback inhibition in stabilizing the classical operon. , 1982, Journal of theoretical biology.

[30]  J Ovádi,et al.  Metabolic consequences of enzyme interactions , 1996, Cell biochemistry and function.

[31]  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.

[32]  Daniel E. Atkinson,et al.  Limitation of Metabolite Concentrations and the Conservation of Solvent Capacity in the Living Cell , 1969 .

[33]  H S Colburn,et al.  Lateral position and interaural discrimination. , 1977, The Journal of the Acoustical Society of America.

[34]  M A Savageau,et al.  Model assessment and refinement using strategies from biochemical systems theory: application to metabolism in human red blood cells. , 1996, Journal of theoretical biology.

[35]  A Sorribas,et al.  Strategies for representing metabolic pathways within biochemical systems theory: reversible pathways. , 1989, Mathematical biosciences.

[36]  M A Savageau,et al.  The tricarboxylic acid cycle in Dictyostelium discoideum. Systemic effects of including protein turnover in the current model. , 1993, The Journal of biological chemistry.

[37]  Rui Alves,et al.  Comparing systemic properties of ensembles of biological networks by graphical and statistical methods , 2000, Bioinform..

[38]  P A Quant.  Experimental application of top-down control analysis to metabolic systems. , 1993, Trends in biochemical sciences.

[39]  Timothy R. Ginn,et al.  A continuous‐time inverse operator for groundwater and contaminant transport modeling: Model identifiability , 1992 .

[40]  M A Savageau.  Optimal design of feedback control by inhibition: dynamic considerations. , 1975, Journal of molecular evolution.

[41]  S. Liang,et al.  Median attractor and transients in random Boolean nets , 1996 .

[42]  B E Wright,et al.  Systems analysis of the tricarboxylic acid cycle in Dictyostelium discoideum. II. Control analysis. , 1992, The Journal of biological chemistry.

[43]  G. Brown,et al.  A 'top-down' approach to the determination of control coefficients in metabolic control theory. , 1990, European journal of biochemistry.

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

[45]  Eberhard O. Voit,et al.  Power-Law Approach to Modeling Biological Systems : I. Theory , 1982 .

[46]  A. Salvador,et al.  Synergism analysis of biochemical systems. I. Conceptual framework. , 2000, Mathematical biosciences.

[47]  Michael L. Shuler,et al.  Generalized differential specific rate equation for microbial growth , 1979 .