Dataset Information

Hierarchical Bayesian inference for ion channel screening dose-response data.

ABSTRACT: Dose-response (or 'concentration-effect') relationships commonly occur in biological and pharmacological systems and are well characterised by Hill curves. These curves are described by an equation with two parameters: the inhibitory concentration 50% (IC50); and the Hill coefficient. Typically just the 'best fit' parameter values are reported in the literature. Here we introduce a Python-based software tool, PyHillFit , and describe the underlying Bayesian inference methods that it uses, to infer probability distributions for these parameters as well as the level of experimental observation noise. The tool also allows for hierarchical fitting, characterising the effect of inter-experiment variability. We demonstrate the use of the tool on a recently published dataset on multiple ion channel inhibition by multiple drug compounds. We compare the maximum likelihood, Bayesian and hierarchical Bayesian approaches. We then show how uncertainty in dose-response inputs can be characterised and propagated into a cardiac action potential simulation to give a probability distribution on model outputs.

SUBMITTER: Johnstone RH

PROVIDER: S-EPMC5134333 | biostudies-literature | 2016

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Hierarchical Bayesian inference for ion channel screening dose-response data.

Johnstone Ross H RH Bardenet Rémi R Gavaghan David J DJ Mirams Gary R GR

Wellcome open research 20161115

Dose-response (or 'concentration-effect') relationships commonly occur in biological and pharmacological systems and are well characterised by Hill curves. These curves are described by an equation with two parameters: the inhibitory concentration 50% (IC50); and the Hill coefficient. Typically just the 'best fit' parameter values are reported in the literature. Here we introduce a Python-based software tool, <b><i>PyHillFit</i></b> , and describe the underlying Bayesian inference methods that i ...[more]

PMID: 27918599

Similar Datasets

Project description:The stochastic behavior of single ion channels is most often described as an aggregated continuous-time Markov process with discrete states. For ligand-gated channels each state can represent a different conformation of the channel protein or a different number of bound ligands. Single-channel recordings show only whether the channel is open or shut: states of equal conductance are aggregated, so transitions between them have to be inferred indirectly. The requirement to filter noise from the raw signal further complicates the modeling process, as it limits the time resolution of the data. The consequence of the reduced bandwidth is that openings or shuttings that are shorter than the resolution cannot be observed; these are known as missed events. Postulated models fitted using filtered data must therefore explicitly account for missed events to avoid bias in the estimation of rate parameters and therefore assess parameter identifiability accurately. In this article, we present the first, to our knowledge, Bayesian modeling of ion-channels with exact missed events correction. Bayesian analysis represents uncertain knowledge of the true value of model parameters by considering these parameters as random variables. This allows us to gain a full appreciation of parameter identifiability and uncertainty when estimating values for model parameters. However, Bayesian inference is particularly challenging in this context as the correction for missed events increases the computational complexity of the model likelihood. Nonetheless, we successfully implemented a two-step Markov chain Monte Carlo method that we called "BICME", which performs Bayesian inference in models of realistic complexity. The method is demonstrated on synthetic and real single-channel data from muscle nicotinic acetylcholine channels. We show that parameter uncertainty can be characterized more accurately than with maximum-likelihood methods. Our code for performing inference in these ion channel models is publicly available.

Project description:When assessing the human risks due to exposure to environmental chemicals, traditional dose-response analyses are not straightforward when there are numerous high-quality epidemiological studies of priority cancer and non-cancer health outcomes. Given this wealth of information, selecting a single "best" study on which to base dose-response analyses is difficult and would potentially ignore much of the available data. Therefore, systematic approaches are necessary for the analysis of these rich databases. Examples are meta-analysis (and further, meta-regression), which are well established methods that consider and incorporate information from multiple studies into the estimation of risks due to exposure to environmental contaminants. In this paper, we propose a hierarchical, Bayesian meta-analysis approach for the dose-response analysis of multiple epidemiological studies. This paper is the second of two papers detailing this approach; the first covered "pre-analysis" steps necessary to prepare the data for dose-response modeling. This paper focuses on the hierarchical Bayesian approach to dose-response modeling and extrapolation of risk to populations of interest using the association between bladder cancer and oral inorganic arsenic (iAs) exposure as an illustrative case study. In particular, this paper addresses the modeling of both case-control and cohort studies with a flexible, logistic model in a hierarchical Bayesian framework that estimates study-specific slopes, as well as a pooled slope across all studies. This approach is akin to a random effects model in which no assumption is made a priori that there is a single, common slope for all included studies. Further, this paper also details extrapolation of the estimates of logistic slope to extra risk in a target population using a lifetable analysis and basic assumptions about background iAs exposure levels. In this case, the target population was the general United States population and information on all-cause mortality and incidence and mortality from bladder cancer was used to perform the lifetable analysis. The methods herein were developed for general use in investigating the association between any pollutant and observed health-effects in epidemiological studies. In order to demonstrate these methods, inorganic arsenic was chosen as a case study given the large epidemiological database that exists for this contaminant.

Dataset Information

Hierarchical Bayesian inference for ion channel screening dose-response data.

Publications

Hierarchical Bayesian inference for ion channel screening dose-response data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure