-
In distributive phosphorylation catalytic constants enable non-trivial dynamics
Authors:
Carsten Conradi,
Maya Mincheva
Abstract:
Ordered distributive double phosphorylation is a recurrent motif in intracellular signaling and control. It is either sequential (where the site phosphorylated last is dephosphorylated first) or cyclic (where the site phosphorylated first is dephosphorylated first). Sequential distributive double phosphorylation has been extensively studied and an inequality involving only the catalytic constants…
▽ More
Ordered distributive double phosphorylation is a recurrent motif in intracellular signaling and control. It is either sequential (where the site phosphorylated last is dephosphorylated first) or cyclic (where the site phosphorylated first is dephosphorylated first). Sequential distributive double phosphorylation has been extensively studied and an inequality involving only the catalytic constants of kinase and phosphatase is known to be sufficient for multistationarity. As multistationarity is necessary for bistability it has been argued that these constants enable bistability.
Here we show for cyclic distributive double phosphorylation that if its catalytic constants satisfy the very same inequality, then Hopf bifurcations and hence sustained oscillations can occur. Hence we argue that in distributive double phosphorylation (sequential or distributive) the catalytic constants enable non-trivial dynamics.
In fact, if the rate constant values in a network of cyclic distributive double phosphorylation are such that Hopf bifurcations and sustained oscillations can occur, then a network of sequential distributive double phosphorylation with the same rate constant values will show multistationarity -- albeit for different values of the total concentrations. For cyclic distributive double phosphorylation we further describe a procedure to generate rate constant values where Hopf bifurcations and hence sustained oscillations can occur. This may, for example, allow for an efficient sampling of oscillatory regions in parameter space.
Our analysis is greatly simplified by the fact that it is possible to reduce the network of cyclic distributive double phosphorylation to what we call a network with a single extreme ray. We summarize key properties of these networks.
△ Less
Submitted 8 February, 2024; v1 submitted 16 August, 2023;
originally announced August 2023.
-
Dynamics of ERK regulation in the processive limit
Authors:
Carsten Conradi,
Nida Obatake,
Anne Shiu,
Xiaoxian Tang
Abstract:
We consider a model of extracellular signal-regulated kinase (ERK) regulation by dual-site phosphorylation and dephosphorylation, which exhibits bistability and oscillations, but loses these properties in the limit in which the mechanisms underlying phosphorylation and dephosphorylation become processive. Our results suggest that anywhere along the way to becoming processive, the model remains bis…
▽ More
We consider a model of extracellular signal-regulated kinase (ERK) regulation by dual-site phosphorylation and dephosphorylation, which exhibits bistability and oscillations, but loses these properties in the limit in which the mechanisms underlying phosphorylation and dephosphorylation become processive. Our results suggest that anywhere along the way to becoming processive, the model remains bistable and oscillatory. More precisely, in simplified versions of the model, precursors to bistability and oscillations (specifically, multistationarity and Hopf bifurcations, respectively) exist at all "processivity levels". Finally, we investigate whether bistability and oscillations can exist together.
△ Less
Submitted 3 September, 2020; v1 submitted 31 October, 2019;
originally announced October 2019.
-
On the existence of Hopf bifurcations in the sequential and distributive double phosphorylation cycle
Authors:
Carsten Conradi,
Elisenda Feliu,
Maya Mincheva
Abstract:
Protein phosphorylation cycles are important mechanisms of the post translational modification of a protein and as such an integral part of intracellular signaling and control. We consider the sequential phosphorylation and dephosphorylation of a protein at two binding sites. While it is known that proteins where phosphorylation is processive and dephosphorylation is distributive admit oscillation…
▽ More
Protein phosphorylation cycles are important mechanisms of the post translational modification of a protein and as such an integral part of intracellular signaling and control. We consider the sequential phosphorylation and dephosphorylation of a protein at two binding sites. While it is known that proteins where phosphorylation is processive and dephosphorylation is distributive admit oscillations (for some value of the rate constants and total concentrations) it is not known whether or not this is the case if both phosphorylation and dephosphorylation are distributive. We study four simplified mass action models of sequential and distributive phosphorylation and show that for each of those there do not exist rate constants and total concentrations where a Hopf bifurcation occurs. To arrive at this result we use convex parameters to parameterize the steady state and Hurwitz matrices.
△ Less
Submitted 4 September, 2019; v1 submitted 20 May, 2019;
originally announced May 2019.
-
Multistationarity in the space of total concentrations for systems that admit a monomial parametrization
Authors:
Carsten Conradi,
Alexandru Iosif,
Thomas Kahle
Abstract:
We apply tools from real algebraic geometry to the problem of multistationarity of chemical reaction networks. A particular focus is on the case of reaction networks whose steady states admit a monomial parametrization. For such systems we show that in the space of total concentrations multistationarity is scale invariant: if there is multistationarity for some value of the total concentrations, t…
▽ More
We apply tools from real algebraic geometry to the problem of multistationarity of chemical reaction networks. A particular focus is on the case of reaction networks whose steady states admit a monomial parametrization. For such systems we show that in the space of total concentrations multistationarity is scale invariant: if there is multistationarity for some value of the total concentrations, then there is multistationarity on the entire ray containing this value (possibly for different rate constants) -- and vice versa. Moreover, for these networks it is possible to decide about multistationarity independent of the rate constants by formulating semi-algebraic conditions that involve only concentration variables. These conditions can easily be extended to include total concentrations. Hence quantifier elimination may give new insights into multistationarity regions in the space of total concentrations. To demonstrate this, we show that for the distributive phosphorylation of a protein at two binding sites multistationarity is only possible if the total concentration of the substrate is larger than either the total concentration of the kinase or the total concentration of the phosphatase. This result is enabled by the chamber decomposition of the space of total concentrations from polyhedral geometry. Together with the corresponding sufficiency result of Bihan et al. this yields a characterization of multistationarity up to lower dimensional regions.
△ Less
Submitted 13 August, 2019; v1 submitted 18 October, 2018;
originally announced October 2018.
-
Emergence of oscillations in a mixed-mechanism phosphorylation system
Authors:
Carsten Conradi,
Maya Mincheva,
Anne Shiu
Abstract:
This work investigates the emergence of oscillations in one of the simplest cellular signaling networks exhibiting oscillations, namely, the dual-site phosphorylation and dephosphorylation network (futile cycle), in which the mechanism for phosphorylation is processive while the one for dephosphorylation is distributive (or vice-versa). The fact that this network yields oscillations was shown rece…
▽ More
This work investigates the emergence of oscillations in one of the simplest cellular signaling networks exhibiting oscillations, namely, the dual-site phosphorylation and dephosphorylation network (futile cycle), in which the mechanism for phosphorylation is processive while the one for dephosphorylation is distributive (or vice-versa). The fact that this network yields oscillations was shown recently by Suwanmajo and Krishnan. Our results, which significantly extend their analyses, are as follows. First, in the three-dimensional space of total amounts, the border between systems with a stable versus unstable steady state is a surface defined by the vanishing of a single Hurwitz determinant. Second, this surface consists generically of simple Hopf bifurcations. Next, simulations suggest that when the steady state is unstable, oscillations are the norm. Finally, the emergence of oscillations via a Hopf bifurcation is enabled by the catalytic and association constants of the distributive part of the mechanism: if these rate constants satisfy two inequalities, then the system generically admits a Hopf bifurcation. Our proofs are enabled by the Routh-Hurwitz criterion, a Hopf-bifurcation criterion due to Yang, and a monomial parametrization of steady states.
△ Less
Submitted 1 February, 2019; v1 submitted 8 September, 2018;
originally announced September 2018.
-
Dynamics of post-translational modification systems: recent progress and future directions
Authors:
Carsten Conradi,
Anne Shiu
Abstract:
Post-translational modification (PTM) of proteins plays a key role in signal transduction, and hence significant effort has gone toward understanding how PTM 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 enabl…
▽ More
Post-translational modification (PTM) of proteins plays a key role in signal transduction, and hence significant effort has gone toward understanding how PTM 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 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 parametrizations, and harnessing results on how incorporating intermediates affects dynamics.
△ Less
Submitted 24 October, 2017; v1 submitted 30 May, 2017;
originally announced May 2017.
-
Identifying parameter regions for multistationarity
Authors:
Carsten Conradi,
Elisenda Feliu,
Maya Mincheva,
Carsten Wiuf
Abstract:
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 giv…
▽ More
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.
△ Less
Submitted 6 May, 2017; v1 submitted 13 August, 2016;
originally announced August 2016.
-
Detecting Binomiality
Authors:
Carsten Conradi,
Thomas Kahle
Abstract:
Binomial ideals are special polynomial ideals with many algorithmically and theoretically nice properties. We discuss the problem of deciding if a given polynomial ideal is binomial. While the methods are general, our main motivation and source of examples is the simplification of steady state equations of chemical reaction networks. For homogeneous ideals we give an efficient, Gröbner-free algori…
▽ More
Binomial ideals are special polynomial ideals with many algorithmically and theoretically nice properties. We discuss the problem of deciding if a given polynomial ideal is binomial. While the methods are general, our main motivation and source of examples is the simplification of steady state equations of chemical reaction networks. For homogeneous ideals we give an efficient, Gröbner-free algorithm for binomiality detection, based on linear algebra only. On inhomogeneous input the algorithm can only give a sufficient condition for binomiality. As a remedy we construct a heuristic toolbox that can lead to simplifications even if the given ideal is not binomial.
△ Less
Submitted 10 September, 2015; v1 submitted 17 February, 2015;
originally announced February 2015.
-
Graph-theoretic analysis of multistationarity using degree theory
Authors:
Carsten Conradi,
Maya Mincheva
Abstract:
Biochemical mechanisms with mass action kinetics are often modeled by systems of polynomial differential equations (DE). Determining directly if the DE system has multiple equilibria (multistationarity) is difficult for realistic systems, since they are large, nonlinear and contain many unknown parameters. Mass action biochemical mechanisms can be represented by a directed bipartite graph with spe…
▽ More
Biochemical mechanisms with mass action kinetics are often modeled by systems of polynomial differential equations (DE). Determining directly if the DE system has multiple equilibria (multistationarity) is difficult for realistic systems, since they are large, nonlinear and contain many unknown parameters. Mass action biochemical mechanisms can be represented by a directed bipartite graph with species and reaction nodes. Graph-theoretic methods can then be used to assess the potential of a given biochemical mechanism for multistationarity by identifying structures in the bipartite graph referred to as critical fragments. In this article we present a graph-theoretic method for conservative biochemical mechanisms characterized by bounded species concentrations, which makes the use of degree theory arguments possible. We illustrate the results with an example of a MAPK network.
△ Less
Submitted 11 November, 2014;
originally announced November 2014.
-
A global convergence result for processive multisite phosphorylation systems
Authors:
Carsten Conradi,
Anne Shiu
Abstract:
Multisite phosphorylation plays an important role in intracellular signaling. There has been much recent work aimed at understanding the dynamics of such systems when the phosphorylation/dephosphorylation mechanism is distributive, that is, when the binding of a substrate and an enzyme molecule results in addition or removal of a single phosphate group and repeated binding therefore is required fo…
▽ More
Multisite phosphorylation plays an important role in intracellular signaling. There has been much recent work aimed at understanding the dynamics of such systems when the phosphorylation/dephosphorylation mechanism is distributive, that is, when the binding of a substrate and an enzyme molecule results in addition or removal of a single phosphate group and repeated binding therefore is required for multisite phosphorylation. In particular, such systems admit bistability. Here we analyze a different class of multisite systems, in which the binding of a substrate and an enzyme molecule results in addition or removal of phosphate groups at all phosphorylation sites. That is, we consider systems in which the mechanism is processive, rather than distributive. We show that in contrast with distributive systems, processive systems modeled with mass-action kinetics do not admit bistability and, moreover, exhibit rigid dynamics: each invariant set contains a unique equilibrium, which is a global attractor. Additionally, we obtain a monomial parametrization of the steady states. Our proofs rely on a technique of Johnston for using "translated" networks to study systems with "toric steady states", recently given sign conditions for injectivity of polynomial maps, and a result from monotone systems theory due to Angeli and Sontag.
△ Less
Submitted 21 October, 2014; v1 submitted 22 April, 2014;
originally announced April 2014.
-
N-site phosphorylation systems with 2N-1 steady states
Authors:
Dietrich Flockerzi,
Katharina Holstein,
Carsten Conradi
Abstract:
Multisite protein phosphorylation plays a prominent role in intracellular processes like signal transduction, cell-cycle control and nuclear signal integration. Many proteins are phosphorylated in a sequential and distributive way at more than one phosphorylation site. Mathematical models of $n$-site sequential distributive phosphorylation are therefore studied frequently. In particular, in {\em W…
▽ More
Multisite protein phosphorylation plays a prominent role in intracellular processes like signal transduction, cell-cycle control and nuclear signal integration. Many proteins are phosphorylated in a sequential and distributive way at more than one phosphorylation site. Mathematical models of $n$-site sequential distributive phosphorylation are therefore studied frequently. In particular, in {\em Wang and Sontag, 2008,} it is shown that models of $n$-site sequential distributive phosphorylation admit at most $2n-1$ steady states. Wang and Sontag furthermore conjecture that for odd $n$, there are at most $n$ and that, for even $n$, there are at most $n+1$ steady states. This, however, is not true: building on earlier work in {\em Holstein et.al., 2013}, we present a scalar determining equation for multistationarity which will lead to parameter values where a $3$-site system has $5$ steady states and parameter values where a $4$-site system has $7$ steady states. Our results therefore are counterexamples to the conjecture of Wang and Sontag. We furthermore study the inherent geometric properties of multistationarity in $n$-site sequential distributive phosphorylation: the complete vector of steady state ratios is determined by the steady state ratios of free enzymes and unphosphorylated protein and there exists a linear relationship between steady state ratios of phosphorylated protein.
△ Less
Submitted 23 April, 2014; v1 submitted 17 December, 2013;
originally announced December 2013.
-
Sign conditions for injectivity of generalized polynomial maps with applications to chemical reaction networks and real algebraic geometry
Authors:
Stefan Müller,
Elisenda Feliu,
Georg Regensburger,
Carsten Conradi,
Anne Shiu,
Alicia Dickenstein
Abstract:
We give necessary and sufficient conditions in terms of sign vectors for the injectivity of families of polynomial maps with arbitrary real exponents defined on the positive orthant. Our work relates and extends existing injectivity conditions expressed in terms of Jacobian matrices and determinants. In the context of chemical reaction networks with power-law kinetics, our results can be used to p…
▽ More
We give necessary and sufficient conditions in terms of sign vectors for the injectivity of families of polynomial maps with arbitrary real exponents defined on the positive orthant. Our work relates and extends existing injectivity conditions expressed in terms of Jacobian matrices and determinants. In the context of chemical reaction networks with power-law kinetics, our results can be used to preclude as well as to guarantee multiple positive steady states. In the context of real algebraic geometry,our work recognizes a prior result of Craciun, Garcia-Puente, and Sottile, together with work of two of the authors, as the first partial multivariate generalization of the classical Descartes' rule, which bounds the number of positive real roots of a univariate real polynomial in terms of the number of sign variations of its coefficients.
△ Less
Submitted 30 October, 2014; v1 submitted 21 November, 2013;
originally announced November 2013.
-
Multistationarity in sequential distributed multisite phosphorylation networks
Authors:
Katharina Holstein,
Dietrich Flockerzi,
Carsten Conradi
Abstract:
Multisite phosphorylation networks are encountered in many intracellular processes like signal transduction, cell-cycle control or nuclear signal integration. In this contribution networks describing the phosphorylation and dephosphorylation of a protein at $n$ sites in a sequential distributive mechanism are considered. Multistationarity (i.e.\ the existence of at least two positive steady state…
▽ More
Multisite phosphorylation networks are encountered in many intracellular processes like signal transduction, cell-cycle control or nuclear signal integration. In this contribution networks describing the phosphorylation and dephosphorylation of a protein at $n$ sites in a sequential distributive mechanism are considered. Multistationarity (i.e.\ the existence of at least two positive steady state solutions of the associated polynomial dynamical system) has been analyzed and established in several contributions. It is, for example, known that there exist values for he rate constants where multistationarity occurs. However, nothing else is known about these rate constants.
Here we present a sign condition that is necessary and sufficient for multistationarity in $n$-site sequential, distributive phosphorylation. We express this sign condition in terms of linear systems and show that solutions of these systems define rate constants where multistationarity is possible. We then present, for $n\geq 2$, a collection of {\em feasible} linear systems and hence give a new and independent proof that multistationarity is possible for $n\geq 2$. Moreover, our results allow to explicitly obtain values for the rate constants where multistationarity is possible. Hence we believe that, for the first time, a systematic exploration of the region in parameter space where multistationarity occurs has become possible.One consequence of our work is that, for any pair of steady states, the ratio of the steady state concentrations of kinase-substrate complexes equals that of phosphatase-substrate complexes.
△ Less
Submitted 2 July, 2013; v1 submitted 24 April, 2013;
originally announced April 2013.
-
Chemical reaction systems with toric steady states
Authors:
Mercedes Perez Millan,
Alicia Dickenstein,
Anne Shiu,
Carsten Conradi
Abstract:
Mass-action chemical reaction systems are frequently used in Computational Biology. The corresponding polynomial dynamical systems are often large (consisting of tens or even hundreds of ordinary differential equations) and poorly parametrized (due to noisy measurement data and a small number of data points and repetitions). Therefore, it is often difficult to establish the existence of (positive)…
▽ More
Mass-action chemical reaction systems are frequently used in Computational Biology. The corresponding polynomial dynamical systems are often large (consisting of tens or even hundreds of ordinary differential equations) and poorly parametrized (due to noisy measurement data and a small number of data points and repetitions). Therefore, it is often difficult to establish the existence of (positive) steady states or to determine whether more complicated phenomena such as multistationarity exist. If, however, the steady state ideal of the system is a binomial ideal, then we show that these questions can be answered easily. The focus of this work is on systems with this property, and we say that such systems have toric steady states. Our main result gives sufficient conditions for a chemical reaction system to have toric steady states. Furthermore, we analyze the capacity of such a system to exhibit positive steady states and multistationarity. Examples of systems with toric steady states include weakly-reversible zero-deficiency chemical reaction systems. An important application of our work concerns the networks that describe the multisite phosphorylation of a protein by a kinase/phosphatase pair in a sequential and distributive mechanism.
△ Less
Submitted 12 July, 2011; v1 submitted 8 February, 2011;
originally announced February 2011.
-
Switching in mass action networks based on linear inequalities
Authors:
Carsten Conradi,
Dietrich Flockerzi
Abstract:
Many biochemical processes can successfully be described by dynamical systems allowing some form of switching when, depending on their initial conditions, solutions of the dynamical system end up in different regions of state space (associated with different biochemical functions). Switching is often realized by a bistable system (i.e. a dynamical system allowing two stable steady state solutions)…
▽ More
Many biochemical processes can successfully be described by dynamical systems allowing some form of switching when, depending on their initial conditions, solutions of the dynamical system end up in different regions of state space (associated with different biochemical functions). Switching is often realized by a bistable system (i.e. a dynamical system allowing two stable steady state solutions) and, in the majority of cases, bistability is established numerically. In our point of view this approach is too restrictive, as, one the one hand, due to predominant parameter uncertainty numerical methods are generally difficult to apply to realistic models originating in Systems Biology. And on the other hand switching already arises with the occurrence of a saddle type steady state (characterized by a Jacobian where exactly one Eigenvalue is positive and the remaining eigenvalues have negative real part). Consequently we derive conditions based on linear inequalities that allow the analytic computation of states and parameters where the Jacobian derived from a mass action network has a defective zero eigenvalue so that -- under certain genericity conditions -- a saddle-node bifurcation occurs. Our conditions are applicable to general mass action networks involving at least one conservation relation, however, they are only sufficient (as infeasibility of linear inequalities does not exclude defective zero eigenvalues).
△ Less
Submitted 26 September, 2011; v1 submitted 4 February, 2010;
originally announced February 2010.
-
Multiple positive steady states in subnetworks defined by stoichiometric generators
Authors:
Carsten Conradi
Abstract:
In Systems Biology there is a growing interest in the question, whether or not a given mathematical model can admit more than one steady state. As parameter values are often unknown or subject to a very high uncertainty, one is often interested in the question, whether or not a given mathematical model can, for some conceivable parameter vector, exhibit multistationarity at all. A partial answer…
▽ More
In Systems Biology there is a growing interest in the question, whether or not a given mathematical model can admit more than one steady state. As parameter values are often unknown or subject to a very high uncertainty, one is often interested in the question, whether or not a given mathematical model can, for some conceivable parameter vector, exhibit multistationarity at all. A partial answer to this question is given in Feinberg's deficiency one algorithm. This algorithm can decide about multistationarity by analyzing systems of linear inequalities that are independent of parameter values. However, the deficiency one algorithm is limited to what its author calls regular deficiency one networks. Many realistic networks have a deficiency higher than one, thus the algorithm cannot be applied directly. In a previous publication it was suggested to analyze certain well defined subnetworks that are guaranteed to be of deficiency one. Realistic reaction networks, however, often lead to subnetworks that are irregular, especially if metabolic networks are considered. Here the special structure of the subnetworks is used to derive conditions for multistationarity. These conditions are independent of the regularity conditions required by the deficiency one algorithm. Thus, in particular, these conditions are applicable to irregular subnetworks.
△ Less
Submitted 22 September, 2009;
originally announced September 2009.