A systematic strategy for estimating hERG block potency and its implications in a new cardiac safety paradigm

Introduction

hERG block potency is widely used to calculate a drug's safety margin against its torsadogenic potential. Previous studies are confounded by use of different patch clamp electrophysiology protocols and a lack of statistical quantification of experimental variability. Since the new cardiac safety paradigm being discussed by the International Council for Harmonisation promotes a tighter integration of nonclinical and clinical data for torsadogenic risk assessment, a more systematic approach to estimate the hERG block potency and safety margin is needed.

Methods

A cross-industry study was performed to collect hERG data on 28 drugs with known torsadogenic risk using a standardized experimental protocol. A Bayesian hierarchical modeling (BHM) approach was used to assess the hERG block potency of these drugs by quantifying both the inter-site and intra-site variability. A modeling and simulation study was also done to evaluate protocol-dependent changes in hERG potency estimates.

Results

A systematic approach to estimate hERG block potency is established. The impact of choosing a safety margin threshold on torsadogenic risk evaluation is explored based on the posterior distributions of hERG potency estimated by this method. The modeling and simulation results suggest any potency estimate is specific to the protocol used.

Discussion

This methodology can estimate hERG block potency specific to a given voltage protocol. The relationship between safety margin thresholds and torsadogenic risk predictivity suggests the threshold should be tailored to each specific context of use, and safety margin evaluation may need to be integrated with other information to form a more comprehensive risk assessment.

•

hERG potency/safety margin is a widely used nonclinical cardiac safety strategy.

•

A new regulatory paradigm promotes the integration of nonclinical and clinical data.

•

Lack of uncertainty quantification hindered using hERG potency in the new paradigm.

•

A systematic method was established to address this limitation.

•

Analysis supports using different safety margin thresholds in different context.

Sinusoidal voltage protocols for rapid characterisation of ion channel kinetics

Key points

Ion current kinetics are commonly represented by current–voltage relationships, time constant–voltage relationships and subsequently mathematical models fitted to these. These experiments take substantial time, which means they are rarely performed in the same cell.

Rather than traditional square‐wave voltage clamps, we fitted a model to the current evoked by a novel sum‐of‐sinusoids voltage clamp that was only 8 s long.

Short protocols that can be performed multiple times within a single cell will offer many new opportunities to measure how ion current kinetics are affected by changing conditions.

The new model predicts the current under traditional square‐wave protocols well, with better predictions of underlying currents than literature models. The current under a novel physiologically relevant series of action potential clamps is predicted extremely well.

The short sinusoidal protocols allow a model to be fully fitted to individual cells, allowing us to examine cell–cell variability in current kinetics for the first time.

Abstract

Understanding the roles of ion currents is crucial to predict the action of pharmaceuticals and mutations in different scenarios, and thereby to guide clinical interventions in the heart, brain and other electrophysiological systems. Our ability to predict how ion currents contribute to cellular electrophysiology is in turn critically dependent on our characterisation of ion channel kinetics – the voltage‐dependent rates of transition between open, closed and inactivated channel states. We present a new method for rapidly exploring and characterising ion channel kinetics, applying it to the hERG potassium channel as an example, with the aim of generating a quantitatively predictive representation of the ion current. We fitted a mathematical model to currents evoked by a novel 8 second sinusoidal voltage clamp in CHO cells overexpressing hERG1a. The model was then used to predict over 5 minutes of recordings in the same cell in response to further protocols: a series of traditional square step voltage clamps, and also a novel voltage clamp comprising a collection of physiologically relevant action potentials. We demonstrate that we can make predictive cell‐specific models that outperform the use of averaged data from a number of different cells, and thereby examine which changes in gating are responsible for cell–cell variability in current kinetics. Our technique allows rapid collection of consistent and high quality data, from single cells, and produces more predictive mathematical ion channel models than traditional approaches.

Ion current kinetics are commonly represented by current–voltage relationships, time constant–voltage relationships and subsequently mathematical models fitted to these. These experiments take substantial time, which means they are rarely performed in the same cell.

Rather than traditional square‐wave voltage clamps, we fitted a model to the current evoked by a novel sum‐of‐sinusoids voltage clamp that was only 8 s long.

Short protocols that can be performed multiple times within a single cell will offer many new opportunities to measure how ion current kinetics are affected by changing conditions.

The new model predicts the current under traditional square‐wave protocols well, with better predictions of underlying currents than literature models. The current under a novel physiologically relevant series of action potential clamps is predicted extremely well.

The short sinusoidal protocols allow a model to be fully fitted to individual cells, allowing us to examine cell–cell variability in current kinetics for the first time.

Hierarchical Bayesian inference for ion channel screening dose-response data

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.

Hierarchical Bayesian inference for ion channel screening dose-response data

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.

Robust cell tracking in epithelial tissues through identification of maximum common subgraphs

Tracking of cells in live-imaging microscopy videos of epithelial sheets is a powerful tool for investigating fundamental processes in embryonic development. Characterizing cell growth, proliferation, intercalation and apoptosis in epithelia helps us to understand how morphogenetic processes such as tissue invagination and extension are locally regulated and controlled. Accurate cell tracking requires correctly resolving cells entering or leaving the field of view between frames, cell neighbour exchanges, cell removals and cell divisions. However, current tracking methods for epithelial sheets are not robust to large morphogenetic deformations and require significant manual interventions. Here, we present a novel algorithm for epithelial cell tracking, exploiting the graph-theoretic concept of a 'maximum common subgraph' to track cells between frames of a video. Our algorithm does not require the adjustment of tissue-specific parameters, and scales in sub-quadratic time with tissue size. It does not rely on precise positional information, permitting large cell movements between frames and enabling tracking in datasets acquired at low temporal resolution due to experimental constraints such as phototoxicity. To demonstrate the method, we perform tracking on the Drosophila embryonic epidermis and compare cell-cell rearrangements to previous studies in other tissues. Our implementation is open source and generally applicable to epithelial tissues.

Hierarchical Bayesian Modelling of Variability and Uncertainty in Synthetic Action Potential Traces

There are many sources of uncertainty in the recording and modelling of membrane action potentials (APs) from cardiomyocytes. For example, there are measurement, parameter, and model uncertainties. There is also extrinsic variability between cells, and intrinsic beat-to-beat variability within a single cell. These combined uncertainties and variability make it very difficult to extrapolate predictions from these models, since current AP models have single parameter values and thus produce a single AP trace. We aim to re-parameterise existing AP models to fit experimental data, and to quantify uncertainty associated with ion current densities when measuring and modelling these APs. We then wish to propagate this uncertainty into model predictions, such as ion channel block effected by a pharmaceutical compound. We perform an in silico study using synthetic data generated from different sets of parameters. We then 'forget' these parameter values and re-infer their distributions using hierarchical Markov chain Monte Carlo methods. We find that we can successfully infer the 'correct' distributions for most ion current densities for each AP trace, along with an approximation to the higher-level distribution from which these parameter values were sampled. There is, however, some level of unidentifiability amongst some of the current densities.

Uncertainty and variability in models of the cardiac action potential: Can we build trustworthy models?

Cardiac electrophysiology models have been developed for over 50 years, and now include detailed descriptions of individual ion currents and sub-cellular calcium handling. It is commonly accepted that there are many uncertainties in these systems, with quantities such as ion channel kinetics or expression levels being difficult to measure or variable between samples. Until recently, the original approach of describing model parameters using single values has been retained, and consequently the majority of mathematical models in use today provide point predictions, with no associated uncertainty.

In recent years, statistical techniques have been developed and applied in many scientific areas to capture uncertainties in the quantities that determine model behaviour, and to provide a distribution of predictions which accounts for this uncertainty. In this paper we discuss this concept, which is termed uncertainty quantification, and consider how it might be applied to cardiac electrophysiology models.

We present two case studies in which probability distributions, instead of individual numbers, are inferred from data to describe quantities such as maximal current densities. Then we show how these probabilistic representations of model parameters enable probabilities to be placed on predicted behaviours. We demonstrate how changes in these probability distributions across data sets offer insight into which currents cause beat-to-beat variability in canine APs. We conclude with a discussion of the challenges that this approach entails, and how it provides opportunities to improve our understanding of electrophysiology.

•

Uncertainty and variability in action potential models can be quantified.

•

A probabilistic method for inferring maximal current densities is developed and applied.

•

We use this to infer the currents responsible for canine beat-to-beat variability.

•

Emulation of mathematical models provides rich information at low computational cost.

•

The importance of considering uncertainty and variability in future is discussed.

Measurement of the proton-air cross section at √s=57 TeV with the Pierre Auger Observatory

We report a measurement of the proton-air cross section for particle production at the center-of-mass energy per nucleon of 57 TeV. This is derived from the distribution of the depths of shower maxima observed with the Pierre Auger Observatory: systematic uncertainties are studied in detail. Analyzing the tail of the distribution of the shower maxima, a proton-air cross section of [505±22(stat)(-36)(+28)(syst)] mb is found.