Mechanism of C-type inactivation in the hERG potassium channel

Honeybee CaV4 has distinct permeation, inactivation, and pharmacology from homologous NaV channels

DSC1, a Drosophila channel with sequence similarity to the voltage-gated sodium channel (NaV), was identified over 20 years ago. This channel was suspected to function as a non-specific cation channel with the ability to facilitate the permeation of calcium ions (Ca2+). A honeybee channel homologous to DSC1 was recently cloned and shown to exhibit strict selectivity for Ca2+, while excluding sodium ions (Na+), thus defining a new family of Ca2+ channels, known as CaV4. In this study, we characterize CaV4, showing that it exhibits an unprecedented type of inactivation, which depends on both an IFM motif and on the permeating divalent cation, like NaV and CaV1 channels, respectively. CaV4 displays a specific pharmacology with an unusual response to the alkaloid veratrine. It also possesses an inactivation mechanism that uses the same structural domains as NaV but permeates Ca2+ ions instead. This distinctive feature may provide valuable insights into how voltage- and calcium-dependent modulation of voltage-gated Ca2+ and Na+ channels occur under conditions involving local changes in intracellular calcium concentrations. Our study underscores the unique profile of CaV4 and defines this channel as a novel class of voltage-gated Ca2+ channels.

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

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

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

BACKGROUND

Cyclic Nucleotide-Binding Domain (CNBD)-family channels display distinct voltage-sensing properties despite sharing sequence and structural similarity. For example, the human Ether-a-go-go Related Gene (hERG) channel and the Hyperpolarization-activated Cyclic Nucleotide-gated (HCN) channel share high amino acid sequence similarity and identical domain structures. hERG conducts outward current and is activated by positive membrane potentials (depolarization), whereas HCN conducts inward current and is activated by negative membrane potentials (hyperpolarization). The structural basis for the "opposite" voltage-sensing properties of hERG and HCN remains unknown.

RESULTS

We found the voltage-sensing domain (VSD) involves in modulating the gating polarity of hERG. We identified that a long-QT syndrome type 2-related mutation within the VSD, K525N, mediated an inwardly rectifying non-deactivating current, perturbing the channel closure, but sparing the open state and inactivated state. K525N rescued the current of a non-functional mutation in the pore helix region (F627Y) of hERG. K525N&F627Y switched hERG into a hyperpolarization-activated channel. The reactivated inward current induced by hyperpolarization mediated by K525N&F627Y can be inhibited by E-4031 and dofetilide quite well. Moreover, we report an extracellular interaction between the S1 helix and the S5-P region is crucial for modulating the gating polarity. The alanine substitution of several residues in this region (F431A, C566A, I607A, and Y611A) impaired the inward current of K525N&F627Y.

CONCLUSIONS

Our data provide evidence that a potential cooperation mechanism in the extracellular vestibule of the VSD and the PD would determine the gating polarity in hERG.

Harnessing AlphaFold to reveal state secrets: Prediction of hERG closed and inactivated states

To design safe, selective, and effective new therapies, there must be a deep understanding of the structure and function of the drug target. One of the most difficult problems to solve has been resolution of discrete conformational states of transmembrane ion channel proteins. An example is KV11.1 (hERG), comprising the primary cardiac repolarizing current, IKr. hERG is a notorious drug anti-target against which all promising drugs are screened to determine potential for arrhythmia. Drug interactions with the hERG inactivated state are linked to elevated arrhythmia risk, and drugs may become trapped during channel closure. However, the structural details of multiple conformational states have remained elusive. Here, we guided AlphaFold2 to predict plausible hERG inactivated and closed conformations, obtaining results consistent with myriad available experimental data. Drug docking simulations demonstrated hERG state-specific drug interactions aligning well with experimental results, revealing that most drugs bind more effectively in the inactivated state and are trapped in the closed state. Molecular dynamics simulations demonstrated ion conduction that aligned with earlier studies. Finally, we identified key molecular determinants of state transitions by analyzing interaction networks across closed, open, and inactivated states in agreement with earlier mutagenesis studies. Here, we demonstrate a readily generalizable application of AlphaFold2 as a novel method to predict discrete protein conformations and novel linkages from structure to function.

Anionic omega currents from single countercharge mutants in the voltage-sensing domain of Ci-VSP

The S4 segment of voltage-sensing domains (VSDs) directly responds to voltage changes by reorienting within the electric field as a permion. A narrow hydrophobic "gasket" or charge transfer center at the core of most VSDs focuses the electric field into a narrow region and catalyzes the sequential and reversible translocation of S4 positive gating charge residues across the electric field while preventing the permeation of physiological ions. Mutating specific S4 gating charges can cause ionic leak currents through the VSDs. These gating pores or omega currents play important pathophysiological roles in many diseases of excitability. Here, we show that mutating D129, a key countercharge residue in the Ciona intestinalis voltage-sensing phosphatase (Ci-VSP), leads to the generation of unique anionic omega currents. Neutralizing D129 causes a dramatic positive shift of activation, facilitates the formation of a continuous water path through the VSD, and creates a positive electrostatic potential landscape inside the VSD that contributes to its unique anionic selectivity. Increasing the population or dwell time of the conducting state by a high external pH or an engineered Cd2+ bridge markedly increases the current magnitude. Our findings uncover a new role of countercharge residues in the impermeable VSD of Ci-VSP and offer insights into mechanisms of the conduction of anionic omega currents linked to countercharge residue mutations.

Integrated Approach Including Docking, MD Simulations, and Network Analysis Highlights the Action Mechanism of the Cardiac hERG Activator RPR260243

hERG is a voltage-gated potassium channel involved in the heart contraction whose defections are associated with the cardiac arrhythmia Long QT Syndrome type 2. The activator RPR260243 (RPR) represents a possible candidate to pharmacologically treat LQTS2 because it enhances the opening of the channel. However, the molecular detail of its action mechanism remains quite elusive. Here, we address the problem using a combination of docking, molecular dynamics simulations, and network analysis. We show that the drug preferably binds at the interface between the voltage sensor and the pore, enhancing the canonical activation path and determining a whole-structure rearrangement of the channel that slightly impairs inactivation.

Mechanism underlying delayed rectifying in human voltage-mediated activation Eag2 channel

The transmembrane voltage gradient is a general physico-chemical cue that regulates diverse biological function through voltage-gated ion channels. How voltage sensing mediates ion flows remains unknown at the molecular level. Here, we report six conformations of the human Eag2 (hEag2) ranging from closed, pre-open, open, and pore dilation but non-conducting states captured by cryo-electron microscopy (cryo-EM). These multiple states illuminate dynamics of the selectivity filter and ion permeation pathway with delayed rectifier properties and Cole-Moore effect at the atomic level. Mechanistically, a short S4-S5 linker is coupled with the constrict sites to mediate voltage transducing in a non-domain-swapped configuration, resulting transitions for constrict sites of F464 and Q472 from gating to open state stabilizing for voltage energy transduction. Meanwhile, an additional potassium ion occupied at positions S6 confers the delayed rectifier property and Cole-Moore effects. These results provide insight into voltage transducing and potassium current across membrane, and shed light on the long-sought Cole-Moore effects.

Noncanonical electromechanical coupling paths in cardiac hERG potassium channel

Voltage-gated potassium channels are involved in many physiological processes such as nerve impulse transmission, the heartbeat, and muscle contraction. However, for many of them the molecular determinants of the gating mechanism remain elusive. Here, using a combination of theoretical and experimental approaches, we address this problem focusing on the cardiac hERG potassium channel. Network analysis of molecular dynamics trajectories reveals the presence of a kinematic chain of residues that couples the voltage sensor domain to the pore domain and involves the S4/S1 and S1/S5 subunit interfaces. Mutagenesis experiments confirm the role of these residues and interfaces in the activation and inactivation mechanisms. Our findings demonstrate the presence of an electromechanical transduction path crucial for the non-domain-swapped hERG channel gating that resembles the noncanonical path identified in domain-swapped K+ channels.

Early Steps in C-Type Inactivation of the hERG Potassium Channel

Fast C-type inactivation confers distinctive functional properties to the hERG potassium channel, and its association to inherited and acquired cardiac arrythmias makes the study of the inactivation mechanism of hERG at the atomic detail of paramount importance. At present, two models have been proposed to describe C-type inactivation in K⁺-channels. Experimental data and computational work on the bacterial KcsA channel support the hypothesis that C-type inactivation results from a closure of the selectivity filter that sterically impedes ion conduction. Alternatively, recent experimental structures of a mutated Shaker channel revealed a widening of the extracellular portion of the selectivity filter, which might diminish conductance by interfering with the mechanism of ion permeation. Here, we performed molecular dynamics simulations of the wild-type hERG, a non-inactivating mutant (hERG-N629D), and a mutant that inactivates faster than the wild-type channel (hERG-F627Y) to find out which and if any of the two reported C-type inactivation mechanisms applies to hERG. Closure events of the selectivity filter were not observed in any of the simulated trajectories but instead, the extracellular section of the selectivity filter deviated from the canonical conductive structure of potassium channels. The degree of widening of the potassium binding sites at the extracellular entrance of the channel was directly related to the degree of inactivation with hERG-F627Y > wild-type hERG > hERG-N629D. These findings support the hypothesis that C-type inactivation in hERG entails a widening of the extracellular entrance of the channel rather than a closure of the selectivity filter.

Structural modeling of the hERG potassium channel and associated drug interactions

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

Mechanisms Underlying C-type Inactivation in Kv Channels: Lessons From Structures of Human Kv1.3 and Fly Shaker-IR Channels

Voltage-gated potassium (Kv) channels modulate the function of electrically-excitable and non-excitable cells by using several types of “gates” to regulate ion flow through the channels. An important gating mechanism, C-type inactivation, limits ion flow by transitioning Kv channels into a non-conducting inactivated state. Here, we highlight two recent papers, one on the human Kv1.3 channel and the second on the Drosophila Shaker Kv channel, that combined cryogenic electron microscopy and molecular dynamics simulation to define mechanisms underlying C-type inactivation. In both channels, the transition to the non-conducting inactivated conformation begins with the rupture of an intra-subunit hydrogen bond that fastens the selectivity filter to the pore helix. The freed filter swings outwards and gets tethered to an external residue. As a result, the extracellular end of the selectivity filter dilates and K⁺ permeation through the pore is impaired. Recovery from inactivation may entail a reversal of this process. Such a reversal, at least partially, is induced by the peptide dalazatide. Binding of dalazatide to external residues in Kv1.3 frees the filter to swing inwards. The extracellular end of the selectivity filter narrows allowing K⁺ to move in single file through the pore typical of conventional knock-on conduction. Inter-subunit hydrogen bonds that stabilize the outer pore in the dalazatide-bound structure are equivalent to those in open-conducting conformations of Kv channels. However, the intra-subunit bond that fastens the filter to the pore-helix is absent, suggesting an incomplete reversal of the process. These mechanisms define how Kv channels self-regulate the flow of K⁺ by changing the conformation of the selectivity filter.

A distinct mechanism of C-type inactivation in the Kv-like KcsA mutant E71V

C-type inactivation is of great physiological importance in voltage-activated K⁺ channels (Kv), but its structural basis remains unresolved. Knowledge about C-type inactivation has been largely deduced from the bacterial K⁺ channel KcsA, whose selectivity filter constricts under inactivating conditions. However, the filter is highly sensitive to its molecular environment, which is different in Kv channels than in KcsA. In particular, a glutamic acid residue at position 71 along the pore helix in KcsA is substituted by a valine conserved in most Kv channels, suggesting that this side chain is a molecular determinant of function. Here, a combination of X-ray crystallography, solid-state NMR and MD simulations of the E71V KcsA mutant is undertaken to explore inactivation in this Kv-like construct. X-ray and ssNMR data show that the filter of the Kv-like mutant does not constrict under inactivating conditions. Rather, the filter adopts a conformation that is slightly narrowed and rigidified. On the other hand, MD simulations indicate that the constricted conformation can nonetheless be stably established in the mutant channel. Together, these findings suggest that the Kv-like KcsA mutant may be associated with different modes of C-type inactivation, showing that distinct filter environments entail distinct C-type inactivation mechanisms.

Molecular dynamics simulations suggest possible activation and deactivation pathways in the hERG channel

The elusive activation/deactivation mechanism of hERG is investigated, a voltage-gated potassium channel involved in severe inherited and drug-induced cardiac channelopathies, including the Long QT Syndrome. Firstly, the available structural data are integrated by providing a homology model for the closed state of the channel. Secondly, molecular dynamics combined with a network analysis revealed two distinct pathways coupling the voltage sensor domain with the pore domain. Interestingly, some LQTS-related mutations known to impair the activation/deactivation mechanism are distributed along the identified pathways, which thus suggests a microscopic interpretation of their role. Split channels simulations clarify a surprising feature of this channel, which is still able to gate when a cut is introduced between the voltage sensor domain and the neighboring helix S5. In summary, the presented results suggest possible activation/deactivation mechanisms of non-domain-swapped potassium channels that may aid in biomedical applications.

Costa et al. present the electro-mechanical coupling between the voltage sensor and pore domain of the hERG channel using a combination of molecular dynamics simulations and theoretical network analyses.

Computational methods and theory for ion channel research

Ion channels are fundamental biological devices that act as gates in order to ensure selective ion transport across cellular membranes; their operation constitutes the molecular mechanism through which basic biological functions, such as nerve signal transmission and muscle contraction, are carried out. Here, we review recent results in the field of computational research on ion channels, covering theoretical advances, state-of-the-art simulation approaches, and frontline modeling techniques. We also report on few selected applications of continuum and atomistic methods to characterize the mechanisms of permeation, selectivity, and gating in biological and model channels.

Structural Plasticity of the Selectivity Filter in Cation Channels

Ion channels allow for the passage of ions across biological membranes, which is essential for the functioning of a cell. In pore loop channels the selectivity filter (SF) is a conserved sequence that forms a constriction with multiple ion binding sites. It is becoming increasingly clear that there are several conformations and dynamic states of the SF in cation channels. Here we outline specific modes of structural plasticity observed in the SFs of various pore loop channels: disorder, asymmetry, and collapse. We summarize the multiple atomic structures with varying SF conformations as well as asymmetric and more dynamic states that were discovered recently using structural biology, spectroscopic, and computational methods. Overall, we discuss here that structural plasticity within the SF is a key molecular determinant of ion channel gating behavior.

Widening the scope of constriction

JGP modeling study suggests that selectivity filter constriction is a plausible mechanism for C-type inactivation of the Shaker voltage-gated potassium channel.

Computational study of non-conductive selectivity filter conformations and C-type inactivation in a voltage-dependent potassium channel

C-type inactivation is a time-dependent process of great physiological significance that is observed in a large class of K⁺ channels. Experimental and computational studies of the pH-activated KcsA channel show that the functional C-type inactivated state, for this channel, is associated with a structural constriction of the selectivity filter at the level of the central glycine residue in the signature sequence, TTV(G)YGD. The structural constriction is allosterically promoted by the wide opening of the intracellular activation gate. However, whether this is a universal mechanism for C-type inactivation has not been established with certainty because similar constricted structures have not been observed for other K⁺ channels. Seeking to ascertain the general plausibility of the constricted filter conformation, molecular dynamics simulations of a homology model of the pore domain of the voltage-gated potassium channel Shaker were performed. Simulations performed with an open intracellular gate spontaneously resulted in a stable constricted-like filter conformation, providing a plausible nonconductive state responsible for C-type inactivation in the Shaker channel. While there are broad similarities with the constricted structure of KcsA, the hypothetical constricted-like conformation of Shaker also displays some subtle differences. Interestingly, those are recapitulated by the Shaker-like E71V KcsA mutant, suggesting that the residue at this position along the pore helix plays a pivotal role in determining the C-type inactivation behavior. Free energy landscape calculations show that the conductive-to-constricted transition in Shaker is allosterically controlled by the degree of opening of the intracellular activation gate, as observed with the KcsA channel. The behavior of the classic inactivating W434F Shaker mutant is also characterized from a 10-μs MD simulation, revealing that the selectivity filter spontaneously adopts a nonconductive conformation that is constricted at the level of the second glycine in the signature sequence, TTVGY(G)D.

Genetic and Molecular Aspects of Drug-Induced QT Interval Prolongation

Long QT syndromes can be either acquired or congenital. Drugs are one of the many etiologies that may induce acquired long QT syndrome. In fact, many drugs frequently used in the clinical setting are a known risk factor for a prolonged QT interval, thus increasing the chances of developing torsade de pointes. The molecular mechanisms involved in the prolongation of the QT interval are common to most medications. However, there is considerable inter-individual variability in drug response, thus making the application of personalized medicine a relevant aspect in long QT syndrome, in order to evaluate the risk of every individual from a pharmacogenetic standpoint.

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.

Molecular determinants of pro-arrhythmia proclivity of d- and l-sotalol via a multi-scale modeling pipeline

Drug isomers may differ in their proarrhythmia risk. An interesting example is the drug sotalol, an antiarrhythmic drug comprising d- and l- enantiomers that both block the hERG cardiac potassium channel and confer differing degrees of proarrhythmic risk. We developed a multi-scale in silico pipeline focusing on hERG channel – drug interactions and used it to probe and predict the mechanisms of pro-arrhythmia risks of the two enantiomers of sotalol. Molecular dynamics (MD) simulations predicted comparable hERG channel binding affinities for d- and l-sotalol, which were validated with electrophysiology experiments. MD derived thermodynamic and kinetic parameters were used to build multi-scale functional computational models of cardiac electrophysiology at the cell and tissue scales. Functional models were used to predict inactivated state binding affinities to recapitulate electrocardiogram (ECG) QT interval prolongation observed in clinical data. Our study demonstrates how modeling and simulation can be applied to predict drug effects from the atom to the rhythm for dl-sotalol and also increased proarrhythmia proclivity of d- vs. l-sotalol when accounting for stereospecific beta-adrenergic receptor blocking.