Intermediates, catalysts, persistence, and boundary steady states

Oscillations and bistability in a model of ERK regulation

This work concerns the question of how two important dynamical properties, oscillations and bistability, emerge in an important biological signaling network. Specifically, we consider a model for dual-site phosphorylation and dephosphorylation of extracellular signal-regulated kinase (ERK). We prove that oscillations persist even as the model is greatly simplified (reactions are made irreversible and intermediates are removed). Bistability, however, is much less robust-this property is lost when intermediates are removed or even when all reactions are made irreversible. Moreover, bistability is characterized by the presence of two reversible, catalytic reactions: as other reactions are made irreversible, bistability persists as long as one or both of the specified reactions is preserved. Finally, we investigate the maximum number of steady states, aided by a network's "mixed volume" (a concept from convex geometry). Taken together, our results shed light on the question of how oscillations and bistability emerge from a limiting network of the ERK network-namely, the fully processive dual-site network-which is known to be globally stable and therefore lack both oscillations and bistability. Our proofs are enabled by a Hopf bifurcation criterion due to Yang, analyses of Newton polytopes arising from Hurwitz determinants, and recent characterizations of multistationarity for networks having a steady-state parametrization.

Dynamics of Posttranslational Modification Systems: Recent Progress and Future Directions

Posttranslational modification of proteins is important for signal transduction, and hence significant effort has gone toward understanding how posttranslational modification networks process information. This involves, on the theory side, analyzing the dynamical systems arising from such networks. Which networks are, for instance, bistable? Which networks admit sustained oscillations? Which parameter values enable such behaviors? In this Biophysical Perspective, we highlight recent progress in this area and point out some important future directions. Along the way, we summarize several techniques for analyzing general networks, such as eliminating variables to obtain steady-state parameterizations, and harnessing results on how incorporating intermediates affects dynamics.

Identifying parameter regions for multistationarity

Mathematical modelling has become an established tool for studying the dynamics of biological systems. Current applications range from building models that reproduce quantitative data to identifying systems with predefined qualitative features, such as switching behaviour, bistability or oscillations. Mathematically, the latter question amounts to identifying parameter values associated with a given qualitative feature. We introduce a procedure to partition the parameter space of a parameterized system of ordinary differential equations into regions for which the system has a unique or multiple equilibria. The procedure is based on the computation of the Brouwer degree, and it creates a multivariate polynomial with parameter depending coefficients. The signs of the coefficients determine parameter regions with and without multistationarity. A particular strength of the procedure is the avoidance of numerical analysis and parameter sampling. The procedure consists of a number of steps. Each of these steps might be addressed algorithmically using various computer programs and available software, or manually. We demonstrate our procedure on several models of gene transcription and cell signalling, and show that in many cases we obtain a complete partitioning of the parameter space with respect to multistationarity.

Intermediates and Generic Convergence to Equilibria

Known graphical conditions for the generic and global convergence to equilibria of the dynamical system arising from a reaction network are shown to be invariant under the so-called successive removal of intermediates, a systematic procedure to simplify the network, making the graphical conditions considerably easier to check.