Simplifying Biochemical Models With Intermediate Species

Abstract

Mathematical models are increasingly being used to understand complex biochemical systems, to analyze experimental data and make predictions about unobserved quantities. However, we rarely know how robust our conclusions are with respect to the choice and uncertainties of the model. Using algebraic techniques we study systematically the effects of intermediate, or transient, species in biochemical systems and provide a simple, yet rigorous mathematical classification of all models obtained from a core model by including intermediates. Main examples include enzymatic and post-translational modification systems, where intermediates often are considered insignificant and neglected in a model, or they are not included because we are unaware of their existence. All possible models obtained from the core model are classified into a finite number of classes. Each class is defined by a mathematically simple canonical model that characterizes crucial dynamical properties, such as mono- and multistationarity and stability of steady states, of all models in the class. We show that if the core model does not have conservation laws, then the introduction of intermediates does not change the steady-state concentrations of the species in the core model, after suitable matching of parameters. Importantly, our results provide guidelines to the modeler in choosing between models and in distinguishing their properties. Further, our work provides a formal way of comparing models that share a common skeleton.

Published in the Journal of the Royal Society Interace. The proof of Proposition 2 in the originally published file was erroneus. The result is though correct and the proof has now been fixed. We thank Magali Giaroli from the University of Buenos Aires for pointing out the error in the proof

References in corpus (3)

• A parameter-free model discrimination criterion based on steady-state coplanarity
• Preclusion of switch behavior in reaction networks with mass-action kinetics
• Cellular compartments cause multistability in biochemical reaction networks and allow cells to process more information

Cited by in corpus (21)

• Identifying parameter regions for multistationarity
• Dynamics of post-translational modification systems: recent progress and future directions
• Memory functions reveal structural properties of gene regulatory networks
• Lower bounds for positive roots and regions of multistationarity in chemical reaction networks
• Elimination of Intermediate Species in Multiscale Stochastic Reaction Networks
• The multistationarity structure of networks with intermediates and a binomial core network
• A proof of unlimited multistability for phosphorylation cycles
• Singular perturbations and scaling
• Topological descriptors of the parameter region of multistationarity: deciding upon connectivity
• Adding species to chemical reaction networks: preserving rank preserves nondegenerate behaviours
• Multistationarity and Bistability for Fewnomial Chemical Reaction Networks
• Building oscillatory chemical reaction networks by adding reversible reactions
• On the sum of chemical reactions
• Parameter regions that give rise to 2[n/2]+1 positive steady states in the n-site phosphorylation system
• Regions of multistationarity in cascades of Goldbeter-Koshland loops
• Groebner bases of reaction networks with intermediate species
• A global convergence result for processive multisite phosphorylation systems
• Graphical Requirements for Multistationarity in Reaction Networks and their Verification in BioModels
• Multistationarity in Structured Reaction Networks
• The inheritance of local bifurcations in mass action networks
• Understanding Multistationarity of Fully Open Reaction Networks