Mechanism of activation at the selectivity filter of the KcsA K+ channel

Do selectivity filter carbonyls in K+ channels flip away from the pore? Two-dimensional infrared spectroscopy study

Molecular dynamics simulations revealed that the carbonyls of the Val residue in the conserved selectivity filter sequence TVGTG of potassium ion channels can flip away from the pore to form hydrogen bonds with the network of water molecules residing behind the selectivity filter. Such a configuration has been proposed to be relevant for C-type inactivation. Experimentally, X-ray crystallography of the KcsA channel admits the possibility that the Val carbonyls can flip, but it cannot decisively confirm the existence of such a configuration. In this study, we combined molecular dynamics simulations and line shape theory to design two-dimensional infrared spectroscopy experiments that can corroborate the existence of the selectivity filter configuration with flipped Val carbonyls. This ability to distinguish between flipped and unflipped carbonyls is based on the varying strength of the electric field inside and outside the pore, which is directly linked to carbonyl stretching frequencies that can be resolved using infrared spectroscopy.

Ionic Coulomb blockade controls the current in a short narrow carbon nanotube

We use all-atom molecular dynamics simulations to investigate ionic conduction in a short, charged, single-wall carbon nanotube. They reveal ionic Coulomb blockade (ICB) oscillations in the current as a function of the fixed charge on the wall, and an associated occupancy staircase. Current peaks related to fluctuations around the 2 → 1 and 1 → 0 steps in occupancy are clearly resolved, in agreement with ICB theory. Current peaks were also observed at constant occupancy. These unpredicted secondary peaks are attributed to edge effects involving a remote knock-on mechanism; they are attenuated, or absent, for certain choices of model parameters. The key parameters of the system that underlie the current oscillations are estimated using ICB theory and the potential of the mean force. Future perspectives opened up by these observations are discussed.

Insight into the K channel's selectivity from binding of K+, Na+ and water to N-methylacetamide

In potassium channels that conduct K+ selectively over Na+, which sites are occupied by K+ or water and the mechanism of selectivity are unresolved questions. The combination of the energetics and the constraints imposed by the protein structure yield the selective permeation and occupancy. To gain insight into the combination of structure and energetics, we performed density functional theory (DFT) calculations of multiple N-methyl acetamide (NMA) ligands binding to K+ and Na+, relative to hydrated K+ and Na+. NMA is an analogue of the amino acid backbone and provides the carbonyl binding to the ions that occurs in most binding sites of the K+ channel. Unconstrained optimal structures are obtained through geometry optimization calculations of the NMA ligand binding. The complexes formed by 8 NMA binding to the cations have the O atoms positioned in nearly identical locations as the O atoms in the selectivity filter. The transfer free energies between bulk water and K+ or Na+ bound to 8 NMA are almost identical, implying there is no selectivity by a single site. For water optimized with 8 NMA, binding is weak and O atoms are not positioned as in the K+ channel selectivity filter, suggesting that the ions are much more favored than water. Optimal structures of 8 NMA binding with two cations (K+ or Na+) are stable and have lower binding free energy than the optimal structures with just one cation. However, in the Na+ case, the optimal structure deforms and does not match the K+ channel; that is, two bound Na+ are destabilizing. In contrast, the two K+ structure is stabilized and the selectivity free energy favors K+. Overall, this study shows that binding site occupancy and the mechanism for K+ selectivity involves multiple K+ binding in multiple neighboring layers or sites of the K+ channel selectivity filter.

Asymmetric contribution of a selectivity filter gate in triggering inactivation of CaV1.3 channels

Voltage-dependent and Ca2+-dependent inactivation (VDI and CDI, respectively) of CaV channels are two biologically consequential feedback mechanisms that fine-tune Ca2+ entry into neurons and cardiomyocytes. Although known to be initiated by distinct molecular events, how these processes obstruct conduction through the channel pore remains poorly defined. Here, focusing on ultrahighly conserved tryptophan residues in the interdomain interfaces near the selectivity filter of CaV1.3, we demonstrate a critical role for asymmetric conformational changes in mediating VDI and CDI. Specifically, mutagenesis of the domain III-IV interface, but not others, enhanced VDI. Molecular dynamics simulations demonstrate that mutations in distinct selectivity filter interfaces differentially impact conformational flexibility. Furthermore, mutations in distinct domains preferentially disrupt CDI mediated by the N- versus C-lobes of CaM, thus uncovering a scheme of structural bifurcation of CaM signaling. These findings highlight the fundamental importance of the asymmetric arrangement of the pseudotetrameric CaV pore domain for feedback inhibition.

Asymmetric Contribution of a Selectivity Filter Gate in Triggering Inactivation of CaV1.3 Channels

Voltage-dependent and Ca2+-dependent inactivation (VDI and CDI, respectively) of CaV channels are two biologically consequential feedback mechanisms that fine-tune Ca2+ entry into neurons and cardiomyocytes. Although known to be initiated by distinct molecular events, how these processes obstruct conduction through the channel pore remains poorly defined. Here, focusing on ultra-highly conserved tryptophan residues in the inter-domain interfaces near the selectivity filter of CaV1.3, we demonstrate a critical role for asymmetric conformational changes in mediating VDI and CDI. Specifically, mutagenesis of the domain III-IV interface, but not others, enhanced VDI. Molecular dynamics simulations demonstrate that mutations in distinct selectivity filter interfaces differentially impact conformational flexibility. Furthermore, mutations in distinct domains preferentially disrupt CDI mediated by the N- versus C-lobes of CaM, thus uncovering a scheme of structural bifurcation of CaM signaling. These findings highlight the fundamental importance of the asymmetric arrangement of the pseudo-tetrameric CaV pore domain for feedback inhibition.

Probing Ion Configurations in the KcsA Selectivity Filter with Single-Isotope Labels and 2D IR Spectroscopy

The potassium ion (K+) configurations of the selectivity filter of the KcsA ion channel protein are investigated with two-dimensional infrared (2D IR) spectroscopy of amide I vibrations. Single 13C-18O isotope labels are used, for the first time, to selectively probe the S1/S2 or S2/S3 binding sites in the selectivity filter. These binding sites have the largest differences in ion occupancy in two competing K+ transport mechanisms: soft-knock and hard-knock. According to the former, water molecules alternate between K+ ions in the selectivity filter while the latter assumes that K+ ions occupy the adjacent sites. Molecular dynamics simulations and computational spectroscopy are employed to interpret experimental 2D IR spectra. We find that in the closed conductive state of the KcsA channel, K+ ions do not occupy adjacent binding sites. The experimental data is consistent with simulated 2D IR spectra of soft-knock ion configurations. In contrast, the simulated spectra for the hard-knock ion configurations do not reproduce the experimental results. 2D IR spectra of the hard-knock mechanism have lower frequencies, homogeneous 2D lineshapes, and multiple peaks. In contrast, ion configurations of the soft-knock model produce 2D IR spectra with a single peak at a higher frequency and inhomogeneous lineshape. We conclude that under equilibrium conditions, in the absence of transmembrane voltage, both water and K+ ions occupy the selectivity filter of the KcsA channel in the closed conductive state. The ion configuration is central to the mechanism of ion transport through potassium channels.

Non-trivial dynamics in a model of glial membrane voltage driven by open potassium pores

Despite the molecular evidence that a nearly linear steady-state current-voltage relationship in mammalian astrocytes reflects a total current resulting from more than one differentially regulated K+ conductance, detailed ordinary differential equation (ODE) models of membrane voltage Vm are still lacking. Various experimental results reporting altered rectification of the major Kir currents in glia, dominated by Kir4.1, have motivated us to develop a detailed model of Vm dynamics incorporating the weaker potassium K2P-TREK1 current in addition to Kir4.1, and study the stability of the resting state Vr. The main question is whether, with the loss of monotonicity in glial I-V curve resulting from altered Kir rectification, the nominal resting state Vr remains stable, and the cell retains the trivial, potassium electrode behavior with Vm after EK. The minimal two-dimensional model of Vm near Vr showed that an N-shape deformed Kir I-V curve induces multistability of Vm in a model that incorporates K2P activation kinetics, and nonspecific K+ leak currents. More specifically, an asymmetrical, nonlinear decrease of outward Kir4.1 conductance, turning the channels into inward rectifiers, introduces instability of Vr. That happens through a robust bifurcation giving birth to a second, more depolarized stable resting state Vdr > -10 mV. Realistic recordings from electrographic seizures were used to perturb the model. Simulations of the model perturbed by constant current through gap junctions and seizure-like discharges as local field potentials led to depolarization and switching of Vm between the two stable states, in a downstate-upstate manner. In the event of prolonged depolarizations near Vdr, such catastrophic instability would affect all aspects of the glial function, from metabolic support to membrane transport, and practically all neuromodulatory roles assigned to glia.

Two gates mediate NMDA receptor activity and are under subunit-specific regulation

Kinetics of NMDA receptor (NMDAR) ion channel opening and closing contribute to their unique role in synaptic signaling. Agonist binding generates free energy to open a canonical gate at the M3 helix bundle crossing. Single channel activity is characterized by clusters, or periods of rapid opening and closing, that are separated by long silent periods. A conserved glycine in the outer most transmembrane helices, the M4 helices, regulates NMDAR function. Here we find that the GluN1 glycine mainly regulates single channel events within a cluster, whereas the GluN2 glycine mainly regulates entry and exit from clusters. Molecular dynamics simulations suggest that, whereas the GluN2 M4 (along with GluN2 pre-M1) regulates the gate at the M3 helix bundle crossing, the GluN1 glycine regulates a 'gate' at the M2 loop. Subsequent functional experiments support this interpretation. Thus, the distinct kinetics of NMDARs are mediated by two gates that are under subunit-specific regulation.

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.

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.

Towards the Idea of Molecular Brains

How can single cells without nervous systems perform complex behaviours such as habituation, associative learning and decision making, which are considered the hallmark of animals with a brain? Are there molecular systems that underlie cognitive properties equivalent to those of the brain? This review follows the development of the idea of molecular brains from Darwin’s “root brain hypothesis”, through bacterial chemotaxis, to the recent discovery of neuron-like r-protein networks in the ribosome. By combining a structural biology view with a Bayesian brain approach, this review explores the evolutionary labyrinth of information processing systems across scales. Ribosomal protein networks open a window into what were probably the earliest signalling systems to emerge before the radiation of the three kingdoms. While ribosomal networks are characterised by long-lasting interactions between their protein nodes, cell signalling networks are essentially based on transient interactions. As a corollary, while signals propagated in persistent networks may be ephemeral, networks whose interactions are transient constrain signals diffusing into the cytoplasm to be durable in time, such as post-translational modifications of proteins or second messenger synthesis. The duration and nature of the signals, in turn, implies different mechanisms for the integration of multiple signals and decision making. Evolution then reinvented networks with persistent interactions with the development of nervous systems in metazoans. Ribosomal protein networks and simple nervous systems display architectural and functional analogies whose comparison could suggest scale invariance in information processing. At the molecular level, the significant complexification of eukaryotic ribosomal protein networks is associated with a burst in the acquisition of new conserved aromatic amino acids. Knowing that aromatic residues play a critical role in allosteric receptors and channels, this observation suggests a general role of π systems and their interactions with charged amino acids in multiple signal integration and information processing. We think that these findings may provide the molecular basis for designing future computers with organic processors.

Ion-dependent structure, dynamics, and allosteric coupling in a non-selective cation channel

The selectivity filter (SF) determines which ions are efficiently conducted through ion channel pores. NaK is a non-selective cation channel that conducts Na⁺ and K⁺ with equal efficiency. Crystal structures of NaK suggested a rigid SF structure, but later solid-state NMR and MD simulations questioned this interpretation. Here, we use solution NMR to characterize how bound Na⁺ vs. K⁺ affects NaK SF structure and dynamics. We find that the extracellular end of the SF is flexible on the ps-ns timescale regardless of bound ion. On a slower timescale, we observe a structural change between the Na⁺ and K⁺-bound states, accompanied by increased structural heterogeneity in Na⁺. We also show direct evidence that the SF structure is communicated to the pore via I88 on the M2 helix. These results support a dynamic SF with multiple conformations involved in non-selective conduction. Our data also demonstrate allosteric coupling between the SF and pore-lining helices in a non-selective cation channel that is analogous to the allosteric coupling previously demonstrated for K⁺-selective channels, supporting the generality of this model.

NaK is a bacterial non-selective cation channel. Here, the authors use solution NMR to show that selectivity filter (SF) in NaK is dynamic, with structural differences between the Na+ and K + -bound states. The conformation of the SF is communicated to the pore-lining helices similarly as in the K + -selective channels.

Thermodynamics of ion binding and occupancy in potassium channels

Potassium channels modulate various cellular functions through efficient and selective conduction of K⁺ ions. The mechanism of ion conduction in potassium channels has recently emerged as a topic of debate. Crystal structures of potassium channels show four K⁺ ions bound to adjacent binding sites in the selectivity filter, while chemical intuition and molecular modeling suggest that the direct ion contacts are unstable. Molecular dynamics (MD) simulations have been instrumental in the study of conduction and gating mechanisms of ion channels. Based on MD simulations, two hypotheses have been proposed, in which the four-ion configuration is an artifact due to either averaged structures or low temperature in crystallographic experiments. The two hypotheses have been supported or challenged by different experiments. Here, MD simulations with polarizable force fields validated by ab initio calculations were used to investigate the ion binding thermodynamics. Contrary to previous beliefs, the four-ion configuration was predicted to be thermodynamically stable after accounting for the complex electrostatic interactions and dielectric screening. Polarization plays a critical role in the thermodynamic stabilities. As a result, the ion conduction likely operates through a simple single-vacancy and water-free mechanism. The simulations explained crystal structures, ion binding experiments and recent controversial mutagenesis experiments. This work provides a clear view of the mechanism underlying the efficient ion conduction and demonstrates the importance of polarization in ion channel simulations.

Polarization shapes the energy landscape of ion conduction in potassium channels.

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.

Shaker-IR K+ channel gating in heavy water: Role of structural water molecules in inactivation

It has been reported earlier that the slow (C-type) inactivated conformation in K_v channels is stabilized by a multipoint hydrogen-bond network behind the selectivity filter. Furthermore, MD simulations revealed that structural water molecules are also involved in the formation of this network locking the selectivity filter in its inactive conformation. We found that the application of an extracellular, but not intracellular, solution based on heavy water (D₂O) dramatically slowed entry into the slow inactivated state in Shaker-IR mutants (T449A, T449A/I470A, and T449K/I470C, displaying a wide range of inactivation kinetics), consistent with the proposed effect of the dynamics of structural water molecules on the conformational stability of the selectivity filter. Alternative hypotheses capable of explaining the observed effects of D₂O were examined. Increased viscosity of the external solution mimicked by the addition of glycerol had a negligible effect on the rate of inactivation. In addition, the inactivation time constants of K⁺ currents in the outward and the inward directions in asymmetric solutions were not affected by a H₂O/D₂O exchange, negating an indirect effect of D₂O on the rate of K⁺ rehydration. The elimination of the nonspecific effects of D₂O on our macroscopic current measurements supports the hypothesis that the rate of structural water exchange at the region behind the selectivity filter determines the rate of slow inactivation, as proposed by molecular modeling.

The Persistent Question of Potassium Channel Permeation Mechanisms

Potassium channels play critical roles in many physiological processes, providing a selective permeation route for K⁺ ions in and out of a cell, by employing a carefully designed selectivity filter, evolutionarily conserved from viruses to mammals. The structure of the selectivity filter was determined at atomic resolution by x-ray crystallography, showing a tight coordination of desolvated K⁺ ions by the channel. However, the molecular mechanism of K⁺ ions permeation through potassium channels remains unclear, with structural, functional and computational studies often providing conflicting data and interpretations. In this review, we will present the proposed mechanisms, discuss their origins, and will critically assess them against all available data. General properties shared by all potassium channels are introduced first, followed by the introduction of two main mechanisms of ion permeation: soft and direct knock-on. Then, we will discuss critical computational and experimental studies that shaped the field. We will especially focus on molecular dynamics (MD) simulations, that provided mechanistic and energetic aspects of K⁺ permeation, but at the same time created long-standing controversies. Further challenges and possible solutions are presented as well.

Unraveling of a Strongly Correlated Dynamical Network of Residues Controlling the Permeation of Potassium in KcsA Ion Channel

The complicated patterns of the single-channel currents in potassium ion channel KcsA are governed by the structural variability of the selectivity filter. A comparative analysis of the dynamics of the wild type KcsA channel and several of its mutants showing different conducting patterns was performed. A strongly correlated dynamical network of interacting residues is found to play a key role in regulating the state of the wild type channel. The network is centered on the aspartate D80 which plays the role of a hub by strong interacting via hydrogen bonds with residues E71, R64, R89, and W67. Residue D80 also affects the selectivity filter via its backbones. This network further compromises ions and water molecules located inside the channel that results in the mutual influence: the permeation depends on the configuration of residues in the network, and the dynamics of network's residues depends on locations of ions and water molecules inside the selectivity filter. Some features of the network provide a further understanding of experimental results describing the KcsA activity. In particular, the necessity of anionic lipids to be present for functioning the channel is explained by the interaction between the lipids and the arginine residues R64 and R89 that prevents destabilizing the structure of the selectivity filter.

Conduction and Gating Properties of the TRAAK Channel from Molecular Dynamics Simulations with Different Force Fields

In recent years, the K2P family of potassium channels has been the subject of intense research activity. Owing to the complex function and regulation of this family of ion channels, it is common practice to complement experimental findings with the atomistic description provided by computational approaches such as molecular dynamics (MD) simulations, especially, in light of the unprecedented timescales accessible at present. However, despite recent substantial improvements, the accuracy of MD simulations is still undermined by the intrinsic limitations of force fields. Here, we systematically assessed the performance of the most popular force fields employed to study ion channels at timescales that are orders of magnitude greater than the ones accessible when these energy functions were first developed. Using 32 μs of trajectories, we investigated the dynamics of a member of the K2P ion channel family, the TRAAK channel, using two established force fields in simulations of biological systems: AMBER and CHARMM. We found that while results are comparable on the nanosecond timescales, significant inconsistencies arise at microsecond timescales.

Selectivity filter ion binding affinity determines inactivation in a potassium channel

Potassium channels can become nonconducting via inactivation at a gate inside the highly conserved selectivity filter (SF) region near the extracellular side of the membrane. In certain ligand-gated channels, such as BK channels and MthK, a Ca²⁺-activated K⁺ channel from Methanobacterium thermoautotrophicum, the SF has been proposed to play a role in opening and closing rather than inactivation, although the underlying conformational changes are unknown. Using X-ray crystallography, identical conductive MthK structures were obtained in wide-ranging K⁺ concentrations (6 to 150 mM), unlike KcsA, whose SF collapses at low permeant ion concentrations. Surprisingly, three of the SF's four binding sites remained almost fully occupied throughout this range, indicating high affinities (likely submillimolar), while only the central S2 site titrated, losing its ion at 6 mM, indicating low K⁺ affinity (∼50 mM). Molecular simulations showed that the MthK SF can also collapse in the absence of K⁺, similar to KcsA, but that even a single K⁺ binding at any of the SF sites, except S4, can rescue the conductive state. The uneven titration across binding sites differs from KcsA, where SF sites display a uniform decrease in occupancy with K⁺ concentration, in the low millimolar range, leading to SF collapse. We found that ions were disfavored in MthK's S2 site due to weaker coordination by carbonyl groups, arising from different interactions with the pore helix and water behind the SF. We conclude that these differences in interactions endow the seemingly identical SFs of KcsA and MthK with strikingly different inactivating phenotypes.

The Voltage-Dependent Deactivation of the KvAP Channel Involves the Breakage of Its S4 Helix

Voltage-gated potassium channels (Kv) allow ion permeation upon changes of the membrane electrostatic potential (Vm). Each subunit of these tetrameric channels is composed of six transmembrane helices, of which the anti-parallel helix bundle S1-S4 constitutes the voltage-sensor domain (VSD) and S5-S6 forms the pore domain. Here, using 82 molecular dynamics (MD) simulations involving 266 replicated VSDs, we report novel responses of the archaebacterial potassium channel KvAP to membrane polarization. We show that the S4 α-helix, which is straight in the experimental crystal structure solved under depolarized conditions (Vm ∼ 0), breaks into two segments when the cell membrane is hyperpolarized (Vm << 0), and reversibly forms a single straight helix following depolarization (Vm = 0). The outermost segment of S4 translates along the normal to the membrane, bringing new perspective to previously paradoxical accessibility experiments that were initially thought to imply the displacement of the whole VSD across the membrane. The novel model is applied through steered and unbiased MD simulations to the recently solved whole structure of KvAP. The simulations show that the resting state involves a re-orientation of the S5 α-helix by ∼ 5-6 degrees in respect to the normal of the bilayer, which could result in the constriction and closure of the selectivity filter. Our findings support the idea that the breakage of S4 under (hyper)polarization is a general feature of Kv channels with a non-swapped topology.

TOK channels use the two gates in classical K+ channels to achieve outward rectification

TOKs are outwardly rectifying K⁺ channels in fungi with two pore-loops and eight transmembrane spans. Here, we describe the TOKs from four pathogens that cause the majority of life-threatening fungal infections in humans. These TOKs pass large currents only in the outward direction like the canonical isolate from Saccharomyces cerevisiae (ScTOK), and distinct from other K⁺ channels. ScTOK, AfTOK1 (Aspergillus fumigatus), and H99TOK (Cryptococcus neoformans grubii) are K⁺ -selective and pass current above the K⁺ reversal potential. CaTOK (Candida albicans) and CnTOK (Cryptococcus neoformans neoformans) pass both K⁺ and Na⁺ and conduct above a reversal potential reflecting the mixed permeability of their selectivity filter. Mutations in CaTOK and ScTOK at sites homologous to those that open the internal gates in classical K⁺ channels are shown to produce inward TOK currents. A favored model for outward rectification is proposed whereby the reversal potential determines ion occupancy, and thus, conductivity, of the selectivity filter gate that is coupled to an imperfectly restrictive internal gate, permitting the filter to sample ion concentrations on both sides of the membrane.

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.

Molecular mechanism of a potassium channel gating through activation gate-selectivity filter coupling

Potassium channels are presumed to have two allosterically coupled gates, the activation gate and the selectivity filter gate, that control channel opening, closing, and inactivation. However, the molecular mechanism of how these gates regulate K+ ion flow through the channel remains poorly understood. An activation process, occurring at the selectivity filter, has been recently proposed for several potassium channels. Here, we use X-ray crystallography and extensive molecular dynamics simulations, to study ion permeation through a potassium channel MthK, for various opening levels of both gates. We find that the channel conductance is controlled at the selectivity filter, whose conformation depends on the activation gate. The crosstalk between the gates is mediated through a collective motion of channel helices, involving hydrophobic contacts between an isoleucine and a conserved threonine in the selectivity filter. We propose a gating model of selectivity filter-activated potassium channels, including pharmacologically relevant two-pore domain (K2P) and big potassium (BK) channels.

Changing perspectives on how the permeation pathway through potassium channels is regulated

The primary means by which ion permeation through potassium channels is controlled, and the key to selective intervention in a range of pathophysiological conditions, is the process by which channels switch between non-conducting and conducting states. Conventionally, this has been explained by a steric mechanism in which the pore alternates between two conformations: a 'closed' state in which the conduction pathway is occluded and an 'open' state in which the pathway is sufficiently wide to accommodate fully hydrated ions. Recently, however, 'non-canonical' mechanisms have been proposed for some classes of K⁺ channels. The purpose of this review is to illuminate structural and dynamic relationships underpinning permeation control in K⁺ channels, indicating where additional data might resolve some of the remaining issues.

Conduction through a narrow inward-rectifier K+ channel pore

G-protein–gated inwardly rectifying potassium channels are important mediators of inhibitory neurotransmission. Based on microsecond-scale molecular dynamics simulations, Bernsteiner et al. propose novel gating details that may enable K⁺ flux via a direct knock-on mechanism.

Inwardly rectifying potassium (Kir) channels play a key role in controlling membrane potentials in excitable and unexcitable cells, thereby regulating a plethora of physiological processes. G-protein–gated Kir channels control heart rate and neuronal excitability via small hyperpolarizing outward K⁺ currents near the resting membrane potential. Despite recent breakthroughs in x-ray crystallography and cryo-EM, the gating and conduction mechanisms of these channels are poorly understood. MD simulations have provided unprecedented details concerning the gating and conduction mechanisms of voltage-gated K⁺ and Na⁺ channels. Here, we use multi-microsecond–timescale MD simulations based on the crystal structures of GIRK2 (Kir3.2) bound to phosphatidylinositol-4,5-bisphosphate to provide detailed insights into the channel’s gating dynamics, including insights into the behavior of the G-loop gate. The simulations also elucidate the elementary steps that underlie the movement of K⁺ ions through an inward-rectifier K⁺ channel under an applied electric field. Our simulations suggest that K⁺ permeation might occur via direct knock-on, similar to the mechanism recently shown for K_v channels.

Inactivation in the potassium channel KcsA

•

C-type inactivation in potassium channels is a nearly universal regulatory mechanism.

•

A major hypothesis states that C-type inactivation involves ion loss at the selectivity filter as an allosteric response to activation.

•

NMR is used to probe protein conformational changes in response to pH and [K⁺], demonstrating that H⁺ and K⁺ binding are allosterically coupled in KcsA.

•

The lipids are integrated parts of potassium channels in terms of structure, energetics and function.

Inactivation, the slow cessation of transmission after activation, is a general feature of potassium channels. It is essential for their function, and malfunctions in inactivation leads to numerous pathologies. The detailed mechanism for the C-type inactivation, distinct from the N-type inactivation, remains an active area of investigation. Crystallography, computational simulations, and NMR have greatly enriched our understanding of the process. Here we review the major hypotheses regarding C-type inactivation, particularly focusing on the key role played by NMR studies of the prokaryotic potassium channel KcsA, which serves as a good model for voltage gated mammalian channels.

PACKMOL-Memgen: A Simple-To-Use, Generalized Workflow for Membrane-Protein-Lipid-Bilayer System Building

We present PACKMOL-Memgen, a simple-to-use, generalized workflow for automated building of membrane-protein-lipid-bilayer systems based on open-source tools including Packmol, memembed, pdbremix, and AmberTools. Compared with web-interface-based related tools, PACKMOL-Memgen allows setup of multiple configurations of a system in a user-friendly and efficient manner within minutes. The generated systems are well-packed and thus well-suited as starting configurations in MD simulations under periodic boundary conditions, requiring only moderate equilibration times. PACKMOL-Memgen is distributed with AmberTools and runs on most computing platforms, and its output can also be used for CHARMM or adapted to other molecular-simulation packages.

Investigation of a KcsA Cytoplasmic pH Gate in Lipoprotein Nanodiscs

The bacterial potassium channel KcsA is gated by pH, opening for conduction under acidic conditions. Molecular determinants responsible for this effect have been identified at the extracellular selectivity filter, at the membrane-cytoplasm interface (TM2 gate), and in the cytoplasmic C-terminal domain (CTD), an amphiphilic four-helix bundle mediated by hydrophobic and electrostatic interactions. Here we have employed NMR and EPR to provide a structural view of the pH-induced open-to-closed CTD transition. KcsA was embedded in lipoprotein nanodiscs (LPNs), selectively methyl-protonated at Leu/Val residues to allow observation of both states by NMR, and spin-labeled for the purposes of EPR studies. We observed a pHinduced structural change between an associated structured CTD at neutral pH and a dissociated flexible CTD at acidic pH, with a transition in the 5.0-5.5 range, consistent with a stabilization of the CTD by channel architecture. A double mutant constitutively open at the TM2 gate exhibited reduced stability of associated CTD, as indicated by weaker spin-spin interactions, a shift to higher transition pH values, and a tenfold reduction in the population of the associated "closed" channels. We extended these findings for isolated CTD-derived peptides to full-length KcsA and have established a contribution of the CTD to KcsA pH-controlled gating, which exhibits a strong correlation with the state of the proximal TM2 gate.

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.

Opening leads to closing: Allosteric crosstalk between the activation and inactivation gates in KcsA

Delemotte appraises new computational work revealing that the intracellular activation gate must open for C-type inactivation to occur in K+ channels.

Rapid constriction of the selectivity filter underlies C-type inactivation in the KcsA potassium channel

C-type inactivation in K⁺ channels is thought to be a result of constriction of the selectivity filter. By using MD simulations, Li et al. show that rapid constriction occurs within 1–2 s when the intracellular activation gate is fully open, but not when the gate is closed or partially open.

C-type inactivation is a time-dependent process observed in many K⁺ channels whereby prolonged activation by an external stimulus leads to a reduction in ionic conduction. While C-type inactivation is thought to be a result of a constriction of the selectivity filter, the local dynamics of the process remain elusive. Here, we use molecular dynamics (MD) simulations of the KcsA channel to elucidate the nature of kinetically delayed activation/inactivation gating coupling. Microsecond-scale MD simulations based on the truncated form of the KcsA channel (C-terminal domain deleted) provide a first glimpse of the onset of C-type inactivation. We observe over multiple trajectories that the selectivity filter consistently undergoes a spontaneous and rapid (within 1–2 µs) transition to a constricted conformation when the intracellular activation gate is fully open, but remains in the conductive conformation when the activation gate is closed or partially open. Multidimensional umbrella sampling potential of mean force calculations and nonequilibrium voltage-driven simulations further confirm these observations. Electrophysiological measurements show that the truncated form of the KcsA channel inactivates faster and greater than full-length KcsA, which is consistent with truncated KcsA opening to a greater degree because of the absence of the C-terminal domain restraint. Together, these results imply that the observed kinetics underlying activation/inactivation gating reflect a rapid conductive-to-constricted transition of the selectivity filter that is allosterically controlled by the slow opening of the intracellular gate.

Ligand discrimination and gating in cyclic nucleotide-gated ion channels from apo and partial agonist-bound cryo-EM structures

Cyclic nucleotide-modulated channels have important roles in visual signal transduction and pacemaking. Binding of cyclic nucleotides (cAMP/cGMP) elicits diverse functional responses in different channels within the family despite their high sequence and structure homology. The molecular mechanisms responsible for ligand discrimination and gating are unknown due to lack of correspondence between structural information and functional states. Using single particle cryo-electron microscopy and single-channel recording, we assigned functional states to high-resolution structures of SthK, a prokaryotic cyclic nucleotide-gated channel. The structures for apo, cAMP-bound, and cGMP-bound SthK in lipid nanodiscs, correspond to no, moderate, and low single-channel activity, respectively, consistent with the observation that all structures are in resting, closed states. The similarity between apo and ligand-bound structures indicates that ligand-binding domains are strongly coupled to pore and SthK gates in an allosteric, concerted fashion. The different orientations of cAMP and cGMP in the ‘resting’ and ‘activated’ structures suggest a mechanism for ligand discrimination.

Ion channels are essential for transmitting signals in the nervous system and brain. One large group of ion channels includes members that are activated by cyclic nucleotides, small molecules used to transmit signals within cells. These cyclic nucleotide-gated channels play an important role in regulating our ability to see and smell.

The activity of these ion channels has been studied for years, but scientists have only recently been able to look into their structure. Since structural biology methods require purified, well-behaved proteins, the members of this ion channel family selected for structural studies do not necessarily match those whose activity has been well established. There is a need for a good model that would allow both the structure and activity of a cyclic nucleotide-gated ion channel to be characterized.

The cyclic nucleotide-gated ion channel, SthK, from bacteria called Spirochaeta thermophila, was identified as such model because both its activity and its structure are accessible. Rheinberger et al. have used cryo electron microscopy to solve several high-resolution structures of SthK channels. In two of the structures, SthK was bound to either one of two types of activating cyclic nucleotides – cAMP or cGMP – and in another structure, no cyclic nucleotides were bound. Separately recording the activity of individual channels allowed the activity states likely to be represented by these structures to be identified.

Combining the results of the experiments revealed no activity from channels in an unbound state, low levels of activity for channels bound to cGMP, and moderate activity for channels bound to cAMP. Rheinberger et al. show that the channel, under the conditions experienced in cryo electron microscopy, is closed in all of the states studied. Unexpectedly, the binding of cyclic nucleotides produced no structural change even in the cyclic nucleotide-binding pocket of the channel, a region that was previously observed to undergo such changes when this region alone was crystallized. Rheinberger et al. deduce from this that the four subunits that make up the channel likely undergo the conformational change towards an open state all at once, rather than one by one.

The structures and the basic functional characterization of SthK channels provide a strong starting point for future research into determining the entire opening and closing cycle for a cyclic nucleotide-gated channel. Human equivalents of the channel are likely to work in similar ways. The results presented by Rheinberger et al. could therefore be built upon to help address diseases that result from deficiencies in cyclic nucleotide-gated channels, such as loss of vision due to retinal degradation (retinitis pigmentosa or progressive cone dystrophy) and achromatopsia.

Inverted allosteric coupling between activation and inactivation gates in K+ channels

The selectivity filter and the activation gate in potassium channels are functionally and structurally coupled. An allosteric coupling underlies C-type inactivation coupled to activation gating in this ion-channel family (i.e., opening of the activation gate triggers the collapse of the channel's selectivity filter). We have identified the second Threonine residue within the TTVGYGD signature sequence of K⁺ channels as a crucial residue for this allosteric communication. A Threonine to Alanine substitution at this position was studied in three representative members of the K⁺-channel family. Interestingly, all of the mutant channels exhibited lack of C-type inactivation gating and an inversion of their allosteric coupling (i.e., closing of the activation gate collapses the channel's selectivity filter). A state-dependent crystallographic study of KcsA-T75A proves that, on activation, the selectivity filter transitions from a nonconductive and deep C-type inactivated conformation to a conductive one. Finally, we provide a crystallographic demonstration that closed-state inactivation can be achieved by the structural collapse of the channel's selectivity filter.