-
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.
-
A graph-theoretic condition for delay stability of reaction systems
Authors:
Gheorghe Craciun,
Maya Mincheva,
Casian Pantea,
Polly Y. Yu
Abstract:
Delay mass-action systems provide a model of chemical kinetics when past states influence the current dynamics. In this work, we provide a graph-theoretic condition for delay stability, i.e., linear stability independent of both rate constants and delay parameters. In particular, the result applies when the system has no delay, implying asymptotic stability for the ODE system. The graph-theoretic…
▽ More
Delay mass-action systems provide a model of chemical kinetics when past states influence the current dynamics. In this work, we provide a graph-theoretic condition for delay stability, i.e., linear stability independent of both rate constants and delay parameters. In particular, the result applies when the system has no delay, implying asymptotic stability for the ODE system. The graph-theoretic condition is about cycles in the directed species-reaction graph of the network, which encodes how different species in the system interact.
△ Less
Submitted 15 May, 2021;
originally announced May 2021.
-
Delay stability of reaction systems
Authors:
Gheorghe Craciun,
Maya Mincheva,
Casian Pantea,
Polly Y. Yu
Abstract:
Delay differential equations are used as a model when the effect of past states has to be taken into account. In this work we consider delay models of chemical reaction networks with mass action kinetics. We obtain a sufficient condition for absolute delay stability of equilibrium concentrations, i.e., local asymptotic stability independent of the delay parameters. Several interesting examples on…
▽ More
Delay differential equations are used as a model when the effect of past states has to be taken into account. In this work we consider delay models of chemical reaction networks with mass action kinetics. We obtain a sufficient condition for absolute delay stability of equilibrium concentrations, i.e., local asymptotic stability independent of the delay parameters. Several interesting examples on sequestration networks with delays are presented.
△ Less
Submitted 4 June, 2020; v1 submitted 10 March, 2020;
originally announced March 2020.
-
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.
-
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.
-
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.
-
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.
-
Large Covariance Estimation by Thresholding Principal Orthogonal Complements
Authors:
Jianqing Fan,
Yuan Liao,
Martina Mincheva
Abstract:
This paper deals with the estimation of a high-dimensional covariance with a conditional sparsity structure and fast-diverging eigenvalues. By assuming sparse error covariance matrix in an approximate factor model, we allow for the presence of some cross-sectional correlation even after taking out common but unobservable factors. We introduce the Principal Orthogonal complEment Thresholding (POET)…
▽ More
This paper deals with the estimation of a high-dimensional covariance with a conditional sparsity structure and fast-diverging eigenvalues. By assuming sparse error covariance matrix in an approximate factor model, we allow for the presence of some cross-sectional correlation even after taking out common but unobservable factors. We introduce the Principal Orthogonal complEment Thresholding (POET) method to explore such an approximate factor structure with sparsity. The POET estimator includes the sample covariance matrix, the factor-based covariance matrix (Fan, Fan, and Lv, 2008), the thresholding estimator (Bickel and Levina, 2008) and the adaptive thresholding estimator (Cai and Liu, 2011) as specific examples. We provide mathematical insights when the factor analysis is approximately the same as the principal component analysis for high-dimensional data. The rates of convergence of the sparse residual covariance matrix and the conditional sparse covariance matrix are studied under various norms. It is shown that the impact of estimating the unknown factors vanishes as the dimensionality increases. The uniform rates of convergence for the unobserved factors and their factor loadings are derived. The asymptotic results are also verified by extensive simulation studies. Finally, a real data application on portfolio allocation is presented.
△ Less
Submitted 13 January, 2013; v1 submitted 30 December, 2011;
originally announced January 2012.
-
High-dimensional covariance matrix estimation in approximate factor models
Authors:
Jianqing Fan,
Yuan Liao,
Martina Mincheva
Abstract:
The variance--covariance matrix plays a central role in the inferential theories of high-dimensional factor models in finance and economics. Popular regularization methods of directly exploiting sparsity are not directly applicable to many financial problems. Classical methods of estimating the covariance matrices are based on the strict factor models, assuming independent idiosyncratic components…
▽ More
The variance--covariance matrix plays a central role in the inferential theories of high-dimensional factor models in finance and economics. Popular regularization methods of directly exploiting sparsity are not directly applicable to many financial problems. Classical methods of estimating the covariance matrices are based on the strict factor models, assuming independent idiosyncratic components. This assumption, however, is restrictive in practical applications. By assuming sparse error covariance matrix, we allow the presence of the cross-sectional correlation even after taking out common factors, and it enables us to combine the merits of both methods. We estimate the sparse covariance using the adaptive thresholding technique as in Cai and Liu [J. Amer. Statist. Assoc. 106 (2011) 672--684], taking into account the fact that direct observations of the idiosyncratic components are unavailable. The impact of high dimensionality on the covariance matrix estimation based on the factor structure is then studied.
△ Less
Submitted 14 March, 2012; v1 submitted 21 May, 2011;
originally announced May 2011.