Bayesian uncertainty analysis for complex systems biology models: emulation, global parameter searches and evaluation of gene functions

Abstract

BackgroundMany mathematical models have now been employed across every area of systems biology. These models increasingly involve large numbers of unknown parameters, have complex structure which can result in substantial evaluation time relative to the needs of the analysis, and need to be compared to observed data of various forms. The correct analysis of such models usually requires a global parameter search, over a high dimensional parameter space, that incorporates and respects the most important sources of uncertainty. This can be an extremely difficult task, but it is essential for any meaningful inference or prediction to be made about any biological system. It hence represents a fundamental challenge for the whole of systems biology.MethodsBayesian statistical methodology for the uncertainty analysis of complex models is introduced, which is designed to address the high dimensional global parameter search problem. Bayesian emulators that mimic the systems biology model but which are extremely fast to evaluate are embeded within an iterative history match: an efficient method to search high dimensional spaces within a more formal statistical setting, while incorporating major sources of uncertainty.ResultsThe approach is demonstrated via application to a model of hormonal crosstalk in Arabidopsis root development, which has 32 rate parameters, for which we identify the sets of rate parameter values that lead to acceptable matches between model output and observed trend data. The multiple insights into the model’s structure that this analysis provides are discussed. The methodology is applied to a second related model, and the biological consequences of the resulting comparison, including the evaluation of gene functions, are described.ConclusionsBayesian uncertainty analysis for complex models using both emulators and history matching is shown to be a powerful technique that can greatly aid the study of a large class of systems biology models. It both provides insight into model behaviour and identifies the sets of rate parameters of interest.Electronic supplementary materialThe online version of this article (doi:10.1186/s12918-017-0484-3) contains supplementary material, which is available to authorized users.

Model dispersion with PRISM; an alternative to MCMC for rapid analysis of models

“…Note that the adjusted values for any known model evaluation are equal to E D i ( f i (x (k) )) = f i (x (k) ) and Var D i ( f i (x (k) )) = 0. For further discussion, we point interested readers to Vernon et al (2010); Vernon et al (2018).…”

“…The process of history matching can be compared to model calibration (discussed in more detail in Vernon et al 2010;Vernon et al 2018), where we assume that there is a single true but unknown parameter set x * and our goal is to make probabilistic statements about x * based on a prior specification, the collection of model evaluations and the observed history (Kennedy & O'Hagan 2001;Goldstein & Rougier 2006). Although history matching and model calibration look alike and are certainly related, they are very different in terms of approach.…”

Bayesian Correlation is a robust similarity measure for single cell RNA-seq data

“…In order to address limitations dependent on the noise, Bayesian statistics have been used to study biological processes [75,74]. Recently, we have proposed a Bayesian correlation scheme to asses similarity between two entities in high-throughput sequencing (HTS) experiments [60,56].…”

How to deal with parameters for whole-cell modelling

“…It is therefore important to consider both structural and parametric uncertainty in every biological modelling study. Scaling such approaches up to WCMs will be a technical challenge and require the development of suitable approximations, or potentially emulation or surrogate modelling approaches [28,103,104].…”