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.

Advances in ion channel high throughput screening: where are we in 2023?

INTRODUCTION

Automated Patch Clamp (APC) technology has become an integral element in ion channel research, drug discovery and development pipelines to overcome the use of the highly time-consuming manual patch clamp (MPC) procedures. This automated technology offers increased throughput and promises a new model in obtaining ion channel recordings, which has significant relevance to the development of novel therapies and safety profiling of candidate therapeutic compounds.

AREAS COVERED

This article reviews the recent innovations in APC technology, including platforms, and highlights how they have facilitated usage in both industry and academia. The review also provides an overview of the ion channel research endeavors and how APC platforms have contributed to the understanding of ion channel research, pharmacological tools and therapeutics. Furthermore, the authors provide their opinion on the challenges and goals for APC technology going forward to accelerate academic research and drug discovery across a host of therapeutic areas.

EXPERT OPINION

It is clear that APC technology has progressed drug discovery programs, specifically in the field of neuroscience and cardiovascular research. The challenge for the future is to keep pace with fundamental research and improve translation of the large datasets obtained.

The Strength of hERG Inhibition by Erythromycin at Different Temperatures Might Be Due to Its Interacting Features with the Channels

Erythromycin is one of the few compounds that remarkably increase ether-a-go-go-related gene (hERG) inhibition from room temperature (RT) to physiological temperature (PT). Understanding how erythromycin inhibits the hERG could help us to decide which compounds are needed for further studies. The whole-cell patch clamp technique was used to investigate the effects of erythromycin on hERG channels at different temperatures. While erythromycin caused a concentration-dependent inhibition of cardiac hERG channels, it also shifted the steady-state activation and steady-state inactivation of the channel to the left and significantly accelerated the onset of inactivation at both temperatures, although temperature itself caused a profound change in the dynamics of hERG channels. Our data also suggest that the binding pattern to S6 of the channels changes at PT. In contrast, cisapride, a well-known hERG blocker whose inhibition is not affected by temperature, does not change its critical binding sites after the temperature is raised to PT. Our data suggest that erythromycin is unique and that the shift in hERG inhibition may not apply to other compounds.

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.

Metabolic and electrolyte abnormalities as risk factors in drug-induced long QT syndrome

Drug-induced long QT syndrome (diLQTS) is the phenomenon by which the administration of drugs causes prolongation of cardiac repolarisation and leads to an increased risk of the ventricular tachycardia known as torsades de pointes (TdP). In most cases of diLQTS, the primary molecular target is the human ether-à-go-go-related gene protein (hERG) potassium channel, which carries the rapid delayed rectifier current (I_Kr) in the heart. However, the proarrhythmic risk associated with drugs that block hERG can be modified in patients by a range of environmental- and disease-related factors, such as febrile temperatures, alterations in pH, dyselectrolytaemias such as hypokalaemia and hypomagnesemia and coadministration with other drugs. In this review, we will discuss the clinical occurrence of drug-induced LQTS in the context of these modifying factors as well as the mechanisms by which they contribute to altered hERG potency and proarrhythmic risk.

Pathophysiological metabolic changes associated with disease modify the proarrhythmic risk profile of drugs with potential to prolong repolarisation

BACKGROUND AND PURPOSE

Hydroxychloroquine, chloroquine and azithromycin are three drugs that were proposed to treat COVID-19. While concern already existed around their proarrhythmic potential there is little data regarding how altered physiological states encountered in patients such as febrile state, electrolyte imbalances or acidosis might change their risk profiles.

EXPERIMENTAL APPROACH

Potency of hERG block was measured using high-throughput electrophysiology in the presence of variable environmental factors. These potencies informed simulations to predict population risk profiles. Effects on cardiac repolarisation were verified in human induced pluripotent stem cell-derived cardiomyocytes from multiple individuals.

KEY RESULTS

Chloroquine and hydroxychloroquine blocked hERG with IC₅₀ of 1.47±0.07 μM and 3.78±0.17 μM respectively, indicating proarrhythmic risk at concentrations effective against SARS-CoV-2 in vitro. Hypokalaemia and hypermagnesemia increased potency of chloroquine and hydroxychloroquine, indicating increased proarrhythmic risk. Acidosis significantly reduced potency of all drugs, whereas increased temperature decreased potency of chloroquine and hydroxychloroquine against hERG but increased potency for azithromycin. In silico simulations demonstrated that proarrhythmic risk was increased by female sex, hypokalaemia and heart failure, and identified specific genetic backgrounds associated with emergence of arrhythmia.

CONCLUSION AND IMPLICATIONS

Our study demonstrates how proarrhythmic risk can be exacerbated by metabolic changes and pre-existing disease. More broadly, the study acts as a blueprint for how high-throughput in vitro screening, combined with in silico simulations can help guide both preclinical screening and clinical management of patients in relation to drugs with potential to prolong repolarisation.

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.

Chloroquine and hydroxychloroquine provoke arrhythmias at concentrations higher than those clinically used to treat COVID-19: A simulation study

Abstract

The risk of fatal arrhythmias is the major concern for using chloroquine (CQ) or hydroxychloroquine (HCQ) to treat coronavirus disease 2019 (COVID‐19), but the reported number of life‐threatening arrhythmic events or deaths is relatively small. The objective of this study was to assess the arrhythmogenic risk of these two drugs using a multiscale heart simulation, which allows testing even at high concentrations, including those that cause fatal arrhythmias. We measured the inhibitory action of CQ, HCQ, and HCQ with 30 μM azithromycin (AZ) on six ion currents (fast [INa] and late [INa,L] components of the sodium current, L‐type calcium current [ICa,L], rapid [IKr/hERG], and slow [IKs] components of delayed rectifier potassium, and inward rectifier potassium [IK1]) over a wide range of concentrations using the automated patch‐clamp system. Using the concentration–inhibition relationship that was thus obtained, we simulated the drug effects while increasing the concentration until the life‐threatening arrhythmia, torsade de pointes (TdP), was observed. The obtained threshold concentrations for TdP were 12.5, 35, and 22.5 μM for CQ, HCQ, and HCQ with AZ, respectively. Adding therapeutic concentrations of mexiletine or verapamil successfully prevented the occurrence of TdP, and verapamil was more effective. CQ, HCQ, and HCQ with AZ thresholds for TdP were larger than both antiviral concentrations that were reported by in vitro experiments and free plasma concentrations that were attained by the clinically used dosage. The current simulation data provided a safety margin to the currently used clinical dose for CQ and HCQ/AZ.

Toward in vivo-relevant hERG safety assessment and mitigation strategies based on relationships between non-equilibrium blocker binding, three-dimensional channel-blocker interactions, dynamic occupancy, dynamic exposure, and cellular arrhythmia

The human ether-a-go-go-related voltage-gated cardiac ion channel (commonly known as hERG) conducts the rapid outward repolarizing potassium current in cardiomyocytes (IKr). Inadvertent blockade of this channel by drug-like molecules represents a key challenge in pharmaceutical R&D due to frequent overlap between the structure-activity relationships of hERG and many primary targets. Building on our previous work, together with recent cryo-EM structures of hERG, we set about to better understand the energetic and structural basis of promiscuous blocker-hERG binding in the context of Biodynamics theory. We propose a two-step blocker binding process consisting of: The initial capture step: diffusion of a single fully solvated blocker copy into a large cavity lined by the intra-cellular cyclic nucleotide binding homology domain (CNBHD). Occupation of this cavity is a necessary but insufficient condition for ion current disruption.The IKr disruption step: translocation of the captured blocker along the channel axis, such that: The head group, consisting of a quasi-rod-shaped moiety, projects into the open pore, accompanied by partial de-solvation of the binding interface.One tail moiety packs along a kink between the S6 helix and proximal C-linker helix adjacent to the intra-cellular entrance of the pore, likewise accompanied by mutual de-solvation of the binding interface (noting that the association barrier is comprised largely of the total head + tail group de-solvation cost).Blockers containing a highly planar moiety that projects into a putative constriction zone within the closed channel become trapped upon closing, as do blockers terminating prior to this region.A single captured blocker copy may conceivably associate and dissociate to/from the pore many times before exiting the CNBHD cavity. Lastly, we highlight possible flaws in the current hERG safety index (SI), and propose an alternate in vivo-relevant strategy factoring in: Benefit/risk.The predicted arrhythmogenic fractional hERG occupancy (based on action potential (AP) simulations of the undiseased human ventricular cardiomyocyte).Alteration of the safety threshold due to underlying disease.Risk of exposure escalation toward the predicted arrhythmic limit due to patient-to-patient pharmacokinetic (PK) variability, drug-drug interactions, overdose, and use for off-label indications in which the hERG safety parameters may differ from their on-label counterparts.

Cross-site and cross-platform variability of automated patch clamp assessments of drug effects on human cardiac currents in recombinant cells

Automated patch clamp (APC) instruments enable efficient evaluation of electrophysiologic effects of drugs on human cardiac currents in heterologous expression systems. Differences in experimental protocols, instruments, and dissimilar site procedures affect the variability of IC50 values characterizing drug block potency. This impacts the utility of APC platforms for assessing a drug's cardiac safety margin. We determined variability of APC data from multiple sites that measured blocking potency of 12 blinded drugs (with different levels of proarrhythmic risk) against four human cardiac currents (hERG [IKr], hCav1.2 [L-Type ICa], peak hNav1.5, [Peak INa], late hNav1.5 [Late INa]) with recommended protocols (to minimize variance) using five APC platforms across 17 sites. IC50 variability (25/75 percentiles) differed for drugs and currents (e.g., 10.4-fold for dofetilide block of hERG current and 4-fold for mexiletine block of hNav1.5 current). Within-platform variance predominated for 4 of 12 hERG blocking drugs and 4 of 6 hNav1.5 blocking drugs. hERG and hNav1.5 block. Bland-Altman plots depicted varying agreement across APC platforms. A follow-up survey suggested multiple sources of experimental variability that could be further minimized by stricter adherence to standard protocols. Adoption of best practices would ensure less variable APC datasets and improved safety margins and proarrhythmic risk assessments.

Evaluation of Possible Proarrhythmic Potency: Comparison of the Effect of Dofetilide, Cisapride, Sotalol, Terfenadine, and Verapamil on hERG and Native IKr Currents and on Cardiac Action Potential

The proarrhythmic potency of drugs is usually attributed to the IKr current block. During safety pharmacology testing analysis of IKr in cardiomyocytes was replaced by human ether-a-go-go-related gene (hERG) test using automated patch-clamp systems in stable transfected cell lines. Aim of this study was to compare the effect of proarrhythmic compounds on hERG and IKr currents and on cardiac action potential. The hERG current was measured by using both automated and manual patch-clamp methods on HEK293 cells. The native ion currents (IKr, INaL, ICaL) were recorded from rabbit ventricular myocytes by manual patch-clamp technique. Action potentials in rabbit ventricular muscle and undiseased human donor hearts were studied by conventional microelectrode technique. Dofetilide, cisapride, sotalol, terfenadine, and verapamil blocked hERG channels at 37°C with an IC50 of 7 nM, 18 nM, 343 μM, 165 nM, and 214 nM, respectively. Using manual patch-clamp, the IC50 values of sotalol and terfenadine were 78 µM and 31 nM, respectively. The IC50 values calculated from IKr measurements at 37°C were 13 nM, 26 nM, 52 μM, 54 nM, and 268 nM, respectively. Cisapride, dofetilide, and sotalol excessively lengthened, terfenadine, and verapamil did not influence the action potential duration. Terfenadine significantly inhibited INaL and moderately ICaL, verapamil blocked only ICaL. Automated hERG assays may over/underestimate proarrhythmic risk. Manual patch-clamp has substantially higher sensitivity to certain drugs. Action potential studies are also required to analyze complex multichannel effects. Therefore, manual patch-clamp and action potential experiments should be a part of preclinical safety tests.

General Principles for the Validation of Proarrhythmia Risk Prediction Models: An Extension of the CiPA In Silico Strategy

This white paper presents principles for validating proarrhythmia risk prediction models for regulatory use as discussed at the In Silico Breakout Session of a Cardiac Safety Research Consortium/Health and Environmental Sciences Institute/US Food and Drug Administration-sponsored Think Tank Meeting on May 22, 2018. The meeting was convened to evaluate the progress in the development of a new cardiac safety paradigm, the Comprehensive in Vitro Proarrhythmia Assay (CiPA). The opinions regarding these principles reflect the collective views of those who participated in the discussion of this topic both at and after the breakout session. Although primarily discussed in the context of in silico models, these principles describe the interface between experimental input and model-based interpretation and are intended to be general enough to be applied to other types of nonclinical models for proarrhythmia assessment. This document was developed with the intention of providing a foundation for more consistency and harmonization in developing and validating different models for proarrhythmia risk prediction using the example of the CiPA paradigm.

Rapid Characterization of hERG Channel Kinetics II: Temperature Dependence

Ion channel behavior can depend strongly on temperature, with faster kinetics at physiological temperatures leading to considerable changes in currents relative to room temperature. These temperature-dependent changes in voltage-dependent ion channel kinetics (rates of opening, closing, inactivating, and recovery) are commonly represented with Q₁₀ coefficients or an Eyring relationship. In this article, we assess the validity of these representations by characterizing channel kinetics at multiple temperatures. We focus on the human Ether-à-go-go-Related Gene (hERG) channel, which is important in drug safety assessment and commonly screened at room temperature so that results require extrapolation to physiological temperature. In Part I of this study, we established a reliable method for high-throughput characterization of hERG1a (Kv11.1) kinetics, using a 15-second information-rich optimized protocol. In this Part II, we use this protocol to study the temperature dependence of hERG kinetics using Chinese hamster ovary cells overexpressing hERG1a on the Nanion SyncroPatch 384PE, a 384-well automated patch-clamp platform, with temperature control. We characterize the temperature dependence of hERG gating by fitting the parameters of a mathematical model of hERG kinetics to data obtained at five distinct temperatures between 25 and 37°C and validate the models using different protocols. Our models reveal that activation is far more temperature sensitive than inactivation, and we observe that the temperature dependency of the kinetic parameters is not represented well by Q₁₀ coefficients; it broadly follows a generalized, but not the standardly-used, Eyring relationship. We also demonstrate that experimental estimations of Q₁₀ coefficients are protocol dependent. Our results show that a direct fit using our 15-s protocol best represents hERG kinetics at any given temperature and suggests that using the Generalized Eyring theory is preferable if no experimental data are available to derive model parameters at a given temperature.

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.

Experimentally Validated Pharmacoinformatics Approach to Predict hERG Inhibition Potential of New Chemical Entities

The hERG (human ether-a-go-go-related gene) encoded potassium ion (K⁺) channel plays a major role in cardiac repolarization. Drug-induced blockade of hERG has been a major cause of potentially lethal ventricular tachycardia termed Torsades de Pointes (TdPs). Therefore, we presented a pharmacoinformatics strategy using combined ligand and structure based models for the prediction of hERG inhibition potential (IC₅₀) of new chemical entities (NCEs) during early stages of drug design and development. Integrated GRid-INdependent Descriptor (GRIND) models, and lipophilic efficiency (LipE), ligand efficiency (LE) guided template selection for the structure based pharmacophore models have been used for virtual screening and subsequent hERG activity (pIC₅₀) prediction of identified hits. Finally selected two hits were experimentally evaluated for hERG inhibition potential (pIC₅₀) using whole cell patch clamp assay. Overall, our results demonstrate a difference of less than ±1.6 log unit between experimentally determined and predicted hERG inhibition potential (IC₅₀) of the selected hits. This revealed predictive ability and robustness of our models and could help in correctly rank the potency order (lower μM to higher nM range) against hERG.