Graphical Requirements for Multistationarity in Reaction Networks and their Verification in BioModels

Thomas' necessary conditions for the existence of multiple steady states in gene networks have been proved by Soulé with high generality for dynamical systems defined by differential equations. When applied to (protein) reaction networks however, those conditions do not provide information since they are trivially satisfied as soon as there is a bimolecular or a reversible reaction. Refined graphical requirements have been proposed to deal with such cases. In this paper, we present for the first time a graph rewriting algorithm for checking the refined conditions given by Soliman, and evaluate its practical performance by applying it systematically to the curated branch of the BioModels repository. This algorithm analyzes all reaction networks (of size up to 430 species) in less than 0.05 second per network, and permits to conclude to the absence of multistationarity in 160 networks over 506. The short computation times obtained in this graphical approach are in sharp contrast to the Jacobian-based symbolic computation approach. We also discuss the case of one extra graphical condition by arc rewiring that allows us to conclude on 20 more networks of this benchmark but with a high computational cost. Finally, we study with some details the case of phosphorylation cycles and MAPK signalling models which show the importance of modelling the intermediate complexations with the enzymes in order to correctly analyze the multistationarity capabilities of such biochemical reaction networks.

[1]  B. Kholodenko,et al.  Signaling switches and bistability arising from multisite phosphorylation in protein kinase cascades , 2004, The Journal of cell biology.

[2]  Alejandra C. Ventura,et al.  A Hidden Feedback in Signaling Cascades Is Revealed , 2008, PLoS Comput. Biol..

[3]  Panos M. Pardalos,et al.  A Greedy Randomized Adaptive Search Procedure for the Feedback Vertex Set Problem , 1998, J. Comb. Optim..

[4]  Liang Qiao,et al.  Bistability and Oscillations in the Huang-Ferrell Model of MAPK Signaling , 2007, PLoS Comput. Biol..

[5]  Nicolas Le Novère,et al.  BioModels Database: a repository of mathematical models of biological processes. , 2013, Methods in molecular biology.

[6]  Martin Feinberg,et al.  Multiple Equilibria in Complex Chemical Reaction Networks: Ii. the Species-reactions Graph , 2022 .

[7]  B. Kholodenko,et al.  Negative feedback and ultrasensitivity can bring about oscillations in the mitogen-activated protein kinase cascades. , 2000, European journal of biochemistry.

[8]  Sylvain Soliman,et al.  A Stronger Necessary Condition for the Multistationarity of Chemical Reaction Networks , 2013, Bulletin of mathematical biology.

[9]  François Fages,et al.  Abstract interpretation and types for systems biology , 2008, Theor. Comput. Sci..

[10]  Sylvain Soliman,et al.  Influence Networks Compared with Reaction Networks: Semantics, Expressivity and Attractors , 2018, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[11]  François Fages,et al.  BIOCHAM: an environment for modeling biological systems and formalizing experimental knowledge , 2006, Bioinform..

[12]  Elisenda Feliu,et al.  A computational method to preclude multistationarity in networks of interacting species , 2013, Bioinform..

[13]  Robert E. Tarjan,et al.  Depth-First Search and Linear Graph Algorithms , 1972, SIAM J. Comput..

[14]  René Thomas.  On the Relation Between the Logical Structure of Systems and Their Ability to Generate Multiple Steady States or Sustained Oscillations , 1981 .

[15]  François Fages,et al.  From Reaction Models to Influence Graphs and Back: A Theorem , 2008, FMSB.

[16]  Elisenda Feliu,et al.  Simplifying biochemical models with intermediate species , 2013, Journal of The Royal Society Interface.

[17]  Hans-Michael Kaltenbach,et al.  A unified view on bipartite species-reaction and interaction graphs for chemical reaction networks , 2012, 1210.0320.

[18]  C. Soulé.  Graphic Requirements for Multistationarity , 2004, Complexus.

[19]  Murad Banaji,et al.  Graph-theoretic criteria for injectivity and unique equilibria in general chemical reaction systems , 2008, Adv. Appl. Math..

[20]  Hanoch Levy,et al.  A Contraction Algorithm for Finding Small Cycle Cutsets , 1988, J. Algorithms.

[21]  El Houssine Snoussi.  Necessary Conditions for Multistationarity and Stable Periodicity , 1998 .

[22]  Sylvain Soliman,et al.  Invariants and Other Structural Properties of Biochemical Models as a Constraint Satisfaction Problem , 2012, Algorithms for Molecular Biology.

[23]  Ursula Klingmüller,et al.  Theoretical and experimental analysis links isoform- specific ERK signalling to cell fate decisions , 2009, Molecular systems biology.

[24]  François Fages,et al.  A constraint programming approach to cutset problems , 2006, Comput. Oper. Res..

[25]  François Fages,et al.  Inferring reaction systems from ordinary differential equations , 2015, Theor. Comput. Sci..

[26]  Murad Banaji,et al.  Graph-theoretic approaches to injectivity and multiple equilibria in systems of interacting elements , 2009, 0903.1190.

[27]  François Fages,et al.  A graphical method for reducing and relating models in systems biology , 2010, Bioinform..

[28]  François Fages,et al.  On enumerating minimal siphons in Petri nets using CLP and SAT solvers: theoretical and practical complexity , 2016, Constraints.

[29]  Michael Hucka,et al.  Systems Biology Markup Language (SBML) Level 2: Structures and Facilities for Model Definitions , 2007, WWW 2007.

[30]  Eduardo Sontag,et al.  On the number of steady states in a multiple futile cycle , 2008, Journal of mathematical biology.