A Numerical Method for Analyzing the Stability of Bi-Parametric Biological Systems

For a biological system modeled by a continuous dynamical system defined by rational functions with two parameters, we propose a numerical method to compute the fold and Hopf bifurcation boundaries of the system restricted in a finite region in the parametric space under certain assumptions. The bifurcation boundaries divide their complement in the region into connected subsets, called cells, such that above each of them the number of equilibria is constant and the stability of each equilibrium remains unchanged. The boundaries are generated by first tracing the fold and Hopf bifurcation curves in a higher dimensional space and then projecting them onto the parameter plane. One advantage of this method is that it can exploit global information of real varieties and generate complete boundaries based on homotopy continuation methods and critical point techniques. The bistability properties of several biological systems are successfully analyzed by our method.

[1]  Yong Feng,et al.  Computing real witness points of positive dimensional polynomial systems , 2017, Theor. Comput. Sci..

[2]  J. R. Pomerening,et al.  Uncovering mechanisms of bistability in biological systems. , 2008, Current opinion in biotechnology.

[3]  M. Gatto,et al.  Quantifying the dynamics of prion infection: a bifurcation analysis of Laurent's model. , 2000, Journal of theoretical biology.

[4]  Changbo Chen,et al.  A Numerical Method for Computing Border Curves of Bi-parametric Real Polynomial Systems and Applications , 2016, CASC.

[5]  Dongming Wang,et al.  Stability analysis of biological systems with real solution classification , 2005, ISSAC.

[6]  Wei Niu,et al.  Algebraic Approaches to Stability Analysis of Biological Systems , 2008, Math. Comput. Sci..

[7]  J R King,et al.  Bistability in a model of mesoderm and anterior mesendoderm specification in Xenopus laevis. , 2009, Journal of theoretical biology.

[8]  M Laurent,et al.  Prion diseases and the 'protein only' hypothesis: a theoretical dynamic study. , 1996, The Biochemical journal.

[9]  Eduardo Sontag,et al.  Untangling the wires: A strategy to trace functional interactions in signaling and gene networks , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[10]  Willy Govaerts,et al.  Numerical methods for bifurcations of dynamical equilibria , 1987 .

[11]  Terence Hwa,et al.  The Innate Growth Bistability and Fitness Landscapes of Antibiotic-Resistant Bacteria , 2013, Science.

[12]  J. A. Kuznecov.  Elements of applied bifurcation theory , 1998 .

[13]  James E. Ferrell,et al.  A positive-feedback-based bistable ‘memory module’ that governs a cell fate decision , 2007, Nature.

[14]  Dima Grigoriev,et al.  Detection of Hopf bifurcations in chemical reaction networks using convex coordinates , 2015, J. Comput. Phys..

[15]  Marc Moreno Maza,et al.  On Solving Parametric Polynomial Systems , 2012, Mathematics in Computer Science.

[16]  C Schimek,et al.  Model of the synthesis of trisporic acid in Mucorales showing bistability. , 2012, IET systems biology.

[17]  Bican Xia,et al.  Special algorithm for stability analysis of multistable biological regulatory systems , 2013, J. Symb. Comput..

[18]  Wenyuan Wu,et al.  Finding points on real solution components and applications to differential polynomial systems , 2013, ISSAC '13.

[19]  Fabrice Rouillier,et al.  Finding at Least One Point in Each Connected Component of a Real Algebraic Set Defined by a Single Equation , 2000, J. Complex..