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.     

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.     

Plants do it differently. A new basis for potassium/sodium selectivity in the pore of an ion channel

Understanding of the molecular architecture necessary for selective K(+) permeation through the pore of ion channels is based primarily on analysis of the crystal structure of the bacterial K(+) channel KcsA, and structure:function studies of cloned animal K(+) channels. Little is known about the conduction properties of a large family of plant proteins with structural similarities to cloned animal cyclic nucleotide-gated channels (CNGCs). Animal CNGCs are nonselective cation channels that do not discriminate between Na(+) and K(+) permeation. These channels all have the same triplet of amino acids in the channel pore ion selectivity filter, and this sequence is different from that of the selectivity filter found in K(+)-selective channels. Plant CNGCs have unique pore selectivity filters; unlike those found in any other family of channels. At present, the significance of the unique pore selectivity filters of plant CNGCs, with regard to discrimination between Na(+) and K(+) permeation is unresolved. Here, we present an electrophysiological analysis of several members of this protein family; identifying the first cloned plant channel (AtCNGC1) that conducts Na(+). Another member of this ion channel family (AtCNGC2) is shown to have a selectivity filter that provides a heretofore unknown molecular basis for discrimination between K(+) and Na(+) permeation. Specific amino acids within the AtCNGC2 pore selectivity filter (Asn-416, Asp-417) are demonstrated to facilitate K(+) over Na(+) conductance. The selectivity filter of AtCNGC2 represents an alternative mechanism to the well-known GYG amino acid triplet of K(+) channels that has been identified as the critical basis for K(+) over Na(+) permeation through the pore of ion channels.     

Contribution of ion binding affinity to ion selectivity and permeation in KcsA, a model potassium channel

Ion permeation and selectivity, key features in ion channel function, are believed to arise from a complex ensemble of energetic and kinetic variables. Here we evaluate the contribution of pore cation binding to ion permeation and selectivity features of KcsA, a model potassium channel. For this, we used E71A and M96V KcsA mutants in which the equilibrium between conductive and nonconductive conformations of the channel is differently shifted. E71A KcsA is a noninactivating channel mutant. Binding of K(+) to this mutant reveals a single set of low-affinity K(+) binding sites, similar to that seen in the binding of K(+) to wild-type KcsA that produces a conductive, low-affinity complex. This seems consistent with the observed K(+) permeation in E71A. Nonetheless, the E71A mutant retains K(+) selectivity, which cannot be explained on the basis of just its low affinity for this ion. At variance, M96V KcsA is a rapidly inactivating mutant that has lost selectivity for K(+) and also conducts Na(+). Here, low-affinity binding and high-affinity binding of both cations are detected, seemingly in agreement with both being permeating species in this mutant channel. In conclusion, binding of the ion to the channel protein seemingly explains certain gating, ion selectivity, and permeation properties. Ion binding stabilizes greatly the channel and, depending upon ion type and concentration, leads to different conformations and ion binding affinities. High-affinity states guarantee binding of specific ions and mediate ion selectivity but are nonconductive. Conversely, low-affinity states would not discriminate well among different ions but allow permeation to occur.     

Ion binding to KcsA: implications in ion selectivity and channel gating

Binding of K+ and Na+ to the potassium channel KcsA has been characterized from the stabilization observed in the heat-induced denaturation of the protein as the ion concentration is increased. KcsA thermal denaturation is known to include (i) dissociation of the homotetrameric channel into its constituent subunits and (ii) protein unfolding. The ion concentration-dependent changes in the thermal stability of the protein, evaluated as the Tm value for thermal-induced denaturation of the protein, may suggest the existence of both high- and low-affinity K+ binding sites of KcsA, which lend support to the tenet that channel gating may be governed by K+ concentration-dependent transitions between different affinity states of the channel selectivity filter. We also found that Na+ binds to KcsA with a KD similar to that estimated electrophysiologically from channel blockade. Therefore, our findings on ion binding to KcsA partly account for K+ over Na+ selectivity and Na+ blockade and argue against the strict “snug fit” hypothesis used initially to explain ion selectivity from the X-ray channel structure. Furthermore, the remarkable effects of increasing the ion concentration, K+ in particular, on the Tm of the denaturation process evidence that synergistic effects of the metal-mediated intersubunit interactions at the channel selectivity filter are a major contributor to the stability of the tetrameric protein. This observation substantiates the notion of a role for ions as structural “effectors” of ion channels.     

Modeling permeation energetics in the KcsA potassium channel

The thermodynamics of cation permeation through the KcsA K(+) channel selectivity filter is studied from the perspective of a physically transparent semimicroscopic model using Monte Carlo free energy integration. The computational approach chosen permits dissection of the separate contributions to ionic stabilization arising from different parts of the channel (selectivity filter carbonyls, single-file water, cavity water, reaction field of bulk water, inner helices, ionizable residues). All features play important roles; their relative significance varies with the ion's position in the filter. The cavity appears to act as an electrostatic buffer, shielding filter ions from structural changes in the inner pore. The model exhibits K(+) vs. Na(+) selectivity, and roughly isoenergetic profiles for K(+) and Rb(+), and discriminates against Cs(+), all in agreement with experimental data. It also indicates that Ba(2+) and Na(+) compete effectively with permeant ions at a site near the boundary between the filter and the cavity, in the vicinity of the barium blocker site.     

Influence of pore residues on permeation properties in the Kv2.1 potassium channel. Evidence for a selective functional interaction of K+ with the outer vestibule

The Kv2.1 potassium channel contains a lysine in the outer vestibule (position 356) that markedly reduces open channel sensitivity to changes in external [K(+)]. To investigate the mechanism underlying this effect, we examined the influence of this outer vestibule lysine on three measures of K(+) and Na(+) permeation. Permeability ratio measurements, measurements of the lowest [K(+)] required for interaction with the selectivity filter, and measurements of macroscopic K(+) and Na(+) conductance, were all consistent with the same conclusion: that the outer vestibule lysine in Kv2.1 interferes with the ability of K(+) to enter or exit the extracellular side of the selectivity filter. In contrast to its influence on K(+) permeation properties, Lys 356 appeared to be without effect on Na(+) permeation. This suggests that Lys 356 limited K(+) flux by interfering with a selective K(+) binding site. Combined with permeation studies, results from additional mutagenesis near the external entrance to the selectivity filter indicated that this site was located external to, and independent from, the selectivity filter. Protonation of a naturally occurring histidine in the same outer vestibule location in the Kv1.5 potassium channel produced similar effects on K(+) permeation properties. Together, these results indicate that a selective, functional K(+) binding site (e.g., local energy minimum) exists in the outer vestibule of voltage-gated K(+) channels. We suggest that this site is the location of K(+) hydration/dehydration postulated to exist based on the structural studies of KcsA. Finally, neutralization of position 356 enhanced outward K(+) current magnitude, but did not influence the ability of internal K(+) to enter the pore. These data indicate that in Kv2.1, exit of K(+) from the selectivity filter, rather than entry of internal K(+) into the channel, limits outward current magnitude. We discuss the implications of these findings in relation to the structural basis of channel conductance in different K(+) channels.     

Absence of ion-binding affinity in the putatively inactivated low-[K+] structure of the KcsA potassium channel

Potassium channels are membrane proteins that selectively conduct K(+) across cellular membranes. The narrowest part of their pore, the selectivity filter, is responsible for distinguishing K(+) from Na(+), and can also act as a gate through a mechanism known as C-type inactivation. It has been proposed that a conformation of the KcsA channel obtained by crystallization in presence of low concentration of K(+) (PDB 1K4D) could correspond to the C-type inactivated state. Here, we show using molecular mechanics simulations that such conformation has little ion-binding affinity and that ions do not contribute to its stability. The simulations suggest that, in this conformation, the selectivity filter is mostly occupied by water molecules. Whether such ion-free state of the KcsA channel is physiologically accessible and representative of the inactivated state of eukaryotic channels remains unclear.     

Two stable, conducting conformations of the selectivity filter in Shaker K+ channels

We have examined the voltage dependence of external TEA block of Shaker K(+) channels over a range of internal K(+) concentrations from 2 to 135 mM. We found that the concentration dependence of external TEA block in low internal K(+) solutions could not be described by a single TEA binding affinity. The deviation from a single TEA binding isotherm was increased at more depolarized membrane voltages. The data were well described by a two-component binding scheme representing two, relatively stable populations of conducting channels that differ in their affinity for external TEA. The relative proportion of these two populations was not much affected by membrane voltage but did depend on the internal K(+) concentration. Low internal K(+) promoted an increase in the fraction of channels with a low TEA affinity. The voltage dependence of the apparent high-affinity TEA binding constant depended on the internal K(+) concentration, becoming almost voltage independent in 5 mM. The K(+) sensitivity of these low- and high-affinity TEA states suggests that they may represent one- and two-ion occupancy states of the selectivity filter, consistent with recent crystallographic results from the bacterial KcsA K(+) channel. We therefore analyzed these data in terms of such a model and found a large (almost 14-fold) difference between the intrinsic TEA affinity of the one-ion and two-ion modes. According to this analysis, the single ion in the one-ion mode (at 0 mV) prefers the inner end of the selectivity filter twofold more than the outer end. This distribution does not change with internal K(+). The two ions in the two-ion mode prefer to occupy the inner end of the selectivity filter at low K(+), but high internal K(+) promotes increased occupancy of the outer sites. Our analysis further suggests that the four K(+) sites in the selectivity filter are spaced between 20 and 25% of the membrane electric field.     

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.     

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.     

Hierarchical approach to predicting permeation in ion channels

A hierarchical computational strategy combining molecular modeling, electrostatics calculations, molecular dynamics, and Brownian dynamics simulations is developed and implemented to compute electrophysiologically measurable properties of the KcsA potassium channel. Models for a series of channels with different pore sizes are developed from the known x-ray structure, using insights into the gating conformational changes as suggested by a variety of published experiments. Information on the pH dependence of the channel gating is incorporated into the calculation of potential profiles for K(+) ions inside the channel, which are then combined with K(+) ion mobilities inside the channel, as computed by molecular dynamics simulations, to provide inputs into Brownian dynamics simulations for computing ion fluxes. The open model structure has a conductance of approximately 110 pS under symmetric 250 mM K(+) conditions, in reasonable agreement with experiments for the largest conducting substate. The dimensions of this channel are consistent with electrophysiologically determined size dependence of quaternary ammonium ion blocking from the intracellular end of this channel as well as with direct structural evidence that tetrabutylammonium ions can enter into the interior cavity of the channel. Realistic values of Ussing flux ratio exponents, distribution of ions within the channel, and shapes of the current-voltage and current-concentration curves are obtained. The Brownian dynamics calculations suggest passage of ions through the selectivity filter proceeds by a "knock-off" mechanism involving three ions, as has been previously inferred from functional and structural studies of barium ion blocking. These results suggest that the present calculations capture the essential nature of K(+) ion permeation in the KcsA channel and provide a proof-of-concept for the integrated microscopic/mesoscopic multitiered approach for predicting ion channel function from structure, which can be applied to other channel structures.     

Effects of conducting and blocking ions on the structure and stability of the potassium channel KcsA

This article reports on the interaction of conducting (K(+)) and blocking (Na(+)) monovalent metal ions with detergent-solubilized and lipid-reconstituted forms of the K(+) channel KcsA. Monitoring of the protein intrinsic fluorescence reveals that the two ions bind competitively to KcsA with distinct affinities (dissociation constants for the KcsA.K(+) and KcsA.Na(+) complexes of approximately 8 and 190 mm, respectively) and induce different conformations of the ion-bound protein. The differences in binding affinity as well as the higher K(+) concentration bathing the intracellular mouth of the channel, through which the cations gain access to the protein binding sites, should favor that only KcsA.K(+) complexes are formed under physiological-like conditions. Nevertheless, despite such prediction, it was also found that concentrations of Na(+) well below its dissociation constant and even in the presence of higher K(+) concentrations, cause a remarkable decrease in the protein thermal stability and facilitate thermal dissociation into subunits of the tetrameric KcsA, as concluded from the temperature dependence of the protein infrared spectra and from gel electrophoresis, respectively. These latter observations cannot be explained based on the occupancy of the binding sites from above and suggest that there must be additional ion binding sites, whose occupancy could not be detected by fluorescence and in which the affinity for Na(+) must be higher or at least similar to that of K(+). Moreover, cation binding as reported by means of fluorescence does not suffice to explain the large differences in free energy of stabilization involved in the formation of the KcsA.Na(+) and KcsA.K(+) complexes, which for the most part should arise from synergistic effects of the ion-mediated intersubunit interactions.     

Homology modeling and molecular dynamics simulation studies of an inward rectifier potassium channel

A homology model has been generated for the pore-forming domain of Kir6.2, a component of an ATP-sensitive K channel, based on the x-ray structure of the bacterial channel KcsA. Analysis of the lipid-exposed and pore-lining surfaces of the model reveals them to be compatible with the known features of membrane proteins and Kir channels, respectively. The Kir6.2 homology model was used as the starting point for nanosecond-duration molecular dynamics simulations in a solvated phospholipid bilayer. The overall drift from the model structure was comparable to that seen for KcsA in previous similar simulations. Preliminary analysis of the interactions of the Kir6.2 channel model with K(+) ions and water molecules during these simulations suggests that concerted single-file motion of K(+) ions and water through the selectivity filter occurs. This is similar to such motion observed in simulations of KcsA. This suggests that a single-filing mechanism is conserved between different K channel structures and may be robust to changes in simulation details. Comparison of Kir6.2 and KcsA suggests some degree of flexibility in the filter, thus complicating models of ion selectivity based upon a rigid 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.     

Exploring the origin of the ion selectivity of the KcsA potassium channel

The availability of structural information about biological ion channels provides an opportunity to gain a detailed understanding of the control of ion selectivity by biological systems. However, accomplishing this task by computer simulation approaches is very challenging. First, although the activation barriers for ion transport can be evaluated by microscopic simulations, it is hard to obtain accurate results by such approaches. Second, the selectivity is related to the actual ion current and not directly to the individual activation barriers. Thus, it is essential to simulate the ion currents and this cannot be accomplished at present by microscopic MD approaches. In order to address this challenge, we developed and refined an approach capable of evaluating ion current while still reflecting the realistic features of the given channel. Our method involves generation of semimacroscopic free energy surfaces for the channel/ions system and Brownian dynamics (BD) simulations of the corresponding ion current. In contrast to most alternative macroscopic models, our approach is able to reproduce the difference between the free energy surfaces of different ions and thus to address the selectivity problem. Our method is used in a study of the selectivity of the KcsA channel toward the K+ and Na+ ions. The BD simulations with the calculated free energy profiles produce an appreciable selectivity. To the best of our knowledge, this is the first time that the trend in the selectivity in the ion current is produced by a computer simulation approach. Nevertheless, the calculated selectivity is still smaller than its experimental estimate. Recognizing that the calculated profiles are not perfect, we examine how changes in these profiles can account for the observed selectivity. It is found that the origin of the selectivity is more complex than generally assumed. The observed selectivity can be reproduced by increasing the barrier at the exit and the entrance of the selectivity filter, but the necessary changes in the barrier approach the limit of the error in the PDLD/S-LRA calculations. Other options that can increase the selectivity are also considered, including the difference between the Na+...Na+ and K+...K+ interaction. However, this interesting effect does not appear to lead to a major difference in selectivity since the Na+ ions at the limit of strong interaction tend to move in a less concerted way than the K+ ions. Changes in the relative binding energies at the different binding sites are also not so effective in changing the selectivity. Finally, it is pointed out that using the calculated profiles as a starting point and forcing the model to satisfy different experimentally based constraints, should eventually provide more detailed understanding of the different complex factors involved in ion selectivity of biological channels.     

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.     

KcsA: it's a potassium channel

Ion conduction and selectivity properties of KcsA, a bacterial ion channel of known structure, were studied in a planar lipid bilayer system at the single-channel level. Selectivity sequences for permeant ions were determined by symmetrical solution conductance (K(+) > Rb(+), NH(4)(+), Tl(+) > Cs(+), Na(+), Li(+)) and by reversal potentials under bi-ionic or mixed-ion conditions (Tl(+) > K(+) > Rb(+) > NH(4)(+) > Na(+), Li(+)). Determination of reversal potentials with submillivolt accuracy shows that K(+) is over 150-fold more permeant than Na(+). Variation of conductance with concentration under symmetrical salt conditions is complex, with at least two ion-binding processes revealing themselves: a high affinity process below 20 mM and a low affinity process over the range 100-1,000 mM. These properties are analogous to those seen in many eukaryotic K(+) channels, and they establish KcsA as a faithful structural model for ion permeation in eukaryotic K(+) channels.     

Ion solvation and structural stability in a sodium channel investigated by molecular dynamics calculations

The stability and ion binding properties of the homo-tetrameric pore domain of a prokaryotic, voltage-gated sodium channel are studied by extensive all-atom molecular dynamics simulations, with the channel protein being embedded in a fully hydrated lipid bilayer. It is found that Na(+) ion presents in a mostly hydrated state inside the wide pore of the selectivity filter of the sodium channel, in sharp contrast to the nearly fully dehydrated state for K(+) ions in potassium channels. Our results also indicate that Na(+) ions make contact with only one or two out of the four polypeptide chains forming the selectivity filter, and surprisingly, the selectivity filter exhibits robust stability for various initial ion configurations even in the absence of ions. These findings are quite different from those in potassium channels. Furthermore, an electric field above 0.5V/nm is suggested to be able to induce Na(+) permeation through the selectivity filter.     

Ion-Ion interactions at the selectivity filter. Evidence from K(+)-dependent modulation of tetraethylammonium efficacy in Kv2.1 potassium channels

In the Kv2.1 potassium channel, binding of K(+) to a high-affinity site associated with the selectivity filter modulates channel sensitivity to external TEA. In channels carrying Na(+) current, K(+) interacts with the TEA modulation site at concentrations 相似文献