Algorithmic Reduction of Biological Networks With Multiple Time Scales
arXiv:2010.10129 · doi:10.1007/s11786-021-00515-2
Abstract
We present a symbolic algorithmic approach that allows to compute invariant manifolds and corresponding reduced systems for differential equations modeling biological networks which comprise chemical reaction networks for cellular biochemistry, and compartmental models for pharmacology, epidemiology and ecology. Multiple time scales of a given network are obtained by scaling, based on tropical geometry. Our reduction is mathematically justified within a singular perturbation setting. The existence of invariant manifolds is subject to hyperbolicity conditions, for which we propose an algorithmic test based on Hurwitz criteria. We finally obtain a sequence of nested invariant manifolds and respective reduced systems on those manifolds. Our theoretical results are generally accompanied by rigorous algorithmic descriptions suitable for direct implementation based on existing off-the-shelf software systems, specifically symbolic computation libraries and Satisfiability Modulo Theories solvers. We present computational examples taken from the well-known BioModels database using our own prototypical implementations.
References in corpus (11)
- Computing Tropical Varieties
- Signal processing in the TGF-beta superfamily ligand-receptor network
- Reduction of dynamical biochemical reaction networks in computational biology
- A Case Study on the Parametric Occurrence of Multiple Steady States
- Identifying the Parametric Occurrence of Multiple Steady States for some Biological Networks
- Symbolic Versus Numerical Computation and Visualization of Parameter Regions for Multistationarity of Biological Networks
- Coordinate-independent singular perturbation reduction for systems with three time scales
- Singular perturbations and scaling
- Efficiently and Effectively Recognizing Toricity of Steady State Varieties
- Computing Tropical Prevarieties with Satisfiability Modulo Theories (SMT) Solvers
- A Logic Based Approach to Finding Real Singularities of Implicit Ordinary Differential Equations