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. 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 Gröbner bases.