Proper lumping in systems biology models.

An algorithm for automatic order reduction of models defined by large systems of differential equations is presented. The algorithm was developed with systems biology models in mind and the motivation behind it is to develop mechanistic pharmacokinetic/pharmacodynamic models from already available systems biology models. The approach used for model reduction is proper lumping of the system's states and is based on a search through the possible combinations of lumps. To avoid combinatorial explosion, a heuristic, greedy search strategy is employed and comparison with the full exhaustive search provides evidence that it performs well. The method takes advantage of an apparent property of this kind of systems that lumps remain consistent over different levels of order reduction. Advantages of the method presented include: the variables and parameters of the reduced model retain a specific physiological meaning; the algorithm is automatic and easy to use; it can be used for nonlinear models and can handle parameter uncertainty and constraints. The algorithm was applied to a model of NF-B signalling pathways in order to demonstrate its use and performance. Significant reduction was achieved for the example model, while agreement with the original model was proportional to the size of the reduced model, as expected. The results of the model reduction were compared with a published, intuitively reduced model of NF-B signalling pathways and were found to be in agreement, in terms of the identified key species for the system's kinetic behaviour. The method may provide useful insights which are complementary to the intuitive reduction approach, especially in large systems.

[1]  L B Sheiner,et al.  Forecasting individual pharmacokinetics , 1979, Clinical pharmacology and therapeutics.

[2]  B. Moore Principal component analysis in linear systems: Controllability, observability, and model reduction , 1981 .

[3]  Raffaella Ocone,et al.  Lumping nonlinear kinetics , 1988 .

[4]  Herschel Rabitz,et al.  A general analysis of exact lumping in chemical kinetics , 1989 .

[5]  Herschel Rabitz,et al.  A general analysis of approximate lumping in chemical kinetics , 1990 .

[6]  Herschel Rabitz,et al.  A GENERAL LUMPING ANALYSIS OF A REACTION SYSTEM COUPLED WITH DIFFUSION , 1991 .

[7]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[8]  Vasilios Manousiouthakis,et al.  KINETIC MODEL REDUCTION USING GENETIC ALGORITHMS , 1998 .

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

[10]  M. D. McKay,et al.  A comparison of three methods for selecting values of input variables in the analysis of output from a computer code , 2000 .

[11]  L B Sheiner,et al.  Pharmacokinetic/pharmacodynamic modeling in drug development. , 2000, Annual review of pharmacology and toxicology.

[12]  Michal Linial,et al.  Using Bayesian Networks to Analyze Expression Data , 2000, J. Comput. Biol..

[13]  A. Hoffmann,et al.  The I (cid:1) B –NF-(cid:1) B Signaling Module: Temporal Control and Selective Gene Activation , 2022 .

[14]  Paul I. Barton,et al.  Optimally-reduced kinetic models: reaction elimination in large-scale kinetic mechanisms , 2003 .

[15]  Malcolm Rowland,et al.  Lumping of Whole-Body Physiologically Based Pharmacokinetic Models , 1998, Journal of Pharmacokinetics and Biopharmaceutics.

[16]  James R. Johnson,et al.  Oscillations in NF-κB Signaling Control the Dynamics of Gene Expression , 2004, Science.

[17]  János Tóth,et al.  Lumping in Pharmacokinetics , 2005, Journal of Pharmacokinetics and Pharmacodynamics.

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

[19]  K. Sneppen,et al.  Minimal model of spiky oscillations in NF-κB signaling , 2006 .

[20]  Douglas B Kell,et al.  Proximate parameter tuning for biochemical networks with uncertain kinetic parameters. , 2008, Molecular bioSystems.