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.     

Comparative study of the energetics of ion permeation in Kv1.2 and KcsA potassium channels

Biological ion channels rely on a multi-ion transport mechanism for fast yet selective permeation of ions. The crystal structure of the KcsA potassium channel provided the first microscopic picture of this process. A similar mechanism is assumed to operate in all potassium channels, but the validity of this assumption has not been well investigated. Here, we examine the energetics of ion permeation in Shaker Kv1.2 and KcsA channels, which exemplify the six-transmembrane voltage-gated and two-transmembrane inward-rectifier channels. We study the feasibility of binding a third ion to the filter and the concerted motion of ions in the channel by constructing the potential of mean force for K⁺ ions in various configurations. For both channels, we find that a pair of K⁺ ions can move almost freely within the filter, but a relatively large free-energy barrier hinders the K⁺ ion from stepping outside the filter. We discuss the effect of the CMAP dihedral energy correction that was recently incorporated into the CHARMM force field on ion permeation dynamics.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Ion behavior in the selectivity filter of HCN1 channels

Hyperpolarization-activated cyclic-nucleotide gated channels (HCNs) are responsible for the generation of pacemaker currents (I_f or I_h) in cardiac and neuronal cells. Despite the overall structural similarity to voltage-gated potassium (Kv) channels, HCNs show much lower selectivity for K⁺ over Na⁺ ions. This increased permeability to Na⁺ is critical to their role in membrane depolarization. HCNs can also select between Na⁺ and Li⁺ ions. Here, we investigate the unique ion selectivity properties of HCNs using molecular-dynamics simulations. Our simulations suggest that the HCN1 pore is flexible and dilated compared with Kv channels with only one stable ion binding site within the selectivity filter. We also observe that ion coordination and hydration differ within the HCN1 selectivity filter compared with those in Kv and cyclic-nucleotide gated channels. Additionally, the C358T mutation further stabilizes the symmetry of the binding site and provides a more fit space for ion coordination, particularly for Li⁺.     

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.     

K binding in the G-loop and water cavity facilitates Ba movement in the Kir2.1 channel

K⁺ are selectively coordinated in the selectivity filter and concerted K⁺ and water movements in this region ensure high conduction rates in K⁺ channels. In channels with long pores many K⁺ binding sites are located intracellular to the selectivity filter (inner vestibule), but their contribution to permeation has not been well studied. We investigated this phenomenon by slowing the ion permeation process via blocking inwardly rectifying Kir2.1 channels with Ba²⁺ in the selectivity filter and observing the effect of K⁺ in the inner vestibule on Ba²⁺ exit. The dose-response effect of the intracellular K⁺ concentration ([K⁺]_i) on Ba²⁺ exit was recorded with and without intracellular polyamines, which compete with K⁺ for binding sites. Ba²⁺ exit was facilitated by the cooperative binding of at least three K⁺. Site-directed mutagenesis studies suggest that K⁺ interacting with Ba²⁺ bound in the selectivity filter were located in the region between selectivity filter and cytoplasmic pore, i.e. the water cavity and G-loop. One of the K⁺ binding sites was located at residue D172 and another was possibly at M301. This study provides functional evidence for the three K⁺ binding sites in the inner vestibule previously identified by crystal structure study.     

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.     

Molecular basis of ion permeability in a voltage‐gated sodium channel

Voltage‐gated sodium channels are essential for electrical signalling across cell membranes. They exhibit strong selectivities for sodium ions over other cations, enabling the finely tuned cascade of events associated with action potentials. This paper describes the ion permeability characteristics and the crystal structure of a prokaryotic sodium channel, showing for the first time the detailed locations of sodium ions in the selectivity filter of a sodium channel. Electrostatic calculations based on the structure are consistent with the relative cation permeability ratios (Na⁺ ≈ Li⁺ ≫ K⁺, Ca²⁺, Mg²⁺) measured for these channels. In an E178D selectivity filter mutant constructed to have altered ion selectivities, the sodium ion binding site nearest the extracellular side is missing. Unlike potassium ions in potassium channels, the sodium ions in these channels appear to be hydrated and are associated with side chains of the selectivity filter residues, rather than polypeptide backbones.     

Ion Conduction through the hERG Potassium Channel

The inward rectifier voltage-gated potassium channel hERG is of primary importance for the regulation of the membrane potential of cardiomyocytes. Unlike most voltage-gated K⁺-channels, hERG shows a low elementary conductance at physiological voltage and potassium concentration. To investigate the molecular features underlying this unusual behavior, we simulated the ion conduction through the selectivity filter at a fully atomistic level by means of molecular dynamics-based methods, using a homology-derived model. According to our calculations, permeation of potassium ions can occur along two pathways, one involving site vacancies inside the filter (showing an energy barrier of about 6 kcal mol⁻¹), and the other characterized by the presence of a knock-on intermediate (about 8 kcal mol⁻¹). These barriers are indeed in accordance with a low conductance behavior, and can be explained in terms of a series of distinctive structural features displayed by the hERG ion permeation pathway.     

Selecting Ions by Size in a Calcium Channel: The Ryanodine Receptor Case Study

Many calcium channels can distinguish between ions of the same charge but different size. For example, when cations are in direct competition with each other, the ryanodine receptor (RyR) calcium channel preferentially conducts smaller cations such as Li⁺ and Na⁺ over larger ones such as K⁺ and Cs⁺. Here, we analyze the physical basis for this preference using a previously established model of RyR permeation and selectivity. Like other calcium channels, RyR has four aspartate residues in its GGGIGDE selectivity filter. These aspartates have their terminal carboxyl group in the pore lumen, which take up much of the available space for permeating ions. We find that small ions are preferred by RyR because they can fit into this crowded environment more easily.     

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.     

Determinants of cation transport selectivity: Equilibrium binding and transport kinetics

ConclusionThe equilibrium ion-binding properties of ion channels and transporters can be difficult to discern from crystal structures alone, as proteins often adopt different lowest energy states depending on the ions bound. In cases where transport is slow, their inherent ion-binding preferences can be used to infer their transport preferences. However, in cases where transport is fast, the transport selectivity can hide their equilibrium preferences by accentuating the kinetics of ions hopping through a channel over its inherent ion-binding preferences. Thus, depending on the arrangement of ion-binding sites in a channel’s selectivity filter, one can achieve either selective or nonselective ion transport.The equilibrium K⁺ selectivity of some nonselective channels suggests a potential mechanism whereby they could evolve into a fast K⁺-selective channel. K⁺ channels and nonselective channels like CNG and HCN are related to one another in both sequence and structure, suggesting an evolutionary link between them. Swap experiments show that only a few mutations separate a nonselective channel from a K⁺-selective channel. One might imagine an evolutionary path between these channels in which the equilibrium preference for a K⁺ ion in a nonselective channel evolves into a K⁺-selective channel through these few mutations to create the selective ion queue. Alternatively, a slow single-ion channel with an equilibrium and transport preference for K⁺ ions could be transformed into a fast multi-ion channel through mutations that create a queue of K⁺-selective ion-binding sites, as is seen in most K⁺ channels studied to date.In the case of multi-ion selectivity filters, such as those found in K⁺ channels, the selectivity filter can be viewed as the active site that interacts with different queues of ions and water molecules. At least three properties emerge from multi-ion queues: (1) high conductance by reducing the affinity of multiple bound ions versus single ions; (2) high selectivity by allowing disfavored ions time to dissociate back into solution; and, consequently, (3) robust selectivity in an environment where ion concentrations can change. For transporters and carriers, the equilibrium preference and slow transport naturally create robust selectivity. In all these cases, equilibrium-based ion selectivity is achieved by slowing transport enough so that the disfavored ion is able to dissociate back into solution before transport takes place.