K(+) versus Na(+) ions in a K channel selectivity filter: a simulation study

Molecular dynamics simulations of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal significant differences in interactions of the selectivity filter with K(+) compared with Na(+) ions. K(+) ions and water molecules within the filter undergo concerted single-file motion in which they translocate between adjacent sites within the filter on a nanosecond timescale. In contrast, Na(+) ions remain bound to sites within the filter and do not exhibit translocation on a nanosecond timescale. Furthermore, entry of a K(+) ion into the filter from the extracellular mouth is observed, whereas this does not occur for a Na(+) ion. Whereas K(+) ions prefer to sit within a cage of eight oxygen atoms of the filter, Na(+) ions prefer to interact with a ring of four oxygen atoms plus two water molecules. These differences in interactions in the selectivity filter may contribute to the selectivity of KcsA for K(+) ions (in addition to the differences in dehydration energy between K(+) and Na(+)) and the block of KcsA by internal Na(+) ions. In our simulations the selectivity filter exhibits significant flexibility in response to changes in ion/protein interactions, with a somewhat greater distortion induced by Na(+) than by K(+) ions.     

Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer

Potassium channels enable K(+) ions to move passively across biological membranes. Multiple nanosecond-duration molecular dynamics simulations (total simulation time 5 ns) of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal motions of ions, water, and protein. Comparison of simulations with and without K(+) ions indicate that the absence of ions destabilizes the structure of the selectivity filter. Within the selectivity filter, K(+) ions interact with the backbone (carbonyl) oxygens, and with the side-chain oxygen of T75. Concerted single-file motions of water molecules and K(+) ions within the selectivity filter of the channel occur on a 100-ps time scale. In a simulation with three K(+) ions (initially two in the filter and one in the cavity), the ion within the central cavity leaves the channel via its intracellular mouth after approximately 900 ps; within the cavity this ion interacts with the Ogamma atoms of two T107 side chains, revealing a favorable site within the otherwise hydrophobically lined cavity. Exit of this ion from the channel is enabled by a transient increase in the diameter of the intracellular mouth. Such "breathing" motions may form the molecular basis of channel gating.     

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.     

Evidence for sequential ion-binding loci along the inner pore of the IRK1 inward-rectifier K+ channel

Steep rectification in IRK1 (Kir2.1) inward-rectifier K(+) channels reflects strong voltage dependence (valence of approximately 5) of channel block by intracellular cationic blockers such as the polyamine spermine. The observed voltage dependence primarily results from displacement, by spermine, of up to five K(+) ions across the narrow K(+) selectivity filter, along which the transmembrane voltage drops steeply. Spermine first binds, with modest voltage dependence, at a shallow site where it encounters the innermost K(+) ion and impedes conduction. From there, spermine can proceed to a deeper site, displacing several more K(+) ions and thereby producing most of the observed voltage dependence. Since in the deeper blocked state the leading amine group of spermine reaches into the cavity region (internal to the selectivity filter) and interacts with residue D172, its trailing end is expected to be near M183. Here, we found that mutation M183A indeed affected the deeper blocked state, which supports the idea that spermine is located in the region lined by the M2 and not deep in the narrow K(+) selectivity filter. As to the shallower site whose location has been unknown, we note that in the crystal structure of homologous GIRK1 (Kir3.1), four aromatic side chains of F255, one from each of the four subunits, constrict the intracellular end of the pore to approximately 10 A. For technical simplicity, we used tetraethylammonium (TEA) as an initial probe to test whether the corresponding residue in IRK1, F254, forms the shallower site. We found that replacing the aromatic side chain with an aliphatic one not only lowered TEA affinity of the shallower site approximately 100-fold but also eliminated the associated voltage dependence and, furthermore, confirmed that similar effects occurred also for spermine. These results establish the evidence for physically separate, sequential ion-binding loci along the long inner pore of IRK1, and strongly suggest that the aromatic side chains of F254 underlie the likely innermost binding locus for both blocker and K(+) ions in the cytoplasmic pore.     

The occupancy of ions in the K+ selectivity filter: charge balance and coupling of ion binding to a protein conformational change underlie high conduction rates

Potassium ions diffuse across the cell membrane in a single file through the narrow selectivity filter of potassium channels. The crystal structure of the KcsA K+ channel revealed the chemical structure of the selectivity filter, which contains four binding sites for K+. In this study, we used Tl+ in place of K+ to address the question of how many ions bind within the filter at a given time, i.e. what is the absolute ion occupancy? By refining the Tl+ structure against data to 1.9A resolution with an anomalous signal, we determined the absolute occupancy of Tl+. Then, by comparing the electron density of Tl+ with that of K+, Rb+ and Cs+, we estimated the absolute occupancy of these three ions. We further analyzed how the ion occupancy affects the conformation of the selectivity filter by analyzing the structure of KcsA at different concentrations of Tl+. Our results indicate that the average occupancy for each site in the selectivity filter is about 0.63 for Tl+ and 0.53 for K+. For K+, Rb+ and Cs+, the total number of ions contained within four sites in the selectivity filter is about two. At low concentrations of permeant ion, the number of ions drops to one in association with a conformational change in the selectivity filter. We conclude that electrostatic balance and coupling of ion binding to a protein conformational change underlie high conduction rates in the setting of high selectivity.     

Insight into the mechanism of inactivation and pH sensitivity in potassium channels from molecular dynamics simulations

Potassium (K (+)) channels can regulate ionic conduction through their pore by a mechanism, involving the selectivity filter, known as C-type inactivation. This process is rapid in the hERG K (+) channel and is fundamental to its physiological role. Although mutations within hERG are known to remove this process, a structural basis for the inactivation mechanism has yet to be characterized. Using MD simulations based on homology modeling, we observe that the carbonyl of the filter aromatic, Phe627, forming the S 0 K (+) binding site, swiftly rotates away from the conduction axis in the wild-type channel. In contrast, in well-characterized non-inactivating mutant channels, this conformational change occurs less frequently. In the non-inactivating channels, interactions with a water molecule located behind the selectivity filter are critical to the enhanced stability of the conducting state. We observe comparable conformational changes in the acid sensitive TASK-1 channel and propose a common mechanism in these channels for regulating efflux of K (+) ions through the selectivity filter.     

Permeability properties of ENaC selectivity filter mutants.

The epithelial Na(+) channel (ENaC), located in the apical membrane of tight epithelia, allows vectorial Na(+) absorption. The amiloride-sensitive ENaC is highly selective for Na(+) and Li(+) ions. There is growing evidence that the short stretch of amino acid residues (preM2) preceding the putative second transmembrane domain M2 forms the outer channel pore with the amiloride binding site and the narrow ion-selective region of the pore. We have shown previously that mutations of the alphaS589 residue in the preM2 segment change the ion selectivity, making the channel permeant to K(+) ions. To understand the molecular basis of this important change in ionic selectivity, we have substituted alphaS589 with amino acids of different sizes and physicochemical properties. Here, we show that the molecular cutoff of the channel pore for inorganic and organic cations increases with the size of the amino acid residue at position alpha589, indicating that alphaS589 mutations enlarge the pore at the selectivity filter. Mutants with an increased permeability to large cations show a decrease in the ENaC unitary conductance of small cations such as Na(+) and Li(+). These findings demonstrate the critical role of the pore size at the alphaS589 residue for the selectivity properties of ENaC. Our data are consistent with the main chain carbonyl oxygens of the alphaS589 residues lining the channel pore at the selectivity filter with their side chain pointing away from the pore lumen. We propose that the alphaS589 side chain is oriented toward the subunit-subunit interface and that substitution of alphaS589 by larger residues increases the pore diameter by adding extra volume at the subunit-subunit interface.     

Conduction properties of KcsA measured using brownian dynamics with flexible carbonyl groups in the selectivity filter

In the narrow segment of an ion conducting pathway, it is likely that a permeating ion influences the positions of the nearby atoms that carry partial or full electronic charges. Here we introduce a method of incorporating the motion of charged atoms lining the pore into Brownian dynamics simulations of ion conduction. The movements of the carbonyl groups in the selectivity filter of the KcsA channel are calculated explicitly, allowing their bond lengths, bond angles, and dihedral angels to change in response to the forces acting upon them. By systematically changing the coefficients of bond stretching and of angle bending, the carbon and oxygen atoms can be made to fluctuate from their fixed positions by varying mean distances. We show that incorporating carbonyl motion in this way does not alter the mechanism of ion conduction and only has a small influence on the computed current. The slope conductance of the channel increases by approximately 25% when the root mean-square fluctuations of the carbonyl groups are increased from 0.01 to 0.61 A. The energy profiles and the number of resident ions in the channel remain unchanged. The method we utilized here can be extended to allow the movement of glutamate or aspartate side chains lining the selectivity filters of other ionic channels.     

The barium site in a potassium channel by x-ray crystallography

X-ray diffraction data were collected from frozen crystals (100 degrees K) of the KcsA K(+) channel equilibrated with solutions containing barium chloride. Difference electron density maps (F(barium) - F(native), 5.0 A resolution) show that Ba(2+) resides at a single location within the selectivity filter. The Ba(2+) blocking site corresponds to the internal aspect (adjacent to the central cavity) of the "inner ion" position where an alkali metal cation is found in the absence of the blocking Ba(2+) ion. The location of Ba(2+) with respect to Rb(+) ions in the pore is in good agreement with the findings on the functional interaction of Ba(2+) with K(+) (and Rb(+)) in Ca(2+)-activated K(+) channels (Neyton, J., and C. Miller. 1988. J. Gen. Physiol. 92:549-567). Taken together, these structural and functional data imply that at physiological ion concentrations a third ion may interact with two ions in the selectivity filter, perhaps by entering from one side and displacing an ion on the opposite side.     

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.     

Ion selectivity in potassium channels

Potassium channels are tetrameric membrane-spanning proteins that provide a selective pore for the conduction of K(+) across the cell membranes. One of the main physiological functions of potassium channels is efficient and very selective transport of K(+) ions through the membrane to the cell. Classical views of ion selectivity are summarized within a historical perspective, and contrasted with the molecular dynamics (MD) simulations free energy perturbation (FEP) performed on the basis of the crystallographic structure of the KcsA phospholipid membrane. The results show that the KcsA channel does not select for K(+) ions by providing a binding site of an appropriate (fixed) cavity size. Rather, selectivity for K(+) arises directly from the intrinsic local physical properties of the ligands coordinating the cation in the binding site, and is a robust feature of a pore symmetrically lined by backbone carbonyl groups. Further analysis reveals that it is the interplay between the attractive ion-ligand (favoring smaller cation) and repulsive ligand-ligand interactions (favoring larger cations) that is the basic element governing Na(+)/K(+) selectivity in flexible protein binding sites. Because the number and the type of ligands coordinating an ion directly modulate such local interactions, this provides a potent molecular mechanism to achieve and maintain a high selectivity in protein binding sites despite a significant conformational flexibility.     

Molecular dynamics of the KcsA K(+) channel in a bilayer membrane

Molecular dynamics (MD) simulations of an atomic model of the KcsA K(+) channel embedded in an explicit dipalmitoylphosphatidylcholine (DPPC) phospholipid bilayer solvated by a 150 mM KCl aqueous salt solution are performed and analyzed. The model includes the KcsA K(+) channel, based on the recent crystallographic structure of, Science. 280:69-77), 112 DPPC, K(+) and Cl(-) ions, as well as over 6500 water molecules for a total of more than 40,000 atoms. Three K(+) ions are explicitly included in the pore. Two are positioned in the selectivity filter on the extracellular side and one in the large water-filled cavity. Different starting configurations of the ions and water molecules in the selectivity filter are considered, and MD trajectories are generated for more than 4 ns. The conformation of KcsA is very stable in all of the trajectories, with a global backbone root mean square (RMS) deviation of less than 1.9 A with respect to the crystallographic structure. The RMS atomic fluctuations of the residues surrounding the selectivity filter on the extracellular side of the channel are significantly lower than those on the intracellular side. The motion of the residues with aromatic side chains surrounding the selectivity filter (Trp(67), Trp(68), Tyr(78), and Tyr(82)) is anisotropic with the smallest RMS fluctuations in the direction parallel to the membrane plane. A concerted dynamic transition of the three K(+) ions in the pore is observed, during which the K(+) ion located initially in the cavity moves into the narrow part of the selectivity filter, while the other two K(+) ions move toward the extracellular side. A single water molecule is stabilized between each pair of ions during the transition, suggesting that each K(+) cation translocating through the narrow pore is accompanied by exactly one water molecule, in accord with streaming potential measurements (, Biophys. J. 55:367-371). The displacement of the ions is coupled with the structural fluctuations of Val(76) and Gly(77), in the selectivity filter, as well as the side chains of Glu(71), Asp(80), and Arg(89), near the extracellular side. Thus the mechanical response of the channel structure at distances as large as 10-20 A from the ions in the selectivity filter appears to play an important role in the concerted transition.     

A computational study of ion binding and protonation states in the KcsA potassium channel

We report results from microscopic molecular dynamics and free energy perturbation simulations of the KcsA potassium channel based on its experimental atomic structure. Conformational properties of selected amino acid residues as well as equilibrium positions of K(+) ions inside the selectivity filter and the internal water cavity are examined. Positions three and four (counting from the extracellular site) in the experimental structure correspond to distinctly separate binding sites for K(+) ions inside the selectivity filter. The protonation states of Glu71 and Asp80, which are close to each other and to the selectivity filter, as well as K(+) binding energies are determined using free energy perturbation calculations. The Glu71 residue which is buried inside a protein cavity is found to be most stable in the neutral form while the solvent exposed Asp80 is ionized. The channel altogether exothermically binds up to three ions, where two of them are located inside the selectivity filter and one in the internal water cavity. Ion permeation mechanisms are discussed in relation to these results.     

Potassium channel, ions, and water: simulation studies based on the high resolution X-ray structure of KcsA

Interactions of Na(+), K(+), Rb(+), and Cs(+) ions within the selectivity filter of a potassium channel have been investigated via multiple molecular dynamics simulations (total simulation time, 48 ns) based on the high resolution structure of KcsA, embedded in a phospholipid bilayer. As in simulations based on a lower resolution structure of KcsA, concerted motions of ions and water within the filter are seen. Despite the use of a higher resolution structure and the inclusion of four buried water molecules thought to stabilize the filter, this region exhibits a significant degree of flexibility. In particular, pronounced distortion of filter occurs if no ions are present within it. The two most readily permeant ions, K(+) and Rb(+), are similar in their interactions with the selectivity filter. In contrast, Na(+) ions tend to distort the filter by binding to a ring of four carbonyl oxygens. The larger Cs(+) ions result in a small degree of expansion of the filter relative to the x-ray structure. Cs(+) ions also appear to interact differently with the gate region of the channel, showing some tendency to bind within a predominantly hydrophobic pocket. The four water molecules buried between the back of the selectivity filter and the remainder of the protein show comparable mobility to the surrounding protein and do not exchange with water molecules within the filter or the central cavity. A preliminary comparison of the use of particle mesh Ewald versus cutoff protocols for the treatment of long-range electrostatics suggests some difference in the kinetics of ion translocation within the filter.     

Inhibition of the collapse of the Shaker K+ conductance by specific scorpion toxins

The Shaker B K(+) conductance (G(K)) collapses when the channels are closed (deactivated) in Na(+) solutions that lack K(+) ions. Also, it is known that external TEA (TEA(o)) impedes the collapse of G(K), and that channel block by TEA(o) and scorpion toxins are two mutually exclusive events. Therefore, we tested the ability of scorpion toxins to inhibit the collapse of G(K) in 0 K(+). We have found that these toxins are not uniform regarding the capacity to protect G(K). Those toxins, whose binding to the channels is destabilized by external K(+), are also effective inhibitors of the collapse of G(K). In addition to K(+), other externally added cations also destabilize toxin block, with an effectiveness that does not match the selectivity sequence of K(+) channels. The inhibition of the drop of G(K) follows a saturation relationship with [toxin], which is fitted well by the Michaelis-Menten equation, with an apparent Kd bigger than that of block of the K(+) current. However, another plausible model is also presented and compared with the Michaelis-Menten model. The observations suggest that those toxins that protect G(K) in 0 K(+) do so by interacting either with the most external K(+) binding site of the selectivity filter (suggesting that the K(+) occupancy of only that site of the pore may be enough to preserve G(K)) or with sites capable of binding K(+) located in the outer vestibule of the pore, above the selectivity filter.     

Molecular dynamics simulations of the TrkH membrane protein

TrkH is a transmembrane protein that mediates uptake of K(+) through the cell membrane. Despite the recent determination of its crystallographic structure, the nature of the permeation mechanism is still unknown, that is, whether K(+) ions move across TrkH by active transport or passive diffusion. Here, molecular dynamics simulations and the umbrella sampling technique have been employed to shed light on this question. The existence of binding site S3 and two alternative binding sites have been characterized. Analysis of the coordination number renders values that are almost constant, with a full contribution from the carbonyls of the protein only at S3. This observation contrasts with observations of K(+) channels, where the contribution of the protein to the coordination number is roughly constant in all four binding sites. An intramembrane loop is found immediately after the selectivity filter at the intracellular side of the protein, which obstructs the permeation pathway, and this is reflected in the magnitude of the energy barriers.     

Ion selectivity and current saturation in inward-rectifier K+ channels

We investigated the features of the inward-rectifier K channel Kir1.1 (ROMK) that underlie the saturation of currents through these channels as a function of permeant ion concentration. We compared values of maximal currents and apparent K(m) for three permeant ions: K(+), Rb(+), and NH(4)(+). Compared with K(+) (i(max) = 4.6 pA and K(m) = 10 mM at -100 mV), Rb(+) had a lower permeability, a lower i(max) (1.8 pA), and a higher K(m) (26 mM). For NH(4)(+), the permeability was reduced more with smaller changes in i(max) (3.7 pA) and K(m) (16 mM). We assessed the role of a site near the outer mouth of channel in the saturation process. This site could be occupied by either permeant ions or low-affinity blocking ions such as Na(+), Li(+), Mg(2+), and Ca(2+) with similar voltage dependence (apparent valence, 0.15-0.20). It prefers Mg(2+) over Ca(2+) and has a monovalent cation selectivity, based on the ability to displace Mg(2+), of K(+) > Li(+) ～ Na(+) > Rb(+) ～ NH(4)(+). Conversely, in the presence of Mg(2+), the K(m) for K(+) conductance was substantially increased. The ability of Mg(2+) to block the channels was reduced when four negatively charged amino acids in the extracellular domain of the channel were mutated to neutral residues. The apparent K(m) for K(+) conduction was unchanged by these mutations under control conditions but became sensitive to the presence of external negative charges when residual divalent cations were chelated with EDTA. The results suggest that a binding site in the outer mouth of the pore controls current saturation. Permeability is more affected by interactions with other sites within the selectivity filter. Most features of permeation (and block) could be simulated by a five-state kinetic model of ion movement through the channel.     

Na+ block and permeation in a K+ channel of known structure

The effects of intracellular Na(+) were studied on K(+) and Rb(+) currents through single KcsA channels. At low voltage, Na(+) produces voltage-dependent block, which becomes relieved at high voltage by a "punchthrough" mechanism representing Na(+) escaping from its blocking site through the selectivity filter. The Na(+) blocking site is located in the wide, hydrated vestibule, and it displays unexpected selectivity for K(+) and Rb(+) against Na(+). The voltage dependence of Na(+) block reflects coordinated movements of the blocker with permeant ions in the selectivity filter.     

Voltage-dependent inactivation gating at the selectivity filter of the MthK K+ channel

Voltage-dependent K(+) channels can undergo a gating process known as C-type inactivation, which involves entry into a nonconducting state through conformational changes near the channel's selectivity filter. C-type inactivation may involve movements of transmembrane voltage sensor domains, although the mechanisms underlying this form of inactivation may be heterogeneous and are often unclear. Here, we report on a form of voltage-dependent inactivation gating observed in MthK, a prokaryotic K(+) channel that lacks a canonical voltage sensor and may thus provide a reduced system to inform on mechanism. In single-channel recordings, we observe that Po decreases with depolarization, with a half-maximal voltage of 96 ± 3 mV. This gating is kinetically distinct from blockade by internal Ca(2+) or Ba(2+), suggesting that it may arise from an intrinsic inactivation mechanism. Inactivation gating was shifted toward more positive voltages by increasing external [K(+)] (47 mV per 10-fold increase in [K(+)]), suggesting that K(+) binding at the extracellular side of the channel stabilizes the open-conductive state. The open-conductive state was stabilized by other external cations, and selectivity of the stabilizing site followed the sequence: K(+) ≈ Rb(+) > Cs(+) > Na(+) > Li(+) ≈ NMG(+). Selectivity of the stabilizing site is weaker than that of sites that determine permeability of these ions, suggesting that the site may lie toward the external end of the MthK selectivity filter. We could describe MthK gating over a wide range of positive voltages and external [K(+)] using kinetic schemes in which the open-conductive state is stabilized by K(+) binding to a site that is not deep within the electric field, with the voltage dependence of inactivation arising from both voltage-dependent K(+) dissociation and transitions between nonconducting (inactivated) states. These results provide a quantitative working hypothesis for voltage-dependent, K(+)-sensitive inactivation gating, a property that may be common to other K(+) channels.