The selectivity of K+ ion channels: testing the hypotheses

How K⁺ channels are able to conduct certain cations yet not others remains an important but unresolved question. The recent elucidation of the structure of NaK, an ion channel that conducts both Na⁺ and K⁺ ions, offers an opportunity to test the various hypotheses that have been put forward to explain the selectivity of K⁺ ion channels. We test the snug-fit, field-strength, and over-coordination hypotheses by comparing their predictions to the results of classical molecular dynamics simulations of the K⁺ selective channel KcsA and the less selective channel NaK embedded in lipid bilayers. Our results are incompatible with the so-called strong variant of the snug-fit hypothesis but are consistent with the over-coordination hypothesis and neither confirm nor refute the field-strength hypothesis. We also find that the ions and waters in the NaK selectivity filter unexpectedly move to a new conformation in seven K⁺ simulations: the two K⁺ ions rapidly move from site S4 to S2 and from the cavity to S4. At the same time, the selectivity filter narrows around sites S1 and S2 and the carbonyl oxygen atoms rotate 20°−40° inwards toward the ion. These motions diminish the large structural differences between the crystallographic structures of the selectivity filters of NaK and KcsA and appear to allow the binding of ions to S2 of NaK at physiological temperature.     

Ion binding properties and structure stability of the NaK channel

Ion distribution in the selectivity filter and ion-water and ion-protein interactions of NaK channel are systematically investigated by all-atom molecular dynamics simulations, with the tetramer channel protein being embedded in a solvated phospholipid bilayer. Analysis of the simulation results indicates that K⁺ ions prefer to bind within the sites formed by two adjacent planes of oxygen atoms from the selectivity filter, while Na⁺ ions are inclined to bind to a single plane of four oxygen atoms. At the same time, both K⁺ and Na⁺ ions can diffuse in the vestibule, accompanying with movements of the water molecules confined in a complex formed by the vestibule together with four small grottos connecting to it. As a result, K⁺ ions show a wide range of coordination numbers (6-8), while Na⁺ ions display a constant coordination number of ∼ 6 in the selectivity filter, which may result in the loss of selectivity of NaK. It is also found that a Ca²⁺ can bind at the extracellular site as reported in the crystal structure in a partially hydrated state, or at a higher site in a full hydration state. Furthermore, the carbonyl group of Asp66 can reorient to point towards the center pore when an ion exists in the vestibule, while that of Gly65 always aligns tangentially to the channel axis, as in the crystallographic structures.     

Sodium Ions Do Not Stabilize the Selectivity Filter of a Potassium Channel

Ion conduction is an essential function for electrical activity in all organisms. The non-selective ion channel NaK was previously shown to adopt two stable conformations of the selectivity filter. Here, we present solid-state NMR measurements of NaK demonstrating a population shift between these conformations induced by changing the ions in the sample while the overall structure of NaK is not affected. We show that two K⁺-selective mutants (NaK2K and NaK2K-Y66F) suffer a complete loss of selectivity filter stability under Na⁺ conditions, but do not collapse into a defined structure. Widespread chemical shift perturbations are seen between the Na⁺ and K⁺ states of the K⁺-selective mutants in the region of the pore helix indicating structural changes. We conclude that the stronger link between the selectivity filter and the pore helix in the K⁺-selective mutants, compared to the non-selective wild-type NaK channel, reduces the ion-dependent conformational flexibility of the selectivity filter.     

Gating at the selectivity filter of ion channels that conduct Na+ and K+ ions

The NaK channel is a cation selective channel with similar permeability for K⁺ and Na⁺. The available crystallographic structure of wild-type (WT) NaK is usually associated with a conductive state of the channel. Here, potential of mean force for complete conduction events of Na⁺ and K⁺ ions through NaK show that: i), large energy barriers prevent the passage of ions through the WT NaK structure, ii), the barriers are correlated to the presence of a hydrogen bond between Asp-66 and Asn-68, and iii), the structure of NaK mutated to mimic cyclic nucleotide-gated channels conducts Na⁺ and K⁺. These results support the hypothesis that the filter of cation selective channels can adopt at least two different structures: a conductive one, represented by the x-ray structures of the NaK-CNG chimeras, and a closed one, represented by the x-ray structures of the WT NaK.     

Hydration valve controlled non-selective conduction of Na and K in the NaK channel

The Na⁺ and K⁺ channels are essential to neural signaling, but our current knowledge at the atomic level is mainly limited to the conducting mechanism of K⁺. Unlike a K⁺ channel having four equivalent K⁺-binding sites in its selectivity filter, a NaK channel has a vestibule in the middle part of its selectivity filter, and can conduct both Na⁺ and K⁺ ions. However, the underlying mechanism for non-selective ion conduction in NaK remains elusive. Here we find four small grottos connecting with the vestibule of the NaK selectivity filter, which form a vestibule-grotto complex perpendicular to the filter pore with a few water molecules within it. It is shown that two or more of the water molecules coming to the vestibule to coordinate the cation are necessary for conducting both Na⁺ and K⁺ ions, while only one water molecule in the vestibule will obstruct ion permeation. Thus, the complex with the aid of interior water movement forms a dynamic hydration valve which is flexible in conveying different cations through the vestibule. Similar exquisite hydration valve mechanisms are expected to be utilized by other non-selective cation channels, and the results should shed new light on the importance of water in neural signaling.     

Preferential binding of K+ ions in the selectivity filter at equilibrium explains high selectivity of K+ channels

K⁺ channels exhibit strong selectivity for K⁺ ions over Na⁺ ions based on electrophysiology experiments that measure ions competing for passage through the channel. During this conduction process, multiple ions interact within the region of the channel called the selectivity filter. Ion selectivity may arise from an equilibrium preference for K⁺ ions within the selectivity filter or from a kinetic mechanism whereby Na⁺ ions are precluded from entering the selectivity filter. Here, we measure the equilibrium affinity and selectivity of K⁺ and Na⁺ ions binding to two different K⁺ channels, KcsA and MthK, using isothermal titration calorimetry. Both channels exhibit a large preference for K⁺ over Na⁺ ions at equilibrium, in line with electrophysiology recordings of reversal potentials and Ba²⁺ block experiments used to measure the selectivity of the external-most ion-binding sites. These results suggest that the high selectivity observed during ion conduction can originate from a strong equilibrium preference for K⁺ ions in the selectivity filter, and that K⁺ selectivity is an intrinsic property of the filter. We hypothesize that the equilibrium preference for K⁺ ions originates in part through the optimal spacing between sites to accommodate multiple K⁺ ions within the selectivity filter.     

Effects of the Protonation State of the EEEE Motif of a Bacterial Na+-channel on Conduction and Pore Structure

A distinctive feature of prokaryotic Na⁺-channels is the presence of four glutamate residues in their selectivity filter. In this study, how the structure of the selectivity filter, and the free-energy profile of permeating Na⁺ ions are altered by the protonation state of Glu177 are analyzed. It was found that protonation of a single glutamate residue was enough to modify the conformation of the selectivity filter and its conduction properties. Molecular dynamics simulations revealed that Glu177 residues may adopt two conformations, with the side chain directed toward the extracellular entrance of the channel or the intracellular cavity. The likelihood of the inwardly directed arrangement increases when Glu177 residues are protonated. The presence of one glutamate residue with its chain directed toward the intracellular cavity increases the energy barrier for translocation of Na⁺ ions. These higher-energy barriers preclude Na⁺ ions to permeate the selectivity filter of prokaryotic Na⁺-channels when one or more Glu177 residues are protonated.     

Effects of the Protonation State of the EEEE Motif of a Bacterial Na-channel on Conduction and Pore Structure

A distinctive feature of prokaryotic Na⁺-channels is the presence of four glutamate residues in their selectivity filter. In this study, how the structure of the selectivity filter, and the free-energy profile of permeating Na⁺ ions are altered by the protonation state of Glu177 are analyzed. It was found that protonation of a single glutamate residue was enough to modify the conformation of the selectivity filter and its conduction properties. Molecular dynamics simulations revealed that Glu177 residues may adopt two conformations, with the side chain directed toward the extracellular entrance of the channel or the intracellular cavity. The likelihood of the inwardly directed arrangement increases when Glu177 residues are protonated. The presence of one glutamate residue with its chain directed toward the intracellular cavity increases the energy barrier for translocation of Na⁺ ions. These higher-energy barriers preclude Na⁺ ions to permeate the selectivity filter of prokaryotic Na⁺-channels when one or more Glu177 residues are protonated.     

Non-Equilibrium Dynamics Contribute to Ion Selectivity in the KcsA Channel

The ability of biological ion channels to conduct selected ions across cell membranes is critical for the survival of both animal and bacterial cells. Numerous investigations of ion selectivity have been conducted over more than 50 years, yet the mechanisms whereby the channels select certain ions and reject others are not well understood. Here we report a new application of Jarzynski’s Equality to investigate the mechanism of ion selectivity using non-equilibrium molecular dynamics simulations of Na⁺ and K⁺ ions moving through the KcsA channel. The simulations show that the selectivity filter of KcsA adapts and responds to the presence of the ions with structural rearrangements that are different for Na⁺ and K⁺. These structural rearrangements facilitate entry of K⁺ ions into the selectivity filter and permeation through the channel, and rejection of Na⁺ ions. A mechanistic model of ion selectivity by this channel based on the results of the simulations relates the structural rearrangement of the selectivity filter to the differential dehydration of ions and multiple-ion occupancy and describes a mechanism to efficiently select and conduct K⁺. Estimates of the K⁺/Na⁺ selectivity ratio and steady state ion conductance for KcsA from the simulations are in good quantitative agreement with experimental measurements. This model also accurately describes experimental observations of channel block by cytoplasmic Na⁺ ions, the “punch through” relief of channel block by cytoplasmic positive voltages, and is consistent with the knock-on mechanism of ion permeation.     

Potassium and sodium ions in a potassium channel studied by molecular dynamics simulations

We have performed simulations of both a single potassium ion and a single sodium ion within the pore of the bacterial potassium channel KcsA. For both ions there is a dehydration energy barrier at the cytoplasmic mouth suggesting that the crystal structure is a closed conformation of the channel. There is a potential energy barrier for a sodium ion in the selectivity filter that is not seen for potassium. Radial distribution functions for both ions with the carbonyl oxygens of the selectivity filter indicate that sodium may interact more tightly with the filter than does potassium. This suggests that the key to the ion selectivity of KcsA is the greater dehydration energy of Na⁺ ions, and helps to explain the block of KcsA by internal Na⁺ ions.     

Coordination Numbers of K and Na Ions Inside the Selectivity Filter of the KcsA Potassium Channel: Insights from First Principles Molecular Dynamics

Quantum mechanics/molecular mechanics (QM/MM) Car-Parrinello simulations were performed to estimate the coordination numbers of K⁺ and Na⁺ ions in the selectivity filter of the KcsA channel, and in water. At the DFT/BLYP level, K⁺ ions were found to display an average coordination number of 6.6 in the filter, and 6.2 in water. Na⁺ ions displayed an average coordination number of 5.2 in the filter, and 5.0 in water. A comparison was made with the average coordination numbers obtained from using classical molecular dynamics (6.7 for K⁺ in the filter, 6.6 for K⁺ in water, 6.0 for Na⁺ in the filter, and 5.2 for Na⁺ in water). The observation that different coordination numbers were displayed by the ions in QM/MM simulations and in classical molecular dynamics is relevant to the discussion of selectivity in K-channels.     

K+ Block Is the Mechanism of Functional Asymmetry in Bacterial Nav Channels

Crystal structures of several bacterial Na_v channels have been recently published and molecular dynamics simulations of ion permeation through these channels are consistent with many electrophysiological properties of eukaryotic channels. Bacterial Na_v channels have been characterized as functionally asymmetric, and the mechanism of this asymmetry has not been clearly understood. To address this question, we combined non-equilibrium simulation data with two-dimensional equilibrium unperturbed landscapes generated by umbrella sampling and Weighted Histogram Analysis Methods for multiple ions traversing the selectivity filter of bacterial Na_vAb channel. This approach provided new insight into the mechanism of selective ion permeation in bacterial Na_v channels. The non-equilibrium simulations indicate that two or three extracellular K⁺ ions can block the entrance to the selectivity filter of Na_vAb in the presence of applied forces in the inward direction, but not in the outward direction. The block state occurs in an unstable local minimum of the equilibrium unperturbed free-energy landscape of two K⁺ ions that can be ‘locked’ in place by modest applied forces. In contrast to K⁺, three Na⁺ ions move favorably through the selectivity filter together as a unit in a loose “knock-on” mechanism of permeation in both inward and outward directions, and there is no similar local minimum in the two-dimensional free-energy landscape of two Na⁺ ions for a block state. The useful work predicted by the non-equilibrium simulations that is required to break the K⁺ block is equivalent to large applied potentials experimentally measured for two bacterial Na_v channels to induce inward currents of K⁺ ions. These results illustrate how inclusion of non-equilibrium factors in the simulations can provide detailed information about mechanisms of ion selectivity that is missing from mechanisms derived from either crystal structures or equilibrium unperturbed free-energy landscapes.     

Nonselective Conduction in a Mutated NaK Channel with Three Cation-Binding Sites

The NaK channel is a cation-selective protein with similar permeability for K⁺ and Na⁺ ions. Crystallographic structures are available for the wild-type and mutated NaK channels with different numbers of cation-binding sites. We have performed a comparison between the potentials of mean force governing the translocation of K⁺ ions and mixtures of one Na⁺ and three K⁺ ions in a mutated NaK channel with only three cation-binding sites (NaK-CNG). Since NaK-CNG is not selective for K⁺ over Na⁺, analysis of its multi-ion potential energy surfaces can provide clues about how selectivity originates. Comparison of the potentials of mean force of NaK-CNG and K⁺-selective channels yields observations that strongly suggest that the number of contiguous ion binding sites in a single-file mechanism is the key determinant of the channel’s selectivity properties, as already proposed by experimental studies. We conclude that the presence of four binding sites in K⁺-selective channels is essential for highly selective and efficient permeation of K⁺ ions, and that a key difference between K⁺-selective and nonselective channels is the absence/presence of a binding site for Na⁺ ions at the boundary between S2 and S3 in the context of multi-ion permeation events.     

Dynamics of K Ion Conduction through Kv1.2

The crystallographic structure of a potassium channel, Kv1.2, in an open state makes it feasible to simulate entire K⁺ ion permeation events driven by a voltage bias and, thereby, elucidate the mechanism underlying ion conduction and selectivity of this type of channel. This Letter demonstrates that molecular dynamics simulations can provide movies of the overall conduction of K⁺ ions through Kv1.2. As suggested earlier, the conduction is concerted in the selectivity filter, involving 2-3 ions residing mainly at sites identified previously by crystallography and modeling. The simulations reveal, however, the jumps of ions between these sites and identify the sequence of multi-ion configurations involved in permeation.     

K+-Dependent Selectivity and External Ca2+ Block of Shab K+ Channels

Potassium channels allow the selective flux of K⁺ excluding the smaller, and more abundant in the extracellular solution, Na⁺ ions. Here we show that Shab is a typical K⁺ channel that excludes Na⁺ under bi-ionic, Na_o/K_i or Na_o/Rb_i, conditions. However, when internal K⁺ is replaced by Cs⁺ (Na_o/Cs_i), stable inward Na⁺ and outward Cs⁺ currents are observed. These currents show that Shab selectivity is not accounted for by protein structural elements alone, as implicit in the snug-fit model of selectivity. Additionally, here we report the block of Shab channels by external Ca²⁺ ions, and compare the effect that internal K⁺ replacement exerts on both Ca²⁺ and TEA block. Our observations indicate that Ca²⁺ blocks the channels at a site located near the external TEA binding site, and that this pore region changes conformation under conditions that allow Na⁺ permeation. In contrast, the latter ion conditions do not significantly affect the binding of quinidine to the pore central cavity. Based on our observations and the structural information derived from the NaK bacterial channel, we hypothesize that Ca²⁺ is probably coordinated by main chain carbonyls of the pore´s first K⁺-binding site.     

Analysis of the selectivity filter of the voltage-gated sodium channel NavRh

NaChBac is a bacterial voltage-gated sodium (Na_v) channel that shows sequence similarity to voltage-gated calcium channels. To understand the ion-permeation mechanism of Na_v channels, we combined molecular dynamics simulation, structural biology and electrophysiological approaches to investigate the recently determined structure of Na_vRh, a marine bacterial NaChBac ortholog. Two Na⁺ binding sites are identified in the selectivity filter (SF) in our simulations: The extracellular Na⁺ ion first approaches site 1 constituted by the side groups of Ser181 and Glu183, and then spontaneously arrives at the energetically more favorable site 2 formed by the carbonyl oxygens of Leu179 and Thr178. In contrast, Ca²⁺ ions are prone to being trapped by Glu183 at site 1, which then blocks the entrance of both Na⁺ and Ca²⁺ to the vestibule of the SF. In addition, Na⁺ permeates through the selective filter in an asymmetrical manner, a feature that resembles that of the mammalian Na_v orthologs. The study reported here provides insights into the mechanism of ion selectivity on Na⁺ over Ca²⁺ in mammalian Na_v channels.     

Ionic interactions of Ba2+ blockades in the MthK K+ channel

The movement and interaction of multiple ions passing through in single file underlie various fundamental K⁺ channel properties, from the effective conduction of K⁺ ions to channel blockade by Ba²⁺ ions. In this study, we used single-channel electrophysiology and x-ray crystallography to probe the interactions of Ba²⁺ with permeant ions within the ion conduction pathway of the MthK K⁺ channel. We found that, as typical of K⁺ channels, the MthK channel was blocked by Ba²⁺ at the internal side, and the Ba²⁺-blocking effect was enhanced by external K⁺. We also obtained crystal structures of the MthK K⁺ channel pore in both Ba²⁺–Na⁺ and Ba²⁺–K⁺ environments. In the Ba²⁺–Na⁺ environment, we found that a single Ba²⁺ ion remained bound in the selectivity filter, preferably at site 2, whereas in the Ba²⁺–K⁺ environment, Ba²⁺ ions were predominantly distributed between sites 3 and 4. These ionic configurations are remarkably consistent with the functional studies and identify a molecular basis for Ba²⁺ blockade of K⁺ channels.     

Dynamics, Energetics, and Selectivity of the Low-K KcsA Channel Structure

Potassium channels are a diverse family of integral membrane proteins through which K⁺ can pass selectively. There is ongoing debate about the nature of conformational changes associated with the opening/closing and conductive/nonconductive states of potassium channels. The channels partly exert their function by varying their conductance through a mechanism known as C-type inactivation. Shortly after the activation of K⁺ channels, their selectivity filter stops conducting ions at a rate that depends on various stimuli. The molecular mechanism of C-type inactivation has not been fully understood yet. However, the X-ray structure of the KcsA channel obtained in the presence of low K⁺ concentration is thought to be representative of a K⁺ channel in the C-type inactivated state. Here, extensive, fully atomistic molecular dynamics and free-energy simulations of the low-K⁺ KcsA structure in an explicit lipid bilayer are performed to evaluate the stability of this structure and the selectivity of its binding sites. We find that the low-K⁺ KcsA structure is stable on the timescale of the molecular dynamics simulations performed, and that ions preferably remain in S1 and S4. In the absence of ions, the selectivity filter evolves toward an asymmetric architecture, as already observed in other computations of the high-K⁺ structure of KcsA and KirBac. The low-K⁺ KcsA structure is not permeable by Na⁺, K⁺, or Rb⁺, and the selectivity of its binding sites is different from that of the high-K⁺ structure.     

Pore-forming transmembrane domains control ion selectivity and selectivity filter conformation in the KirBac1.1 potassium channel

Potassium (K⁺) channels are membrane proteins with the remarkable ability to very selectively conduct K⁺ ions across the membrane. High-resolution structures have revealed that dehydrated K⁺ ions permeate through the narrowest region of the pore, formed by the backbone carbonyls of the signature selectivity filter (SF) sequence TxGYG. However, the existence of nonselective channels with similar SF sequences, as well as effects of mutations in other regions on selectivity, suggest that the SF is not the sole determinant of selectivity. We changed the selectivity of the KirBac1.1 channel by introducing mutations at residue I131 in transmembrane helix 2 (TM2). These mutations increase Na⁺ flux in the absence of K⁺ and introduce significant proton conductance. Consistent with K⁺ channel crystal structures, single-molecule FRET experiments show that the SF is conformationally constrained and stable in high-K⁺ conditions but undergoes transitions to dilated low-FRET states in high-Na⁺/low-K⁺ conditions. Relative to wild-type channels, I131M mutants exhibit marked shifts in the K⁺ and Na⁺ dependence of SF dynamics to higher K⁺ and lower Na⁺ concentrations. These results illuminate the role of I131, and potentially other structural elements outside the SF, in controlling ion selectivity, by suggesting that the physical interaction of these elements with the SF contributes to the relative stability of the constrained K⁺-induced SF configuration versus nonselective dilated conformations.     

The Mechanism of Na+/K+ Selectivity in Mammalian Voltage-Gated Sodium Channels Based on Molecular Dynamics Simulation

Voltage-gated sodium (Na_v) channels and their Na⁺/K⁺ selectivity are of great importance in the mammalian neuronal signaling. According to mutational analysis, the Na⁺/K⁺ selectivity in mammalian Na_v channels is mainly determined by the Lys and Asp/Glu residues located at the constriction site within the selectivity filter. Despite successful molecular dynamics simulations conducted on the prokaryotic Na_v channels, the lack of Lys at the constriction site of prokaryotic Na_v channels limits how much can be learned about the Na⁺/K⁺ selectivity in mammalian Na_v channels. In this work, we modeled the mammalian Na_v channel by mutating the key residues at the constriction site in a prokaryotic Na_v channel (Na_vRh) to its mammalian counterpart. By simulating the mutant structure, we found that the Na⁺ preference in mammalian Na_v channels is collaboratively achieved by the deselection from Lys and the selection from Asp/Glu within the constriction site.