Eshan D Mitra1, William S Hlavacek1. 1. Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545, USA.
Abstract
MOTIVATION: Recent work has demonstrated the feasibility of using non-numerical, qualitative data to parameterize mathematical models. However, uncertainty quantification (UQ) of such parameterized models has remained challenging because of a lack of a statistical interpretation of the objective functions used in optimization. RESULTS: We formulated likelihood functions suitable for performing Bayesian UQ using qualitative observations of underlying continuous variables or a combination of qualitative and quantitative data. To demonstrate the resulting UQ capabilities, we analyzed a published model for immunoglobulin E (IgE) receptor signaling using synthetic qualitative and quantitative datasets. Remarkably, estimates of parameter values derived from the qualitative data were nearly as consistent with the assumed ground-truth parameter values as estimates derived from the lower throughput quantitative data. These results provide further motivation for leveraging qualitative data in biological modeling. AVAILABILITY AND IMPLEMENTATION: The likelihood functions presented here are implemented in a new release of PyBioNetFit, an open-source application for analyzing Systems Biology Markup Language- and BioNetGen Language-formatted models, available online at www.github.com/lanl/PyBNF. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
MOTIVATION: Recent work has demonstrated the feasibility of using non-numerical, qualitative data to parameterize mathematical models. However, uncertainty quantification (UQ) of such parameterized models has remained challenging because of a lack of a statistical interpretation of the objective functions used in optimization. RESULTS: We formulated likelihood functions suitable for performing Bayesian UQ using qualitative observations of underlying continuous variables or a combination of qualitative and quantitative data. To demonstrate the resulting UQ capabilities, we analyzed a published model for immunoglobulin E (IgE) receptor signaling using synthetic qualitative and quantitative datasets. Remarkably, estimates of parameter values derived from the qualitative data were nearly as consistent with the assumed ground-truth parameter values as estimates derived from the lower throughput quantitative data. These results provide further motivation for leveraging qualitative data in biological modeling. AVAILABILITY AND IMPLEMENTATION: The likelihood functions presented here are implemented in a new release of PyBioNetFit, an open-source application for analyzing Systems Biology Markup Language- and BioNetGen Language-formatted models, available online at www.github.com/lanl/PyBNF. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Authors: Katherine C Chen; Laurence Calzone; Attila Csikasz-Nagy; Frederick R Cross; Bela Novak; John J Tyson Journal: Mol Biol Cell Date: 2004-05-28 Impact factor: 4.138
Authors: William S Hlavacek; Jennifer A Csicsery-Ronay; Lewis R Baker; María Del Carmen Ramos Álamo; Alexander Ionkov; Eshan D Mitra; Ryan Suderman; Keesha E Erickson; Raquel Dias; Joshua Colvin; Brandon R Thomas; Richard G Posner Journal: Methods Mol Biol Date: 2019
Authors: Noga Kozer; Dipak Barua; Suzanne Orchard; Eduoard C Nice; Antony W Burgess; William S Hlavacek; Andrew H A Clayton Journal: Mol Biosyst Date: 2013-04-30
Authors: Minghan Chen; Brandon D Amos; Layne T Watson; John J Tyson; Young Cao; Clifford A Shaffer; Michael W Trosset; Cihan Oguz; Gisella Kakoti Journal: IEEE/ACM Trans Comput Biol Bioinform Date: 2017-11-14 Impact factor: 3.710
Authors: Neil R Adames; P Logan Schuck; Katherine C Chen; T M Murali; John J Tyson; Jean Peccoud Journal: Mol Biol Cell Date: 2015-08-26 Impact factor: 4.138