Numerical Matrices Method for nonlinear system identification and description of dynamics of biochemical reaction networks.

A flexible Numerical Matrices Method (NMM) for nonlinear system identification has been developed based on a description of the dynamics of the system in terms of kinetic complexes. A set of related methods are presented that include increasing amounts of prior information about the reaction network structure, resulting in increased accuracy of the reconstructed rate constants. The NMM is based on an analytical least squares solution for a set of linear equations to determine the rate parameters. In the absence of prior information, all possible unimolecular and bimolecular reactions among the species in the system are considered, and the elements of a general kinetic matrix are determined. Inclusion of prior information is facilitated by formulation of the kinetic matrix in terms of a stoichiometry matrix or a more general set of representation matrices. A method for determination of the stoichiometry matrix beginning only with time-dependent concentration data is presented. In addition, we demonstrate that singularities that arise from linear dependencies among the species can be avoided by inclusion of data collected from a number of different initial states. The NMM provides a flexible set of tools for analysis of complex kinetic data, in particular for analysis of chemical and biochemical reaction networks.

[1]  R. Jackson,et al.  General mass action kinetics , 1972 .

[2]  C. R. Jones,et al.  Determination of Rate Constants for Complex Kinetics Models , 1967 .

[3]  M. Feinberg,et al.  Chemical mechanism structure and the coincidence of the stoichiometric and kinetic subspaces , 1977 .

[4]  Piet Hemker,et al.  Numerical methods for differential equations in system simulation and in parameter estimation : (analysis and simulation of biochemical systems; proc. of the 8th febs meeting, amsterdam, 1972, p 59-80) , 1972 .

[5]  Bernhard O Palsson,et al.  The convex basis of the left null space of the stoichiometric matrix leads to the definition of metabolically meaningful pools. , 2003, Biophysical journal.

[6]  O. Levenspiel.  Chemical Reaction Engineering , 1972 .

[7]  Adam P. Arkin,et al.  Statistical Construction of Chemical Reaction Mechanisms from Measured Time-Series , 1995 .

[8]  Piet Hemker,et al.  Mathematical modelling in blood coagulation ; Simulation and parameter estimation , 1997 .

[9]  P. Traub,et al.  Structure and function of E. coli ribosomes. V. Reconstitution of functionally active 30S ribosomal particles from RNA and proteins. , 1968, Proceedings of the National Academy of Sciences of the United States of America.

[10]  S. Vajda,et al.  Direct integral method via splne-approximation for estimating rate constants , 1982 .

[11]  J. Ross,et al.  A Test Case of Correlation Metric Construction of a Reaction Pathway from Measurements , 1997 .

[12]  Peter P. Valko,et al.  A direct-indirect procedure for estimation of kinetic parameters , 1986 .

[13]  P. McSharry,et al.  Mathematical and computational techniques to deduce complex biochemical reaction mechanisms. , 2004, Progress in biophysics and molecular biology.

[14]  Ossi Renkonen,et al.  Numerical Differentiation and Integration in Chemical Kinetics. , 1963 .

[15]  F. Schluenzen,et al.  Structure of Functionally Activated Small Ribosomal Subunit at 3.3 Å Resolution , 2000, Cell.

[16]  Yonathan Bard,et al.  Comparison of Gradient Methods for the Solution of Nonlinear Parameter Estimation Problems , 1970 .

[17]  J. Ross,et al.  Determination of causal connectivities of species in reaction networks , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[18]  A. V. Karnaukhov,et al.  Application of a New Method of Nonlinear Dynamical System Identification to Biochemical Problems , 2003, Biochemistry (Moscow).

[19]  Gheorghe Maria,et al.  A Review of Algorithms and Trends in Kinetic Model Identification for Chemical and Biochemical Systems , 2004 .

[20]  Carmen G. Moles,et al.  Parameter estimation in biochemical pathways: a comparison of global optimization methods. , 2003, Genome research.

[21]  T. Brubaker,et al.  Nonlinear Parameter Estimation , 1979 .

[22]  G. Siuzdak,et al.  An assembly landscape for the 30S ribosomal subunit , 2005, Nature.

[23]  D. Fell,et al.  A general definition of metabolic pathways useful for systematic organization and analysis of complex metabolic networks , 2000, Nature Biotechnology.

[24]  Benno Hess,et al.  Analysis and simulation of biochemical systems , 1972 .

[25]  C. Vonrhein,et al.  Structure of the 30S ribosomal subunit , 2000, Nature.

[26]  B. Palsson,et al.  Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. , 2003, Genome research.

[27]  D. Fell,et al.  Reaction routes in biochemical reaction systems: Algebraic properties, validated calculation procedure and example from nucleotide metabolism , 2002, Journal of mathematical biology.

[28]  S. Schuster,et al.  Metabolic network structure determines key aspects of functionality and regulation , 2002, Nature.

[29]  Bernhard O Palsson,et al.  Systemic metabolic reactions are obtained by singular value decomposition of genome-scale stoichiometric matrices. , 2003, Journal of theoretical biology.

[30]  W. J. H Stortelder.  Parameter estimation in nonlinear dynamic systems , 1998 .

[31]  L. Hosten,et al.  A comparative study of short cut procedures for parameter estimation in differential equations , 1979 .