Dynamics of hERG Closure Allow Novel Insights into hERG Blocking by Small Molecules

Machine-learning technique, QSAR and molecular dynamics for hERG-drug interactions

One of the most well-known anti-targets defining medication cardiotoxicity is the voltage-dependent hERG K + channel, which is well-known for its crucial involvement in cardiac action potential repolarization. Torsades de Pointes, QT prolongation, and sudden death are all caused by hERG (the human Ether-à-go-go-Related Gene) inhibition. There is great interest in creating predictive computational (in silico) tools to identify and weed out potential hERG blockers early in the drug discovery process because testing for hERG liability and the traditional experimental screening are complicated, expensive and time-consuming. This study used 2D descriptors of a large curated dataset of 6766 compounds and machine learning approaches to build robust descriptor-based QSAR and predictive classification models for KCNH2 liability. Decision Tree, Random Forest, Logistic Regression, Ada Boosting, kNN, SVM, Naïve Bayes, neural network and stochastic gradient classification classifier algorithms were used to build classification models. If a compound's IC50 value was between 10 μM and less, it was classified as a blocker (hERG-positive), and if it was more, it was classified as a non-blocker (hERG-negative). Matthew's correlation coefficient formula and F1score were applied to compare and track the developed models' performance. Molecular docking and dynamics studies were performed to understand the cardiotoxicity relating to the hERG-gene. The hERG residues interacting after 100 ns are LEU:697, THR:708, PHE:656, HIS:674, HIS:703, TRP:705 and ASN:709 and the hERG-ligand-16 complex trajectory showed stable behaviour with lesser fluctuations in the entire simulation of 200 ns.Communicated by Ramaswamy H. Sarma.

Structural modeling of the hERG potassium channel and associated drug interactions

The voltage-gated potassium channel, KV11.1, encoded by the human Ether-à-go-go-Related Gene (hERG), is expressed in cardiac myocytes, where it is crucial for the membrane repolarization of the action potential. Gating of the hERG channel is characterized by rapid, voltage-dependent, C-type inactivation, which blocks ion conduction and is suggested to involve constriction of the selectivity filter. Mutations S620T and S641A/T within the selectivity filter region of hERG have been shown to alter the voltage dependence of channel inactivation. Because hERG channel blockade is implicated in drug-induced arrhythmias associated with both the open and inactivated states, we used Rosetta to simulate the effects of hERG S620T and S641A/T mutations to elucidate conformational changes associated with hERG channel inactivation and differences in drug binding between the two states. Rosetta modeling of the S641A fast-inactivating mutation revealed a lateral shift of the F627 side chain in the selectivity filter into the central channel axis along the ion conduction pathway and the formation of four lateral fenestrations in the pore. Rosetta modeling of the non-inactivating mutations S620T and S641T suggested a potential molecular mechanism preventing F627 side chain from shifting into the ion conduction pathway during the proposed inactivation process. Furthermore, we used Rosetta docking to explore the binding mechanism of highly selective and potent hERG blockers - dofetilide, terfenadine, and E4031. Our structural modeling correlates well with much, but not all, existing experimental evidence involving interactions of hERG blockers with key residues in hERG pore and reveals potential molecular mechanisms of ligand interactions with hERG in an inactivated state.

Molecular dynamics of hERG channel: insights into understanding the binding of small molecules for detuning cardiotoxicity

Evaluation of cardiotoxicity potential of new chemical entities (NCEs) has lately become one of the stringent filters in the drug discovery and development process. Cardiotoxicity is caused mainly by the inhibition of human ether-a-go-go related gene (hERG) channel protein. Inhibition of the hERG channel leads to a life-threatening condition known as cardiac arrhythmia. Knowledge of the structural behaviour of the hERG would aid greatly in the design of new drug molecules that do not interact with the protein and add to the safety index. In this study, a computational model for the active-state of hERG was developed. This model was equilibrated by performing the molecular dynamics simulations for 100 ns followed by clustering and selection of a representative structure based on the largest populated cluster. To study the changes in the protein structure on inhibition, three inhibitory ligands, namely, dofetilide, cisapride and terfenadine were docked, followed by molecular dynamics simulations of 200 ns for the apo and each ligand-bound structure. It was observed that docking and simulation studies of the hERG model exhibited noticeable conformational changes in the protein upon ligand-binding. A significant change in the kink of the S6-transmembrane helix was observed. Inter-chain distances between the crucial residues Y652 and F656 (present below the ion-selectivity filter), their side-chain orientation and hydrogen bonding indicated a probable collapse of the pore. These changes may infer the initiation in transition of hERG from an open to an inactive state. Hence, these findings would help in designing compounds devoid of hERG inhibition with reduced cardiotoxicity.Communicated by Ramaswamy H. Sarma.

hERG toxicity assessment: Useful guidelines for drug design

All along the drug development process, one of the most frequent adverse side effects, leading to the failure of drugs, is the cardiac arrhythmias. Such failure is mostly related to the capacity of the drug to inhibit the human ether-à-go-go-related gene (hERG) cardiac potassium channel. The early identification of hERG inhibition properties of biological active compounds has focused most of attention over the years. In order to prevent the cardiac side effects, a great number of in silico, in vitro and in vivo assays have been performed. The main goal of these studies is to understand the reasons of these effects, and then to give information or instructions to scientists involved in drug development to avoid the cardiac side effects. To evaluate anticipated cardiovascular effects, early evaluation of hERG toxicity has been strongly recommended for instance by the regulatory agencies such as U.S. Food and Drug Administration (FDA) and European Medicines Agency (EMA). Thus, following an initial screening of a collection of compounds to find hits, a great number of pharmacomodulation studies on the novel identified chemical series need to be performed including activity evaluation towards hERG. We provide in this concise review clear guidelines, based on described examples, illustrating successful optimization process to avoid hERG interactions as cases studies and to spur scientists to develop safe drugs.

Molecular Docking Guided Grid-Independent Descriptor Analysis to Probe the Impact of Water Molecules on Conformational Changes of hERG Inhibitors in Drug Trapping Phenomenon

Human ether a-go-go related gene (hERG) or KV11.1 potassium channels mediate the rapid delayed rectifier current (I_Kr) in cardiac myocytes. Drug-induced inhibition of hERG channels has been implicated in the development of acquired long QT syndrome type (aLQTS) and fatal arrhythmias. Several marketed drugs have been withdrawn for this reason. Therefore, there is considerable interest in developing better tests for predicting drugs which can block the hERG channel. The drug-binding pocket in hERG channels, which lies below the selectivity filter, normally contains K⁺ ions and water molecules. In this study, we test the hypothesis that these water molecules impact drug binding to hERG. We developed 3D QSAR models based on alignment independent descriptors (GRIND) using docked ligands in open and closed conformations of hERG in the presence (solvated) and absence (non-solvated) of water molecules. The ligand–protein interaction fingerprints (PLIF) scheme was used to summarize and compare the interactions. All models delineated similar 3D hERG binding features, however, small deviations of about ~0.4 Å were observed between important hotspots of molecular interaction fields (MIFs) between solvated and non-solvated hERG models. These small changes in conformations do not affect the performance and predictive power of the model to any significant extent. The model that exhibits the best statistical values was attained with a cryo_EM structure of the hERG channel in open state without water. This model also showed the best R² of 0.58 and 0.51 for the internal and external validation test sets respectively. Our results suggest that the inclusion of water molecules during the docking process has little effect on conformations and this conformational change does not impact the predictive ability of the 3D QSAR models.

Drug-likeness of linear pentamidine analogues and their impact on the hERG K+channel – correlation with structural features

The pentamidines with S atoms or sulfanilide groups in the linker have favorable drug-likeness parameters and low toxicity.

Open-Access Activity Prediction Tools for Natural Products. Case Study: hERG Blockers

Interference with the hERG potassium ion channel may cause cardiac arrhythmia and can even lead to death. Over the last few decades, several drugs, already on the market, and many more investigational drugs in various development stages, have had to be discontinued because of their hERG-associated toxicity. To recognize potential hERG activity in the early stages of drug development, a wide array of computational tools, based on different principles, such as 3D QSAR, 2D and 3D similarity, and machine learning, have been developed and are reviewed in this chapter. The various available prediction tools Similarity Ensemble Approach, SuperPred, SwissTargetPrediction, HitPick, admetSAR, PASSonline, Pred-hERG, and VirtualToxLab™ were used to screen a dataset of known hERG synthetic and natural product actives and inactives to quantify and compare their predictive power. This contribution will allow the reader to evaluate the suitability of these computational methods for their own related projects. There is an unmet need for natural product-specific prediction tools in this field.

The current limits in virtual screening and property prediction

Beyond finding inhibitors that show high binding affinity to the respective target, there is the challenge of optimizing their properties with respect to metabolic and toxicological issues, as well as further off-target effects. To reduce the experimental effort of synthesizing and testing actual substances in corresponding assays, virtual screening has become an indispensable toolbox in preclinical development. The scope of application covers the prediction of molecular properties including solubility, metabolic liability and binding to antitargets, such as the hERG channel. Furthermore, prediction of binding sites and drugable targets are emerging aspects of virtual screening. Issues involved with the currently applied computational models including machine learning algorithms are outlined, such as limitations to the accuracy of prediction and overfitting.

New potential binding determinant for hERG channel inhibitors

Human ether-à-go-go related gene (hERG) 1 channels conduct the rapid delayed rectifier K⁺ current (I_Kr) and are essential for the repolarization of the cardiac action potential. hERG1 inhibition by structurally diverse drugs may lead to life threatening arrhythmia. Putative binding determinants of hERG1 channel blockers include T623, S624 and V625 on the pore helix, and residues G648, Y652 and F656, located on segment S6. We and others have previously hypothesized that additional binding determinants may be located on helix S5, which is in close contact with the S6 segments. In order to test this hypothesis, we performed a detailed investigation combining ionic current measurements with two-microelectrode voltage clamp and molecular modeling techniques. We identified a novel aromatic high affinity binding determinant for blockers located in helix S5, F557, which is equally potent as Y652. Modeling supports a direct interaction with the outer pore helix.

Astemizole Derivatives as Fluorescent Probes for hERG Potassium Channel Imaging

The detection and imaging of hERG potassium channels in living cells can provide useful information for hERG-correlation studies. Herein, three small-molecule fluorescent probes, based on the potent hERG channel inhibitor astemizole, for the imaging of hERG channels in hERG-transfected HEK293 cells (hERG-HEK293) and human colorectal cancer cells (HT-29), are described. These probes are expected to be applied in the physiological and pathological studies of hERG channels.

Drug trapping in hERG K+ channels: (not) a matter of drug size?

Inhibition of hERG K+ channels by structurally diverse drugs prolongs the ventricular action potential and increases the risk of torsade de pointes arrhythmias and sudden cardiac death. The capture of drugs behind closed channel gates, so-called drug trapping, is suggested to harbor an increased pro-arrhythmic risk. In this study, the trapping mechanisms of a trapped hERG blocker propafenone and a bulky derivative (MW: 647.24 g mol-1) were studied by making use of electrophysiological measurements in combination with molecular dynamics simulations. Our study suggests that the hERG cavity is able to accommodate very bulky compounds without disturbing gate closure.

Computational investigations of hERG channel blockers: New insights and current predictive models

Identification of potential human Ether-a-go-go Related-Gene (hERG) potassium channel blockers is an essential part of the drug development and drug safety process in pharmaceutical industries or academic drug discovery centers, as they may lead to drug-induced QT prolongation, arrhythmia and Torsade de Pointes. Recent reports also suggest starting to address such issues at the hit selection stage. In order to prioritize molecules during the early drug discovery phase and to reduce the risk of drug attrition due to cardiotoxicity during pre-clinical and clinical stages, computational approaches have been developed to predict the potential hERG blockage of new drug candidates. In this review, we will describe the current in silico methods developed and applied to predict and to understand the mechanism of actions of hERG blockers, including ligand-based and structure-based approaches. We then discuss ongoing research on other ion channels and hERG polymorphism susceptible to be involved in LQTS and how systemic approaches can help in the drug safety decision.

In silico prediction of hERG inhibition

The voltage-gated potassium channel encoded by hERG carries a delayed rectifying potassium current (IKr) underlying repolarization of the cardiac action potential. Pharmacological blockade of the hERG channel results in slowed repolarization and therefore prolongation of action potential duration and an increase in the QT interval as measured on an electrocardiogram. Those are possible to cause sudden death, leading to the withdrawals of many drugs, which is the reason for hERG screening. Computational in silico prediction models provide a rapid, economic way to screen compounds during early drug discovery. In this review, hERG prediction models are classified as 2D and 3D quantitative structure–activity relationship models, pharmacophore models, classification models, and structure based models (using homology models of hERG).