A method for zooming of nonlinear models of biochemical systems

BackgroundModels of biochemical systems are typically complex, which may complicate the discovery of cardinal biochemical principles. It is therefore important to single out the parts of a model that are essential for the function of the system, so that the remaining non-essential parts can be eliminated. However, each component of a mechanistic model has a clear biochemical interpretation, and it is desirable to conserve as much of this interpretability as possible in the reduction process. Furthermore, it is of great advantage if we can translate predictions from the reduced model to the original model.ResultsIn this paper we present a novel method for model reduction that generates reduced models with a clear biochemical interpretation. Unlike conventional methods for model reduction our method enables the mapping of predictions by the reduced model to the corresponding detailed predictions by the original model. The method is based on proper lumping of state variables interacting on short time scales and on the computation of fraction parameters, which serve as the link between the reduced model and the original model. We illustrate the advantages of the proposed method by applying it to two biochemical models. The first model is of modest size and is commonly occurring as a part of larger models. The second model describes glucose transport across the cell membrane in baker's yeast. Both models can be significantly reduced with the proposed method, at the same time as the interpretability is conserved.ConclusionsWe introduce a novel method for reduction of biochemical models that is compatible with the concept of zooming. Zooming allows the modeler to work on different levels of model granularity, and enables a direct interpretation of how modifications to the model on one level affect the model on other levels in the hierarchy. The method extends the applicability of the method that was previously developed for zooming of linear biochemical models to nonlinear models.

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

[2]  F. A. Davidson,et al.  Effects of periodic input on the quasi-steady state assumptions for enzyme-catalysed reactions , 2005, Journal of mathematical biology.

[3]  Edda Klipp,et al.  Biochemical network models simplified by balanced truncation , 2005, The FEBS journal.

[4]  A. Bulychev,et al.  Photo-electrochemical control of photosystem II chlorophyll fluorescence in vivo. , 2002, Bioelectrochemistry.

[5]  Herschel Rabitz,et al.  The Effect of Lumping and Expanding on Kinetic Differential Equations , 1997, SIAM J. Appl. Math..

[6]  Michael B. Yaffe,et al.  Data-driven modelling of signal-transduction networks , 2006, Nature Reviews Molecular Cell Biology.

[7]  T. Maiwald,et al.  Materials and Methods SOM Text Figs. S1 to S16 References Materials and Methods , 2022 .

[8]  Katherine C. Chen,et al.  Integrative analysis of cell cycle control in budding yeast. , 2004, Molecular biology of the cell.

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

[10]  Mikael Sunnåker,et al.  Zooming of states and parameters using a lumping approach including back-translation , 2010, BMC Systems Biology.

[11]  Jacob Roll,et al.  Systems biology: model based evaluation and comparison of potential explanations for given biological data , 2009, The FEBS journal.

[12]  S. Schnell,et al.  Use and abuse of the quasi-steady-state approximation. , 2006, Systems biology.

[13]  John J Tyson,et al.  Bifurcation analysis of a model of the budding yeast cell cycle. , 2004, Chaos.

[14]  R. M. Murray,et al.  Model reduction of interconnected linear systems , 2009 .

[15]  E. Gilles,et al.  Computational modeling of the dynamics of the MAP kinase cascade activated by surface and internalized EGF receptors , 2002, Nature Biotechnology.

[16]  L. Ljung,et al.  Control theory : multivariable and nonlinear methods , 2000 .

[17]  Julio Saez-Rodriguez,et al.  A domain-oriented approach to the reduction of combinatorial complexity in signal transduction networks , 2006, BMC Bioinformatics.

[18]  E D Gilles,et al.  Reduction of mathematical models of signal transduction networks: simulation-based approach applied to EGF receptor signalling. , 2004, Systems biology.

[19]  Olivia Eriksson.  Simplicity within Complexity : Understanding dynamics of cellular networks by model reduction , 2008 .

[20]  Xingming Zhao,et al.  Computational Systems Biology , 2013, TheScientificWorldJournal.

[21]  Moshe Y. Vardi,et al.  Dynamic and static limitation in multiscale reaction networks, revisited , 2007, physics/0703278.

[22]  Alexander Lanzon.  Directions, applications and methods in robust control , 2009 .

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

[24]  L. Aarons,et al.  Proper lumping in systems biology models. , 2009, IET systems biology.

[25]  Iliya V. Karlin,et al.  Method of invariant manifold for chemical kinetics , 2003 .

[27]  D. Lauffenburger,et al.  Physicochemical modelling of cell signalling pathways , 2006, Nature Cell Biology.

[28]  Henning Schmidt,et al.  Complexity reduction of biochemical rate expressions , 2008, Bioinform..

[29]  Ernst Dieter Gilles,et al.  Exact model reduction of combinatorial reaction networks , 2008, BMC Systems Biology.

[30]  E. Querfurth,et al.  In vivo investigations of glucose transport in Saccharomyces cerevisiae , 2000, Biotechnology and bioengineering.

[31]  A. Tzafriri,et al.  The total quasi-steady-state approximation is valid for reversible enzyme kinetics. , 2004, Journal of theoretical biology.

[32]  Jake V. Bouvrie,et al.  Model Reduction for Nonlinear Control Systems using Kernel Subspace Methods , 2011, ArXiv.

[33]  Stephen Gilmore,et al.  Modelling the Influence of RKIP on the ERK Signalling Pathway Using the Stochastic Process Algebra PEPA , 2006, Trans. Comp. Sys. Biology.

[34]  E W Jacobsen,et al.  Structural robustness of biochemical network models-with application to the oscillatory metabolism of activated neutrophils. , 2008, IET systems biology.

[35]  Michael Ederer,et al.  Reduced modeling of signal transduction – a modular approach , 2007, BMC Bioinformatics.

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

[37]  Thomas F. Edgar,et al.  An improved method for nonlinear model reduction using balancing of empirical gramians , 2002 .

[38]  Gunnar Cedersund,et al.  Mass and Information Feedbacks through Receptor Endocytosis Govern Insulin Signaling as Revealed Using a Parameter-free Modeling Framework* , 2010, The Journal of Biological Chemistry.

[39]  Alexander N Gorban,et al.  Dynamical robustness of biological networks with hierarchical distribution of time scales. , 2007, IET systems biology.

[40]  Mats Jirstrand,et al.  Systems biology Systems Biology Toolbox for MATLAB : a computational platform for research in systems biology , 2006 .

[41]  Barbara M. Bakker,et al.  Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry. , 2000, European journal of biochemistry.

[42]  J. Doyle,et al.  Robust and optimal control , 1995, Proceedings of 35th IEEE Conference on Decision and Control.

[43]  Eric Bullinger,et al.  Approximations and their consequences for dynamic modelling of signal transduction pathways. , 2007, Mathematical biosciences.

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

[45]  O. Radulescu,et al.  Dynamic and static limitation in reaction networks , revisited , 2007 .