Measuring kinetics and potency of hERG block for CiPA

Experimentally validated modeling of dynamic drug-hERG channel interactions reproducing the binding mechanisms and its importance in action potential duration

BACKGROUND AND OBJECTIVE

Assessment of drug cardiotoxicity is critical in the development of new compounds and modeling of drug-binding dynamics to hERG can improve early cardiotoxicity assessment. We previously developed a methodology to generate Markovian models reproducing preferential state-dependent binding properties, trapping dynamics and the onset of IKr block using simple voltage clamp protocols. Here, we test this methodology with real IKr blockers and investigate the impact of drug dynamics on action potential prolongation.

METHODS

Experiments were performed on HEK cells stably transfected with hERG and using the Nanion SyncroPatch 384i. Three protocols, P-80, P0 and P 40, were applied to obtain the experimental data from the drugs and the Markovian models were generated using our pipeline. The corresponding static models were also generated and a modified version of the O´Hara-Rudy action potential model was used to simulate the action potential duration.

RESULTS

The experimental Hill plots and the onset of IKr block of ten compounds were obtained using our voltage clamp protocols and the models generated successfully mimicked these experimental data, unlike the CiPA dynamic models. Marked differences in APD prolongation were observed when drug effects were simulated using the dynamic models and the static models.

CONCLUSIONS

These new dynamic models of ten well-known IKr blockers constitute a validation of our methodology to model dynamic drug-hERG channel interactions and highlight the importance of state-dependent binding, trapping dynamics and the time-course of IKr block to assess drug effects even at the steady-state.

In Silico Human Cardiomyocyte Action Potential Modeling: Exploring Ion Channel Input Combinations

In silico modeling offers an opportunity to supplement and accelerate cardiac safety testing. With in silico modeling, computational simulation methods are used to predict electrophysiological interactions and pharmacological effects of novel drugs on critical physiological processes. The O'Hara-Rudy's model was developed to predict the response to different ion channel inhibition levels on cardiac action potential duration (APD) which is known to directly correlate with the QT interval. APD data at 30% 60% and 90% inhibition were derived from the model to delineate possible ventricular arrhythmia scenarios and the marginal contribution of each ion channel to the model. Action potential values were calculated for epicardial, myocardial, and endocardial cells, with action potential curve modeling. This study assessed cardiac ion channel inhibition data combinations to consider when undertaking in silico modeling of proarrhythmic effects as stipulated in the Comprehensive in Vitro Proarrhythmia Assay (CiPA). As expected, our data highlight the importance of the delayed rectifier potassium channel (IKr) as the most impactful channel for APD prolongation. The impact of the transient outward potassium channel (Ito) inhibition on APD was minimal while the inward rectifier (IK1) and slow component of the delayed rectifier potassium channel (IKs) also had limited APD effects. In contrast, the contribution of fast sodium channel (INa) and/or L-type calcium channel (ICa) inhibition resulted in substantial APD alterations supporting the pharmacological relevance of in silico modeling using input from a limited number of cardiac ion channels including IKr, INa, and ICa, at least at an early stage of drug development.

The impact of uncertainty in hERG binding mechanism on in silico predictions of drug-induced proarrhythmic risk

BACKGROUND AND PURPOSE

Drug-induced reduction of the rapid delayed rectifier potassium current carried by the human Ether-à-go-go-Related Gene (hERG) channel is associated with increased risk of arrhythmias. Recent updates to drug safety regulatory guidelines attempt to capture each drug's hERG binding mechanism by combining in vitro assays with in silico simulations. In this study, we investigate the impact on in silico proarrhythmic risk predictions due to uncertainty in the hERG binding mechanism and physiological hERG current model.

EXPERIMENTAL APPROACH

Possible pharmacological binding models were designed for the hERG channel to account for known and postulated small molecule binding mechanisms. After selecting a subset of plausible binding models for each compound through calibration to available voltage-clamp electrophysiology data, we assessed their effects, and the effects of different physiological models, on proarrhythmic risk predictions.

KEY RESULTS

For some compounds, multiple binding mechanisms can explain the same data produced under the safety testing guidelines, which results in different inferred binding rates. This can result in substantial uncertainty in the predicted torsade risk, which often spans more than one risk category. By comparison, we found that the effect of a different hERG physiological current model on risk classification was subtle.

CONCLUSION AND IMPLICATIONS

The approach developed in this study assesses the impact of uncertainty in hERG binding mechanisms on predictions of drug-induced proarrhythmic risk. For some compounds, these results imply the need for additional binding data to decrease uncertainty in safety-critical applications.

Harnessing AlphaFold to reveal state secrets: Prediction of hERG closed and inactivated states

To design safe, selective, and effective new therapies, there must be a deep understanding of the structure and function of the drug target. One of the most difficult problems to solve has been resolution of discrete conformational states of transmembrane ion channel proteins. An example is KV11.1 (hERG), comprising the primary cardiac repolarizing current, IKr. hERG is a notorious drug anti-target against which all promising drugs are screened to determine potential for arrhythmia. Drug interactions with the hERG inactivated state are linked to elevated arrhythmia risk, and drugs may become trapped during channel closure. However, the structural details of multiple conformational states have remained elusive. Here, we guided AlphaFold2 to predict plausible hERG inactivated and closed conformations, obtaining results consistent with myriad available experimental data. Drug docking simulations demonstrated hERG state-specific drug interactions aligning well with experimental results, revealing that most drugs bind more effectively in the inactivated state and are trapped in the closed state. Molecular dynamics simulations demonstrated ion conduction that aligned with earlier studies. Finally, we identified key molecular determinants of state transitions by analyzing interaction networks across closed, open, and inactivated states in agreement with earlier mutagenesis studies. Here, we demonstrate a readily generalizable application of AlphaFold2 as a novel method to predict discrete protein conformations and novel linkages from structure to function.

Automatic modeling of dynamic drug-hERG channel interactions using three voltage protocols and machine learning techniques: A simulation study

BACKGROUND

Assessment of drug cardiac safety is critical in the development of new compounds and is commonly addressed by evaluating the half-maximal blocking concentration of the potassium human ether-à-go-go related gene (hERG) channels. However, recent works have evidenced that the modelling of drug-binding dynamics to hERG can help to improve early cardiac safety assessment. Our goal is to develop a methodology to automatically generate Markovian models of the drug-hERG channel interactions.

METHODS

The training and the test sets consisted of 20800 and 5200 virtual drugs, respectively, distributed into 104 groups with different affinities and kinetics to the conformational states of the channel. In our system, drugs may bind to any state (individually or simultaneously), with different degrees of preference for a conformational state and the change of the conformational state of the drug bound channels may be restricted or allowed. To model such a wide range of possibilities, 12 Markovian chains are considered. Our approach uses the response of the drugs to our three previously developed voltage clamp protocols, which enhance the differences in the probabilities of occupying a certain conformational state of the channel (open, closed and inactivated). The computing tool is comprised of a classifier and a parameter optimizer and uses linear interpolation, support vector machines and a simplex method for function minimization.

RESULTS

We propose a novel methodology that automatically generates dynamic drug models using Markov model formulations and that elucidates the states where the drug binds and unbinds and the preferential binding state using data obtained from simple voltage clamp protocols that captures the preferential state-dependent binding properties, the relative affinities, trapping and non-trapping dynamics and the onset of I_Kr block. Overall, the tool correctly predicted the class of 92.04% of the drugs and the model provided by the tool accurately fitted the response of the target compound, the mean accuracy being 97.53%. Moreover, generation of the dynamic model of an I_Kr blocker from its response to our voltage clamp protocols usually takes less than an hour on a common desktop computer.

CONCLUSION

Our methodology could be very useful to model and simulate dynamic drug-hERG channel interactions. It would contribute to the improvement of the preclinical assessment of the proarrhythmic risk of drugs that inhibit I_Kr and the efficacy of antiarrhythmic I_Kr blockers.

Safety pharmacology in 2022: Taking one small step for cardiovascular safety assay development but one giant leap for regulatory drug safety assessment

The 2021 Annual Safety Pharmacology (SP) Society (SPS) meeting was held virtually October 4–8, 2021 due to the continuing COVID-19 global pandemic. This themed issue of J Pharmacol Toxicol Methods comprises articles arising from the meeting. As in previous years the manuscripts reflect various areas of innovation in SP including a perspective on aging and its impact on drug attrition during safety assessments, an integrated assessment of respiratory, cardiovascular and animal activity of in vivo nonclinical studies, development of a dynamic QT-rate correction method in primates, evaluation of the “comprehensive in vitro proarrhythmia assay” (CiPA) ion channel protocol to the automated patch clamp, and best practices regarding the conduct of hERG electrophysiology studies and an analysis of secondary pharmacology assays by the FDA. The meeting also generated 85 abstracts (reproduced in the current volume of J Pharmacol Toxicol Methods). It appears that the validation of methods remains a challenge in SP. Nevertheless, the continued efforts to mine approaches to detection of proarrhythmia liability remains a baffling obsession given the ability of Industry to completely prevent drugs entering into clinical study only to be found to have proarrhythmic properties, with no reports of such for at least ten years. Perhaps it is time to move on from CiPA and find genuine problems to solve?

Translating the measurement of hERG kinetics and drug block for CiPA to a high throughput platform

The Comprehensive in vitro Proarrhythmic Assay (CiPA) has promoted use of in silico models of drug effects on cardiac repolarization to improve proarrhythmic risk prediction. These models contain a pharmacodynamic component describing drug binding to hERG channels that required in vitro data for kinetics of block, in addition to potency, to constrain them. To date, development and validation has been undertaken using data from manual patch-clamp. The application of this approach at scale requires the development of a high-throughput, automated patch-clamp (APC) implementation. Here, we present a comprehensive analysis of the implementation of the Milnes, or CiPA dynamic protocol, on an APC platform, including quality control and data analysis. Kinetics and potency of block were assessed for bepridil, cisapride, terfenadine and verapamil with data retention/QC pass rate of 21.8% overall, or as high as 50.4% when only appropriate sweep lengths were considered for drugs with faster kinetics. The variability in IC₅₀ and kinetics between manual and APC was comparable to that seen between sites/platforms in previous APC studies of potency. Whilst the experimental success is less than observed in screens of potency alone, it is still significantly greater than manual patch. With the modifications to protocol design, including sweep length, number of repetitions, and leak correction recommended in this study, this protocol can be applied on APC to acquire data comparable to manual patch clamp.

Assessment of proarrhythmogenic risk for chloroquine and hydroxychloroquine using the CiPA concept

Chloroquine and hydroxychloroquine have been proposed recently as therapy for SARS-CoV-2-infected patients, but during 3 months of extensive use concerns were raised related to their clinical effectiveness and arrhythmogenic risk. Therefore, we estimated for these compounds several proarrhythmogenic risk predictors according to the Comprehensive in vitro Proarrhythmia Assay (CiPA) paradigm. Experiments were performed with either CytoPatch™2 automated or manual patch-clamp setups on HEK293T cells stably or transiently transfected with hERG1, hNav1.5, hKir2.1, hKv7.1+hMinK, and on Pluricyte® cardiomyocytes (Ncardia), using physiological solutions. Dose-response plots of hERG1 inhibition fitted with Hill functions yielded IC50 values in the low micromolar range for both compounds. We found hyperpolarizing shifts of tens of mV, larger for chloroquine, in the voltage-dependent activation but not inactivation, as well as a voltage-dependent block of hERG current, larger at positive potentials. We also found inhibitory effects on peak and late I_Na and on I_K1, with IC50 of tens of μM and larger for chloroquine. The two compounds, tested on Pluricyte® cardiomyocytes using the β-escin-perforated method, inhibited I_Kr, I_CaL, I_{Na peak}, but had no effect on I_f. In current-clamp they caused action potential prolongation. Our data and those from literature for I_to were used to compute proarrhythmogenic risk predictors B_net (Mistry HB, 2018) and Q_net (Dutta S et al., 2017), with hERG1 blocking/unblocking rates estimated from time constants of fractional block. Although the two antimalarials are successfully used in autoimmune diseases, and chloroquine may be effective in atrial fibrillation, assays place these drugs in the intermediate proarrhythmogenic risk group.

Rutaecarpine targets hERG channels and participates in regulating electrophysiological properties leading to ventricular arrhythmia

Drug-mediated or medical condition-mediated disruption of hERG function accounts for the main cause of acquired long-QT syndrome (acLQTs), which predisposes affected individuals to ventricular arrhythmias (VA) and sudden death. Many Chinese herbal medicines, especially alkaloids, have risks of arrhythmia in clinical application. The characterized mechanisms behind this adverse effect are frequently associated with inhibition of cardiac hERG channels. The present study aimed to assess the potent effect of Rutaecarpine (Rut) on hERG channels. hERG-HEK293 cell was applied for evaluating the effect of Rut on hERG channels and the underlying mechanism. hERG current (I_hERG ) was measured by patch-clamp technique. Protein levels were analysed by Western blot, and the phosphorylation of Sp1 was determined by immunoprecipitation. Optical mapping and programmed electrical stimulation were used to evaluate cardiac electrophysiological activities, such as APD, QT/QTc, occurrence of arrhythmia, phase singularities (PSs), and dominant frequency (DF). Our results demonstrated that Rut reduced the I_hERG by binding to F656 and Y652 amino acid residues of hERG channel instantaneously, subsequently accelerating the channel inactivation, and being trapped in the channel. The level of hERG channels was reduced by incubating with Rut for 24 hours, and Sp1 in nucleus was inhibited simultaneously. Mechanismly, Rut reduced threonine (Thr)/ tyrosine (Tyr) phosphorylation of Sp1 through PI3K/Akt pathway to regulate hERG channels expression. Cell-based model unables to fully reveal the pathological process of arrhythmia. In vivo study, we found that Rut prolonged QT/QTc intervals and increased induction rate of ventricular fibrillation (VF) in guinea pig heart after being dosed Rut for 2 weeks. The critical reasons led to increased incidence of arrhythmias eventually were prolonged APD₉₀ and APD₅₀ and the increase of DF, numbers of PSs, incidence of early after-depolarizations (EADs). Collectively, the results of this study suggest that Rut could reduce the I_hERG by binding to hERG channels through F656 and Y652 instantaneously. While, the PI3K/Akt/Sp1 axis may play an essential role in the regulation of hERG channels, from the perspective of the long-term effects of Rut (incubating for 24 hours). Importantly, the changes of electrophysiological properties by Rut were the main cause of VA.

Cardiac hERG K+ Channel as Safety and Pharmacological Target

The human ether-á-go-go related gene (hERG, KCNH2) encodes the pore-forming subunit of the potassium channel responsible for a fast component of the cardiac delayed rectifier potassium current (I_Kr). Outward I_Kr is an important determinant of cardiac action potential (AP) repolarization and effectively controls the duration of the QT interval in humans. Dysfunction of hERG channel can cause severe ventricular arrhythmias and thus modulators of the channel, including hERG inhibitors and activators, continue to attract intense pharmacological interest. Certain inhibitors of hERG channel prolong the action potential duration (APD) and effective refractory period (ERP) to suppress premature ventricular contraction and are used as class III antiarrhythmic agents. However, a reduction of the hERG/I_Kr current has been recognized as a predominant mechanism responsible for the drug-induced delayed repolarization known as acquired long QT syndromes (LQTS), which is linked to an increased risk for "torsades de pointes" (TdP) ventricular arrhythmias and sudden cardiac death. Many drugs of different classes and structures have been identified to carry TdP risk. Hence, assessing hERG/I_Kr blockade of new drug candidates is mandatory in the drug development process according to the regulatory agencies. In contrast, several hERG channel activators have been shown to enhance I_Kr and shorten the APD and thus might have potential antiarrhythmic effects against pathological LQTS. However, these activators may also be proarrhythmic due to excessive shortening of APD and the ERP.

When Does the IC50 Accurately Assess the Blocking Potency of a Drug?

Preclinical assessment of drug-induced proarrhythmicity is typically evaluated by the potency of the drug to block the potassium human ether-à-go-go-related gene (hERG) channels, which is currently quantified by the IC₅₀. However, channel block depends on the experimental conditions. Our aim is to improve the evaluation of the blocking potency of drugs by designing experimental stimulation protocols to measure the IC₅₀ that will help to decide whether the IC₅₀ is representative enough. We used the state-of-the-art mathematical models of the cardiac electrophysiological activity to design three stimulation protocols that enhance the differences in the probabilities to occupy a certain conformational state of the channel and, therefore, the potential differences in the blocking effects of a compound. We simulated an extensive set of 144 in silico I_Kr blockers with different kinetics and affinities to conformational states of the channel and we also experimentally validated our key predictions. Our results show that the IC₅₀ protocol dependency relied on the tested compounds. Some of them showed no differences or small differences on the IC₅₀ value, which suggests that the IC₅₀ could be a good indicator of the blocking potency in these cases. However, others provided highly protocol dependent IC₅₀ values, which could differ by even 2 orders of magnitude. Moreover, the protocols yielding the maximum IC₅₀ and minimum IC₅₀ depended on the drug, which complicates the definition of a "standard" protocol to minimize the influence of the stimulation protocol on the IC₅₀ measurement in safety pharmacology. As a conclusion, we propose the adoption of our three-protocol IC₅₀ assay to estimate the potency to block hERG in vitro. If the IC₅₀ values obtained for a compound are similar, then the IC₅₀ could be used as an indicator of its blocking potency, otherwise kinetics and state-dependent binding properties should be accounted.

Human Heart Cardiomyocytes in Drug Discovery and Research: New Opportunities in Translational Sciences

In preclinical drug development, accurate prediction of drug effects on the human heart is critically important, whether in the context of cardiovascular safety or for the purpose of modulating cardiac function to treat heart disease. Current strategies have significant limitations, whereby, cardiotoxic drugs can escape detection or potential life-saving therapies are abandoned due to false positive toxicity signals. Thus, new and more reliable translational approaches are urgently needed to help accelerate the rate of new therapy development. Renewed efforts in the recovery of human donor hearts for research and in cardiomyocyte isolation methods, are providing new opportunities for preclinical studies in adult primary cardiomyocytes. These cells exhibit the native physiological and pharmacological properties, overcoming the limitations presented by artificial cellular models, animal models and have great potential for providing an excellent tool for preclinical drug testing. Adult human primary cardiomyocytes have already shown utility in assessing drug-induced cardiotoxicity risk and helping in the identification of new treatments for cardiac diseases, such as heart failure and atrial fibrillation. Finally, strategies with actionable decision-making trees that rely on data derived from adult human primary cardiomyocytes will provide the holistic insights necessary to accurately predict human heart effects of drugs.

hERG Function in Light of Structure

The human ether-a-go-go-related gene1 (hERG) ion channel has been the subject of fascination since it was identified as a target of long QT syndrome more than 20 years ago. In this Biophysical Perspective, we look at what makes hERG intriguing and vexingly unique. By probing recent high-resolution structures in the context of functional and biochemical data, we attempt to summarize new insights into hERG-specific function and articulate important unanswered questions. X-ray crystallography and cryo-electron microscopy have revealed features not previously on the radar-the "nonswapped" transmembrane architecture, an "intrinsic ligand," and hydrophobic pockets off a pore cavity that is surprisingly small. Advances in our understanding of drug block and inactivation mechanisms are noted, but a full picture will require more investigation.

Novel method for action potential measurements from intact cardiac monolayers with multiwell microelectrode array technology

The cardiac action potential (AP) is vital for understanding healthy and diseased cardiac biology and drug safety testing. However, techniques for high throughput cardiac AP measurements have been limited. Here, we introduce a novel technique for reliably increasing the coupling of cardiomyocyte syncytium to planar multiwell microelectrode arrays, resulting in a stable, label-free local extracellular action potential (LEAP). We characterized the reliability and stability of LEAP, its relationship to the field potential, and its efficacy for quantifying AP morphology of human induced pluripotent stem cell derived and primary rodent cardiomyocytes. Rise time, action potential duration, beat period, and triangulation were used to quantify compound responses and AP morphology changes induced by genetic modification. LEAP is the first high throughput, non-invasive, label-free, stable method to capture AP morphology from an intact cardiomyocyte syncytium. LEAP can accelerate our understanding of stem cell models, while improving the automation and accuracy of drug testing.

Recording of multiple ion current components and action potentials in human induced pluripotent stem cell-derived cardiomyocytes via automated patch-clamp

INTRODUCTION

The Comprehensive in vitro Proarrhythmia Assay (CiPA) initiative proposes a three-step approach to evaluate proarrhythmogenic liability of drug candidates: effects on individual ion channels in heterologous expression systems, integrating these data into in-silico models of the electrical activity of human cardiomyocytes, and comparison with experiments on human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM). Here we introduce patch-clamp electrophysiology techniques on hiPSC-CM to combine two of the CiPA steps in one assay.

METHODS

We performed automated patch-clamp experiments on hiPSC-CM (Cor.4U^®, Ncardia) using the CytoPatch™2 platform in ruptured whole-cell and β-escin-perforated-patch configurations. A combination of three voltage-clamp protocols allowed recording of five distinct ion current components (voltage-gated Na⁺ current, L-type Ca²⁺ current, transient outward K⁺ current, delayed rectifier K⁺ current, and "funny" hyperpolarization-activated current) from the same cell. We proved their molecular identity by either Na⁺ replacement with choline or by applying specific blockers: nifedipine, cisapride, chromanol 293B, phrixotoxin-1, ZD7288. We developed a C++ script for automated analysis of voltage-clamp recordings and computation of ion current/conductance surface density for these five cardiac ion currents.

RESULTS

The distributions from n = 54 hiPSC-CM in "ruptured" patch-clamp vs. n = 35 hiPSC-CM in β-escin-perforated patch-clamp were similar for membrane capacitance, access resistance, and ion current/conductance surface densities. The β-escin-perforated configuration resulted in improved stability of action potential (AP) shape and duration over a 10-min interval, with APD90 decay rate 0.7 ± 1.6%/min (mean ± SD, n = 4) vs. 4.6 ± 1.1%/min. (n = 3) for "ruptured" approach (p = 0.0286, one-tailed Mann-Whitney test).

DISCUSSION

The improved stability obtained here will allow development of CiPA-compliant automated patch-clamp assays on hiPSC-CM. Future applications include the study of multi ion-channel blocking properties of drugs using dynamic-clamp protocols, adding a valuable new tool to the arsenal of safety-pharmacology.

The influence of hERG1a and hERG1b isoforms on drug safety screening in iPSC-CMs

The human Ether-à-go-go Related Gene (hERG) encodes the pore forming subunit of the channel that conducts the rapid delayed rectifier potassium current I_Kr. I_Kr drives repolarization in the heart and when I_Kr is dysfunctional, cardiac repolarization delays, the QT interval on the electrocardiogram (ECG) prolongs and the risk of developing lethal arrhythmias such as Torsade de Pointes (TdP) increases. TdP risk is incorporated in drug safety screening for cardiotoxicity where hERG is the main target since the I_Kr channels appear highly sensitive to blockage. hERG block is also included as an important read-out in the Comprehensive in Vitro Proarrhythmia Assay (CiPA) initiative which aims to combine in vitro and in silico experiments on induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) to screen for cardiotoxicity. However, the hERG channel has some unique features to consider for drug safety screening, which we will discuss in this study. The hERG channel consists of different isoforms, hERG1a and hERG1b, which individually influence the kinetics of the channel and the drug response in the human heart and in iPSC-CMs. hERG1b is often underappreciated in iPSC-CM studies, drug screening assays and in silico models, and the fact that its contribution might substantially differ between iPSC-CM and healthy but also diseased human heart, adds to this problem. In this study we show that the activation kinetics in iPSC-CMs resemble hERG1b kinetics using Cs⁺ as a charge carrier. Not including hERG1b in drug safety testing might underestimate the actual role of hERG1b in repolarization and drug response, and might lead to inappropriate conclusions. We stress to focus more on including hERG1b in drug safety testing concerning I_Kr.

Protocol-Dependent Differences in IC50 Values Measured in Human Ether-Á-Go-Go-Related Gene Assays Occur in a Predictable Way and Can Be Used to Quantify State Preference of Drug Binding

Current guidelines around preclinical screening for drug-induced arrhythmias require the measurement of the potency of block of voltage-gated potassium channel subtype 11.1 (K_v11.1) as a surrogate for risk. A shortcoming of this approach is that the measured IC₅₀ of K_v11.1 block varies widely depending on the voltage protocol used in electrophysiological assays. In this study, we aimed to investigate the factors that contribute to these differences and to identify whether it is possible to make predictions about protocol-dependent block that might facilitate the comparison of potencies measured using different assays. Our data demonstrate that state preferential binding, together with drug-binding kinetics and trapping, is an important determinant of the protocol dependence of K_v11.1 block. We show for the first time that differences in IC₅₀ measured between protocols occurs in a predictable way, such that machine-learning algorithms trained using a selection of simple voltage protocols can indeed predict protocol-dependent potency. Furthermore, we also show that the preference of a drug for binding to the open versus the inactivated state of K_v11.1 can also be inferred from differences in IC₅₀ values measured between protocols. Our work therefore identifies how state preferential drug binding is a major determinant of the protocol dependence of IC₅₀ values measured in preclinical K_v11.1 assays. It also provides a novel method for quantifying the state dependence of K_v11.1 drug binding that will facilitate the development of more complete models of drug binding to K_v11.1 and improve our understanding of proarrhythmic risk associated with compounds that block K_v11.1.

Drug-induced Proarrhythmia and Torsade de Pointes: A Primer for Students and Practitioners of Medicine and Pharmacy

Multiple marketing withdrawals due to proarrhythmic concerns occurred in the United States, Canada, and the United Kingdom in the late 1980s to early 2000s. This primer reviews the clinical implications of a drug's identified proarrhythmic liability, the issues associated with these safety-related withdrawals, and the actions taken by the International Council for Harmonisation of Technical Requirements for Pharmaceuticals for Human Use (ICH) and by regulatory agencies in terms of changing drug development practices and introducing new nonclinical and clinical tests to asses proarrhythmic liability. ICH Guidelines S7B and E14 were released in 2005. Since then, they have been adopted by many regional regulatory authorities and have guided nonclinical and clinical proarrhythmic cardiac safety assessments during drug development. While this regulatory paradigm has been successful in preventing drugs with unanticipated potential for inducing the rare but potentially fatal polymorphic ventricular arrhythmia torsade de pointes from entering the market, it has led to the termination of drug development programs for other potentially useful medicines because of isolated results from studies with limited predictive value. Research efforts are now exploring alternative approaches to better predict potential proarrhythmic liabilities. For example, in the domain of human electrocardiographic assessments, concentration-response modeling conducted during phase 1 clinical development has recently become an accepted alternate primary methodology to the ICH E14 "thorough QT/QTc" study for defining a drug's corrected QT interval prolongation liability under certain conditions. When a drug's therapeutic benefit is considered important at a public health level but there is also an identified proarrhythmic liability that may result from administration of the single drug in certain individuals and/or drug-drug interactions, marketing approval will be accompanied by appropriate directions in the drug's prescribing information. Health-care professionals in the fields of medicine and pharmacy need to consider the prescribing information in conjunction with individual patients' clinical characteristics and concomitant medications when prescribing and dispensing such drugs.

Translating New Science Into the Drug Review Process: The US FDA's Division of Applied Regulatory Science

In 2011, the US Food and drug Administration (FDA) developed a strategic plan for regulatory science that focuses on developing new tools, standards, and approaches to assess the safety, efficacy, quality, and performance of FDA-regulated products. In line with this, the Division of Applied Regulatory Science was created to move new science into the Center for Drug Evaluation and Research (CDER) review process and close the gap between scientific innovation and drug review. The Division, located in the Office of Clinical Pharmacology, is unique in that it performs mission-critical applied research and review across the translational research spectrum including in vitro and in vivo laboratory research, in silico computational modeling and informatics, and integrated clinical research covering clinical pharmacology, experimental medicine, and postmarket analyses. The Division collaborates with Offices throughout CDER, across the FDA, other government agencies, academia, and industry. The Division is able to rapidly form interdisciplinary teams of pharmacologists, biologists, chemists, computational scientists, and clinicians to respond to challenging regulatory questions for specific review issues and for longer-range projects requiring the development of predictive models, tools, and biomarkers to speed the development and regulatory evaluation of safe and effective drugs. This article reviews the Division's recent work and future directions, highlighting development and validation of biomarkers; novel humanized animal models; translational predictive safety combining in vitro, in silico, and in vivo clinical biomarkers; chemical and biomedical informatics tools for safety predictions; novel approaches to speed the development of complex generic drugs, biosimilars, and antibiotics; and precision medicine.

Digging into Lipid Membrane Permeation for Cardiac Ion Channel Blocker d-Sotalol with All-Atom Simulations

Interactions of drug molecules with lipid membranes play crucial role in their accessibility of cellular targets and can be an important predictor of their therapeutic and safety profiles. Very little is known about spatial localization of various drugs in the lipid bilayers, their active form (ionization state) or translocation rates and therefore potency to bind to different sites in membrane proteins. All-atom molecular simulations may help to map drug partitioning kinetics and thermodynamics, thus providing in-depth assessment of drug lipophilicity. As a proof of principle, we evaluated extensively lipid membrane partitioning of d-sotalol, well-known blocker of a cardiac potassium channel K_v11.1 encoded by the hERG gene, with reported substantial proclivity for arrhythmogenesis. We developed the positively charged (cationic) and neutral d-sotalol models, compatible with the biomolecular CHARMM force field, and subjected them to all-atom molecular dynamics (MD) simulations of drug partitioning through hydrated lipid membranes, aiming to elucidate thermodynamics and kinetics of their translocation and thus putative propensities for hydrophobic and aqueous hERG access. We found that only a neutral form of d-sotalol accumulates in the membrane interior and can move across the bilayer within millisecond time scale, and can be relevant to a lipophilic channel access. The computed water-membrane partitioning coefficient for this form is in good agreement with experiment. There is a large energetic barrier for a cationic form of the drug, dominant in water, to cross the membrane, resulting in slow membrane translocation kinetics. However, this form of the drug can be important for an aqueous access pathway through the intracellular gate of hERG. This route will likely occur after a neutral form of a drug crosses the membrane and subsequently re-protonates. Our study serves to demonstrate a first step toward a framework for multi-scale in silico safety pharmacology, and identifies some of the challenges that lie therein.

Ion Channels in Drug Discovery and Safety Pharmacology

Ion channels are membrane proteins involved in almost all physiological processes, including neurotransmission, muscle contraction, pace-making activity, secretion, electrolyte and water balance, immune response, and cell proliferation. Due to their broad distribution in human body and physiological roles, ion channels are attractive targets for drug discovery and safety pharmacology. Over the years ion channels have been associated to many genetic diseases ("channelopathies"). For most of these diseases the therapy is mainly empirical and symptomatic, often limited by lack of efficacy and tolerability for a number of patients. The search for the development of new and more specific therapeutic approaches is therefore strongly pursued. At the same time acquired channelopathies or dangerous side effects (such as proarrhythmic risk) can develop as a consequence of drugs unexpectedly targeting ion channels. Several noncardiovascular drugs are known to block cardiac ion channels, leading to potentially fatal delayed ventricular repolarization. Thus, the search of reliable preclinical cardiac safety testing in early stage of drug discovery is mandatory. To fulfill these needs, both ion channels drug discovery and toxicology strategies are evolving toward comprehensive research approaches integrating ad hoc designed in silico predictions and experimental studies for a more reliable and quick translation of results to the clinic side.Here we discuss two examples of how the combination of in silico methods and patch clamp experiments can help addressing drug discovery and safety issues regarding ion channels.

Mechanistic Model-Informed Proarrhythmic Risk Assessment of Drugs: Review of the "CiPA" Initiative and Design of a Prospective Clinical Validation Study

The Comprehensive in vitro Proarrhythmia Assay (CiPA) initiative is developing and validating a mechanistic-based assessment of the proarrhythmic risk of drugs. CiPA proposes to assess a drug's effect on multiple ion channels and integrate the effects in a computer model of the human cardiomyocyte to predict proarrhythmic risk. Unanticipated or missed effects will be assessed with human stem cell-derived cardiomyocytes and electrocardiogram (ECG) analysis in early phase I clinical trials. This article provides an overview of CiPA and the rationale and design of the CiPA phase I ECG validation clinical trial, which involves assessing an additional ECG biomarker (J-Tpeak) for QT prolonging drugs. If successful, CiPA will 1) create a pathway for drugs with hERG block / QT prolongation to advance without intensive ECG monitoring in phase III trials if they have low proarrhythmic risk; and 2) enable updating drug labels to be more informative about proarrhythmic risk, not just QT prolongation.

Adult Human Primary Cardiomyocyte-Based Model for the Simultaneous Prediction of Drug-Induced Inotropic and Pro-arrhythmia Risk

Cardiac safety remains the leading cause of drug development discontinuation. We developed a human cardiomyocyte-based model that has the potential to provide a predictive preclinical approach for simultaneously predicting drug-induced inotropic and pro-arrhythmia risk.

Methods: Adult human primary cardiomyocytes from ethically consented organ donors were used to measure contractility transients. We used measures of changes in contractility parameters as markers to infer both drug-induced inotropic effect (sarcomere shortening) and pro-arrhythmia (aftercontraction, AC); contractility escape (CE); time to 90% relaxation (TR90). We addressed the clinical relevance of this approach by evaluating the effects of 23 torsadogenic and 10 non-torsadogenic drugs. Each drug was tested separately at four multiples of the free effective therapeutic plasma concentration (fETPC).

Results: Human cardiomyocyte-based model differentiated between torsadogenic and non-torsadogenic drugs. For example, dofetilide, a torsadogenic drug, caused ACs and increased TR90 starting at 10-fold the fETPC, while CE events were observed at the highest multiple of fETPC (100-fold). Verapamil, a non-torsadogenic drug, did not change TR90 and induced no AC or CE up to the highest multiple of fETPCs tested in this study (222-fold). When drug pro-arrhythmic activity was evaluated at 10-fold of the fETPC, AC parameter had excellent assay sensitivity and specificity values of 96 and 100%, respectively. This high predictivity supports the translational safety potential of this preparation and of the selected marker. The data demonstrate that human cardiomyocytes could also identify drugs associated with inotropic effects. hERG channel blockers, like dofetilide, had no effects on sarcomere shortening, while multi-ion channel blockers, like verapamil, inhibited sarcomere shortening.

Conclusions: Isolated adult human primary cardiomyocytes can simultaneously predict risks associated with inotropic activity and pro-arrhythmia and may enable the generation of reliable and predictive data for assessing human cardiotoxicity at an early stage in drug discovery.

In Vitro and In Silico Risk Assessment in Acquired Long QT Syndrome: The Devil Is in the Details

Acquired long QT syndrome, mostly as a result of drug block of the Kv11. 1 potassium channel in the heart, is characterized by delayed cardiac myocyte repolarization, prolongation of the T interval on the ECG, syncope and sudden cardiac death due to the polymorphic ventricular arrhythmia Torsade de Pointes (TdP). In recent years, efforts are underway through the Comprehensive in vitro proarrhythmic assay (CiPA) initiative, to develop better tests for this drug induced arrhythmia based in part on in silico simulations of pharmacological disruption of repolarization. However, drug binding to Kv11.1 is more complex than a simple binary molecular reaction, meaning simple steady state measures of potency are poor surrogates for risk. As a result, there is a plethora of mechanistic detail describing the drug/Kv11.1 interaction—such as drug binding kinetics, state preference, temperature dependence and trapping—that needs to be considered when developing in silico models for risk prediction. In addition to this, other factors, such as multichannel pharmacological profile and the nature of the ventricular cell models used in simulations also need to be considered in the search for the optimum in silico approach. Here we consider how much of mechanistic detail needs to be included for in silico models to accurately predict risk and further, how much of this detail can be retrieved from protocols that are practical to implement in high throughout screens as part of next generation of preclinical in silico drug screening approaches?

Quantitative Comparison of Effects of Dofetilide, Sotalol, Quinidine, and Verapamil between Human Ex vivo Trabeculae and In silico Ventricular Models Incorporating Inter-Individual Action Potential Variability

Background:In silico modeling could soon become a mainstream method of pro-arrhythmic risk assessment in drug development. However, a lack of human-specific data and appropriate modeling techniques has previously prevented quantitative comparison of drug effects between in silico models and recordings from human cardiac preparations. Here, we directly compare changes in repolarization biomarkers caused by dofetilide, dl-sotalol, quinidine, and verapamil, between in silico populations of human ventricular cell models and ex vivo human ventricular trabeculae. Methods and Results:Ex vivo recordings from human ventricular trabeculae in control conditions were used to develop populations of in silico human ventricular cell models that integrated intra- and inter-individual variability in action potential (AP) biomarker values. Models were based on the O'Hara-Rudy ventricular cardiomyocyte model, but integrated experimental AP variability through variation in underlying ionic conductances. Changes to AP duration, triangulation and early after-depolarization occurrence from application of the four drugs at multiple concentrations and pacing frequencies were compared between simulations and experiments. To assess the impact of variability in IC50 measurements, and the effects of including state-dependent drug binding dynamics, each drug simulation was repeated with two different IC50 datasets, and with both the original O'Hara-Rudy hERG model and a recently published state-dependent model of hERG and hERG block. For the selective hERG blockers dofetilide and sotalol, simulation predictions of AP prolongation and repolarization abnormality occurrence showed overall good agreement with experiments. However, for multichannel blockers quinidine and verapamil, simulations were not in agreement with experiments across all IC50 datasets and IKr block models tested. Quinidine simulations resulted in overprolonged APs and high incidence of repolarization abnormalities, which were not observed in experiments. Verapamil simulations showed substantial AP prolongation while experiments showed mild AP shortening. Conclusions: Results for dofetilide and sotalol show good agreement between experiments and simulations for selective compounds, however lack of agreement from simulations of quinidine and verapamil suggest further work is needed to understand the more complex electrophysiological effects of these multichannel blocking drugs.

Using automated patch clamp electrophysiology platforms in pain-related ion channel research: insights from industry and academia

Automated patch clamp (APC) technology was first developed at the turn of the millennium. The increased throughput it afforded promised a new paradigm in ion channel recordings, offering the potential to overcome the time-consuming, low-throughput bottleneck, arising from manual patch clamp investigations. This has relevance to the fast-paced development of novel therapies for chronic pain. This review highlights the advances in technology, using select examples that have facilitated APC usage in both industry and academia. It covers both first generation and the latest developments in second-generation platforms. In addition, it also provides an overview of the pain research field and how APC platforms have furthered our understanding of ion channel research and the development of pharmacological tools and therapeutics. APC platforms have much to offer to the ion channel research community, and this review highlights areas of best practice for both academia and industry. The impact of APC platforms and the prospects of ion channel research and improved therapeutics for chronic pain will be evaluated.

LINKED ARTICLES

This article is part of a themed section on Recent Advances in Targeting Ion Channels to Treat Chronic Pain. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v175.12/issuetoc.

Computational approaches to understand the adverse drug effect on potassium, sodium and calcium channels for predicting TdP cardiac arrhythmias

Ion channels play a crucial role in the cardiovascular system. Our understanding of cardiac ion channel function has improved since their first discoveries. The flow of potassium, sodium and calcium ions across cardiomyocytes is vital for regular cardiac rhythm. Blockage of these channels, delays cardiac repolarization or tend to shorten repolarization and may induce arrhythmia. Detection of drug risk by channel blockade is considered essential for drug regulators. Advanced computational models can be used as an early screen for torsadogenic potential in drug candidates. New drug candidates that are determined to not cause blockage are more likely to pass successfully through preclinical trials and not be withdrawn later from the marketplace by manufacturer. Several different approved drugs, however, can cause a distinctive polymorphic ventricular arrhythmia known as torsade de pointes (TdP), which may lead to sudden death. The objective of the present study is to review the mechanisms and computational models used to assess the risk that a drug may TdP.

KEY POINTS

There is strong evidence from multiple studies that blockage of the L-type calcium current reduces risk of TdP. Blockage of sodium channels slows cardiac action potential conduction, however, not all sodium channel blocking antiarrhythmic drugs produce a significant effect, while late sodium channel block reduces TdP. Interestingly, there are some drugs that block the hERG potassium channel and therefore cause QT prolongation, but they are not associated with TdP. Recent studies confirmed the necessity of studying multiple distinctionic ion channels which are responsible for cardiac related diseases or TdP, to obtain an improved clinical TdP risk prediction of compound interactions and also for designing drugs.

Safety pharmacology methods and models in an evolving regulatory environment

This editorial prefaces the annual themed issue on safety pharmacology (SP) methods published in the Journal of Pharmacological and Toxicological Methods (JPTM). We highlight here the content derived from the recent 2016 Safety Pharmacology Society (SPS), Canadian Society of Pharmacology and Therapeutics (CSPT), and Japanese Safety Pharmacology Society (JSPS) joint meeting held in Vancouver, B.C., Canada. This issue of JPTM continues the tradition of providing a publication summary of articles primarily presented at the joint meeting with direct bearing on the discipline of SP. As the regulatory landscape is expected to evolve with revision announced for the existing guidance document on non-clinical proarrhythmia risk assessment (ICHS7B) there is also imminent inception of the Comprehensive in vitro Proarrhythmia Assay (CiPA) initiative. Thus, the field of SP is dynamically progressing with characterization and implementation of numerous alternative non-clinical safety models. Novel method development and refinement in all areas of the discipline are reflected in the content.