Machine learning and deep learning approaches for enhanced prediction of hERG blockade: a comprehensive QSAR modeling study.

BACKGROUND: Cardiotoxicity is a major cause of drug withdrawal. The hERG channel, regulating ion flow, is pivotal for heart and nervous system function. Its blockade is a concern in drug development. Predicting hERG blockade is essential for identifying cardiac safety issues. Various QSAR models exist, but their performance varies. Ongoing improvements show promise, necessitating continued efforts to enhance accuracy using emerging deep learning algorithms in predicting potential hERG blockade. STUDY DESIGN AND METHOD: Using a large training dataset, six individual QSAR models were developed. Additionally, three ensemble models were constructed. All models were evaluated using 10-fold cross-validations and two external datasets. RESULTS: The 10-fold cross-validations resulted in Mathews correlation coefficient (MCC) values from 0.682 to 0.730, surpassing the best-reported model on the same dataset (0.689). External validations yielded MCC values from 0.520 to 0.715 for the first dataset, exceeding those of previously reported models (0-0.599). For the second dataset, MCC values fell between 0.025 and 0.215, aligning with those of reported models (0.112-0.220). CONCLUSIONS: The developed models can assist the pharmaceutical industry and regulatory agencies in predicting hERG blockage activity, thereby enhancing safety assessments and reducing the risk of adverse cardiac events associated with new drug candidates.

Rescue of expression and function of long QT syndrome-causing mutant hERG channels by enhancing channel stability in the plasma membrane.

The human ether-g-go-go-related gene (hERG) encodes the Kv11.1 (or hERG) channel that conducts the rapidly activating delayed rectifier potassium current (IKr). Naturally occurring mutations in hERG impair the channel function and cause long QT syndrome type 2 (LQT2). Many missense hERG mutations lead to a lack of channel expression on the cell surface, representing a major mechanism for the loss-of-function of mutant channels. While it is generally thought that a trafficking defect underlies the lack of channel expression on the cell surface, in the present study, we demonstrate that the trafficking defective mutant hERG G601S can reach the plasma membrane but is unstable and quickly degrades, which is akin to Wild Type (WT) hERG channels under low K+ conditions. We previously showed that Serine (S) residue at 624 in the innermost position of the selectivity filter of hERG is involved in hERG membrane stability such that substitution of Serine 624 with Threonine (S624T) enhances hERG stability and renders hERG insensitive to low K+ culture. Here, we report that the intragenic addition of S624T substitution to trafficking defective hERG mutants G601S, N470D and P596R led to a complete rescue of the function of these otherwise loss-of-function mutant channels to a level similar to the WT channel, representing the most effective rescue means for the function of mutant hERG channels. These findings not only provide novel insights into hERG mutation-mediated channel dysfunction, but also point to the critical role of S624 in hERG stability on the plasma membrane.

Analysis of cardiohemodynamic and electrophysiological effects of morphine along with its toxicokinetic profile using the halothane-anesthetized dogs.

Although morphine has been used for treatment-resistant dyspnea in end-stage heart failure patients, information on its cardiovascular safety profile remains limited. Morphine was intravenously administered to halothane-anesthetized dogs (n=4) in doses of 0.1, 1 and 10 mg/kg/10 min with 20 min of observation period. The low and middle doses attained therapeutic (0.13 µg/mL) and supratherapeutic (0.97 µg/mL) plasma concentrations, respectively. The low dose hardly altered any of the cardiovascular variables except that the QT interval was prolonged for 10-15 min after its start of infusion. The middle dose reduced the preload and afterload to the left ventricle for 5-15 min, then decreased the left ventricular contractility and mean blood pressure for 10-30 min, and finally suppressed the heart rate for 15-30 min. Moreover, the middle dose gradually but progressively prolonged the atrioventricular conduction time, QT interval/QTcV, ventricular late repolarization period and ventricular effective refractory period without altering the intraventricular conduction time, ventricular early repolarization period or terminal repolarization period. A reverse-frequency-dependent delay of ventricular repolarization was confirmed. The high dose induced cardiohemodynamic collapse mainly due to vasodilation in the initial 2 animals by 1.9 and 3.3 min after its start of infusion, respectively, which needed circulatory support to treat. The high dose was not tested further in the remaining 2 animals. Thus, intravenously administered morphine exerts a rapidly appearing vasodilator action followed by slowly developing cardiosuppressive effects. Morphine can delay the ventricular repolarization possibly through IKr inhibition in vivo, but its potential to develop torsade de pointes will be small.

Molecular Pathways and Animal Models of Arrhythmias.

Arrhythmias account for over 300,000 annual deaths in the United States, and approximately half of all deaths are associated with heart disease. Mechanisms underlying arrhythmia risk are complex; however, work in humans and animal models over the past 25 years has identified a host of molecular pathways linked with both arrhythmia substrates and triggers. This chapter will focus on select arrhythmia pathways solved by linking human clinical and genetic data with animal models.

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.

Use of Surface Plasmon Resonance Technique for Studies of Inter-domain Interactions in Ion Channels.

Ion channels are transmembrane proteins essential for cellular functions and are important drug targets. Surface plasmon resonance (SPR) is a powerful technique for investigating protein-protein and protein-small molecule ligand interactions. SPR has been underutilized for studies of ion channels, even though it could provide a wealth of information on the mechanisms of ion channel regulation and aid in ion channel drug discovery. Here we provide a detailed description of the use of SPR technology for investigating inter-domain interactions in KCNH potassium-selective and voltage-gated ion channels.

Supporting an integrated QTc risk assessment using the hERG margin distributions for three positive control agents derived from multiple laboratories and on multiple occasions.

BACKGROUND: Determination of a drug's potency in blocking the hERG channel is an established safety pharmacology study. Best practice guidelines have been published for reliable assessment of hERG potency. In addition, a set of plasma concentration and plasma protein binding fraction data were provided as denominators for margin calculations. The aims of the current analysis were five-fold: provide data allowing creation of consistent denominators for the hERG margin distributions of the key reference agents, explore the variation in hERG margins within and across laboratories, provide a hERG margin to 10 ms QTc prolongation based on several newer studies, provide information to use these analyses for reference purposes, and provide recommended hERG margin 'cut-off' values. METHODS: The analyses used 12 hERG IC50 'best practice' data sets (for the 3 reference agents). A group of 5 data sets came from a single laboratory. The other 7 data sets were collected by 6 different laboratories. RESULTS: The denominator exposure distributions were consistent with the ICH E14/S7B Training Materials. The inter-occasion and inter-laboratory variability in hERG IC50 values were comparable. Inter-drug differences were most important in determining the pooled margin variability. The combined data provided a robust hERG margin reference based on best practice guidelines and consistent exposure denominators. The sensitivity of hERG margin thresholds were consistent with the sensitivity described over the course of the last two decades. CONCLUSION: The current data provide further insight into the sensitivity of the 30-fold hERG margin 'cut-off' used for two decades. Using similar hERG assessments and these analyses, a future researcher can use a hERG margin threshold to support a negative QTc integrated risk assessment.

The proteostasis interactomes of trafficking-deficient variants of the voltage-gated potassium channel KV11.1 associated with long QT syndrome.

The voltage-gated potassium ion channel KV11.1 plays a critical role in cardiac repolarization. Genetic variants that render Kv11.1 dysfunctional cause long QT syndrome (LQTS), which is associated with fatal arrhythmias. Approximately 90% of LQTS-associated variants cause intracellular protein transport (trafficking) dysfunction, which pharmacological chaperones like E-4031 can rescue. Protein folding and trafficking decisions are regulated by chaperones, protein quality control factors, and trafficking machinery comprising the cellular proteostasis network. Here, we test whether trafficking dysfunction is associated with alterations in the proteostasis network of pathogenic Kv11.1 variants and whether pharmacological chaperones can normalize the proteostasis network of responsive variants. We used affinity-purification coupled with tandem mass tag-based quantitative mass spectrometry to assess protein interaction changes of WT KV11.1 or trafficking-deficient channel variants in the presence or absence of E-4031. We identified 572 core KV11.1 protein interactors. Trafficking-deficient variants KV11.1-G601S and KV11.1-G601S-G965∗ had significantly increased interactions with proteins responsible for folding, trafficking, and degradation compared to WT. We confirmed previous findings that the proteasome is critical for KV11.1 degradation. Our report provides the first comprehensive characterization of protein quality control mechanisms of KV11.1. We find extensive interactome remodeling associated with trafficking-deficient KV11.1 variants and with pharmacological chaperone rescue of KV11.1 cell surface expression. The identified protein interactions could be targeted therapeutically to improve KV11.1 trafficking and treat LQTS.

Macrocyclization strategy for improving candidate profiles in medicinal chemistry.

Macrocycles are defined as cyclic compounds with 12 or more members. In medicinal chemistry, they are categorized based on their core chemistry into cyclic peptides and macrocycles. Macrocycles are advantageous because of their structural diversity and ability to achieve high affinity and selectivity towards challenging targets that are often not addressable by conventional small molecules. The potential of macrocyclization to optimize drug-like properties while maintaining adequate bioavailability and permeability has been emphasized as a key innovation in medicinal chemistry. This review provides a detailed case study of the application of macrocyclization over the past 5 years, starting from the initial analysis of acyclic active compounds to optimization of the resulting macrocycles for improved efficacy and drug-like properties. Additionally, it illustrates the strategic value of macrocyclization in contemporary drug discovery efforts.

A novel mutation in hERG gene associated with azithromycin-induced acquired long QT syndrome.

BACKGROUND: Mutations in human ether-à-go-go-related gene (hERG) potassium channels are closely associated with long QT syndrome (LQTS). Previous studies have demonstrated that macrolide antibiotics increase the risk of cardiovascular diseases. To date, the mechanisms underlying acquired LQTS remain elusive. METHODS: A novel hERG mutation I1025N was identified in an azithromycin-treated patient with acquired long QT syndrome via Sanger sequencing. The mutant I1025N plasmid was transfected into HEK-293 cells, which were subsequently incubated with azithromycin. The effect of azithromycin and mutant I1025N on the hERG channel was evaluated via western blot, immunofluorescence, and electrophysiology techniques. RESULTS: The protein expression of the mature hERG protein was down-regulated, whereas that of the immature hERG protein was up-regulated in mutant I1025N HEK-293 cells. Azithromycin administration resulted in a negative effect on the maturation of the hERG protein. Additionally, the I1025N mutation exerted an inhibitory effect on hERG channel current. Moreover, azithromycin inhibited hERG channel current in a concentration-dependent manner. The I1025N mutation and azithromycin synergistically decreased hERG channel expression and hERG current. However, the I1025N mutation and azithromycin did not alter channel gating dynamics. CONCLUSIONS: These findings suggest that hERG gene mutations might be involved in the genetic susceptibility mechanism underlying acquired LQTS induced by azithromycin.

The scorpion toxin BeKm-1 blocks hERG cardiac potassium channels using an indispensable arginine residue.

BeKm-1 is a peptide toxin from scorpion venom that blocks the pore of the potassium channel hERG (Kv11.1) in the human heart. Although individual protein structures have been resolved, the structure of the complex between hERG and BeKm-1 is unknown. Here, we used molecular dynamics and ensemble docking, guided by previous double-mutant cycle analysis data, to obtain an in silico model of the hERG-BeKm-1 complex. Adding to the previous mutagenesis study of BeKm-1, our model uncovers the key role of residue Arg20, which forms three interactions (a salt bridge and hydrogen bonds) with the channel vestibule simultaneously. Replacement of this residue even by lysine weakens the interactions significantly. In accordance, the recombinantly produced BeKm-1R20K mutant exhibited dramatically decreased activity on hERG. Our model may be useful for future drug design attempts.

Population-specific variations in KCNH2 predispose patients to delayed ventricular repolarization upon dihydroartemisinin-piperaquine therapy.

Dihydroartemisinin-piperaquine is efficacious for the treatment of uncomplicated malaria and its use is increasing globally. Despite the positive results in fighting malaria, inhibition of the Kv11.1 channel (hERG; encoded by the KCNH2 gene) by piperaquine has raised concerns about cardiac safety. Whether genetic factors could modulate the risk of piperaquine-mediated QT prolongations remained unclear. Here, we first profiled the genetic landscape of KCNH2 variability using data from 141,614 individuals. Overall, we found 1,007 exonic variants distributed over the entire gene body, 555 of which were missense. By optimizing the gene-specific parametrization of 16 partly orthogonal computational algorithms, we developed a KCNH2-specific ensemble classifier that identified a total of 116 putatively deleterious missense variations. To evaluate the clinical relevance of KCNH2 variability, we then sequenced 293 Malian patients with uncomplicated malaria and identified 13 variations within the voltage sensing and pore domains of Kv11.1 that directly interact with channel blockers. Cross-referencing of genetic and electrocardiographic data before and after piperaquine exposure revealed that carriers of two common variants, rs1805121 and rs41314375, experienced significantly higher QT prolongations (ΔQTc of 41.8 ms and 61 ms, respectively, vs 14.4 ms in controls) with more than 50% of carriers having increases in QTc >30 ms. Furthermore, we identified three carriers of rare population-specific variations who experienced clinically relevant delayed ventricular repolarization. Combined, our results map population-scale genetic variability of KCNH2 and identify genetic biomarkers for piperaquine-induced QT prolongation that could help to flag at-risk patients and optimize efficacy and adherence to antimalarial therapy.

Design, synthesis, and biological evaluation of aralkyl piperazine and piperidine derivatives targeting SSRI/5-HT1A/5-HT7.

Serotonin reuptake inhibition combined with the action targeting 5-hydroxytryptamine receptor subtypes can serve as a potential target for the development of antidepressant drugs. Herein a series of new aralkyl piperazines and piperidines were designed and synthesized by the structural modifications of the previously discovered aralkyl piperidine compound 1, targeting SSRI/5-HT1A/5-HT7. The results exhibited that compound 5a showed strong binding to 5-HT1A and 5-HT7 (Ki of 0.46 nM, 2.7 nM, respectively) and a high level of serotonin reuptake inhibition (IC50 of 1.9 nM), all of which were significantly elevated compared to 1. In particular, compound 5a showed weaker inhibitory activity against hERG than 1, and demonstrated good stability in liver microsomes in vitro. The preliminary screening using FST indicated that orally administered 5a, at a high dose, could reduce immobility time in mice markedly, indicating potential antidepressant activity.

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.

An LQT2-related mutation in the voltage-sensing domain is involved in switching the gating polarity of hERG.

BACKGROUND: Cyclic Nucleotide-Binding Domain (CNBD)-family channels display distinct voltage-sensing properties despite sharing sequence and structural similarity. For example, the human Ether-a-go-go Related Gene (hERG) channel and the Hyperpolarization-activated Cyclic Nucleotide-gated (HCN) channel share high amino acid sequence similarity and identical domain structures. hERG conducts outward current and is activated by positive membrane potentials (depolarization), whereas HCN conducts inward current and is activated by negative membrane potentials (hyperpolarization). The structural basis for the "opposite" voltage-sensing properties of hERG and HCN remains unknown. RESULTS: We found the voltage-sensing domain (VSD) involves in modulating the gating polarity of hERG. We identified that a long-QT syndrome type 2-related mutation within the VSD, K525N, mediated an inwardly rectifying non-deactivating current, perturbing the channel closure, but sparing the open state and inactivated state. K525N rescued the current of a non-functional mutation in the pore helix region (F627Y) of hERG. K525N&F627Y switched hERG into a hyperpolarization-activated channel. The reactivated inward current induced by hyperpolarization mediated by K525N&F627Y can be inhibited by E-4031 and dofetilide quite well. Moreover, we report an extracellular interaction between the S1 helix and the S5-P region is crucial for modulating the gating polarity. The alanine substitution of several residues in this region (F431A, C566A, I607A, and Y611A) impaired the inward current of K525N&F627Y. CONCLUSIONS: Our data provide evidence that a potential cooperation mechanism in the extracellular vestibule of the VSD and the PD would determine the gating polarity in hERG.

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.

Recording ten-fold larger IKr conductances with automated patch clamping using equimolar Cs+ solutions.

Background: The rapid delayed rectifier potassium current (IKr) is important for cardiac repolarization and is most often involved in drug-induced arrhythmias. However, accurately measuring this current can be challenging in human-induced pluripotent stem cell (hiPSC)-derived cardiomyocytes because of its small current density. Interestingly, the ion channel conducting IKr, hERG channel, is not only permeable to K+ ions but also to Cs+ ions when present in equimolar concentrations inside and outside of the cell. Methods: In this study, IhERG was measured from Chinese hamster ovary (CHO)-hERG cells and hiPSC-CM using either Cs+ or K+ as the charge carrier. Equimolar Cs+ has been used in the literature in manual patch-clamp experiments, and here, we apply this approach using automated patch-clamp systems. Four different (pre)clinical drugs were tested to compare their effects on Cs+- and K+-based currents. Results: Using equimolar Cs+ solutions gave rise to approximately ten-fold larger hERG conductances. Comparison of Cs+- and K+-mediated currents upon application of dofetilide, desipramine, moxifloxacin, or LUF7244 revealed many similarities in inhibition or activation properties of the drugs studied. Using equimolar Cs+ solutions gave rise to approximately ten-fold larger hERG conductances. In hiPSC-CM, the Cs+-based conductance is larger compared to the known K+-based conductance, and the Cs+ hERG conductance can be inhibited similarly to the K+-based conductance. Conclusion: Using equimolar Cs+ instead of K+ for IhERG measurements in an automated patch-clamp system gives rise to a new method by which, for example, quick scans can be performed on effects of drugs on hERG currents. This application is specifically relevant when such experiments are performed using cells which express small IKr current densities in combination with small membrane capacitances.

Whole-Cell Configuration of the Patch-Clamp Technique in the hERG Channel Assay.

In vitro electrophysiological safety studies have become an integral part of the drug development process because, in many instances, compound-induced QT prolongation has been associated with a direct block of human ether-a-go-go-related gene (hERG) potassium channels or their native current, the rapidly activating delayed rectifier potassium current (IKr ). Therefore, according to the ICH S7B guideline, the in vitro hERG channel patch-clamp assay is commonly used as an early screen to predict the ability of a compound to prolong the QT interval prior to first-in-human testing. The protocols described in this article are designed to assess the effects of acute or long-term exposure to new chemical entities on the amplitude of IKr in HEK293 cells stably transfected with the hERG channel (whole-cell configuration of the patch-clamp technique). Examples of results obtained with moxifloxacin, terfenadine, arsenic, pentamidine, erythromycin, and sotalol are provided for illustrative purposes. © 2024 Wiley Periodicals LLC. Basic Protocol: Measurement of the acute effects of test items in the hERG channel test Alternate Protocol: Measurement of the long-term effects of test items in the hERG channel test.

Toward Digital Twin Technology for Precision Pharmacology.

The authors demonstrate the feasibility of technological innovation for personalized medicine in the context of drug-induced arrhythmia. The authors use atomistic-scale structural models to predict rates of drug interaction with ion channels and make predictions of their effects in digital twins of induced pluripotent stem cell-derived cardiac myocytes. The authors construct a simplified multilayer, 1-dimensional ring model with sufficient path length to enable the prediction of arrhythmogenic dispersion of repolarization. Finally, the authors validate the computational pipeline prediction of drug effects with data and quantify drug-induced propensity to repolarization abnormalities in cardiac tissue. The technology is high throughput, computationally efficient, and low cost toward personalized pharmacologic prediction.