Towards a Structural View of Drug Binding to hERG K+ Channels

CardioDPi: An explainable deep-learning model for identifying cardiotoxic chemicals targeting hERG, Cav1.2, and Nav1.5 channels

The cardiotoxic effects of various pollutants have been a growing concern in environmental and material science. These effects encompass arrhythmias, myocardial injury, cardiac insufficiency, and pericardial inflammation. Compounds such as organic solvents and air pollutants disrupt the potassium, sodium, and calcium ion channels cardiac cell membranes, leading to the dysregulation of cardiac function. However, current cardiotoxicity models have disadvantages of incomplete data, ion channels, interpretability issues, and inability of toxic structure visualization. Herein, an interpretable deep-learning model known as CardioDPi was developed, which is capable of discriminating cardiotoxicity induced by the human Ether-à-go-go-related gene (hERG) channel, sodium channel (Na_v1.5), and calcium channel (Ca_v1.5) blockade. External validation yielded promising area under the ROC curve (AUC) values of 0.89, 0.89, and 0.94 for the hERG, Na_v1.5, and Ca_v1.5 channels, respectively. The CardioDPi can be freely accessed on the web server CardioDPipredictor (http://cardiodpi.sapredictor.cn/). Furthermore, the structural characteristics of cardiotoxic compounds were analyzed and structural alerts (SAs) can be extracted using the user-friendly CardioDPi-SAdetector web service (http://cardiosa.sapredictor.cn/). CardioDPi is a valuable tool for identifying cardiotoxic chemicals that are environmental and health risks. Moreover, the SA system provides essential insights for mode-of-action studies concerning cardiotoxic compounds.

P-Glycoprotein-Mediated Interaction Is a Risk Factor for QT Prolongation in Concomitant Use of Antipsychotics and SSRIs as P-Glycoprotein-Mediated Inhibitors: Analysis of the Japanese Adverse Drug Event Report Database

The inhibition of human ether-a-go-go-related gene (hERG) channels is a known cause of QT prolongation triggered by antipsychotic drugs. Our previous studies suggest that P-glycoprotein (P-gp)-mediated drug interactions may lead to increased gastrointestinal absorption of pimozide and its accumulation in cardiomyocytes, thereby enhancing the inhibitory effect of hERG channels. There is a paucity of epidemiological studies examining the risk of QT prolongation by antipsychotic drugs in terms of P-gp-mediated interactions with concomitant drugs. Therefore, using the Japanese Adverse Event Reporting Database, we investigated whether the risk of QT prolongation triggered by antipsychotic drugs associated with hERG inhibition is affected by the concomitant use of selective serotonin reuptake inhibitors (SSRIs) associated with P-gp inhibition. The results showed that the frequency of QT prolongation increased when the antipsychotic drugs quetiapine and sulpiride, which are P-gp substrates, were combined with SSRIs with P-gp inhibition. In contrast, no association with QT prolongation was observed in patients on non-P-gp-substrate antipsychotics, irrespective of the P-gp inhibitory effect of the concomitant SSRI. These results suggest that P-gp-mediated interactions are a risk factor for antipsychotic-induced QT prolongation. There is a need for further investigation into the risks of specific drug combinations.

Structural modeling of hERG channel-drug interactions using Rosetta

The human ether-a-go-go-related gene (hERG) not only encodes a potassium-selective voltage-gated ion channel essential for normal electrical activity in the heart but is also a major drug anti-target. Genetic hERG mutations and blockage of the channel pore by drugs can cause long QT syndrome, which predisposes individuals to potentially deadly arrhythmias. However, not all hERG-blocking drugs are proarrhythmic, and their differential affinities to discrete channel conformational states have been suggested to contribute to arrhythmogenicity. We used Rosetta electron density refinement and homology modeling to build structural models of open-state hERG channel wild-type and mutant variants (Y652A, F656A, and Y652A/F656 A) and a closed-state wild-type channel based on cryo-electron microscopy structures of hERG and EAG1 channels. These models were used as protein targets for molecular docking of charged and neutral forms of amiodarone, nifekalant, dofetilide, d/l-sotalol, flecainide, and moxifloxacin. We selected these drugs based on their different arrhythmogenic potentials and abilities to facilitate hERG current. Our docking studies and clustering provided atomistic structural insights into state-dependent drug-channel interactions that play a key role in differentiating safe and harmful hERG blockers and can explain hERG channel facilitation through drug interactions with its open-state hydrophobic pockets.

Facilitation of hERG Activation by Its Blocker: A Mechanism to Reduce Drug-Induced Proarrhythmic Risk

Modulation of the human Ether-à-go-go-Related Gene (hERG) channel, a crucial voltage-gated potassium channel in the repolarization of action potentials in ventricular myocytes of the heart, has significant implications on cardiac electrophysiology and can be either antiarrhythmic or proarrhythmic. For example, hERG channel blockade is a leading cause of long QT syndrome and potentially life-threatening arrhythmias, such as torsades de pointes. Conversely, hERG channel blockade is the mechanism of action of Class III antiarrhythmic agents in terminating ventricular tachycardia and fibrillation. In recent years, it has been recognized that less proarrhythmic hERG blockers with clinical potential or Class III antiarrhythmic agents exhibit, in addition to their hERG-blocking activity, a second action that facilitates the voltage-dependent activation of the hERG channel. This facilitation is believed to reduce the proarrhythmic potential by supporting the final repolarizing of action potentials. This review covers the pharmacological characteristics of hERG blockers/facilitators, the molecular mechanisms underlying facilitation, and their clinical significance, as well as unresolved issues and requirements for research in the fields of ion channel pharmacology and drug-induced arrhythmias.

Utilizing Heteroatom Types and Numbers from Extensive Ligand Libraries to Develop Novel hERG Blocker QSAR Models Using Machine Learning-Based Classifiers

The human ether-à-go-go-related gene (hERG) channel plays a crucial role in membrane repolarization. Any disruptions in its function can lead to severe cardiovascular disorders such as long QT syndrome (LQTS), which increases the risk of serious cardiovascular problems such as tachyarrhythmia and sudden cardiac death. Drug-induced LQTS is a significant concern and has resulted in drug withdrawals from the market in the past. The main objective of this study is to pinpoint crucial heteroatoms present in ligands that initiate interactions leading to the effective blocking of the hERG channel. To achieve this aim, ligand-based quantitative structure-activity relationships (QSAR) models were constructed using extensive ligand libraries, considering the heteroatom types and numbers, and their associated hERG channel blockage pIC50 values. Machine learning-assisted QSAR models were developed to analyze the key structural components influencing compound activity. Among the various methods, the KPLS method proved to be the most efficient, allowing the construction of models based on eight distinct fingerprints. The study delved into investigating the influence of heteroatoms on the activity of hERG blockers, revealing their significant role. Furthermore, by quantifying the effect of heteroatom types and numbers on ligand activity at the hERG channel, six compound pairs were selected for molecular docking. Subsequent molecular dynamics simulations and per residue MM/GBSA calculations were performed to comprehensively analyze the interactions of the selected pair compounds.

In silico prediction of hERG blockers using machine learning and deep learning approaches

The human ether-à-go-go-related gene (hERG) is associated with drug cardiotoxicity. If the hERG channel is blocked, it will lead to prolonged QT interval and cause sudden death in severe cases. Therefore, it is important to evaluate the hERG-blocking property of compounds in early drug discovery. In this study, a dataset containing 4556 compounds with IC50 values determined by patch clamp techniques on mammalian lineage cells was collected, and hERG blockers and non-blockers were distinguished according to three single thresholds and two binary thresholds. Four machine learning (ML) algorithms combining four molecular fingerprints and molecular descriptors as well as graph convolutional neural networks (GCNs) were used to construct a series of binary classification models. The results showed that the best models varied for different thresholds. The ML models implemented by support vector machine and random forest performed well based on Morgan fingerprints and molecular descriptors, with AUCs ranging from 0.884 to 0.950. GCN showed superior prediction performance with AUCs above 0.952, which might be related to its direct extraction of molecular features from the original input. Meanwhile, the classification of binary threshold was better than that of single threshold, which could provide us with a more accurate prediction of hERG blockers. At last, the applicability domain for the model was defined, and seven structural alerts that might generate hERG blockage were identified by information gain and substructure frequency analysis. Our work would be beneficial for identifying hERG blockers in chemicals.

Synthetic Routes to Approved Drugs Containing a Spirocycle

The use of spirocycles in drug discovery and medicinal chemistry has been booming in the last two decades. This has clearly translated into the landscape of approved drugs. Among two dozen clinically used medicines containing a spirocycle, 50% have been approved in the 21st century. The present review focuses on the notable synthetic routes to such drugs invented in industry and academia, and is intended to serve as a useful reference source of synthetic as well as general drug information for researchers engaging in the design of new spirocyclic scaffolds for medicinal use or embarking upon analog syntheses inspired by the existing approved drugs.

Virtual screening of DrugBank database for hERG blockers using topological Laplacian-assisted AI models

The human ether-a-go-go (hERG) potassium channel (Kv11.1) plays a critical role in mediating cardiac action potential. The blockade of this ion channel can potentially lead fatal disorder and/or long QT syndrome. Many drugs have been withdrawn because of their serious hERG-cardiotoxicity. It is crucial to assess the hERG blockade activity in the early stage of drug discovery. We are particularly interested in the hERG-cardiotoxicity of compounds collected in the DrugBank database considering that many DrugBank compounds have been approved for therapeutic treatments or have high potential to become drugs. Machine learning-based in silico tools offer a rapid and economical platform to virtually screen DrugBank compounds. We design accurate and robust classifiers for blockers/non-blockers and then build regressors to quantitatively analyze the binding potency of the DrugBank compounds on the hERG channel. Molecular sequences are embedded with two natural language processing (NLP) methods, namely, autoencoder and transformer. Complementary three-dimensional (3D) molecular structures are embedded with two advanced mathematical approaches, i.e., topological Laplacians and algebraic graphs. With our state-of-the-art tools, we reveal that 227 out of the 8641 DrugBank compounds are potential hERG blockers, suggesting serious drug safety problems. Our predictions provide guidance for the further experimental interrogation of DrugBank compounds' hERG-cardiotoxicity.

Intracellular Binding of Terfenadine Competes with Its Access to Pancreatic ß-cell ATP-Sensitive K+ Channels and Human ether-à-go-go-Related Gene Channels

Most blockers of both hERG (human ether-à-go-go-related gene) channels and pancreatic ß-cell ATP-sensitive K⁺ (K_ATP) channels access their binding sites from the cytoplasmic side of the plasma membrane. It is unknown whether binding to intracellular components competes with binding of these substances to K⁺ channels. The whole-cell configuration of the patch-clamp technique, a laser-scanning confocal microscope, and fluorescence correlation spectroscopy (FCS) were used to study hERG channels expressed in HEK (human embryonic kidney) 293 cells and K_ATP channels from the clonal insulinoma cell line RINm5F. When applied via the pipette solution in the whole-cell configuration, terfenadine blocked both hERG and K_ATP currents with much lower potency than after application via the bath solution, which was not due to P-glycoprotein-mediated efflux of terfenadine. Such a difference was not observed with dofetilide and tolbutamide. 37-68% of hERG/EGFP (enhanced green-fluorescent protein) fusion proteins expressed in HEK 293 cells were slowly diffusible as determined by laser-scanning microscopy in the whole-cell configuration and by FCS in intact cells. Bath application of a green-fluorescent sulphonylurea derivative (Bodipy-glibenclamide) induced a diffuse fluorescence in the cytosol of RINm5F cells under whole-cell patch-clamp conditions. These observations demonstrate the presence of intracellular binding sites for hERG and K_ATP channel blockers not dialyzable by the patch-pipette solution. Intracellular binding of terfenadine was not influenced by a mutated hERG (Y652A) channel. In conclusion, substances with high lipophilicity are not freely diffusible inside the cell but steep concentration gradients might exist within the cell and in the sub-membrane space.

General Pharmacological Activation Mechanism of K+ Channels Bypassing Channel Gates

Under the known pharmacological activation mechanisms, activators allosterically or directly open potassium channel gates. However, herein, molecular dynamics simulations on TREK-1, a member of the channel class gated at the filter, suggested that negatively charged activators act with a gate-independent mechanism where compounds increase currents by promoting ions passing through the central cavity. Then, based on studies of KCNQ2, we uncovered that this noncanonical activation mechanism is shared by the other channel class gated at the helix-bundle crossing. Rational drug design found a novel KCNQ2 agonist, CLE030, which stably binds to the central cavity. Functional analysis, molecular dynamics simulations, and calculations of the potential of mean force revealed that the carbonyl oxygen of CLE030 influences permeant ions in the central cavity to contribute to its activation effects. Together, this study discovered a ligand-to-ion activation mechanism for channels that bypasses their gates and thus is conserved across subfamilies with different gates.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Bifurcations and Proarrhythmic Behaviors in Cardiac Electrical Excitations

The heart is a hierarchical dynamic system consisting of molecules, cells, and tissues, and acts as a pump for blood circulation. The pumping function depends critically on the preceding electrical activity, and disturbances in the pattern of excitation propagation lead to cardiac arrhythmia and pump failure. Excitation phenomena in cardiomyocytes have been modeled as a nonlinear dynamical system. Because of the nonlinearity of excitation phenomena, the system dynamics could be complex, and various analyses have been performed to understand the complex dynamics. Understanding the mechanisms underlying proarrhythmic responses in the heart is crucial for developing new ways to prevent and control cardiac arrhythmias and resulting contractile dysfunction. When the heart changes to a pathological state over time, the action potential (AP) in cardiomyocytes may also change to a different state in shape and duration, often undergoing a qualitative change in behavior. Such a dynamic change is called bifurcation. In this review, we first summarize the contribution of ion channels and transporters to AP formation and our knowledge of ion-transport molecules, then briefly describe bifurcation theory for nonlinear dynamical systems, and finally detail its recent progress, focusing on the research that attempts to understand the developing mechanisms of abnormal excitations in cardiomyocytes from the perspective of bifurcation phenomena.

Structure-Based Prediction of hERG-Related Cardiotoxicity: A Benchmark Study

Drug-induced blockade of the human ether-à-go-go-related gene (hERG) channel is today considered the main cause of cardiotoxicity in postmarketing surveillance. Hence, several ligand-based approaches were developed in the last years and are currently employed in the early stages of a drug discovery process for in silico cardiac safety assessment of drug candidates. Herein, we present the first structure-based classifiers able to discern hERG binders from nonbinders. LASSO regularized support vector machines were applied to integrate docking scores and protein–ligand interaction fingerprints. A total of 396 models were trained and validated based on: (i) high-quality experimental bioactivity information returned by 8337 curated compounds extracted from ChEMBL (version 25) and (ii) structural predictor data. Molecular docking simulations were performed using GLIDE and GOLD software programs and four different hERG structural models, namely, the recently published structures obtained by cryoelectron microscopy (PDB codes: 5VA1 and 7CN1) and two published homology models selected for comparison. Interestingly, some classifiers return performances comparable to ligand-based models in terms of area under the ROC curve (AUC_MAX = 0.86 ± 0.01) and negative predictive values (NPV_MAX = 0.81 ± 0.01), thus putting forward the herein proposed computational workflow as a valuable tool for predicting hERG-related cardiotoxicity without the limitations of ligand-based models, typically affected by low interpretability and a limited applicability domain. From a methodological point of view, our study represents the first example of a successful integration of docking scores and protein–ligand interaction fingerprints (IFs) through a support vector machine (SVM) LASSO regularized strategy. Finally, the study highlights the importance of using hERG structural models accounting for ligand-induced fit effects and allowed us to select the best-performing protein conformation (made available in the Supporting Information, SI) to be employed for a reliable structure-based prediction of hERG-related cardiotoxicity.

Molecular dynamics: a powerful tool for studying the medicinal chemistry of ion channel modulators

Molecular dynamics (MD) simulations allow researchers to investigate the behavior of desired biological targets at ever-decreasing costs with ever-increasing precision. Among the biological macromolecules, ion channels are remarkable transmembrane proteins, capable of performing special biological processes and revealing a complex regulatory matrix, including modulation by small molecules, either endogenous or exogenous. Recently, given the developments in ion channel structure determination and accessibility of bio-computational techniques, MD and related tools are becoming increasingly popular in the intense research area regarding ligand-channel interactions. This review synthesizes and presents the most important fields of MD involvement in investigating channel-molecule interactions, including, but not limited to, deciphering the binding modes of ligands to their ion channel targets and the mechanisms through which chemical compounds exert their effect on channel function. Special attention is devoted to the importance of more elaborate methods, such as free energy calculations, while principles regarding drug design and discovery are highlighted. Several technical aspects involving the creation and simulation of channel-molecule MD systems (ligand parameterization, proper membrane setup, system building, etc.) are also presented.

Quinazolinone dimers as a potential new class of safer Kv1 inhibitors: Overcoming hERG, sodium and calcium channel affinities

The discovery of more selective and safer voltage-gated potassium channel blockers is an extremely demanding approach. Designing selective Kv1.5 inhibitors is very challenging as only limited data is available on this target due to a lacking crystal structure for this ion channel receptor. Herein, we synthesized a series of 21 novel quinazolinone dimers 3a-i, 5a-i and 10a-c. We tried to avoid structural features responsible for non-selectivity and for most potassium channel blockers' side effects in our design. In contrast to other works, which lack investigation over wide ranges of potassium and sodium channels, we screened the inhibitory activity of our synthesized compounds over multiple voltage-gated potassium channels, including six different human Kv1 channel subtypes Kv1.1, Kv1.2, Kv1.3, Kv1.4, Kv1.5 and Kv1.6 channels as well as Kv2.1, Kv3.1, Kv4.3, Kv7.2, Kv7.3, Kv10.1, hERG, and Shaker IR. Moreover, these compounds' selectivity was investigated on sodium channels Nav1.2, Nav1.4 and Nav1.5 and calcium channels Cav3.1-Cav3.3. The results revealed two compounds (3a and 3e) with low micromolar Kv1.5 inhibition activity with EC₅₀ values of 5.1 ± 0.9 µM and 12.5 ± 1.1 µM, respectively. However, at higher concentrations, they also showed inhibitory activity on Kv1.3 and Kv1.1 channels. This might be due to structural similarities between these three Kv1 channel isoforms. Compound 3a shows a slight preference for Kv1.5. Interestingly, they lack any activity on other potassium channels (including hERG), sodium channels, and calcium channels. Our findings recommend quinazolinone dimers with ethylene linker as a potential new class of safer Kv1 inhibitors and a good start for designing more selective and potent Kv1.5 inhibitors.

Discovery of KV 1.3 ion channel inhibitors: Medicinal chemistry approaches and challenges

The KV 1.3 voltage-gated potassium ion channel is involved in many physiological processes both at the plasma membrane and in the mitochondria, chiefly in the immune and nervous systems. Therapeutic targeting KV 1.3 with specific peptides and small molecule inhibitors shows great potential for treating cancers and autoimmune diseases, such as multiple sclerosis, type I diabetes mellitus, psoriasis, contact dermatitis, rheumatoid arthritis, and myasthenia gravis. However, no KV 1.3-targeted compounds have been approved for therapeutic use to date. This review focuses on the presentation of approaches for discovering new KV 1.3 peptide and small-molecule inhibitors, and strategies to improve the selectivity of active compounds toward KV 1.3. Selectivity of dalatazide (ShK-186), a synthetic derivate of the sea anemone toxin ShK, was achieved by chemical modification and has successfully reached clinical trials as a potential therapeutic for treating autoimmune diseases. Other peptides and small-molecule inhibitors are critically evaluated for their lead-like characteristics and potential for progression into clinical development. Some small-molecule inhibitors with well-defined structure-activity relationships have been optimized for selective delivery to mitochondria, and these offer therapeutic potential for the treatment of cancers. This overview of KV 1.3 inhibitors and methodologies is designed to provide a good starting point for drug discovery to identify novel effective KV 1.3 modulators against this target in the future.

Proton Pump Inhibitors Directly Block hERG-Potassium Channel and Independently Increase the Risk of QTc Prolongation in a Large Cohort of US Veterans

[Figure: see text].

Molecular Insights Into the Gating Kinetics of the Cardiac hERG Channel, Illuminated by Structure and Molecular Dynamics

The rapidly activating delayed rectifier K⁺ current generated by the cardiac hERG potassium channel encoded by KCNH2 is the most important reserve current for cardiac repolarization. The unique inward rectification characteristics of the hERG channel depend on the gating regulation, which involves crucial structural domains and key single amino acid residues in the full-length hERG channel. Identifying critical molecules involved in the regulation of gating kinetics for the hERG channel requires high-resolution structures and molecular dynamics simulation models. Based on the latest progress in hERG structure and molecular dynamics simulation research, summarizing the molecules involved in the changes in the channel state helps to elucidate the unique gating characteristics of the channel and the reason for its high affinity to cardiotoxic drugs. In this review, we aim to summarize the significant advances in understanding the voltage gating regulation of the hERG channel based on its structure obtained from cryo-electron microscopy and computer simulations, which reveal the critical roles of several specific structural domains and amino acid residues.

Anthracycline Associated Disturbances of Cardiovascular Homeostasis

Despite the dynamic progress of modern medicine, oncological and cardiovascular diseases (CVD) remain a severe economic burden worldwide. Therefore, the study of chemotherapeutic cardiotoxicity appears to be comprehensively demanded. Nowadays, pharmacological therapy in oncology has undoubtedly unprecedented development, but at the same time, the rates of cardiovascular complications of chemotherapy still remain unchanged. The well-established and highly effective, but at the same time, cardiotoxic anthracyclines have not lost their relevance. Furthermore, they remain indispensable components of an immense amount of chemotherapy regimens, such as AC, FAC, etc. Moreover, the anthracycline-containing chemotherapy regimens have become a standard of care in several cancer types. In the context of the above mentioned, the study of the pathophysiological mechanisms, biochemical aspects, and dynamics of the morphological remodeling of doxorubicin-induced cardiovascular homeostasis disturbances will enable finding new targets of pharmacological therapy, which either in the short or long perspectives, will have a beneficial effect, improving both the quality of life and prognosis of oncological patients. This article covers a versatile overview of the molecular mechanisms of doxorubicin-induced cardiotoxicity. The pathogenesis of cardiotoxicity assessment could help to explore specific molecular mechanisms that initiate cardiovascular alteration that may favorably affect the future development of targeted drugs that could prevent cardiovascular events in cancer patients.

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.

Flavonoids and hERG channels: Friends or foes?

The cardiac action potential is regulated by several ion channels. Drugs capable to block these channels, in particular the human ether-à-go-go-related gene (hERG) channel, also known as K_V11.1 channel, may lead to a potentially lethal ventricular tachyarrhythmia called "Torsades de Pointes". Thus, evaluation of the hERG channel off-target activity of novel chemical entities is nowadays required to safeguard patients as well as to avoid attrition in drug development. Flavonoids, a large class of natural compounds abundantly present in food, beverages, herbal medicines, and dietary food supplements, generally escape this assessment, though consumed in consistent amounts. Continuously growing evidence indicates that these compounds may interact with the hERG channel and block it. The present review, by examining numerous studies, summarizes the state-of-the-art in this field, describing the most significant examples of direct and indirect inhibition of the hERG channel current operated by flavonoids. A description of the molecular interactions between a few of these natural molecules and the Rattus norvegicus channel protein, achieved by an in silico approach, is also presented.

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.

Refinement of a cryo-EM structure of hERG: Bridging structure and function

The human-ether-a-go-go-related gene (hERG) encodes the voltage-gated potassium channel (KCNH2 or Kv11.1, commonly known as hERG). This channel plays a pivotal role in the stability of phase 3 repolarization of the cardiac action potential. Although a high-resolution cryo-EM structure is available for its depolarized (open) state, the structure surprisingly did not feature many functionally important interactions established by previous biochemical and electrophysiology experiments. Using molecular dynamics flexible fitting (MDFF), we refined the structure and recovered the missing functionally relevant salt bridges in hERG in its depolarized state. We also performed electrophysiology experiments to confirm the functional relevance of a novel salt bridge predicted by our refinement protocol. Our work shows how refinement of a high-resolution cryo-EM structure helps to bridge the existing gap between the structure and function in the voltage-sensing domain (VSD) of hERG.

Off-label use of chloroquine, hydroxychloroquine, azithromycin and lopinavir/ritonavir in COVID-19 risks prolonging the QT interval by targeting the hERG channel

Coronavirus disease-2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), poses an enormous challenge to the medical system, especially the lack of safe and effective COVID-19 treatment methods, forcing people to look for drugs that may have therapeutic effects as soon as possible. Some old drugs have shown clinical benefits after a few small clinical trials that attracted great attention. Clinically, however, many drugs, including those currently used in COVID-19, such as chloroquine, hydroxychloroquine, azithromycin, and lopinavir/ritonavir, may cause cardiotoxicity by acting on cardiac potassium channels, especially hERG channel through their off-target effects. The blocking of the hERG channel prolongs QT intervals on electrocardiograms; thus, it might induce severe ventricular arrhythmias and even sudden cardiac death. Therefore, while focusing on the efficacy of COVID-19 drugs, the fact that they block hERG channels to cause arrhythmias cannot be ignored. To develop safer and more effective drugs, it is necessary to understand the interactions between drugs and the hERG channel and the molecular mechanism behind this high affinity. In this review, we focus on the biochemical and molecular mechanistic aspects of drug-related blockade of the hERG channel to provide insights into QT prolongation caused by off-label use of related drugs in COVID-19, and hope to weigh the risks and benefits when using these drugs.

Critical Assessment of Artificial Intelligence Methods for Prediction of hERG Channel Inhibition in the "Big Data" Era

The rise of novel artificial intelligence (AI) methods necessitates their benchmarking against classical machine learning for a typical drug-discovery project. Inhibition of the potassium ion channel, whose alpha subunit is encoded by the human ether-à-go-go-related gene (hERG), leads to a prolonged QT interval of the cardiac action potential and is a significant safety pharmacology target for the development of new medicines. Several computational approaches have been employed to develop prediction models for the assessment of hERG liabilities of small molecules including recent work using deep learning methods. Here, we perform a comprehensive comparison of hERG effect prediction models based on classical approaches (random forests and gradient boosting) and modern AI methods [deep neural networks (DNNs) and recurrent neural networks (RNNs)]. The training set (∼9000 compounds) was compiled by integrating the hERG bioactivity data from the ChEMBL database with experimental data generated from an in-house, high-throughput thallium flux assay. We utilized different molecular descriptors including the latent descriptors, which are real-value continuous vectors derived from chemical autoencoders trained on a large chemical space (>1.5 million compounds). The models were prospectively validated on ∼840 in-house compounds screened in the same thallium flux assay. The best results were obtained with the XGBoost method and RDKit descriptors. The comparison of models based only on latent descriptors revealed that the DNNs performed significantly better than the classical methods. The RNNs that operate on SMILES provided the highest model sensitivity. The best models were merged into a consensus model that offered superior performance compared to reference models from academic and commercial domains. Furthermore, we shed light on the potential of AI methods to exploit the big data in chemistry and generate novel chemical representations useful in predictive modeling and tailoring a new chemical space.

A current understanding of drug-induced QT prolongation and its implications for anticancer therapy

The QT interval, a global index of ventricular repolarization, varies among individuals and is influenced by diverse physiologic and pathophysiologic stimuli such as gender, age, heart rate, electrolyte concentrations, concomitant cardiac disease, and other diseases such as diabetes. Many drugs produce a small but reproducible effect on QT interval but in rare instances this is exaggerated and marked QT prolongation can provoke the polymorphic ventricular tachycardia 'torsades de pointes', which can cause syncope or sudden cardiac death. The generally accepted common mechanism whereby drugs prolong QT is block of a key repolarizing potassium current in heart, IKr, generated by expression of KCNH2, also known as HERG. Thus, evaluation of the potential that a new drug entity may cause torsades de pointes has relied on exposure of normal volunteers or patients to drug at usual and high concentrations, and on assessment of IKr block in vitro. More recent work, focusing on anticancer drugs with QT prolonging liability, is defining new pathways whereby drugs can prolong QT. Notably, the in vitro effects of some tyrosine kinase inhibitors to prolong cardiac action potentials (the cellular correlate of QT) can be rescued by intracellular phosphatidylinositol 3,4,5-trisphosphate, the downstream effector of phosphoinositide 3-kinase. This finding supports a role for inhibition of this enzyme, either directly or by inhibition of upstream kinases, to prolong QT through mechanisms that are being worked out, but include enhanced inward 'late' sodium current during the plateau of the action potential. The definition of non-IKr-dependent pathways to QT prolongation will be important for assessing risk, not only with anticancer therapies but also with other QT prolonging drugs and for generating a refined understanding how variable activity of intracellular signalling systems can modulate QT and associated arrhythmia risk.

Structures Illuminate Cardiac Ion Channel Functions in Health and in Long QT Syndrome

The cardiac action potential is critical to the production of a synchronized heartbeat. This electrical impulse is governed by the intricate activity of cardiac ion channels, among them the cardiac voltage-gated potassium (Kv) channels KCNQ1 and hERG as well as the voltage-gated sodium (Nav) channel encoded by SCN5A. Each channel performs a highly distinct function, despite sharing a common topology and structural components. These three channels are also the primary proteins mutated in congenital long QT syndrome (LQTS), a genetic condition that predisposes to cardiac arrhythmia and sudden cardiac death due to impaired repolarization of the action potential and has a particular proclivity for reentrant ventricular arrhythmias. Recent cryo-electron microscopy structures of human KCNQ1 and hERG, along with the rat homolog of SCN5A and other mammalian sodium channels, provide atomic-level insight into the structure and function of these proteins that advance our understanding of their distinct functions in the cardiac action potential, as well as the molecular basis of LQTS. In this review, the gating, regulation, LQTS mechanisms, and pharmacological properties of KCNQ1, hERG, and SCN5A are discussed in light of these recent structural findings.

The EAG Voltage-Dependent K+ Channel Subfamily: Similarities and Differences in Structural Organization and Gating

EAG (ether-à-go-go or KCNH) are a subfamily of the voltage-gated potassium (Kv) channels. Like for all potassium channels, opening of EAG channels drives the membrane potential toward its equilibrium value for potassium, thus setting the resting potential and repolarizing action potentials. As voltage-dependent channels, they switch between open and closed conformations (gating) when changes in membrane potential are sensed by a voltage sensing domain (VSD) which is functionally coupled to a pore domain (PD) containing the permeation pathway, the potassium selectivity filter, and the channel gate. All Kv channels are tetrameric, with four VSDs formed by the S1-S4 transmembrane segments of each subunit, surrounding a central PD with the four S5-S6 sections arranged in a square-shaped structure. Structural information, mutagenesis, and functional experiments, indicated that in "classical/Shaker-type" Kv channels voltage-triggered VSD reorganizations are transmitted to PD gating via the α-helical S4-S5 sequence that links both modules. Importantly, these Shaker-type channels share a domain-swapped VSD/PD organization, with each VSD contacting the PD of the adjacent subunit. In this case, the S4-S5 linker, acting as a rigid mechanical lever (electromechanical lever coupling), would lead to channel gate opening at the cytoplasmic S6 helices bundle. However, new functional data with EAG channels split between the VSD and PD modules indicate that, in some Kv channels, alternative VSD/PD coupling mechanisms do exist. Noticeably, recent elucidation of the architecture of some EAG channels, and other relatives, showed that their VSDs are non-domain swapped. Despite similarities in primary sequence and predicted structural organization for all EAG channels, they show marked kinetic differences whose molecular basis is not completely understood. Thus, while a common general architecture may establish the gating system used by the EAG channels and the physicochemical coupling of voltage sensing to gating, subtle changes in that common structure, and/or allosteric influences of protein domains relatively distant from the central gating machinery, can crucially influence the gating process. We consider here the latest advances on these issues provided by the elucidation of eag1 and erg1 three-dimensional structures, and by both classical and more recent functional studies with different members of the EAG subfamily.

Revealing Molecular Determinants of hERG Blocker and Activator Binding

The Kv11.1 potassium channel, encoded by the human ether-a-go-go-related gene (hERG), plays an essential role in the cardiac action potential. hERG blockade by small molecules can induce "torsade de pointes" arrhythmias and sudden death; as such, it is an important off-target to avoid during drug discovery. Recently, a cryo-EM structure of the open channel state of hERG was reported, opening the door to in silico docking analyses and interpretation of hERG structure-activity relationships, with a view to avoiding blocking activity. Despite this, docking directly to this cryo-EM structure has been reported to yield binding modes that are unable to explain known mutagenesis data. In this work, we use molecular dynamics simulations to sample a range of channel conformations and run ensemble docking campaigns at the known hERG binding site below the selectivity filter, composed of the central cavity and the four deep hydrophobic pockets. We identify a hERG conformational state allowing discrimination of blockers vs nonblockers from docking; furthermore, the binding pocket agrees with mutagenesis data, and blocker binding modes fit the hERG blocker pharmacophore. We then use the same protocol to identify a binding pocket in the hERG channel pore for hERG activators, again agreeing with the reported mutagenesis. Our approach may be useful in drug discovery campaigns to prioritize candidate compounds based on hERG liability via virtual docking screens.

Human ether-à-go-go-related potassium channel: exploring SAR to improve drug design

hERG is best known as a primary anti-target, the inhibition of which is responsible for serious side effects. A renewed interest in hERG as a desired target, especially in oncology, was sparked because of its role in cellular proliferation and apoptosis. In this study, we survey the most recent advances regarding hERG by focusing on SAR in the attempt to elucidate, at a molecular level, off-target and on-target actions of potential hERG binders, which are highly promiscuous and largely varying in structure. Understanding the rationale behind hERG interactions and the molecular determinants of hERG activity is a real challenge and comprehension of this is of the utmost importance to prioritize compounds in early stages of drug discovery and to minimize cardiotoxicity attrition in preclinical and clinical studies.

IK1 channel agonist zacopride suppresses ventricular arrhythmias in conscious rats with healing myocardial infarction

AIMS

Arrhythmogenesis of chronic myocardial infarction (MI) is associated with the prolongation of action potential, reduction of inward rectifier potassium (I_K1, Kir) channels and hyper-activity of Calcium/calmodulin-dependent kinase II (CaMKII) in cardiomyocytes. Zacopride, a selective I_K1 agonist, was applied to clarify the cardioprotection of I_K1 agonism via a CaMKII signaling on arrhythmias post-MI.

METHODS

Male SD rats were implanted wireless transmitter in the abdominal cavity and subjected to left main coronary artery ligation or sham operation. The telemetric ECGs were monitored per day throughout 4 weeks. At the endpoint, isoproterenol (1.28 mg/kg, i.v.) was administered for provocation test. The expressions of Kir2.1 (dominant subunit of I_K1 in ventricle) and CaMKII were detected by Western-blotting.

KEY FINDINGS

In the telemetric rats post-MI, zacopride significantly reduced the episodes of atrioventricular conduction block (AVB), premature ventricular contraction (PVC), ventricular tachycardia (VT) and ventricular fibrillation (VF), without significant effect on superventricular premature contraction (SPVC). In provocation test, zacopride suppressed the onset of ventricular arrhythmias in conscious PMI or sham rats. The expression of Kir2.1 was significantly downregulated and p-CaMKII was upregulated post-MI, whereas both were restored by zacopride treatment.

SIGNIFICANCE

I_K1/Kir2.1 might be an attractive target for pharmacological controlling of lethal arrhythmias post MI.

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.

Cardiac toxicity of Triptergium wilfordii Hook F. may correlate with its inhibition to hERG channel

Tripterygium wilfordii Hook F. (TWHF) is a Chinese traditional medicine with cardiac toxicities. However, the mechanism of acute cardiac toxicity is not very clear. By using patch clamp techniques, we found that 0.05 mg/ml and 0.1 mg/ml of the aqueous crude extract of TWHF inhibit 21.4 ± 1.6% and 86.7 ± 5.7% (n = 5) of hERG current Amplitudes (I_hERG) respectively. We further found that Celastrol, one of main components of TWHF, inhibits hERG with an IC₅₀ of 0.83 μM. Additional mutagenesis studies show that mutations of T623A, S624A and F656A significantly alter the inhibition and S624A has the strongest effect, supported by our docking model. Our data suggest that inhibition of hERG channel activity by Celastrol contributed to TWHF cardiotoxicity.

A pharmacological master key mechanism that unlocks the selectivity filter gate in K+ channels

Potassium (K⁺) channels have been evolutionarily tuned for activation by diverse biological stimuli, and pharmacological activation is thought to target these specific gating mechanisms. Here we report a class of negatively charged activators (NCAs) that bypass the specific mechanisms but act as master keys to open K⁺ channels gated at their selectivity filter (SF), including many two-pore domain K⁺ (K_2P) channels, voltage-gated hERG (human ether-à-go-go-related gene) channels and calcium (Ca²⁺)-activated big-conductance potassium (BK)-type channels. Functional analysis, x-ray crystallography, and molecular dynamics simulations revealed that the NCAs bind to similar sites below the SF, increase pore and SF K⁺ occupancy, and open the filter gate. These results uncover an unrecognized polypharmacology among K⁺ channel activators and highlight a filter gating machinery that is conserved across different families of K⁺ channels with implications for rational drug design.

Cryo-Electron Microscopy: Moving Beyond X-Ray Crystal Structures for Drug Receptors and Drug Development

Electron cryo-microscopy (cryo-EM) has revolutionized structure determination of membrane proteins and holds great potential for structure-based drug discovery. Here we discuss the potential of cryo-EM in the rational design of therapeutics for membrane proteins compared to X-ray crystallography. We also detail recent progress in the field of drug receptors, focusing on cryo-EM of two protein families with established therapeutic value, the γ-aminobutyric acid A receptors (GABA_ARs) and G protein-coupled receptors (GPCRs). GABA_ARs are pentameric ion channels, and cryo-EM structures of physiological heteromeric receptors in a lipid environment have uncovered the molecular basis of receptor modulation by drugs such as diazepam. The structures of ten GPCR-G protein complexes from three different classes of GPCRs have now been determined by cryo-EM. These structures give detailed insights into molecular interactions with drugs, GPCR-G protein selectivity, how accessory membrane proteins alter receptor-ligand pharmacology, and the mechanism by which HIV uses GPCRs to enter host cells.

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.

Importance of the Choice of a Recombinant System to Produce Large Amounts of Functional Membrane Protein hERG

Human ether-a-gogo related gene (hERG) product is the membrane potassium channel Kv11.1, which is involved in the electrical activity of the heart. As such, it is a key player in the toxicity of many drug candidates. Therefore, having this protein at hand during earlier stages of drug discovery is important for preventing later toxicity. Furthermore, having a fair quantity of functional channels may help in the development of the necessary techniques for gaining insight in this channel structure. Thus, we performed a comparative study of methods for over-expressing a mutated but functional, hERG in different orthologous hosts, such as yeast, bacteria, insect and human cell lines. We also engineered the protein to test various constructs of a functional channel. We obtained a significant amount of a functional mutant channel from HEK cells that we thoroughly characterized. The present work paves the way for the expression of large amounts of this protein, with which protein crystallization or cryo-electronic microscopy will be attempted. This will be a way to gain information on the structure of the hERG active site and its modelization to obtain data on the pauses of various reference compounds from the pharmacopeia, as well as to gain information about the thermodynamics of the hERG/ligand relationship.

The Pore-Lipid Interface: Role of Amino-Acid Determinants of Lipophilic Access by Ivabradine to the hERG1 Pore Domain

Abnormal cardiac electrical activity is a common side effect caused by unintended block of the promiscuous drug target human ether-à-go-go-related gene (hERG1), the pore-forming domain of the delayed rectifier K⁺ channel in the heart. hERG1 block leads to a prolongation of the QT interval, a phase of the cardiac cycle that underlies myocyte repolarization detectable on the electrocardiogram. Even newly released drugs such as heart-rate lowering agent ivabradine block the rapid delayed rectifier current I_Kr, prolong action potential duration, and induce potentially lethal arrhythmia known as torsades de pointes. In this study, we describe a critical drug-binding pocket located at the lateral pore surface facing the cellular membrane. Mutations of the conserved M651 residue alter ivabradine-induced block but not by the common hERG1 blocker dofetilide. As revealed by molecular dynamics simulations, binding of ivabradine to a lipophilic pore access site is coupled to a state-dependent reorientation of aromatic residues F557 and F656 in the S5 and S6 helices. We show that the M651 mutation impedes state-dependent dynamics of F557 and F656 aromatic cassettes at the protein-lipid interface, which has a potential to disrupt drug-induced block of the channel. This fundamentally new mechanism coupling the channel dynamics and small-molecule access from the membrane into the hERG1 intracavitary site provides a simple rationale for the well established state-dependence of drug blockade. SIGNIFICANCE STATEMENT: The drug interference with the function of the cardiac hERG channels represents one of the major sources of drug-induced heart disturbances. We found a novel and a critical drug-binding pocket adjacent to a lipid-facing surface of the hERG1 channel, which furthers our molecular understanding of drug-induced QT syndrome.

Antibodies and venom peptides: new modalities for ion channels

Ion channels play fundamental roles in both excitable and non-excitable tissues and therefore constitute attractive drug targets for myriad neurological, cardiovascular and metabolic diseases as well as for cancer and immunomodulation. However, achieving selectivity for specific ion channel subtypes with small-molecule drugs has been challenging, and there currently is a growing trend to target ion channels with biologics. One approach is to improve the pharmacokinetics of existing or novel venom-derived peptides. In parallel, after initial studies with polyclonal antibodies demonstrated the technical feasibility of inhibiting channel function with antibodies, multiple preclinical programmes are now using the full spectrum of available technologies to generate conventional monoclonal and engineered antibodies or nanobodies against extracellular loops of ion channels. After a summary of the current state of ion channel drug discovery, this Review discusses recent developments using the purinergic receptor channel P2X purinoceptor 7 (P2X7), the voltage-gated potassium channel KV1.3 and the voltage-gated sodium channel NaV1.7 as examples of targeting ion channels with biologics.

New Structures and Gating of Voltage-Dependent Potassium (Kv) Channels and Their Relatives: A Multi-Domain and Dynamic Question

Voltage-dependent potassium channels (Kv channels) are crucial regulators of cell excitability that participate in a range of physiological and pathophysiological processes. These channels are molecular machines that display a mechanism (known as gating) for opening and closing a gate located in a pore domain (PD). In Kv channels, this mechanism is triggered and controlled by changes in the magnitude of the transmembrane voltage sensed by a voltage-sensing domain (VSD). In this review, we consider several aspects of the VSD–PD coupling in Kv channels, and in some relatives, that share a common general structure characterized by a single square-shaped ion conduction pore in the center, surrounded by four VSDs located at the periphery. We compile some recent advances in the knowledge of their architecture, based in cryo-electron microscopy (cryo-EM) data for high-resolution determination of their structure, plus some new functional data obtained with channel variants in which the covalent continuity between the VSD and PD modules has been interrupted. These advances and new data bring about some reconsiderations about the use of exclusively a classical electromechanical lever model of VSD–PD coupling by some Kv channels, and open a view of the Kv-type channels as allosteric machines in which gating may be dynamically influenced by some long-range interactional/allosteric mechanisms.

Bipolar switching by HCN voltage sensor underlies hyperpolarization activation

Hyperpolarization-activated, cyclic nucleotide-gated (HCN) channels show an inverted voltage response compared with virtually all other voltage-gated channels, opening on hyperpolarization rather than depolarization. Although the structure of the HCN1 channel was recently solved, the structural element(s) responsible for the inverted gating polarity of HCN is not known. Here, we use a hierarchical approach, by first characterizing the functional contribution of each structural element to channel gating, and then identifying the critical interactions between these elements. Our studies reveal that the HCN voltage sensor can gate the same pore open on both depolarization and hyperpolarization, thereby acting as a bipolar switch. Elements in the pore domain shut off the depolarization-activation pathway in wild-type channels.

Despite sharing a common architecture with archetypal voltage-gated ion channels (VGICs), hyperpolarization- and cAMP-activated ion (HCN) channels open upon hyperpolarization rather than depolarization. The basic motions of the voltage sensor and pore gates are conserved, implying that these domains are inversely coupled in HCN channels. Using structure-guided protein engineering, we systematically assembled an array of mosaic channels that display the full complement of voltage-activation phenotypes observed in the VGIC superfamily. Our studies reveal that the voltage sensor of the HCN channel has an intrinsic ability to drive pore opening in either direction and that the extra length of the HCN S4 is not the primary determinant for hyperpolarization activation. Tight interactions at the HCN voltage sensor–pore interface drive the channel into an hERG-like inactivated state, thereby obscuring its opening upon depolarization. This structural element in synergy with the HCN cyclic nucleotide-binding domain and specific interactions near the pore gate biases the channel toward hyperpolarization-dependent opening. Our findings reveal an unexpected common principle underpinning voltage gating in the VGIC superfamily and identify the essential determinants of gating polarity.

Challenges and advances in atomistic simulations of potassium and sodium ion channel gating and permeation

Ion channels are implicated in many essential physiological events such as electrical signal propagation and cellular communication. The advent of K⁺ and Na⁺ ion channel structure determination has facilitated numerous investigations of molecular determinants of their behaviour. At the same time, rapid development of computer hardware and molecular simulation methodologies has made computational studies of large biological molecules in all-atom representation tractable. The concurrent evolution of experimental structural biology with biomolecular computer modelling has yielded mechanistic details of fundamental processes unavailable through experiments alone, such as ion conduction and ion channel gating. This review is a short survey of the atomistic computational investigations of K⁺ and Na⁺ ion channels, focusing on KcsA and several voltage-gated channels from the K_V and Na_V families, which have garnered many successes and engendered several long-standing controversies regarding the nature of their structure-function relationship. We review the latest advancements and challenges facing the field of molecular modelling and simulation regarding the structural and energetic determinants of ion channel function and their agreement with experimental observations.

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.

Diurnal Profile of the QTc Interval Following Moxifloxacin Administration

Understanding the physiological fluctuations in the corrected QT (QTc) interval is important to accurately interpret the variations in drug-induced prolongation. The present study aimed to define the time course of the effect of moxifloxacin on the QT interval to understand the duration of the responses to moxifloxacin. This retrospective analysis was performed on data taken from a thorough QT 4-way crossover study with 40 subjects. Each period consisted of a baseline electrocardiogram (ECG) day (day -1) and a treatment day (day 1). On both days, ECGs were recorded simultaneously using 2 different systems operating in parallel: a bedside ECG and a continuous Holter recording. The subjects were randomized to 1 of 4 treatments: 5 mg and 40 mg of intravenous amisulpride, a single oral dose of moxifloxacin (400 mg), or placebo. Standardized meals, identical in all 4 periods, with similar nutritional value were served. Bedside ECG results confirmed that the moxifloxacin peak effect was delayed in the fed state and showed that the Fridericia corrected QT prolongation induced by moxifloxacin persisted until the end of the 24-hour measurement period. The use of continuous Holter monitoring provided further insight, as it revealed that the moxifloxacin effect on QTc was influenced by diurnal and nocturnal environmental factors, and hysteresis effects were noticeable. The findings suggested that moxifloxacin prolongs QTc beyond its elimination from the blood circulation. This is of relevance to current concentration-effect modeling approaches, which presume the absence of hysteresis effects.

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.

Never at rest: insights into the conformational dynamics of ion channels from cryo-electron microscopy

The tightly regulated opening and closure of ion channels underlies the electrical signals that are vital for a wide range of physiological processes. Two decades ago the first atomic level view of ion channel structures led to a detailed understanding of ion selectivity and conduction. In recent years, spectacular developments in the field of cryo-electron microscopy have resulted in cryo-EM superseding crystallography as the technique of choice for determining near-atomic resolution structures of ion channels. Here, we will review the recent developments in cryo-EM and its specific application to the study of ion channel gating. We will highlight the advantages and disadvantages of the current technology and where the field is likely to head in the next few years.