Approximation of dynamical systems using s-systems theory: application to biological systems

In this article we propose a new symbolic-numeric algorithm to find positive equilibria of a n-dimensional dynamical system. This algorithm uses a symbolic manipulation of ODE in order to give a local approximation of differential equations with power-law dynamics (S-systems). A numerical calculus is then performed to converge towards an equilibrium, giving at the same time a S-system approximating the initial system around this equilibrium. This algorithm has been applied to a real biological example in 14 dimensions which is a subsystem of a metabolic pathway in Arabidopsis Thaliana.

[1]  H. D. Jong,et al.  Qualitative simulation of genetic regulatory networks using piecewise-linear models , 2004, Bulletin of mathematical biology.

[2]  Michael A Savageau,et al.  Alternative designs for a genetic switch: analysis of switching times using the piecewise power-law representation. , 2002, Mathematical biosciences.

[3]  Reinhard C. Laubenbacher.  A computer algebra approach to biological systems , 2003, ISSAC '03.

[4]  V. Klee,et al.  When is a Matrix Sign Stable? , 1977, Canadian Journal of Mathematics.

[5]  L. Glass.  Combinatorial and topological methods in nonlinear chemical kinetics , 1975 .

[6]  Alberto Policriti,et al.  XS-systems: eXtended S-Systems and Algebraic Differential Automata for Modeling Cellular Behavior , 2002, HiPC.

[7]  R. Thomas,et al.  Multistationarity, the basis of cell differentiation and memory. II. Logical analysis of regulatory networks in terms of feedback circuits. , 2001, Chaos.

[8]  Gilles Curien,et al.  A kinetic model of the branch-point between the methionine and threonine biosynthesis pathways in Arabidopsis thaliana. , 2003, European journal of biochemistry.

[9]  T. Mestl,et al.  A mathematical framework for describing and analysing gene regulatory networks. , 1995, Journal of theoretical biology.

[10]  F R Adler,et al.  How to make a biological switch. , 2000, Journal of theoretical biology.

[11]  E. Voit,et al.  Recasting nonlinear differential equations as S-systems: a canonical nonlinear form , 1987 .

[12]  James Quirk,et al.  Qualitative Problems in Matrix Theory , 1969 .

[13]  A. Goriely,et al.  Universal canonical forms for time-continuous dynamical systems. , 1989, Physical review. A, General physics.

[14]  J. Ortega.  Numerical Analysis: A Second Course , 1974 .

[15]  M A Savageau,et al.  Accuracy of alternative representations for integrated biochemical systems. , 1987, Biochemistry.

[16]  A. Goldbeter,et al.  Biochemical Oscillations And Cellular Rhythms: Contents , 1996 .

[17]  J. Yorke,et al.  Fractal Basin Boundaries, Long-Lived Chaotic Transients, And Unstable-Unstable Pair Bifurcation , 1983 .

[18]  S. Wiggins.  Introduction to Applied Nonlinear Dynamical Systems and Chaos , 1989 .