Simplification of biochemical models: a general approach based on the analysis of the impact of individual species and reactions on the systems dynamics

BackgroundGiven the complex mechanisms underlying biochemical processes systems biology researchers tend to build ever increasing computational models. However, dealing with complex systems entails a variety of problems, e.g. difficult intuitive understanding, variety of time scales or non-identifiable parameters. Therefore, methods are needed that, at least semi-automatically, help to elucidate how the complexity of a model can be reduced such that important behavior is maintained and the predictive capacity of the model is increased. The results should be easily accessible and interpretable. In the best case such methods may also provide insight into fundamental biochemical mechanisms.ResultsWe have developed a strategy based on the Computational Singular Perturbation (CSP) method which can be used to perform a "biochemically-driven" model reduction of even large and complex kinetic ODE systems. We provide an implementation of the original CSP algorithm in COPASI (a COmplex PAthway SImulator) and applied the strategy to two example models of different degree of complexity - a simple one-enzyme system and a full-scale model of yeast glycolysis.ConclusionThe results show the usefulness of the method for model simplification purposes as well as for analyzing fundamental biochemical mechanisms. COPASI is freely available at http://www.copasi.org.

[1]  Ursula Kummer,et al.  A New Time-Dependent Complexity Reduction Method for Biochemical Systems , 2005, Trans. Comp. Sys. Biology.

[2]  Jacky L. Snoep,et al.  Web-based kinetic modelling using JWS Online , 2004, Bioinform..

[3]  Ulrich Maas,et al.  Simplifying chemical kinetics: Intrinsic low-dimensional manifolds in composition space , 1992 .

[4]  J. Angilella,et al.  Dust trapping in vortex pairs , 2010, 1003.3842.

[5]  Gunnar Cedersund,et al.  Reduction of a biochemical model with preservation of its basic dynamic properties , 2006, The FEBS journal.

[6]  Irina Surovtsova,et al.  Accessible methods for the dynamic time-scale decomposition of biochemical systems , 2009, Bioinform..

[7]  Thomas Pfeiffer,et al.  Exploring the pathway structure of metabolism: decomposition into subnetworks and application to Mycoplasma pneumoniae , 2002, Bioinform..

[8]  Robert Roskoski Michaelis-Menten Kinetics , 2007 .

[9]  B O Palsson,et al.  Mathematical modelling of dynamics and control in metabolic networks. I. On Michaelis-Menten kinetics. , 1984, Journal of theoretical biology.

[10]  Jason A. Papin,et al.  Analysis of metabolic capabilities using singular value decomposition of extreme pathway matrices. , 2003, Biophysical journal.

[11]  F. Hynne,et al.  Full-scale model of glycolysis in Saccharomyces cerevisiae. , 2001, Biophysical chemistry.

[12]  Neema Jamshidi,et al.  Description and analysis of metabolic connectivity and dynamics in the human red blood cell. , 2002, Biophysical journal.

[13]  Bernhard O. Palsson,et al.  Top-Down Analysis of Temporal Hierarchy in Biochemical Reaction Networks , 2008, PLoS Comput. Biol..

[14]  J Gorecki,et al.  Derivation of a quantitative minimal model from a detailed elementary-step mechanism supported by mathematical coupling analysis. , 2005, The Journal of chemical physics.

[15]  Mudita Singhal,et al.  COPASI - a COmplex PAthway SImulator , 2006, Bioinform..

[16]  Dimitris A. Goussis,et al.  Physical understanding of complex multiscale biochemical models via algorithmic simplification: Glycolysis in Saccharomyces cerevisiae , 2010 .

[17]  C. Daub,et al.  BMC Systems Biology , 2007 .

[18]  L. Petzold Automatic Selection of Methods for Solving Stiff and Nonstiff Systems of Ordinary Differential Equations , 1983 .

[19]  Habib N. Najm,et al.  Model Reduction and Physical Understanding of Slowly Oscillating Processes: The Circadian Cycle , 2006, Multiscale Model. Simul..

[20]  M. Mavrovouniotis,et al.  Simplification of Mathematical Models of Chemical Reaction Systems. , 1998, Chemical reviews.

[21]  S. H. Lam,et al.  Using CSP to Understand Complex Chemical Kinetics ∗ , 1992 .

[22]  Sven Sahle,et al.  Applications and trends in systems biology in biochemistry , 2011, The FEBS journal.

[23]  Melanie I. Stefan,et al.  BioModels Database: An enhanced, curated and annotated resource for published quantitative kinetic models , 2010, BMC Systems Biology.

[24]  Alexander N. Gorban,et al.  Robust simplifications of multiscale biochemical networks , 2008, BMC Systems Biology.

[25]  Gene H. Golub,et al.  Matrix computations , 1983 .

[26]  K. Mauch,et al.  Tendency modeling: a new approach to obtain simplified kinetic models of metabolism applied to Saccharomyces cerevisiae. , 2000, Metabolic engineering.

[27]  Ned S. Wingreen,et al.  Chemotaxis in Escherichia coli: A Molecular Model for Robust Precise Adaptation , 2007, PLoS Comput. Biol..