Identifying the parametric occurrence of multiple steady states for some biological networks

Russell Bradford, James H. Davenport, Matthew England, Hassan Errami, Vladimir Gerdt, Dima Grigoriev, Charles Hoyt, Marek Košta, Ovidiu Radulescu, Thomas Sturm, Andreas Weber

Research output: Contribution to journalArticle

Abstract

We consider a problem from biological network analysis of determining regions in a parameter space over which there are multiple steady states for positive real values of variables and parameters. We describe multiple approaches to address the problem using tools from Symbolic Computation. We describe how progress was made to achieve semi-algebraic descriptions of the multistationarity regions of parameter space, and compare symbolic and numerical methods. The biological networks studied are models of the mitogen-activated protein kinases (MAPK) network which has already consumed considerable effort using special insights into its structure of corresponding models. Our main example is a model with 11 equations in 11 variables and 19 parameters, 3 of which are of interest for symbolic treatment. The model also imposes positivity conditions on all variables and parameters. We apply combinations of symbolic computation methods designed for mixed equality / inequality systems, specifically virtual substitution, lazy real triangularization and cylindrical algebraic decomposition, as well as a simplification technique adapted from Gaussian elimination and graph theory. We are able to determine semi-algebraic conditions for multistationarity of our main example over a 2-dimensional parameter space. We also study a second MAPK model and a symbolic grid sampling technique which can locate such regions in 3-dimensional parameter space.

Original languageEnglish
Pages (from-to)84-119
Number of pages36
JournalJournal of Symbolic Computation
Volume98
Early online date15 Jul 2019
DOIs
Publication statusE-pub ahead of print - 15 Jul 2019

Keywords

  • Biological networks
  • MAPK
  • Mixed equation/inequality solving
  • Real quantifier elimination
  • Signalling pathways

ASJC Scopus subject areas

  • Algebra and Number Theory
  • Computational Mathematics

Cite this

Identifying the parametric occurrence of multiple steady states for some biological networks. / Bradford, Russell; Davenport, James H.; England, Matthew; Errami, Hassan; Gerdt, Vladimir; Grigoriev, Dima; Hoyt, Charles; Košta, Marek; Radulescu, Ovidiu; Sturm, Thomas; Weber, Andreas.

In: Journal of Symbolic Computation, Vol. 98, 01.05.2020, p. 84-119.

Research output: Contribution to journalArticle

Bradford, R, Davenport, JH, England, M, Errami, H, Gerdt, V, Grigoriev, D, Hoyt, C, Košta, M, Radulescu, O, Sturm, T & Weber, A 2020, 'Identifying the parametric occurrence of multiple steady states for some biological networks', Journal of Symbolic Computation, vol. 98, pp. 84-119. https://doi.org/10.1016/j.jsc.2019.07.008
Bradford, Russell ; Davenport, James H. ; England, Matthew ; Errami, Hassan ; Gerdt, Vladimir ; Grigoriev, Dima ; Hoyt, Charles ; Košta, Marek ; Radulescu, Ovidiu ; Sturm, Thomas ; Weber, Andreas. / Identifying the parametric occurrence of multiple steady states for some biological networks. In: Journal of Symbolic Computation. 2020 ; Vol. 98. pp. 84-119.
@article{c2dbf0ec5df94a81abcd4e6fd0727ba0,
title = "Identifying the parametric occurrence of multiple steady states for some biological networks",
abstract = "We consider a problem from biological network analysis of determining regions in a parameter space over which there are multiple steady states for positive real values of variables and parameters. We describe multiple approaches to address the problem using tools from Symbolic Computation. We describe how progress was made to achieve semi-algebraic descriptions of the multistationarity regions of parameter space, and compare symbolic and numerical methods. The biological networks studied are models of the mitogen-activated protein kinases (MAPK) network which has already consumed considerable effort using special insights into its structure of corresponding models. Our main example is a model with 11 equations in 11 variables and 19 parameters, 3 of which are of interest for symbolic treatment. The model also imposes positivity conditions on all variables and parameters. We apply combinations of symbolic computation methods designed for mixed equality / inequality systems, specifically virtual substitution, lazy real triangularization and cylindrical algebraic decomposition, as well as a simplification technique adapted from Gaussian elimination and graph theory. We are able to determine semi-algebraic conditions for multistationarity of our main example over a 2-dimensional parameter space. We also study a second MAPK model and a symbolic grid sampling technique which can locate such regions in 3-dimensional parameter space.",
keywords = "Biological networks, MAPK, Mixed equation/inequality solving, Real quantifier elimination, Signalling pathways",
author = "Russell Bradford and Davenport, {James H.} and Matthew England and Hassan Errami and Vladimir Gerdt and Dima Grigoriev and Charles Hoyt and Marek Košta and Ovidiu Radulescu and Thomas Sturm and Andreas Weber",
year = "2019",
month = "7",
day = "15",
doi = "10.1016/j.jsc.2019.07.008",
language = "English",
volume = "98",
pages = "84--119",
journal = "Journal of Symbolic Computation",
issn = "0747-7171",
publisher = "Elsevier Academic Press Inc",

}

TY - JOUR

T1 - Identifying the parametric occurrence of multiple steady states for some biological networks

AU - Bradford, Russell

AU - Davenport, James H.

AU - England, Matthew

AU - Errami, Hassan

AU - Gerdt, Vladimir

AU - Grigoriev, Dima

AU - Hoyt, Charles

AU - Košta, Marek

AU - Radulescu, Ovidiu

AU - Sturm, Thomas

AU - Weber, Andreas

PY - 2019/7/15

Y1 - 2019/7/15

N2 - We consider a problem from biological network analysis of determining regions in a parameter space over which there are multiple steady states for positive real values of variables and parameters. We describe multiple approaches to address the problem using tools from Symbolic Computation. We describe how progress was made to achieve semi-algebraic descriptions of the multistationarity regions of parameter space, and compare symbolic and numerical methods. The biological networks studied are models of the mitogen-activated protein kinases (MAPK) network which has already consumed considerable effort using special insights into its structure of corresponding models. Our main example is a model with 11 equations in 11 variables and 19 parameters, 3 of which are of interest for symbolic treatment. The model also imposes positivity conditions on all variables and parameters. We apply combinations of symbolic computation methods designed for mixed equality / inequality systems, specifically virtual substitution, lazy real triangularization and cylindrical algebraic decomposition, as well as a simplification technique adapted from Gaussian elimination and graph theory. We are able to determine semi-algebraic conditions for multistationarity of our main example over a 2-dimensional parameter space. We also study a second MAPK model and a symbolic grid sampling technique which can locate such regions in 3-dimensional parameter space.

AB - We consider a problem from biological network analysis of determining regions in a parameter space over which there are multiple steady states for positive real values of variables and parameters. We describe multiple approaches to address the problem using tools from Symbolic Computation. We describe how progress was made to achieve semi-algebraic descriptions of the multistationarity regions of parameter space, and compare symbolic and numerical methods. The biological networks studied are models of the mitogen-activated protein kinases (MAPK) network which has already consumed considerable effort using special insights into its structure of corresponding models. Our main example is a model with 11 equations in 11 variables and 19 parameters, 3 of which are of interest for symbolic treatment. The model also imposes positivity conditions on all variables and parameters. We apply combinations of symbolic computation methods designed for mixed equality / inequality systems, specifically virtual substitution, lazy real triangularization and cylindrical algebraic decomposition, as well as a simplification technique adapted from Gaussian elimination and graph theory. We are able to determine semi-algebraic conditions for multistationarity of our main example over a 2-dimensional parameter space. We also study a second MAPK model and a symbolic grid sampling technique which can locate such regions in 3-dimensional parameter space.

KW - Biological networks

KW - MAPK

KW - Mixed equation/inequality solving

KW - Real quantifier elimination

KW - Signalling pathways

UR - http://www.scopus.com/inward/record.url?scp=85069739042&partnerID=8YFLogxK

U2 - 10.1016/j.jsc.2019.07.008

DO - 10.1016/j.jsc.2019.07.008

M3 - Article

AN - SCOPUS:85069739042

VL - 98

SP - 84

EP - 119

JO - Journal of Symbolic Computation

JF - Journal of Symbolic Computation

SN - 0747-7171

ER -