Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations

Biological system's dynamics are increasingly studied with nonlinear ordinary differential equations, whose parameters are estimated from input/output experimental data. Structural identifiability analysis addresses the theoretical question whether the inverse problem of recovering the unknown...

Full description

Bibliographic Details
Main Authors: Maria Pia Saccomani, Karl Thomaseth
Format: Article
Language:English
Published: AIMS Press 2019-07-01
Series:Mathematical Biosciences and Engineering
Subjects:
Online Access:https://www.aimspress.com/article/10.3934/mbe.2019322?viewType=HTML
_version_ 1811234571330519040
author Maria Pia Saccomani
Karl Thomaseth
author_facet Maria Pia Saccomani
Karl Thomaseth
author_sort Maria Pia Saccomani
collection DOAJ
description Biological system's dynamics are increasingly studied with nonlinear ordinary differential equations, whose parameters are estimated from input/output experimental data. Structural identifiability analysis addresses the theoretical question whether the inverse problem of recovering the unknown parameters from noise-free data is uniquely solvable (global), or if there is a finite (local), or an infinite number (non identifiable) of parameter values that generate identical input/output trajectories. In contrast, practical identifiability analysis aims to assess whether the experimental data provide information on the parameter estimates in terms of precision and accuracy. A main difference between the two identifiability approaches is that the former is mostly carried out analytically and provides exact results at a cost of increased computational complexity, while the latter is usually numerically tested by calculating statistical confidence regions and relies on decision thresholds. Here we focus on local identifiability, a critical issue in biological modeling. This is the case when a model has multiple parameter solutions which equivalently describe the input/output data, but predict different behaviours of the unmeasured variables, often those of major interest. We present theoretical background and applications to locally identifiable ODE models described by rational functions. We show how structural identifiability analysis completes the practical identifiability results. In particular we propose an algorithmic approach, implemented with our software DAISY, to calculate all numerical parameter solutions and to predict the corresponding behaviour of the unmeasured variables, which otherwise would remain hidden. A case study of a locally identifiable HIV model shows that one should be aware of the presence of multiple parameter solutions to comprehensively describe the biological system and avoid biological misinterpretation of the results.
first_indexed 2024-04-12T11:39:43Z
format Article
id doaj.art-5e0bb63bee3647d3903b3f5e1d87ad87
institution Directory Open Access Journal
issn 1551-0018
language English
last_indexed 2024-04-12T11:39:43Z
publishDate 2019-07-01
publisher AIMS Press
record_format Article
series Mathematical Biosciences and Engineering
spelling doaj.art-5e0bb63bee3647d3903b3f5e1d87ad872022-12-22T03:34:43ZengAIMS PressMathematical Biosciences and Engineering1551-00182019-07-011666438645310.3934/mbe.2019322Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretationsMaria Pia Saccomani0Karl Thomaseth11. Department of Information Engineering, University of Padova, Padova, 35131 Italy2. National Research Council (IEIIT-CNR) c/o Department of Information Engineering, University of Padova, Padova, 35131 ItalyBiological system's dynamics are increasingly studied with nonlinear ordinary differential equations, whose parameters are estimated from input/output experimental data. Structural identifiability analysis addresses the theoretical question whether the inverse problem of recovering the unknown parameters from noise-free data is uniquely solvable (global), or if there is a finite (local), or an infinite number (non identifiable) of parameter values that generate identical input/output trajectories. In contrast, practical identifiability analysis aims to assess whether the experimental data provide information on the parameter estimates in terms of precision and accuracy. A main difference between the two identifiability approaches is that the former is mostly carried out analytically and provides exact results at a cost of increased computational complexity, while the latter is usually numerically tested by calculating statistical confidence regions and relies on decision thresholds. Here we focus on local identifiability, a critical issue in biological modeling. This is the case when a model has multiple parameter solutions which equivalently describe the input/output data, but predict different behaviours of the unmeasured variables, often those of major interest. We present theoretical background and applications to locally identifiable ODE models described by rational functions. We show how structural identifiability analysis completes the practical identifiability results. In particular we propose an algorithmic approach, implemented with our software DAISY, to calculate all numerical parameter solutions and to predict the corresponding behaviour of the unmeasured variables, which otherwise would remain hidden. A case study of a locally identifiable HIV model shows that one should be aware of the presence of multiple parameter solutions to comprehensively describe the biological system and avoid biological misinterpretation of the results.https://www.aimspress.com/article/10.3934/mbe.2019322?viewType=HTMLbiological systemslocal identifiabilitymultiple parameter solutionsdifferential equationparameter estimationhiv model
spellingShingle Maria Pia Saccomani
Karl Thomaseth
Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations
Mathematical Biosciences and Engineering
biological systems
local identifiability
multiple parameter solutions
differential equation
parameter estimation
hiv model
title Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations
title_full Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations
title_fullStr Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations
title_full_unstemmed Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations
title_short Calculating all multiple parameter solutions of ODE models to avoid biological misinterpretations
title_sort calculating all multiple parameter solutions of ode models to avoid biological misinterpretations
topic biological systems
local identifiability
multiple parameter solutions
differential equation
parameter estimation
hiv model
url https://www.aimspress.com/article/10.3934/mbe.2019322?viewType=HTML
work_keys_str_mv AT mariapiasaccomani calculatingallmultipleparametersolutionsofodemodelstoavoidbiologicalmisinterpretations
AT karlthomaseth calculatingallmultipleparametersolutionsofodemodelstoavoidbiologicalmisinterpretations