The evolutionarily conserved residue A653 plays a key role in HERG channel closing

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.

Full opening of helix bundle crossing does not lead to NaK channel activation

A critical part of ion channel function is the ability to open and close in response to stimuli and thus conduct ions in a regulated fashion. While x-ray diffraction studies of ion channels suggested a general steric gating mechanism located at the helix bundle crossing (HBC), recent functional studies on several channels indicate that the helix bundle crossing is wide-open even in functionally nonconductive channels. Two NaK channel variants were crystallized in very different open and closed conformations, which served as important models of the HBC gating hypothesis. However, neither of these NaK variants is conductive in liposomes unless phenylalanine 92 is mutated to alanine (F92A). Here, we use NMR to probe distances at near-atomic resolution of the two NaK variants in lipid bicelles. We demonstrate that in contrast to the crystal structures, both NaK variants are in a fully open conformation, akin to Ca2+-bound MthK channel structure where the HBC is widely open. While we were not able to determine what a conductive NaK structure is like, our further inquiry into the gating mechanism suggests that the selectivity filter and pore helix are coupled to the M2 helix below and undergo changes in the structure when F92 is mutated. Overall, our data show that NaK exhibits coupling between the selectivity filter and HBC, similar to K+ channels, and has a more complex gating mechanism than previously thought, where the full opening of HBC does not lead to channel activation.

The role of site-directed point mutations in protein misfolding

Mutations inducing higher clashing and lower matching residue pairs lead to misfolding.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

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

The hERG K(+) channel is important for establishing normal electrical activity in the human heart. The channel's unique gating response to membrane potential changes indicates specific interactions between voltage sensor and pore domains that are poorly understood. In the absence of a crystal structure we constructed a homology model of the full hERG membrane domain and performed 0.5 μs molecular dynamics (MD) simulations in a hydrated membrane. The simulations identify potential interactions involving residues at the extracellular surface of S1 in the voltage sensor and at the N-terminal end of the pore helix in the hERG model. In addition, a diffuse interface involving hydrophobic residues on S4 (voltage sensor) and pore domain S5 of an adjacent subunit was stable during 0.5 μs of simulation. To assess the ability of the model to give insight into the effects of channel mutation we simulated a hERG mutant that contains a Leu to Pro substitution in the voltage sensor S4 helical segment (hERG L532P). Consistent with the retention of gated K(+) conductance, the L532P mutation was accommodated in the S4 helix with little disruption of helical structure. The mutation reduced the extent of interaction across the S4-S5 interface, suggesting a structural basis for the greatly enhanced deactivation rate in hERG L532P. The study indicates that pairwise comparison of wild-type and mutated channel models is a useful approach to interpreting functional data where uncertainty in model structures exist.

Identification of functionally critical residues in the channel domain of inositol trisphosphate receptors

We have combined alanine mutagenesis and functional assays to identify amino acid residues in the channel domain that are critical for inositol 1,4,5-trisphosphate receptor (IP(3)R) channel function. The residues selected were highly conserved in all three IP(3)R isoforms and were located in the cytosolic end of the S6 pore-lining helix and proximal portion of the C-tail. Two adjacent hydrophobic amino acids (Ile-2588 and Ile-2589) at the putative cytosolic interface of the S6 helix inactivated channel function and could be candidates for the channel gate. Of five negatively charged residues mutated, none completely eliminated channel function. Of five positively charged residues mutated, only one inactivated the channel (Arg-2596). In addition to the previously identified role of a pair of cysteines in the C-tail (Cys-2610 and Cys-2613), a pair of highly conserved histidines (His-2630 and His-2635) were also essential for channel function. Expression of the H2630A and H2635A mutants (but not R2596A) produced receptors with destabilized interactions between the N-terminal fragment and the channel domain. A previously unrecognized association between the cytosolic C-tail and the TM 4,5-loop was demonstrated using GST pulldown assays. However, none of the mutations in the C-tail interfered with this interaction or altered the ability of the C-tail to assemble into dimers. Our present findings and recent information on IP(3)R structure from electron microscopy and crystallography are incorporated into a revised model of channel gating.

Gating of two pore domain potassium channels

Two-pore-domain potassium (K2P) channels are responsible for background leak currents which regulate the membrane potential and excitability of many cell types. Their activity is modulated by a variety of chemical and physical stimuli which act to increase or decrease the open probability of individual K2P channels. Crystallographic data and homology modelling suggest that all K(+) channels possess a highly conserved structure for ion selectivity and gating mechanisms. Like other K(+) channels, K2P channels are thought to have two primary conserved gating mechanisms: an inactivation (or C-type) gate at the selectivity filter close to the extracellular side of the channel and an activation gate at the intracellular entrance to the channel involving key, identified, hinge glycine residues. Zinc and hydrogen ions regulate Drosophila KCNK0 and mammalian TASK channels, respectively, by interacting with the inactivation gate of these channels. In contrast, the voltage dependence of TASK3 channels is mediated through its activation gate. For KCNK0 it has been shown that the gates display positive cooperativity. It is of much interest to determine whether other K2P regulatory compounds interact with either the activation gate or the inactivation gate to alter channel activity or, indeed, whether additional regulatory gating pathways exist.