Multistationarity questions in reduced versus extended biochemical networks.

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, mainly from Feliu and Wiuf (J R Soc Interface 10:20130484, 2013), Sadeghimanesh and Feliu (Bull Math Biol 81:2428-2462, 2019), Pérez Millán and Dickenstein (SIAM J Appl Dyn Syst 17(2):1650-1682, 2018), Dickenstein et al. (Bull Math Biol 81:1527-1581, 2019).

Identifiability from a Few Species for a Class of Biochemical Reaction Networks.

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.

Multistationarity in Structured Reaction Networks.

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.

MAPK's networks and their capacity for multistationarity due to toric steady states.

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.

Implicit dose-response curves.

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.

Complex-linear invariants of biochemical networks.

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.

Chemical reaction systems with toric steady states.

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 parameterized (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.

How far is complex balancing from detailed balancing?

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. These systems have a wide range of applications in chemistry and biology. Their main properties have been set by Horn, Jackson and Feinberg.