Interactions between voltage sensor and pore domains in a hERG K+ channel model from molecular simulations and the effects of a voltage sensor mutation

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.

Deciphering hERG Mutation in Long QT Syndrome Type 2 Using Antisense Oligonucleotide-Mediated Techniques: Lessons from Cystic Fibrosis

Long QT syndrome type 2 (LQT2) is a genetic disorder caused by mutations in the KCNH2 gene, also known as the human ether-a-go-go-related gene (hERG). Over 30% of hERG mutations result in a premature termination codon (PTC) that triggers a process called nonsense-mediated mRNA decay (NMD), where the mRNA transcript is degraded. NMD is a quality control mechanism that removes faulty mRNA to prevent the translation of truncated proteins. Recent advances in antisense oligonucleotide (ASO) technology in the field of cystic fibrosis (CF) have yielded significant progress, including the ASO-mediated comprehensive characterization of key NMD factors and exon-skipping therapy. These advances have contributed to our understanding of the role of PTC-containing mutations in disease phenotypes and have also led to the development of potentially useful therapeutic strategies. Historically, studies of CF have provided valuable insights for the research on LQT2, particularly concerning increasing the expression of hERG. In this article, we outline the current state of knowledge regarding ASO, NMD, and hERG and discuss the introduction of ASO technology in the CF to elucidate the pathogenic mechanisms through targeting NMD. We also discuss the potential clinical therapeutic benefits and limitations of ASO for the management of LQT2. By drawing on lessons learned from CF research, we explore the potential translational values of these advances into LQT2 studies.

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.

In silico Exploration of Interactions Between Potential COVID-19 Antiviral Treatments and the Pore of the hERG Potassium Channel-A Drug Antitarget

Background: In the absence of SARS-CoV-2 specific antiviral treatments, various repurposed pharmaceutical approaches are under investigation for the treatment of COVID-19. Antiviral drugs considered for this condition include atazanavir, remdesivir, lopinavir-ritonavir, and favipiravir. Whilst the combination of lopinavir and ritonavir has been previously linked to prolongation of the QTc interval on the ECG and risk of torsades de pointes arrhythmia, less is known in this regard about atazanavir, remdesivir, and favipiravir. Unwanted abnormalities of drug-induced QTc prolongation by diverse drugs are commonly mediated by a single cardiac anti-target, the hERG potassium channel. This computational modeling study was undertaken in order to explore the ability of these five drugs to interact with known determinants of drug binding to the hERG channel pore. Methods: Atazanavir, remdesivir, ritonavir, lopinavir and favipiravir were docked to in silico models of the pore domain of hERG, derived from cryo-EM structures of hERG and the closely related EAG channel. Results: Atazanavir was readily accommodated in the open hERG channel pore in proximity to the S6 Y652 and F656 residues, consistent with published experimental data implicating these aromatic residues in atazanavir binding to the channel. Lopinavir, ritonavir, and remdesivir were also accommodated in the open channel, making contacts in a model-dependent fashion with S6 aromatic residues and with residues at the base of the selectivity filter/pore helix. The ability of remdesivir (at 30 μM) to inhibit the channel was confirmed using patch-clamp recording. None of these four drugs could be accommodated in the closed channel structure. Favipiravir, a much smaller molecule, was able to fit within the closed channel and could adopt multiple binding poses in the open channel, but with few simultaneous interactions with key binding residues. Only favipiravir and remdesivir showed the potential to interact with lateral pockets below the selectivity filter of the channel. Conclusions: All the antiviral drugs studied here can, in principle, interact with components of the hERG potassium channel canonical binding site, but are likely to differ in their ability to access lateral binding pockets. Favipiravir's small size and relatively paucity of simultaneous interactions may confer reduced hERG liability compared to the other drugs. Experimental structure-function studies are now warranted to validate these observations.

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.

Electrophysiological characterization of the modified hERGT potassium channel used to obtain the first cryo-EM hERG structure

The voltage-gated hERG (human-Ether-à-go-go Related Gene) K+ channel plays a fundamental role in cardiac action potential repolarization. Loss-of-function mutations or pharmacological inhibition of hERG leads to long QT syndrome, whilst gain-of-function mutations lead to short QT syndrome. A recent open channel cryo-EM structure of hERG represents a significant advance in the ability to interrogate hERG channel structure-function. In order to suppress protein aggregation, a truncated channel construct of hERG (hERGT ) was used to obtain this structure. In hERGT cytoplasmic domain residues 141 to 350 and 871 to 1,005 were removed from the full-length channel protein. There are limited data on the electrophysiological properties of hERGT channels. Therefore, this study was undertaken to determine how hERGT influences channel function at physiological temperature. Whole-cell measurements of hERG current (IhERG ) were made at 37°C from HEK 293 cells expressing wild-type (WT) or hERGT channels. With a standard +20 mV activating command protocol, neither end-pulse nor tail IhERG density significantly differed between WT and hERGT . However, the IhERG deactivation rate was significantly slower for hERGT . Half-maximal activation voltage (V0.5 ) was positively shifted for hERGT by ~+8 mV (p < .05 versus WT), without significant change to the activation relation slope factor. Neither the voltage dependence of inactivation, nor time course of development of inactivation significantly differed between WT and hERGT , but recovery of IhERG from inactivation was accelerated for hERGT (p < .05 versus WT). Steady-state "window" current was positively shifted for hERGT with a modest increase in the window current peak. Under action potential (AP) voltage clamp, hERGT IhERG showed modestly increased current throughout the AP plateau phase with a significant increase in current integral during the AP. The observed consequences for hERGT IhERG of deletion of the two cytoplasmic regions may reflect changes to electrostatic interactions influencing the voltage sensor domain.

The molecular determinants of R-roscovitine block of hERG channels

Human ether-à-go-go-related gene (Kv11.1, or hERG) is a potassium channel that conducts the delayed rectifier potassium current (I_Kr) during the repolarization phase of cardiac action potentials. hERG channels have a larger pore than other K⁺channels and can trap many unintended drugs, often resulting in acquired LQTS (aLQTS). R-roscovitine is a cyclin-dependent kinase (CDK) inhibitor that induces apoptosis in colorectal, breast, prostate, multiple myeloma, other cancer cell lines, and tumor xenografts, in micromolar concentrations. It is well tolerated in phase II clinical trials. R-roscovitine inhibits open hERG channels but does not become trapped in the pore. Two-electrode voltage clamp recordings from Xenopus oocytes expressing wild-type (WT) or hERG pore mutant channels (T623A, S624A, Y652A, F656A) demonstrated that compared to WT hERG, T623A, Y652A, and F656A inhibition by 200 μM R-roscovitine was ~ 48%, 29%, and 73% weaker, respectively. In contrast, S624A hERG was inhibited more potently than WT hERG, with a ~ 34% stronger inhibition. These findings were further supported by the IC₅₀ values, which were increased for T623A, Y652A and F656A (by ~5.5, 2.75, and 42 fold respectively) and reduced 1.3 fold for the S624A mutant. Our data suggest that while T623, Y652, and F656 are critical for R-roscovitine-mediated inhibition, S624 may not be. Docking studies further support our findings. Thus, R-roscovitine’s relatively unique features, coupled with its tolerance in clinical trials, could guide future drug screens.

Functional and pharmacological characterization of an S5 domain hERG mutation associated with short QT syndrome

Congenital short QT syndrome (SQTS) is a repolarization disorder characterized by abbreviated QT intervals, atrial and ventricular arrhythmias and a risk of sudden death. This study characterized a missense mutation (I560T) in the S5 domain of the hERG K⁺ channel that has been associated with variant 1 of the SQTS. Whole cell patch clamp recordings of wild-type (WT) and I560T hERG current (I_hERG) were made at 37 °C from hERG expressing HEK 293 cells, and the structural context of the mutation was investigated using a recently reported cryo-EM structure of hERG. Under conventional voltage clamp, the I560T mutation increased I_hERG amplitude without altering the voltage-dependence of activation, although it accelerated activation time-course and also slowed deactivation time-course at some voltages. The voltage dependence of I_hERG inactivation was positively shifted (by ∼24 mV) and the time-course of inactivation was slowed by the I560T mutation. There was also a modest decrease in K⁺ over Na⁺ ion selectivity with the I560T mutation. Under action potential (AP) voltage clamp, the net charge carried by hERG was significantly increased during ventricular, Purkinje fibre and atrial APs, with maximal I_hERG also occurring earlier during the plateau phase of ventricular and Purkinje fibre APs. The I560T mutation exerted only a modest effect on quinidine sensitivity of I_hERG: the IC₅₀ for mutant I_hERG was 2.3 fold that for WT I_hERG under conventional voltage clamp. Under AP voltage clamp the inhibitory effect of 1 μM quinidine was largely retained for I560T hERG and the timing of peak I560T I_hERG was altered towards that of the WT channel. In both the open channel structure and a closed hERG channel model based on the closely-related EAG structure, I560T side-chains were oriented towards membrane lipid and away from adjacent domains of the channel, contrasting with previous predictions based on homology modelling. In summary, the I560T mutation produces multiple effects on hERG channel operation that result in a gain-of-function that is expected to abbreviate ventricular, atrial and Purkinje fibre repolarization. Quinidine is likely to be of value in offsetting the increase in I_hERG and altered I_hERG timing during ventricular APs in SQTS with this mutation.

The S1 helix critically regulates the finely tuned gating of Kv11.1 channels

Congenital mutations in the cardiac Kv11.1 channel can cause long QT syndrome type 2 (LQTS2), a heart rhythm disorder associated with sudden cardiac death. Mutations act either by reducing protein expression at the membrane and/or by perturbing the intricate gating properties of Kv11.1 channels. A number of clinical LQTS2-associated mutations have been reported in the first transmembrane segment (S1) of Kv11.1 channels, but the role of this region of the channel is largely unexplored. In part, this is due to problems defining the extent of the S1 helix, as a consequence of its low sequence homology with other Kv family members. Here, we used NMR spectroscopy and electrophysiological characterization to show that the S1 of Kv11.1 channels extends seven helical turns, from Pro-405 to Phe-431, and is flanked by unstructured loops. Functional analysis suggests that pre-S1 loop residues His-402 and Tyr-403 play an important role in regulating the kinetics and voltage dependence of channel activation and deactivation. Multiple residues within the S1 helix also play an important role in fine-tuning the voltage dependence of activation, regulating slow deactivation, and modulating C-type inactivation of Kv11.1 channels. Analyses of LQTS2-associated mutations in the pre-S1 loop or S1 helix of Kv11.1 channels demonstrate perturbations to both protein expression and most gating transitions. Thus, S1 region mutations would reduce both the action potential repolarizing current passed by Kv11.1 channels in cardiac myocytes, as well as the current passed in response to premature depolarizations that normally helps protect against the formation of ectopic beats.

Interactions between amiodarone and the hERG potassium channel pore determined with mutagenesis and in silico docking

The antiarrhythmic drug amiodarone delays cardiac repolarisation through inhibition of hERG-encoded potassium channels responsible for the rapid delayed rectifier potassium current (I_Kr). This study aimed to elucidate molecular determinants of amiodarone binding to the hERG channel. Whole-cell patch-clamp recordings were made at 37 °C of ionic current (I_hERG) carried by wild-type (WT) or mutant hERG channels expressed in HEK293 cells. Alanine mutagenesis and ligand docking were used to investigate the roles of pore cavity amino-acid residues in amiodarone binding. Amiodarone inhibited WT outward I_hERG tails with a half-maximal inhibitory concentration (IC₅₀) of ∼45 nM, whilst inward I_hERG tails in a high K⁺ external solution ([K⁺]_e) of 94 mM were blocked with an IC₅₀ of 117.8 nM. Amiodarone’s inhibitory action was contingent upon channel gating. Alanine-mutagenesis identified multiple residues directly or indirectly involved in amiodarone binding. The IC₅₀ for the S6 aromatic Y652A mutation was increased to ∼20-fold that of WT I_hERG, similar to the pore helical mutant S624A (∼22-fold WT control). The IC₅₀ for F656A mutant I_hERG was ∼17-fold its corresponding WT control. Computational docking using a MthK-based hERG model differentiated residues likely to interact directly with drug and those whose Ala mutation may affect drug block allosterically. The requirements for amiodarone block of aromatic residues F656 and Y652 within the hERG pore cavity are smaller than for other high affinity I_hERG inhibitors, with relative importance to amiodarone binding of the residues investigated being S624A ∼ Y652A > F656A > V659A > G648A > T623A.

NS1643 interacts around L529 of hERG to alter voltage sensor movement on the path to activation

Activators of hERG1 such as NS1643 are being developed for congenital/acquired long QT syndrome. Previous studies identify the neighborhood of L529 around the voltage-sensor as a putative interacting site for NS1643. With NS1643, the V1/2 of activation of L529I (-34 ± 4 mV) is similar to wild-type (WT) (-37 ± 3 mV; P > 0.05). WT and L529I showed no difference in the slope factor in the absence of NS1643 (8 ± 0 vs. 9 ± 0) but showed a difference in the presence of NS1643 (9 ± 0.3 vs. 22 ± 1; P < 0.01). Voltage-clamp-fluorimetry studies also indicated that in L529I, NS1643 reduces the voltage-sensitivity of S4 movement. To further assess mechanism of NS1643 action, mutations were made in this neighborhood. NS1643 shifts the V1/2 of activation of both K525C and K525C/L529I to hyperpolarized potentials (-131 ± 4 mV for K525C and -120 ± 21 mV for K525C/L529I). Both K525C and K525C/K529I had similar slope factors in the absence of NS1643 (18 ± 2 vs. 34 ± 5, respectively) but with NS1643, the slope factor of K525C/L529I increased from 34 ± 5 to 71 ± 10 (P < 0.01) whereas for K525C the slope factor did not change (18 ± 2 at baseline and 16 ± 2 for NS1643). At baseline, K525R had a slope factor similar to WT (9 vs. 8) but in the presence of NS1643, the slope factor of K525R was increased to 24 ± 4 vs. 9 ± 0 mV for WT (P < 0.01). Molecular modeling indicates that L529I induces a kink in the S4 voltage-sensor helix, altering a salt-bridge involving K525. Moreover, docking studies indicate that NS1643 binds to the kinked structure induced by the mutation with a higher affinity. Combining biophysical, computational, and electrophysiological evidence, a mechanistic principle governing the action of some activators of hERG1 channels is proposed.

Arginine side chain interactions and the role of arginine as a gating charge carrier in voltage sensitive ion channels

Gating charges in voltage-sensing domains (VSD) of voltage-sensitive ion channels and enzymes are carried on arginine side chains rather than lysine. This arginine preference may result from the unique hydration properties of the side chain guanidinium group which facilitates its movement through a hydrophobic plug that seals the center of the VSD, as suggested by molecular dynamics simulations. To test for side chain interactions implicit in this model we inspected interactions of the side chains of arginine and lysine with each of the 19 non-glycine amino acids in proteins in the protein data bank. The arginine guanidinium interacts with non-polar aromatic and aliphatic side chains above and below the guanidinium plane while hydrogen bonding with polar side chains is restricted to in-plane positions. In contrast, non-polar side chains interact largely with the aliphatic part of the lysine side chain. The hydration properties of arginine and lysine are strongly reflected in their respective interactions with non-polar and polar side chains as observed in protein structures and in molecular dynamics simulations, and likely underlie the preference for arginine as a mobile charge carrier in VSD.

Voltage sensor gating charge transfer in a hERG potassium channel model

Relaxation of a hERG K⁺ channel model during molecular-dynamics simulation in a hydrated POPC bilayer was accompanied by transitions of an arginine gating charge across a charge transfer center in two voltage sensor domains. Inspection of the passage of arginine side chains across the charge transfer center suggests that the unique hydration properties of the arginine guanidine cation facilitates charge transfer during voltage sensor responses to changes in membrane potential, and underlies the preference of Arg over Lys as a mobile charge carrier in voltage-sensitive ion channels.

Getting to the heart of hERG K(+) channel gating

Potassium ion channels encoded by the human ether-a-go-go related gene (hERG) form the ion-conducting subunit of the rapid delayed rectifier potassium current (IKr ). Although hERG channels exhibit a widespread tissue distribution they play a particularly important role in the heart. There has been considerable interest in hERG K(+) channels for three main reasons. First, they have very unusual gating kinetics, most notably rapid and voltage-dependent inactivation coupled to slow deactivation, which has led to the suggestion that they may play a specific role in the suppression of arrhythmias. Second, mutations in hERG are the cause of 30-40% of cases of congenital long QT syndrome (LQTS), the commonest inherited primary arrhythmia syndrome. Third, hERG is the molecular target for the vast majority of drugs that cause drug-induced LQTS, the commonest cause of drug-induced arrhythmias and cardiac death. Drug-induced LQTS has now been reported for a large range of both cardiac and non-cardiac drugs, in which this side effect is entirely undesired. In recent years there have been comprehensive reviews published on hERG K(+) channels (Vandenberg et al. 2012) and we will not re-cover this ground. Rather, we focus on more recent work on the structural basis and dynamics of hERG gating with an emphasis on how the latest developments may facilitate translational research in the area of stratifying risk of arrhythmias.

Assessing hERG pore models as templates for drug docking using published experimental constraints: the inactivated state in the context of drug block

Many structurally and therapeutically diverse drugs interact with the human heart K⁺ channel hERG by binding within the K⁺ permeation pathway of the open channel, leading to drug-induced ‘long QT syndrome’. Drug binding to hERG is often stabilized by inactivation gating. In the absence of a crystal structure, hERG pore homology models have been used to characterize drug interactions. Here we assess potentially inactivated states of the bacterial K⁺ channel, KcsA, as templates for inactivated state hERG pore models in the context of drug binding using computational docking. Although Flexidock and GOLD docking produced low energy score poses in the models tested, each method selected a MthK K⁺ channel-based model over models based on the putative inactivated state KcsA structures for each of the 9 drugs tested. The variety of docking poses found indicates that an optimal arrangement for drug binding of aromatic side chains in the hERG pore can be achieved in several different configurations. This plasticity of the drug “binding site” is likely to be a feature of the hERG inactivated state. The results demonstrate that experimental data on specific drug interactions can be used as structural constraints to assess and refine hERG homology models.