-
Beyond Boolean networks, a multi-valued approach
Authors:
J. García Galofre,
M. Pérez Millán,
A. Galarza Rial,
R. Laubenbacher,
A. Dickenstein
Abstract:
Boolean networks can be viewed as functions on the set of binary strings of a given length, described via logical rules. They were introduced as dynamic models into biology, in particular as logical models of intracellular regulatory networks involving genes, proteins, and metabolites. Since genes can have several modes of action, depending on their expression levels, binary variables are often no…
▽ More
Boolean networks can be viewed as functions on the set of binary strings of a given length, described via logical rules. They were introduced as dynamic models into biology, in particular as logical models of intracellular regulatory networks involving genes, proteins, and metabolites. Since genes can have several modes of action, depending on their expression levels, binary variables are often not sufficiently rich, requiring the use of multi-valued networks instead. The steady state analysis of Boolean networks is computationally complex, and increasing the number of variable values beyond $2$ adds substantially to this complexity, and no general methods are available beyond simulation. The main contribution of this paper is to give an algorithm to compute the steady states of a multi-valued network that has a complexity that, in many cases, is essentially the same as that for the case of binary values. Our approach is based on a representation of multi-valued networks using multi-valued logic functions, providing a biologically intuitive representation of the network. Furthermore, it uses tools to compute lattice points in rational polytopes, tapping a rich area of algebraic combinatorics as a source for combinatorial algorithms for Boolean network analysis. An implementation of the algorithm is provided.
△ Less
Submitted 25 April, 2024;
originally announced April 2024.
-
Absolute concentration robustness: Algebra and geometry
Authors:
Luis David García Puente,
Elizabeth Gross,
Heather A Harrington,
Matthew Johnston,
Nicolette Meshkat,
Mercedes Pérez Millán,
Anne Shiu
Abstract:
Motivated by the question of how biological systems maintain homeostasis in changing environments, Shinar and Feinberg introduced in 2010 the concept of absolute concentration robustness (ACR). A biochemical system exhibits ACR in some species if the steady-state value of that species does not depend on initial conditions. Thus, a system with ACR can maintain a constant level of one species even a…
▽ More
Motivated by the question of how biological systems maintain homeostasis in changing environments, Shinar and Feinberg introduced in 2010 the concept of absolute concentration robustness (ACR). A biochemical system exhibits ACR in some species if the steady-state value of that species does not depend on initial conditions. Thus, a system with ACR can maintain a constant level of one species even as the environment changes. Despite a great deal of interest in ACR in recent years, the following basic question remains open: How can we determine quickly whether a given biochemical system has ACR? Although various approaches to this problem have been proposed, we show that they are incomplete. Accordingly, we present new methods for deciding ACR, which harness computational algebra. We illustrate our results on several biochemical signaling networks.
△ Less
Submitted 3 October, 2024; v1 submitted 29 December, 2023;
originally announced January 2024.
-
Multistationarity questions in reduced vs extended biochemical networks
Authors:
Alicia Dickenstein,
Magalí Giaroli,
Mercedes Pérez Millán,
Rick Rischter
Abstract:
We address several questions in reduced versus extended networks via the elimination or addition of intermediate complexes in the framework of chemical reaction networks with mass-action kinetics. We clarify and extend advances in the literature concerning multistationarity in this context. We establish general results about MESSI systems, which we use to compute the circuits of multistationarity…
▽ More
We address several questions in reduced versus extended networks via the elimination or addition of intermediate complexes in the framework of chemical reaction networks with mass-action kinetics. We clarify and extend advances in the literature concerning multistationarity in this context. We establish general results about MESSI systems, which we use to compute the circuits of multistationarity for significant biochemical networks.
△ Less
Submitted 12 October, 2023; v1 submitted 3 October, 2023;
originally announced October 2023.
-
Parameter regions that give rise to 2[n/2]+1 positive steady states in the n-site phosphorylation system
Authors:
Magalí Giaroli,
Rick Rischter,
Mercedes Pérez Millán,
Alicia Dickenstein
Abstract:
The distributive sequential n-site phosphorylation/dephosphorylation system is an important building block in networks of chemical reactions arising in molecular biology, which has been intensively studied. In the nice paper of Wang and Sontag (2008) it is shown that for certain choices of the reaction rate constants and total conservation constants, the system can have 2[n/2]+1 positive steady st…
▽ More
The distributive sequential n-site phosphorylation/dephosphorylation system is an important building block in networks of chemical reactions arising in molecular biology, which has been intensively studied. In the nice paper of Wang and Sontag (2008) it is shown that for certain choices of the reaction rate constants and total conservation constants, the system can have 2[n/2]+1 positive steady states (that is, n+1 positive steady states for n even and n positive steady states for n odd). In this paper we give open parameter regions in the space of reaction rate constants and total conservation constants that ensure these number of positive steady states, while assuming in the modeling that roughly only 1/4 of the intermediates occur in the reaction mechanism. This result is based on the general framework developed by Bihan, Dickenstein, and Giaroli (2018), which can be applied to other networks. We also describe how to implement these tools to search for multistationarity regions in a computer algebra system and present some computer aided results.
△ Less
Submitted 25 April, 2019;
originally announced April 2019.
-
Identifiability from a few species for a class of biochemical reaction networks
Authors:
Gabriela Jeronimo,
Mercedes Pérez Millán,
Pablo Solernó
Abstract:
Under mass-action kinetics, biochemical reaction networks give rise to polynomial autonomous dynamical systems whose parameters are often difficult to estimate. We deal in this paper with the problem of identifying the kinetic parameters of a class of biochemical networks which are abundant, such as multisite phosphorylation systems and phosphorylation cascades (for example, MAPK cascades). For an…
▽ More
Under mass-action kinetics, biochemical reaction networks give rise to polynomial autonomous dynamical systems whose parameters are often difficult to estimate. We deal in this paper with the problem of identifying the kinetic parameters of a class of biochemical networks which are abundant, such as multisite phosphorylation systems and phosphorylation cascades (for example, MAPK cascades). For any system of this class we explicitly exhibit a single species for each connected component of the associated digraph such that the successive total derivatives of its concentration allow us to identify all the parameters occurring in the component. The number of derivatives needed is bounded essentially by the length of the corresponding connected component of the digraph. Moreover, in the particular case of the cascades, we show that the parameters can be identified from a bounded number of successive derivatives of the last product of the last layer. This theoretical result induces also a heuristic interpolation-based identifiability procedure to recover the values of the rate constants from exact measurements.
△ Less
Submitted 25 April, 2019; v1 submitted 20 November, 2018;
originally announced November 2018.
-
Multistationarity in Structured Reaction Networks
Authors:
Alicia Dickenstein,
Mercedes Perez Millan,
Anne Shiu,
Xiaoxian Tang
Abstract:
Many dynamical systems arising in biology and other areas exhibit multistationarity (two or more positive steady states with the same conserved quantities). Although deciding multistationarity for a polynomial dynamical system is an effective question in real algebraic geometry, it is in general difficult to determine whether a given network can give rise to a multistationary system, and if so, to…
▽ More
Many dynamical systems arising in biology and other areas exhibit multistationarity (two or more positive steady states with the same conserved quantities). Although deciding multistationarity for a polynomial dynamical system is an effective question in real algebraic geometry, it is in general difficult to determine whether a given network can give rise to a multistationary system, and if so, to identify witnesses to multistationarity, that is, specific parameter values for which the system exhibits multiple steady states. Here we investigate both problems. First, we build on work of Conradi, Feliu, Mincheva, and Wiuf, who showed that for certain reaction networks whose steady states admit a positive parametrization, multistationarity is characterized by whether a certain "critical function" changes sign. Here, we allow for more general parametrizations, which make it much easier to determine the existence of a sign change. This is particularly simple when the steady-state equations are linearly equivalent to binomials; we give necessary conditions for this to happen, which hold for many networks studied in the literature. We also give a sufficient condition for multistationarity of networks whose steady-state equations can be replaced by equivalent triangular-form equations. Finally, we present methods for finding witnesses to multistationarity, which we show work well for certain structured reaction networks, including those common to biological signaling pathways. Our work relies on results from degree theory, on the existence of explicit rational parametrizations of the steady states, and on the specialization of Groebner bases.
△ Less
Submitted 7 February, 2019; v1 submitted 12 October, 2018;
originally announced October 2018.
-
The structure of MESSI biological systems
Authors:
Mercedes Pérez Millán,
Alicia Dickenstein
Abstract:
We introduce a general framework for biological systems, called MESSI systems, that describe Modifications of type Enzyme-Substrate or Swap with Intermediates, and we prove general results based on the network structure. Many post-translational modification networks are MESSI systems. For example: the motifs in [Feliu and Wiuf (2012a)], sequential distributive and processive multisite phosphorylat…
▽ More
We introduce a general framework for biological systems, called MESSI systems, that describe Modifications of type Enzyme-Substrate or Swap with Intermediates, and we prove general results based on the network structure. Many post-translational modification networks are MESSI systems. For example: the motifs in [Feliu and Wiuf (2012a)], sequential distributive and processive multisite phosphorylation networks, most of the examples in [Angeli et al. (2007)], phosphorylation cascades, two component systems as in [Kothamachu et al. (2015)], the bacterial EnvZ/OmpR network in [Shinar and Feinberg (2010)], and all linear networks. We show that, under mass-action kinetics, MESSI systems are conservative. We simplify the study of steady states of these systems by explicit elimination of intermediate complexes and we give conditions to ensure an explicit rational parametrization of the variety of steady states (inspired by [Feliu and Wiuf (2013a, 2013b), Thomson and Gunawardena (2009)]). We define an important subclass of MESSI systems with toric steady states [Pérez Millán et al. (2012)] and we give for MESSI systems with toric steady states an easy algorithm to determine the capacity for multistationarity. In this case, the algorithm provides rate constants for which multistationarity takes place, based on the theory of oriented matroids.
△ Less
Submitted 5 May, 2018; v1 submitted 27 December, 2016;
originally announced December 2016.
-
MAPK's networks and their capacity for multistationarity due to toric steady states
Authors:
Mercedes Pérez Millán,
Adrián G. Turjanski
Abstract:
Mitogen-activated protein kinase (MAPK) signaling pathways play an essential role in the transduction of environmental stimuli to the nucleus, thereby regulating a variety of cellular processes, including cell proliferation, differentiation and programmed cell death. The components of the MAPK extracellular activated protein kinase (ERK) cascade represent attractive targets for cancer therapy as t…
▽ More
Mitogen-activated protein kinase (MAPK) signaling pathways play an essential role in the transduction of environmental stimuli to the nucleus, thereby regulating a variety of cellular processes, including cell proliferation, differentiation and programmed cell death. The components of the MAPK extracellular activated protein kinase (ERK) cascade represent attractive targets for cancer therapy as their aberrant activation is a frequent event among highly prevalent human cancers. MAPK networks are a model for computational simulation, mostly using Ordinary and Partial Differential Equations. Key results showed that these networks can have switch-like behavior, bistability and oscillations. In this work, we consider three representative ERK networks, one with a negative feedback loop, which present a binomial steady state ideal under mass-action kinetics. We therefore apply the theoretical result present in Pérez Millán et. al (2012) to find a set of rate constants that allow two significantly different stable steady states in the same stoichiometric compatibility class for each network. Our approach makes it possible to study certain aspects of the system, such as multistationarity, without relying on simulation, since we do not assume a priori any constant but the topology of the network. As the performed analysis is general it could be applied to many other important biochemical networks.
△ Less
Submitted 26 March, 2014;
originally announced March 2014.
-
Implicit dose-response curves
Authors:
Mercedes Pérez Millán,
Alicia Dickenstein
Abstract:
We develop tools from computational algebraic geometry for the study of steady state features of autonomous polynomial dynamical systems via elimination of variables. In particular, we obtain nontrivial bounds for the steady state concentration of a given species in biochemical reaction networks with mass-action kinetics. This species is understood as the output of the network and we thus bound th…
▽ More
We develop tools from computational algebraic geometry for the study of steady state features of autonomous polynomial dynamical systems via elimination of variables. In particular, we obtain nontrivial bounds for the steady state concentration of a given species in biochemical reaction networks with mass-action kinetics. This species is understood as the output of the network and we thus bound the maximal response of the system. The improved bounds give smaller starting boxes to launch numerical methods. We apply our results to the sequential enzymatic network studied in Markevich et al.(2004) to find nontrivial upper bounds for the different substrate concentrations at steady state.
Our approach does not require any simulation, analytical expression to describe the output in terms of the input, or the absence of multistationarity. Instead, we show how to extract information from effectively computable implicit dose-response curves with the use of resultants and discriminants. We moreover illustrate in the application to an enzymatic network, the relation between the exact implicit dose-response curve we obtain symbolically and the standard hysteresis diagram provided by a numerical solver.
The setting and tools we propose could yield many other results adapted to any autonomous polynomial dynamical system, beyond those where it is possible to get explicit expressions.
△ Less
Submitted 11 July, 2014; v1 submitted 30 January, 2014;
originally announced January 2014.
-
Complex-linear invariants of biochemical networks
Authors:
R. L. Karp,
M. Pérez Millán,
T. Dasgupta,
A. Dickenstein,
J. Gunawardena
Abstract:
The nonlinearities found in molecular networks usually prevent mathematical analysis of network behaviour, which has largely been studied by numerical simulation. This can lead to difficult problems of parameter determination. However, molecular networks give rise, through mass-action kinetics, to polynomial dynamical systems, whose steady states are zeros of a set of polynomial equations. These e…
▽ More
The nonlinearities found in molecular networks usually prevent mathematical analysis of network behaviour, which has largely been studied by numerical simulation. This can lead to difficult problems of parameter determination. However, molecular networks give rise, through mass-action kinetics, to polynomial dynamical systems, whose steady states are zeros of a set of polynomial equations. These equations may be analysed by algebraic methods, in which parameters are treated as symbolic expressions whose numerical values do not have to be known in advance. For instance, an "invariant" of a network is a polynomial expression on selected state variables that vanishes in any steady state. Invariants have been found that encode key network properties and that discriminate between different network structures. Although invariants may be calculated by computational algebraic methods, such as Gröbner bases, these become computationally infeasible for biologically realistic networks. Here, we exploit Chemical Reaction Network Theory (CRNT) to develop an efficient procedure for calculating invariants that are linear combinations of "complexes", or the monomials coming from mass action. We show how this procedure can be used in proving earlier results of Horn and Jackson and of Shinar and Feinberg for networks of deficiency at most one. We then apply our method to enzyme bifunctionality, including the bacterial EnvZ/OmpR osmolarity regulator and the mammalian 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase glycolytic regulator, whose networks have deficiencies up to four. We show that bifunctionality leads to different forms of concentration control that are robust to changes in initial conditions or total amounts. Finally, we outline a systematic procedure for using complex-linear invariants to analyse molecular networks of any deficiency.
△ Less
Submitted 14 July, 2012;
originally announced July 2012.
-
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.
-
How far is complex balancing from detailed balancing?
Authors:
Alicia Dickenstein,
Mercedes Perez Millan
Abstract:
The aim of this article is to build on the use of tools from computational algebra initiated in Craciun, Dickenstein, Shiu, Sturmfels (JSC, 2009), for the study of general kinetic systems, which have a wide range of applications in chemistry and biology. We clarify the relation between the algebraic conditions that must be satisfied by the reaction constants in general (mass action) kinetics sys…
▽ More
The aim of this article is to build on the use of tools from computational algebra initiated in Craciun, Dickenstein, Shiu, Sturmfels (JSC, 2009), for the study of general kinetic systems, which have a wide range of applications in chemistry and biology. We clarify the relation between the algebraic conditions that must be satisfied by the reaction constants in general (mass action) kinetics systems for the existence of detailed or complex balancing equilibria. The main properties of these systems have been set by Horn, Jackson and Feinberg.
△ Less
Submitted 6 January, 2010;
originally announced January 2010.