PuSH - Publication Server of Helmholtz Zentrum München

Benchmarking of numerical integration methods for ODE models of biological systems.

Sci. Rep. 11:2696 (2021)
Publ. Version/Full Text Research data DOI
Open Access Gold
Creative Commons Lizenzvertrag
Ordinary differential equation (ODE) models are a key tool to understand complex mechanisms in systems biology. These models are studied using various approaches, including stability and bifurcation analysis, but most frequently by numerical simulations. The number of required simulations is often large, e.g., when unknown parameters need to be inferred. This renders efficient and reliable numerical integration methods essential. However, these methods depend on various hyperparameters, which strongly impact the ODE solution. Despite this, and although hundreds of published ODE models are freely available in public databases, a thorough study that quantifies the impact of hyperparameters on the ODE solver in terms of accuracy and computation time is still missing. In this manuscript, we investigate which choices of algorithms and hyperparameters are generally favorable when dealing with ODE models arising from biological processes. To ensure a representative evaluation, we considered 142 published models. Our study provides evidence that most ODEs in computational biology are stiff, and we give guidelines for the choice of algorithms and hyperparameters. We anticipate that our results will help researchers in systems biology to choose appropriate numerical methods when dealing with ODE models.
Additional Metrics?
Edit extra informations Login
Publication type Article: Journal article
Document type Scientific Article
ISSN (print) / ISBN 2045-2322
e-ISSN 2045-2322
Quellenangaben Volume: 11, Issue: 1, Pages: , Article Number: 2696 Supplement: ,
Publisher Nature Publishing Group
Publishing Place London
Reviewing status Peer reviewed
Grants Clusters of Excellence, Germany
Deutsche Forschungsgemeinschaft