A molecular switch driving inactivation in the cardiac K+ channel HERG

Mutation in pore-helix modulates interplay between filter gate and Ba2+ block in a Kcv channel pore

The selectivity filter of K+ channels catalyzes a rapid and highly selective transport of K+ while serving as a gate. To understand the control of this filter gate, we use the pore-only K+ channel KcvNTS in which gating is exclusively determined by the activity of the filter gate. It has been previously shown that a mutation at the C-terminus of the pore-helix (S42T) increases K+ permeability and introduces distinct voltage-dependent and K+-sensitive channel closures at depolarizing voltages. Here, we report that the latter are not generated by intrinsic conformational changes of the filter gate but by a voltage-dependent block caused by nanomolar trace contaminations of Ba2+ in the KCl solution. Channel closures can be alleviated by extreme positive voltages and they can be completely abolished by the high-affinity Ba2+ chelator 18C6TA. By contrast, the same channel closures can be augmented by adding Ba2+ at submicromolar concentrations to the cytosolic buffer. These data suggest that a conservative exchange of Ser for Thr in a crucial position of the filter gate increases the affinity of the filter for Ba2+ by >200-fold at positive voltages. While Ba2+ ions apparently remain only for a short time in the filter-binding sites of the WT channel before passing the pore, they remain much longer in the mutant channel. Our findings suggest that the dwell times of permeating and blocking ions in the filter-binding sites are tightly controlled by interactions between the pore-helix and 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.

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.

Toward a Structural View of hERG Activation by the Small-Molecule Activator ICA-105574

Outward current conducted by human ether-à-go-go-related gene type 1 (hERG1) K+ channels is important for action potential repolarization in the human ventricle. Rapid, voltage-dependent inactivation greatly reduces outward currents conducted by hERG1 channels and involves conformational changes in the ion selectivity filter (SF). Recently, compounds have been found that activate hERG1 channel function by modulating gating mechanisms such as reducing inactivation. Such activating compounds could represent a novel approach to prevent arrhythmias associated with prolonged ventricular repolarization associated with inherited or acquired long QT syndrome. ICA-105574 (ICA), a 3-nitro-n-(4-phenoxyphenyl) benzamide derivative activates hERG1 by strongly attenuating pore-type inactivation. We previously mapped the putative binding site for ICA to a hydrophobic pocket located between two adjacent subunits. Here, we used the recently reported cryoelectron microscopy structures of hERG1 to elucidate the structural mechanisms by which ICA influences the stability of the SF. By combining molecular dynamics simulations, voltage-clamp electrophysiology, and the synthesis of novel ICA derivatives, we provide atomistic insights into SF dynamics and propose a structural link between the SF and S6 segments. Further, our study highlights the importance of the nitro moiety, at the meta position of the benzamide ring, for the activity of ICA and reveals that the (bio)isosteric substitution of this side chain can switch the activity to weak inhibitors. Our findings indicate that ICA increases the stability of the SF to attenuate channel inactivation, and this action requires a fine-tuned compound geometry.

Selectivity filter modalities and rapid inactivation of the hERG1 channel

The human ether-á-go-go-related gene (hERG1) channel conducts small outward K⁺ currents that are critical for cardiomyocyte membrane repolarization. The gain-of-function mutation N629D at the outer mouth of the selectivity filter (SF) disrupts inactivation and K⁺-selective transport in hERG1, leading to arrhythmogenic phenotypes associated with long-QT syndrome. Here, we combined computational electrophysiology with Markov state model analysis to investigate how SF-level gating modalities control selective cation transport in wild-type (WT) and mutant (N629D) hERG1 variants. Starting from the recently reported cryogenic electron microscopy (cryo-EM) open-state channel structure, multiple microseconds-long molecular-dynamics (MD) trajectories were generated using different cation configurations at the filter, voltages, electrolyte concentrations, and force-field parameters. Most of the K⁺ permeation events observed in hERG1-WT simulations occurred at microsecond timescales, influenced by the spontaneous dehydration/rehydration dynamics at the filter. The SF region displayed conductive, constricted, occluded, and dilated states, in qualitative agreement with the well-documented flickering conductance of hERG1. In line with mutagenesis studies, these gating modalities resulted from dynamic interaction networks involving residues from the SF, outer-mouth vestibule, P-helices, and S5-P segments. We found that N629D mutation significantly stabilizes the SF in a state that is permeable to both K⁺ and Na⁺, which is reminiscent of the SF in the nonselective bacterial NaK channel. Increasing the external K⁺ concentration induced "WT-like" SF dynamics in N629D, in qualitative agreement with the recovery of flickering currents in experiments. Overall, our findings provide an understanding of the molecular mechanisms controlling selective transport in K⁺ channels with a nonconventional SF sequence.

Dynamics of the EAG1 K+ channel selectivity filter assessed by molecular dynamics simulations

EAG1 channels belong to the KCNH family of voltage gated potassium channels. They are expressed in several brain regions and increased expression is linked to certain cancer types. Recent cryo-EM structure determination finally revealed the structure of these channels in atomic detail, allowing computational investigations. In this study, we performed molecular dynamics simulations to investigate the ion binding sites and the dynamical behavior of the selectivity filter. Our simulations suggest that sites S2 and S4 form stable ion binding sites, while ions placed at sites S1 and S3 rapidly switched to sites S2 and S4. Further, ions tended to dissociate away from S0 within less than 20 ns, due to increased filter flexibility. This was followed by water influx from the extracellular side, leading to a widening of the filter in this region, and likely non-conductive filter configurations. Simulations with the inactivation-enhancing mutant Y464A or Na⁺ ions lead to trapped water molecules behind the SF, suggesting that these simulations captured early conformational changes linked to C-type inactivation.

New potential binding determinant for hERG channel inhibitors

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

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.

Cooperative subunit interactions mediate fast C-type inactivation of hERG1 K+ channels

At depolarized membrane potentials, the conductance of some voltage-gated K(+) channels is reduced by C-type inactivation. This gating process is voltage independent in Kv1 and involves a conformational change in the selectivity filter that is mediated by cooperative subunit interactions. C-type inactivation in hERG1 K(+) channels is voltage-dependent, much faster in onset and greatly attenuates currents at positive potentials. Here we investigate the potential role of subunit interactions in C-type inactivation of hERG1 channels. Point mutations in hERG1 known to eliminate (G628C/S631C), inhibit (S620T or S631A) or enhance (T618A or M645C) C-type inactivation were introduced into subunits that were combined with wild-type subunits to form concatenated tetrameric channels with defined subunit composition and stoichiometry. Channels were heterologously expressed in Xenopus oocytes and the two-microelectrode voltage clamp was used to measure the kinetics and steady-state properties of inactivation of whole cell currents. The effect of S631A or T618A mutations on inactivation was a graded function of the number of mutant subunits within a concatenated tetramer as predicted by a sequential model of cooperative subunit interactions, whereas M645C subunits increased the rate of inactivation of concatemers, as predicted for subunits that act independently of one another. For mutations located within the inactivation gate proper (S620T or G628C/S631C), the presence of a single subunit in a concatenated hERG1 tetramer disrupted gating to the same extent as that observed for mutant homotetramers. Together, our findings indicate that the final step of C-type inactivation of hERG1 channels involves a concerted, all-or-none cooperative interaction between all four subunits, and that probing the mechanisms of channel gating with concatenated heterotypic channels should be interpreted with care, as conclusions regarding the nature of subunit interactions may depend on the specific mutation used to probe the gating process.

ICA-105574 interacts with a common binding site to elicit opposite effects on inactivation gating of EAG and ERG potassium channels

Rapid and voltage-dependent inactivation greatly attenuates outward currents in ether-a-go-go-related gene (ERG) K(+) channels. In contrast, inactivation of related ether-a-go-go (EAG) K(+) channels is very slow and minimally reduces outward currents. ICA-105574 (ICA, or 3-nitro-N-[4-phenoxyphenyl]-benzamide) has opposite effects on inactivation of these two channel types. Although ICA greatly attenuates ERG inactivation by shifting its voltage dependence to more positive potentials, it enhances the rate and extent of EAG inactivation without altering its voltage dependence. Here, we investigate whether the inverse functional response to ICA in EAG and ERG channels is related to differences in ICA binding site or to intrinsic mechanisms of inactivation. Molecular modeling coupled with site-directed mutagenesis suggests that ICA binds in a channel-specific orientation to a hydrophobic pocket bounded by the S5/pore helix/S6 of one subunit and S6 of an adjacent subunit. ICA is a mixed agonist of mutant EAG and EAG/ERG chimera channels that inactivate by a combination of slow and fast mechanisms. With the exception of three residues, the specific amino acids that form the putative binding pocket for ICA in ERG are conserved in EAG. Mutations introduced into EAG to replicate the ICA binding site in ERG did not alter the functional response to ICA. Together these findings suggest that ICA binds to the same site in EAG and ERG channels to elicit opposite functional effects. The resultant agonist or antagonist activity is determined solely by channel-specific differences in the mechanisms of inactivation gating.