Human ex-vivo action potential model for pro-arrhythmia risk assessment

While current S7B/E14 guidelines have succeeded in protecting patients from QT-prolonging drugs, the absence of a predictive paradigm identifying pro-arrhythmic risks has limited the development of valuable drug programs. We investigated if a human ex-vivo action potential (AP)-based model could provide a more predictive approach for assessing pro-arrhythmic risk in man. Human ventricular trabeculae from ethically consented organ donors were used to evaluate the effects of dofetilide, d,l-sotalol, quinidine, paracetamol and verapamil on AP duration (APD) and recognized pro-arrhythmia predictors (short-term variability of APD at 90% repolarization (STV(APD90)), triangulation (ADP90-APD30) and incidence of early afterdepolarizations at 1 and 2Hz to quantitatively identify the pro-arrhythmic risk. Each drug was blinded and tested separately with 3 concentrations in triplicate trabeculae from 5 hearts, with one vehicle time control per heart. Electrophysiological stability of the model was not affected by sequential applications of vehicle (0.1% dimethyl sulfoxide). Paracetamol and verapamil did not significantly alter anyone of the AP parameters and were classified as devoid of pro-arrhythmic risk. Dofetilide, d,l-sotalol and quinidine exhibited an increase in the manifestation of pro-arrhythmia markers. The model provided quantitative and actionable activity flags and the relatively low total variability in tissue response allowed for the identification of pro-arrhythmic signals. Power analysis indicated that a total of 6 trabeculae derived from 2 hearts are sufficient to identify drug-induced pro-arrhythmia. Thus, the human ex-vivo AP-based model provides an integrative translational assay assisting in shaping clinical development plans that could be used in conjunction with the new CiPA-proposed approach.

Computational cardiology and risk stratification for sudden cardiac death: one of the grand challenges for cardiology in the 21st century

Risk stratification in the context of sudden cardiac death has been acknowledged as one of the major challenges facing cardiology for the past four decades. In recent years, the advent of high performance computing has facilitated organ-level simulation of the heart, meaning we can now examine the causes, mechanisms and impact of cardiac dysfunction in silico. As a result, computational cardiology, largely driven by the Physiome project, now stands at the threshold of clinical utility in regards to risk stratification and treatment of patients at risk of sudden cardiac death. In this white paper, we outline a roadmap of what needs to be done to make this translational step, using the relatively well-developed case of acquired or drug-induced long QT syndrome as an exemplar case.

A New Perspective in the Field of Cardiac Safety Testing through the Comprehensive In Vitro Proarrhythmia Assay Paradigm

For the past decade, cardiac safety screening to evaluate the propensity of drugs to produce QT interval prolongation and Torsades de Pointes (TdP) arrhythmia has been conducted according to ICH S7B and ICH E14 guidelines. Central to the existing approach are hERG channel assays and in vivo QT measurements. Although effective, the present paradigm carries a risk of unnecessary compound attrition and high cost, especially when considering costly thorough QT (TQT) studies conducted later in drug development. The C: omprehensive I: n Vitro P: roarrhythmia A: ssay (CiPA) initiative is a public-private collaboration with the aim of updating the existing cardiac safety testing paradigm to better evaluate arrhythmia risk and remove the need for TQT studies. It is hoped that CiPA will produce a standardized ion channel assay approach, incorporating defined tests against major cardiac ion channels, the results of which then inform evaluation of proarrhythmic actions in silico, using human ventricular action potential reconstructions. Results are then to be confirmed using human (stem cell-derived) cardiomyocytes. This perspective article reviews the rationale, progress of, and challenges for the CiPA initiative, if this new paradigm is to replace existing practice and, in time, lead to improved and widely accepted cardiac safety testing guidelines.

Pharmacological and electrophysiological characterization of AZSMO-23, an activator of the hERG K(+) channel

BACKGROUND AND PURPOSE

We aimed to characterize the pharmacology and electrophysiology of N-[3-(1H-benzimidazol-2-yl)-4-chloro-phenyl]pyridine-3-carboxamide (AZSMO-23), an activator of the human ether-a-go-go-related gene (hERG)-encoded K(+) channel (Kv 11.1).

EXPERIMENTAL APPROACH

Automated electrophysiology was used to study the pharmacology of AZSMO-23 on wild-type (WT), Y652A, F656T or G628C/S631C hERG, and on other cardiac ion channels. Its mechanism of action was characterized with conventional electrophysiology.

KEY RESULTS

AZSMO-23 activated WT hERG pre-pulse and tail current with EC50 values of 28.6 and 11.2 μM respectively. At 100 μM, pre-pulse current at +40 mV was increased by 952 ± 41% and tail current at -30 mV by 238 ± 13% compared with vehicle values. The primary mechanism for this effect was a 74.5 mV depolarizing shift in the voltage dependence of inactivation, without any shift in the voltage dependence of activation. Structure-activity relationships for this effect were remarkably subtle, with close analogues of AZSMO-23 acting as hERG inhibitors. AZSMO-23 blocked the mutant channel, hERG Y652A, but against another mutant channel, hERG F656T, its activator activity was enhanced. It inhibited activity of the G628C/S631C non-inactivating hERG mutant channel. AZSMO-23 was not hERG selective, as it blocked hKv 4.3-hKChIP2.2, hCav 3.2 and hKv 1.5 and activated hCav 1.2/β2/α2δ channels.

CONCLUSION AND IMPLICATIONS

The activity of AZSMO-23 and those of its close analogues suggest these compounds may be of value to elucidate the mechanism of type 2 hERG activators to better understand the pharmacology of this area from both a safety perspective and in relation to treatment of congenital long QT syndrome.

Assessment of cardiomyocyte contraction in human-induced pluripotent stem cell-derived cardiomyocytes

Functional changes to cardiomyocytes are a common cause of attrition in preclinical and clinical drug development. Current approaches to assess cardiomyocyte contractility in vitro are limited to low-throughput methods not amenable to early drug discovery. Human-induced pluripotent stem cell-derived cardiomyocytes (hiPS-CMs) were used to assess their suitability to detect drug-induced changes in cardiomyocyte contraction. Application of field stimulation and measurement of cardiac contraction (IonOptix edge detection) and Ca(2+) transients confirmed hiPS-CMs to be a suitable model to investigate drug-induced changes in cardiomyocyte contractility. Using a live cell, fast kinetic fluorescent assay with a Ca(2+) sensitive dye to test 31 inotropic and 20 non-inotropic compounds in vivo, we report that hiPS-CMs provide a high-throughput experimental model to detect changes in cardiomyocyte contraction that is applicable to early drug discovery with a sensitivity and specificity of 87% and 70%, respectively. Moreover, our data provide evidence of the detection of this liability at therapeutically relevant concentrations with throughput amenable to influencing chemical design in drug discovery. Measurement of multiple parameters of the Ca(2+) transient in addition to the number of Ca(2+) transients offered no insight into the mechanism of cardiomyocyte contraction.

Prediction of Thorough QT study results using action potential simulations based on ion channel screens

INTRODUCTION

Detection of drug-induced pro-arrhythmic risk is a primary concern for pharmaceutical companies and regulators. Increased risk is linked to prolongation of the QT interval on the body surface ECG. Recent studies have shown that multiple ion channel interactions can be required to predict changes in ventricular repolarisation and therefore QT intervals. In this study we attempt to predict the result of the human clinical Thorough QT (TQT) study, using multiple ion channel screening which is available early in drug development.

METHODS

Ion current reduction was measured, in the presence of marketed drugs which have had a TQT study, for channels encoded by hERG, CaV1.2, NaV1.5, KCNQ1/MinK, and Kv4.3/KChIP2.2. The screen was performed on two platforms - IonWorks Quattro (all 5 channels, 34 compounds), and IonWorks Barracuda (hERG & CaV1.2, 26 compounds). Concentration-effect curves were fitted to the resulting data, and used to calculate a percentage reduction in each current at a given concentration. Action potential simulations were then performed using the ten Tusscher and Panfilov (2006), Grandi et al. (2010) and O'Hara et al. (2011) human ventricular action potential models, pacing at 1Hz and running to steady state, for a range of concentrations.

RESULTS

We compared simulated action potential duration predictions with the QT prolongation observed in the TQT studies. At the estimated concentrations, simulations tended to underestimate any observed QT prolongation. When considering a wider range of concentrations, and conventional patch clamp rather than screening data for hERG, prolongation of ≥5ms was predicted with up to 79% sensitivity and 100% specificity.

DISCUSSION

This study provides a proof-of-principle for the prediction of human TQT study results using data available early in drug development. We highlight a number of areas that need refinement to improve the method's predictive power, but the results suggest that such approaches will provide a useful tool in cardiac safety assessment.

An impedance-based cellular assay using human iPSC-derived cardiomyocytes to quantify modulators of cardiac contractility

Cardiovascular toxicity, a prominent reason for late-stage failures in drug development, has resulted in a demand for in vitro assays that can predict this liability in early drug discovery. Current in vitro cardiovascular safety testing primarily focuses on ion channel modulation and low throughput cardiomyocyte (CM) contractility measurements. We evaluated both human induced pluripotent stem cell-derived CMs (hiPSC-CMs) and rat neonatal CMs (rat CMs) on the xCELLigence Cardio system which uses impedance technology to quantify CM beating properties in a 96-well format. Forty-nine compounds were tested in concentration-response mode to determine potency for modulation of CM beating, a surrogate biomarker for contractility. These compounds had previously been tested in vivo and in a low throughput in vitro optical-based contractility assay that measures sarcomere shortening in electrically paced dog CMs. In comparison with in vivo contractility effects, hiPSC-CM impedance had assay sensitivity, specificity, and accuracy values of 90%, 74%, and 82%, respectively. These values compared favorably to values reported for the dog CM optical assay (83%, 84%, and 82%) and were slightly better than impedance using rat CMs (77%, 74%, and 74%). The potency values from the hiPSC-CM and rat CM assays spanned four orders of magnitude and correlated with values from the dog CM optical assay (r(2 )= 0.76 and 0.70, respectively). The Cardio system assay has >5× higher throughput than the optical assay. Thus, hiPSC-CM impedance testing can help detect the human cardiotoxic potential of novel therapeutics early in drug discovery, and if a hazard is identified, has sufficient throughput to support the design-make-test-analyze cycle to mitigate this liability.

Interventricular differences in β-adrenergic responses in the canine heart: role of phosphodiesterases

Background

RV and LV have different embryologic, structural, metabolic, and electrophysiologic characteristics, but whether interventricular differences exist in β‐adrenergic (β‐AR) responsiveness is unknown. In this study, we examine whether β‐AR response and signaling differ in right (RV) versus left (LV) ventricles.

Methods and Results

Sarcomere shortening, Ca²⁺ transients, I_Ca,L and I_Ks currents were recorded in isolated dog LV and RV midmyocytes. Intracellular [cAMP] and PKA activity were measured by live cell imaging using FRET‐based sensors. Isoproterenol increased sarcomere shortening ≈10‐fold and Ca²⁺‐transient amplitude ≈2‐fold in LV midmyocytes (LVMs) versus ≈25‐fold and ≈3‐fold in RVMs. FRET imaging using targeted Epac2camps sensors revealed no change in subsarcolemmal [cAMP], but a 2‐fold higher β‐AR stimulation of cytoplasmic [cAMP] in RVMs versus LVMs. Accordingly, β‐AR regulation of I_Ca,L and I_Ks were similar between LVMs and RVMs, whereas cytoplasmic PKA activity was increased in RVMs. Both PDE3 and PDE4 contributed to the β‐AR regulation of cytoplasmic [cAMP], and the difference between LVMs and RVMs was abolished by PDE3 inhibition and attenuated by PDE4 inhibition. Finally LV and RV intracavitary pressures were recorded in anesthetized beagle dogs. A bolus injection of isoproterenol increased RV dP/dt_max≈5‐fold versus 3‐fold in LV.

Conclusion

Canine RV and LV differ in their β‐AR response due to intrinsic differences in myocyte β‐AR downstream signaling. Enhanced β‐AR responsiveness of the RV results from higher cAMP elevation in the cytoplasm, due to a decreased degradation by PDE3 and PDE4 in the RV compared to the LV.

hERG inhibitors with similar potency but different binding kinetics do not pose the same proarrhythmic risk: implications for drug safety assessment

INTRODUCTION

Since the discovery of the link that exists between drug-induced hERG inhibition and Torsade de Pointes (TdP), extreme attention has been given to avoid new drugs inhibiting this channel. hERG inhibition is routinely screened for in new drugs and, typically, IC50 values are compared to projected plasma concentrations to define a safety margin.

METHODS AND RESULTS

We aimed to show that drugs with similar hERG potency are not uniformly pro-arrhythmic-this depends on the drug binding kinetics and mode of action (trapped or not) rather than the IC50 value only. We used a mathematical model of hERG and its related encoded current IKr to simulate drug binding in different configurations. Expression systems mimicking the screening process were first investigated. hERG model was then incorporated into a canine action potential (AP) and tissue model to study the impact of drug binding configurations on AP and pseudo-ECG (QT interval prolongation). Our data show that: (1) trapped and not trapped configurations and different binding kinetics could be identified during hERG screening; (2) slow binding, not trapped drugs, induced less AP prolongation and minimal QT interval prolongation (4.7%) at a concentration equal to the IC50 whereas maximal pro-arrhythmic risk was observed for trapped drugs at the same concentration (QT interval prolongation, 23.1%).

CONCLUSION

Our study demonstrates the need for screening for hERG binding configurations rather than potency alone. It also demonstrates the potential link between hERG, drug mode of action and TdP, and the need to question the current regulatory guidance.

Variability in high-throughput ion-channel screening data and consequences for cardiac safety assessment

Introduction

Unwanted drug interactions with ionic currents in the heart can lead to an increased proarrhythmic risk to patients in the clinic. It is therefore a priority for safety pharmacology teams to detect block of cardiac ion channels, and new technologies have enabled the development of automated and high-throughput screening assays using cell lines. As a result of screening multiple ion-channels there is a need to integrate information, particularly for compounds affecting more than one current, and mathematical electrophysiology in-silico action potential models are beginning to be used for this.

Methods

We quantified the variability associated with concentration-effect curves fitted to recordings from high-throughput Molecular Devices IonWorks® Quattro™ screens when detecting block of I_Kr (hERG), I_Na (NaV1.5), I_CaL (CaV1.2), I_Ks (KCNQ1/minK) and I_to (Kv4.3/KChIP2.2), and the Molecular Devices FLIPR® Tetra fluorescence screen for I_CaL (CaV1.2), for control compounds used at AstraZeneca and GlaxoSmithKline. We examined how screening variability propagates through in-silico action potential models for whole cell electrical behaviour, and how confidence intervals on model predictions can be estimated with repeated simulations.

Results

There are significant levels of variability associated with high-throughput ion channel electrophysiology screens. This variability is of a similar magnitude for different cardiac ion currents and different compounds. Uncertainty in the Hill coefficients of reported concentration-effect curves is particularly high. Depending on a compound’s ion channel blocking profile, the uncertainty introduced into whole-cell predictions can become significant.

Discussion

Our technique allows confidence intervals to be placed on computational model predictions that are based on high-throughput ion channel screens. This allows us to suggest when repeated screens should be performed to reduce uncertainty in a compound’s action to acceptable levels, to allow a meaningful interpretation of the data.

High-throughput screening of drug-binding dynamics to HERG improves early drug safety assessment

The use of computational models to predict drug-induced changes in the action potential (AP) is a promising approach to reduce drug safety attrition but requires a better representation of more complex drug-target interactions to improve the quantitative prediction. The blockade of the human ether-a-go-go-related gene (HERG) channel is a major concern for QT prolongation and Torsade de Pointes risk. We aim to develop quantitative in-silico AP predictions based on a new electrophysiological protocol (suitable for high-throughput HERG screening) and mathematical modeling of ionic currents. Electrophysiological recordings using the IonWorks device were made from HERG channels stably expressed in Chinese hamster ovary cells. A new protocol that delineates inhibition over time was applied to assess dofetilide, cisapride, and almokalant effects. Dynamic effects displayed distinct profiles for these drugs compared with concentration-effects curves. Binding kinetics to specific states were identified using a new HERG Markov model. The model was then modified to represent the canine rapid delayed rectifier K+ current at 37°C and carry out AP predictions. Predictions were compared with a simpler model based on conductance reduction and were found to be much closer to experimental data. Improved sensitivity to concentration and pacing frequency variables was obtained when including binding kinetics. Our new electrophysiological protocol is suitable for high-throughput screening and is able to distinguish drug-binding kinetics. The association of this protocol with our modeling approach indicates that quantitative predictions of AP modulation can be obtained, which is a significant improvement compared with traditional conductance reduction methods.

Diastolic spontaneous calcium release from the sarcoplasmic reticulum increases beat-to-beat variability of repolarization in canine ventricular myocytes after β-adrenergic stimulation

RATIONALE

Spontaneous Ca(2+) release (SCR) from the sarcoplasmic reticulum can cause delayed afterdepolarizations and triggered activity, contributing to arrhythmogenesis during β-adrenergic stimulation. Excessive beat-to-beat variability of repolarization duration (BVR) is a proarrhythmic marker. Previous research has shown that BVR is increased during intense β-adrenergic stimulation, leading to SCR.

OBJECTIVE

We aimed to determine ionic mechanisms controlling BVR under these conditions.

METHODS AND RESULTS

Membrane potentials and cell shortening or Ca(2+) transients were recorded from isolated canine left ventricular myocytes in the presence of isoproterenol. Action-potential (AP) durations after delayed afterdepolarizations were significantly prolonged. Addition of slowly activating delayed rectifier K(+) current (I(Ks)) blockade led to further AP prolongation after SCR, and this strongly correlated with exaggerated BVR. Suppressing SCR via inhibition of ryanodine receptors, Ca(2+)/calmodulin-dependent protein kinase II inhibition, or by using Mg(2+) or flecainide eliminated delayed afterdepolarizations and decreased BVR independent of effects on AP duration. Computational analyses and voltage-clamp experiments measuring L-type Ca(2+) current (I(CaL)) with and without previous SCR indicated that I(CaL) was increased during Ca(2+)-induced Ca(2+) release after SCR, and this contributes to AP prolongation. Prolongation of QT, T(peak)-T(end) intervals, and left ventricular monophasic AP duration of beats after aftercontractions occurred before torsades de pointes in an in vivo dog model of drug-induced long-QT1 syndrome.

CONCLUSIONS

SCR contributes to increased BVR by interspersed prolongation of AP duration, which is exacerbated during I(Ks) blockade. Attenuation of Ca(2+)-induced Ca(2+) release by SCR underlies AP prolongation via increased I(CaL.) These data provide novel insights into arrhythmogenic mechanisms during β-adrenergic stimulation besides triggered activity and illustrate the importance of I(Ks) function in preventing excessive BVR.

An in silico canine cardiac midmyocardial action potential duration model as a tool for early drug safety assessment

Cell lines expressing ion channels (IC) and the advent of plate-based electrophysiology device have enabled a molecular understanding of the action potential (AP) as a means of early QT assessment. We sought to develop an in silico AP (isAP) model that provides an assessment of the effect of a compound on the myocyte AP duration (APD) using concentration-effect curve data from a panel of five ICs (hNav1.5, hCav1.2, hKv4.3/hKChIP2.2, hKv7.1/hminK, hKv11.1). A test set of 53 compounds was selected to cover a range of selective and mixed IC modulators that were tested for their effects on optically measured APD. A threshold of >10% change in APD at 90% repolarization (APD90) was used to signify an effect at the top test concentration. To capture the variations observed in left ventricular midmyocardial myocyte APD data from 19 different dogs, the isAP model was calibrated to produce an ensemble of 19 model variants that could capture the shape and form of the APs and also quantitatively replicate dofetilide- and diltiazem-induced APD90 changes. Provided with IC panel data only, the isAP model was then used, blinded, to predict APD90 changes greater than 10%. At a simulated concentration of 30 μM and based on a criterion that six of the variants had to agree, isAP prediction was scored as showing greater than 80% predictivity of compound activity. Thus, early in drug discovery, the isAP model allows integrating separate IC data and is amenable to the throughput required for use as a virtual screen.

hERG subunit composition determines differential drug sensitivity

BACKGROUND AND PURPOSE

The majority of human ether-a-go-go-related gene (hERG) screens aiming to minimize the risk of drug-induced long QT syndrome have been conducted using heterologous systems expressing the hERG 1a subunit, although both hERG 1a and 1b subunits contribute to the K+ channels producing the repolarizing current I(Kr) . We tested a range of compounds selected for their diversity to determine whether hERG 1a and 1a/1b channels exhibit different sensitivities that may influence safety margins or contribute to a stratified risk analysis.

EXPERIMENTAL APPROACH

We used the IonWorks™ plate-based electrophysiology device to compare sensitivity of hERG 1a and 1a/1b channels stably expressed in HEK293 cells to 50 compounds previously shown to target hERG channels. Potency was determined as IC₅₀ values (µM) obtained from non-cumulative, eight-point concentration-effect curves of normalized data, fitted to the Hill equation. To minimize possible sources of variability, compound potency was assessed using test plates arranged in alternating columns of cells expressing hERG 1a and 1a/1b.

KEY RESULTS

Although the potency of most compounds was similar for the two targets, some surprising differences were observed. Fluoxetine (Prozac) was more potent at blocking hERG 1a/1b than 1a channels, yielding a corresponding reduction in the safety margin. In contrast, E-4031 was a more potent blocker of hERG 1a compared with 1a/1b channels, as previously reported, as was dofetilide, another high-affinity blocker.

CONCLUSIONS AND IMPLICATIONS

The current assays may underestimate the risk of some drugs to cause torsades de pointes arrhythmia, and overestimate the risk of others.