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.     

K Conduction in the Selectivity Filter of Potassium Channels Is Monitored by the Charge Distribution along Their Sequence

Potassium channels display a high conservation of sequence of the selectivity filter (SF), yet nature has designed a variety of channels that present a wide range of absolute rates of K⁺ permeation. In KcsA, the structural archetype for K channels, under physiological concentrations, two K⁺ ions reside in the SF in configurations 1,3 (up state) and 2,4 (down state) and ion conduction is believed to follow a throughput cycle involving a transition between these states. Using free-energy calculations of KcsA, Kv1.2, and mutant channels, we show that this transition is characterized by a channel-dependent energy barrier. This barrier is strongly influenced by the charges partitioned along the sequence of each channel. These results unveil therefore how, for similar structures of the SF, the rate of K⁺ turnover may be fine-tuned within the family of potassium channels.     

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.     

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.     

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.     

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.     

Interaction of local anesthetics with the K+ channel pore domain

Local anesthetics and related drugs block ionic currents of Na⁺, K⁺ and Ca²⁺ conducted across the cell membrane by voltage-dependent ion channels. Many of these drugs bind in the permeation pathway, occlude the pore and stop ion movement. However channel-blocking drugs have also been associated with decreased membrane stability of certain tetrameric K⁺ channels, similar to the destabilization of channel function observed at low extracellular K⁺ concentration. Such drug-dependent stability may result from electrostatic repulsion of K⁺ from the selectivity filter by a cationic drug molecule bound in the central cavity of the channel. In this study we used the pore domain of the KcsA K⁺ channel protein to test this hypothesis experimentally with a biochemical assay of tetramer stability and theoretically by computational simulation of local anesthetic docking to the central cavity. We find that two common local anesthetics, lidocaine and tetracaine, promote thermal dissociation of the KcsA tetramer in a K⁺-dependent fashion. Docking simulations of these drugs with open, open-inactivated and closed crystal structures of KcsA yield many energetically favorable drug-channel complexes characterized by nonbonded attraction to pore-lining residues and electrostatic repulsion of K⁺. The results suggest that binding of cationic drugs to the inner cavity can reduce tetramer stability of K⁺ channels.     

Ion–protein interactions of a potassium ion channel studied by attenuated total reflection Fourier transform infrared spectroscopy

An understanding of ion–protein interactions is key to a better understanding of the molecular mechanisms of proteins, such as enzymes, ion channels, and ion pumps. A potassium ion channel, KcsA, has been extensively studied in terms of ion selectivity. Alkali metal cations in the selectivity filter were visualized by X-ray crystallography. Infrared spectroscopy has an intrinsically higher structural sensitivity due to frequency changes in molecular vibrations interacting with different ions. In this review article, I attempt to summarize ion-exchange-induced differences in Fourier transform infrared spectroscopy, as applied to KcsA, to explain how this method can be utilized to study ion–protein interactions in the KcsA selectivity filter. A band at 1680 cm^?1 in the amide I region would be a marker band for the ion occupancy of K⁺, Rb⁺, and Cs⁺ in the filter. The band at 1627 cm^?1 observed in both Na⁺ and Li⁺ conditions suggests that the selectivity filter similarly interacts with these ions. In addition to the structural information, the results show that the titration of K⁺ ions provides quantitative information on the ion affinity of the selectivity filter.     

Sodium permeability and sensitivity induced by mutations in the selectivity filter of the KcsA channel towards Kir channels

The bacterial potassium (K⁺) channel KcsA provides an attractive model system to study ion permeation behavior in a selective K⁺-channel. We changed residue at the N-terminal end of the selectivity filter of KcsA (T74V) to its counterpart in inwardly rectifying K⁺-channels (Kir). The tetramer was found to be stable as unmodified KcsA. Under symmetrical and asymmetrical conditions, Na⁺ increased the inward current in the virtual absence of K⁺ however outward currents were nearly abolished which could be recovered upon internal K⁺ addition. Na⁺ also drastically increased the channel open time either in the presence or virtual absence of K⁺. Furthermore, the T74V mutation decreased the internal Ba²⁺ affinity of the channel possibly by binding to a K⁺ site in the pore. In additional experiments, another point mutation V76I in T74V mutant was carried out thus the selectivity filter resembled more the selectivity filter of Kir channels. The mutant tetramer was converted into monomers as determined by conventional gel electrophoresis. However, native like gel electrophoresis, Trp fluorescence and acrylamide quenching experiments indicated that this mutant still formed a tetramer and apparently adopted similar folding properties as unmodified KcsA. Single-channel experiments further demonstrated that the channel was selective for K⁺ over Na⁺ as Na⁺ blocked channel currents. These data suggest that single point mutation T74V alters the selectivity filter and allows simultaneous occupancy and conduction of K⁺ and Na⁺ probably via ion–ion interaction in the pore. In contrast, both mutations (T74V and V76I) in the same molecule seem to reorganize the pore conformation which controls the overall stability of a selective K⁺-channel.     

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.     

Energy distribution and ion selectivity of the bacterial potassium channel

The ion selectivity of the bacterial potassium channel K_CSA is explained upon comparing the energy characteristics of the interaction of cations (Li⁺, Na⁺, K⁺) with atoms of the selectivity filter of the protein pore. Quantum-chemical calculations reveal a deeper potential well for potassium ions, which accounts for preferred K⁺ permeation. It is shown that the conventional methods with AMBER, CHARMM, OPLS force fields in standard parametrization as well as partial re-parametrization give incorrect estimates of ion energy distribution in the channel.     

The conserved potassium channel filter can have distinct ion binding profiles: Structural analysis of rubidium,cesium, and barium binding in NaK2K

Potassium channels are highly selective for K⁺ over the smaller Na⁺. Intriguingly, they are permeable to larger monovalent cations such as Rb⁺ and Cs⁺ but are specifically blocked by the similarly sized Ba²⁺. In this study, we used structural analysis to determine the binding profiles for these permeant and blocking ions in the selectivity filter of the potassium-selective NaK channel mutant NaK2K and also performed permeation experiments using single-channel recordings. Our data revealed that some ion binding properties of NaK2K are distinct from those of the canonical K⁺ channels KcsA and MthK. Rb⁺ bound at sites 1, 3, and 4 in NaK2K, as it does in KcsA. Cs⁺, however, bound predominantly at sites 1 and 3 in NaK2K, whereas it binds at sites 1, 3, and 4 in KcsA. Moreover, Ba²⁺ binding in NaK2K was distinct from that which has been observed in KcsA and MthK, even though all of these channels show similar Ba²⁺ block. In the presence of K⁺, Ba²⁺ bound to the NaK2K channel at site 3 in conjunction with a K⁺ at site 1; this led to a prolonged block of the channel (the external K⁺-dependent Ba²⁺ lock-in state). In the absence of K⁺, however, Ba²⁺ acts as a permeating blocker. We found that, under these conditions, Ba²⁺ bound at sites 1 or 0 as well as site 3, allowing it to enter the filter from the intracellular side and exit from the extracellular side. The difference in the Ba²⁺ binding profile in the presence and absence of K⁺ thus provides a structural explanation for the short and prolonged Ba²⁺ block observed in NaK2K.     

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.     

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.     

Structures of Gating Intermediates in a K+ channel

Regulation of ion conduction through the pore of a K⁺ channel takes place through the coordinated action of the activation gate at the bundle crossing of the inner helices and the inactivation gate located at the selectivity filter. The mechanism of allosteric coupling of these gates is of key interest. Here we report new insights into this allosteric coupling mechanism from studies on a W67F mutant of the KcsA channel. W67 is in the pore helix and is highly conserved in K⁺ channels. The KcsA W67F channel shows severely reduced inactivation and an enhanced rate of activation. We use continuous wave EPR spectroscopy to establish that the KcsA W67F channel shows an altered pH dependence of activation. Structural studies on the W67F channel provide the structures of two intermediate states: a pre- open state and a pre-inactivated state of the KcsA channel. These structures highlight key nodes in the allosteric pathway. The structure of the KcsA W67F channel with the activation gate open shows altered ion occupancy at the second ion binding site (S2) in the selectivity filter. This finding in combination with previous studies strongly support a requirement for ion occupancy at the S2 site for the channel to inactivate.     

Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K⁺ ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K⁺ channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K⁺ to Na⁺ substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na⁺ locally contracts the SF around Na⁺, releasing ligands that do not participate in Na⁺ coordination, and occluding the permeation pathway. K⁺ selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na⁺ ions due to basic structural differences between liquid water and the “quasi-liquid” SF matrix.     

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.     

Mutations in the K+-Channel KcsA Toward Kir Channels Alter Salt-Induced Clusterization and Blockade by Quaternary Alkylammonium Ions

Protein aggregation is a result of malfunction in protein folding, assembly, and transport, caused by protein mutation and/or changes in the cell environment, thus triggering many human diseases. We have shown that bacterial K⁺-channel KcsA, which acts as a representative model for ion channels, forms salt-induced large conductive complexes in a particular environment. In the present study, we investigated the effects of point mutations in the selectivity filter of KcsA on intrinsic stability, aggregation, and channel blocking behavior. First, we found that a low sodium chloride concentration in potassium-containing media induced fast transfer of single channels to a planar lipid bilayer. Second, increasing the sodium chloride concentration drastically increased the total channel current, indicating enhanced vesicle fusion and transfer of multiple channels to a planar lipid bilayer. However, such complexes exhibited high conductance as well as higher open probability compared to the unmodified KcsA behavior shown previously. Interestingly, the affinity of aggregated complexes for larger symmetric quaternary alkylammonium ions (QAs) was found to be much higher than that for tetraethylammonium, a classical blocker of the K⁺ channel. Based on these findings, we propose that mutant channel complexes exhibit larger pore dimensions, thus resembling more the topological properties of voltage-gated and inwardly rectifying K⁺ channels.     

Conduction of Na+ and K+ through the NaK channel: molecular and Brownian dynamics studies

Conduction of ions through the NaK channel, with M0 helix removed, was studied using both Brownian dynamics and molecular dynamics. Brownian dynamics simulations predict that the truncated NaK has approximately a third of the conductance of the related KcsA K⁺ channel, is outwardly rectifying, and has a Michaelis-Menten current-concentration relationship. Current magnitude increases when the glutamine residue located near the intracellular gate is replaced with a glutamate residue. The channel is blocked by extracellular Ca²⁺. Molecular dynamics simulations show that, under the influence of a strong applied potential, both Na⁺ and K⁺ move across the selectivity filter, although conduction rates for Na⁺ ions are somewhat lower. The mechanism of conduction of Na⁺ differs significantly from that of K⁺ in that Na⁺ is preferentially coordinated by single planes of pore-lining carbonyl oxygens, instead of two planes as in the usual K⁺ binding sites. The water-containing filter pocket resulting from a single change in the selectivity filter sequence (compared to potassium channels) disrupts several of the planes of carbonyl oxygens, and thus reduces the filter's ability to discriminate against sodium.