Voltage-dependent gating at the KcsA selectivity filter

The prokaryotic K(+) channel KcsA, although lacking a 'standard' voltage-sensing domain, shows voltage-dependent gating that leads to an increase in steady-state open probability of almost two orders of magnitude between +150 and -150 mV. Here we show that voltage-dependent gating in KcsA is associated with the movement of approximately 0.7 equivalent electronic charges. This charge movement produces an increase in the rate of entry into a long-lived inactivated state and seems to be independent of the proton-activation mechanism. Charge neutralization at position 71 renders the channel essentially voltage-independent by preventing entry into the inactivated state. A mechanism for voltage-dependent gating at the selectivity filter is proposed that is based on the reorientation of the carboxylic moiety of Glu71 and its influence in the conformational dynamics of the selectivity filter.     

Molecular dynamics estimates of ion diffusion in model hydrophobic and KcsA potassium channels

Molecular dynamics simulations are carried out to obtain estimates of diffusion coefficients of biologically important Na+, K+, Ca2+ and Cl- ions in hydrophobic cylindrical channels with varying radii and large reservoirs. Calculations for the cylindrical channels are compared to those for the KcsA potassium channel, for which the protein structure has recently been determined from X-ray diffraction experiments. Our results show that ion diffusion is maintained at reasonably high levels even within narrow channels, and does not support the very small diffusion coefficients used in some continuum models in order to fit experimental data. The present estimates of ion diffusion coefficients are useful in the calculation of channel conductance using the Poisson-Nernst-Planck theory, or Brownian dynamics.     

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.     

Mechanism of K+/Na+ selectivity in potassium channels from the perspective of the non-selective bacterial channel NaK

Studies of the non-selective NaK channel as a model system provided detailed structural and functional insights into different factors that modulate selective ion conductance. The combination of structural and functional studies of WT- and mutant forms of NaK channel show that apparent conservation of the multiple binding sites in the canonical selectivity filter of K-channel is essential component of the observed robust Na (+) discrimination.     

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.     

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.     

Molecular dynamics study of the KcsA potassium channel

The structural, dynamical, and thermodynamic properties of a model potassium channel are studied using molecular dynamics simulations. We use the recently unveiled protein structure for the KcsA potassium channel from Streptomyces lividans. Total and free energy profiles of potassium and sodium ions reveal a considerable preference for the larger potassium ions. The selectivity of the channel arises from its ability to completely solvate the potassium ions, but not the smaller sodium ions. Self-diffusion of water within the narrow selectivity filter is found to be reduced by an order of magnitude from bulk levels, whereas the wider hydrophobic section of the pore maintains near-bulk self-diffusion. Simulations examining multiple ion configurations suggest a two-ion channel. Ion diffusion is found to be reduced to approximately (1)/(3) of bulk diffusion within the selectivity filter. The reduced ion mobility does not hinder the passage of ions, as permeation appears to be driven by Coulomb repulsion within this multiple ion channel.     

Mechanisms of selectivity in channels and enzymes studied with interactive molecular dynamics

Interactive molecular dynamics, a new modeling tool for rapid investigation of the physical mechanisms of biological processes at the atomic level, is applied to study selectivity and regulation of the membrane channel protein GlpF and the enzyme glycerol kinase. These proteins facilitate the first two steps of Escherichia coli glycerol metabolism. Despite their different function and architecture the proteins are found to employ common mechanisms for substrate selectivity: an induced geometrical fit by structurally homologous binding sites and an induced rapid dipole moment reversal. Competition for hydrogen bonding sites with water in both proteins is critical for substrate motion. In glycerol kinase, it is shown that the proposed domain motion prevents competition with water, in turn regulating the binding of glycerol.     

Extracellular blockade of K(+) channels by TEA: results from molecular dynamics simulations of the KcsA channel

TEA is a classical blocker of K(+) channels. From mutagenesis studies, it has been shown that external blockade by TEA is strongly dependent upon the presence of aromatic residue at Shaker position 449 which is located near the extracellular entrance to the pore (Heginbotham, L., and R. MacKinnon. 1992. Neuron. 8:483-491). The data suggest that TEA interacts simultaneously with the aromatic residues of the four monomers. The determination of the 3-D structure of the KcsA channel using X-ray crystallography (Doyle, D.A., J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, and R. MacKinnon. 1998. Science. 280:69-77) has raised some issues that remain currently unresolved concerning the interpretation of these observations. In particular, the center of the Tyr82 side chains in KcsA (corresponding to position 449 in Shaker) forms a square of 11.8-A side, a distance which is too large to allow simultaneous interactions of a TEA molecule with the four aromatic side chains. In this paper, the external blockade by TEA is explored by molecular dynamics simulations of an atomic model of KcsA in an explicit phospholipid bilayer with aqueous salt solution. It is observed, in qualitative accord with the experimental results, that TEA is stable when bound to the external side of the wild-type KcsA channel (with Tyr82), but is unstable when bound to a mutant channel in which the tyrosine residue has been substituted by a threonine. The free energy profile of TEA relative to the pore is calculated using umbrella sampling simulations to characterize quantitatively the extracellular blockade. It is found, in remarkable agreement with the experiment, that the TEA is more stably bound by 2.3 kcal/mol to the channel with four tyrosine residues. In the case of the wild-type KcsA channel, TEA (which has the shape of a flattened oblate spheroid) acts as an ideal plug blocking the pore. In contrast, it is considerably more off-centered and tilted in the case of the mutant channel. The enhanced stability conferred by the tyrosine residues does not arise from Pi-cation interactions, but appears to be due to differences in the hydration structure of the TEA. Finally, it is shown that the experimentally observed voltage dependence of TEA block, which is traditionally interpreted in terms of the physical position of the TEA along the axis of the pore, must arise indirectly via coupling with the ions in the pore.     

Molecular dynamics simulations and KcsA channel gating

The gating mechanism of a bacterial potassium channel, KcsA, has been investigated via multi-nanosecond molecular dynamic simulations of the channel molecules embedded in a fully solvated palmitoyloleoylphosphatidylcholine bilayer. Four events are seen in which a cation (K(+) or, in one case, Na(+)) initially present in the central cavity exits through the intracellular mouth (the presumed gate) of the channel. Whilst in the cavity a cation interacts with the sidechain T107 O gamma atom of one of the subunits prior to its exit from the channel. Secondary structure analysis as a function of time reveals a break in the helicity of one of the M2 helices. This break is expected to lend flexibility to the helices, enabling them to "open" (minimum pore radius >0.13 nm) and "close" (minimum pore radius <0.13 nm) the channel. Fluctuations in the pore radius at the intracellular gate region are of the order of 0.05 nm, with an average radius in the region of the gate of ca. 0.1 nm. However, around the time of exit of a cation, the pore widens to about 0.15 nm. The distances between the C alpha atoms of the inner helices M2 reveal a coupled increase and decrease between the opposite pair of helices at about the time of exit of the ion. This suggests a breathing motion of the M2 helices that may form the basis for a gating mechanism.     

Hydropathic analysis and comparison of KcsA and Shaker potassium channels

The similarity in structure of potassium (K(+)) channels from different families has been revealed by only recently available crystallographic 3D structural data. The hydropathic analysis presented in this work illuminates whether homologous residues perform the same functions in channels that use different gating mechanisms. We calculated and compared the hydropathic profiles of two K(+) channels, KcsA and Kv1.2 (the latter a member of the Shaker family), at their pore-forming domain. Quantitative information describing important interactions stabilizing the protein beyond obvious secondary-structure elements was extracted from the analysis and applied as a template for subsequent molecular-dynamics (MD) analyses. For example, two key groups of interactions, defining the turns that connect the transmembrane helices and responsible for the orientation of the pore helix, were identified. Our results also indicate that Asp(80) and Asp(379) play a similar role in stabilizing the P-loop of KcsA and Kv1.2, respectively, but to significantly different extents.     

Potassium-selective block of barium permeation through single KcsA channels

Ba(2+), a doubly charged analogue of K(+), specifically blocks K(+) channels by virtue of electrostatic stabilization in the permeation pathway. Ba(2+) block is used here as a tool to determine the equilibrium binding affinity for various monovalent cations at specific sites in the selectivity filter of a noninactivating mutant of KcsA. At high concentrations of external K(+), the block-time distribution is double exponential, marking at least two Ba(2+) sites in the selectivity filter, in accord with a Ba(2+)-containing crystal structure of KcsA. By analyzing block as a function of extracellular K(+), we determined the equilibrium dissociation constant of K(+) and of other monovalent cations at an extracellular site, presumably S1, to arrive at a selectivity sequence for binding at this site: Rb(+) (3 μM) > Cs(+) (23 μM) > K(+) (29 μM) > NH(4)(+) (440 μM) > Na(+) and Li(+) (>1 M). This represents an unusually high selectivity for K(+) over Na(+), with |ΔΔG(0)| of at least 7 kcal mol(-1). These results fit well with other kinetic measurements of selectivity as well as with the many crystal structures of KcsA in various ionic conditions.     

Water and potassium dynamics inside the KcsA K(+) channel

Molecular dynamics simulations and electrostatic modeling are used to investigate structural and dynamical properties of the potassium ions and of water molecules inside the KcsA channel immersed in a membrane-mimetic environment. Two potassium ions, initially located in the selectivity filter binding sites, maintain their position during 2 ns of dynamics. A third potassium ion is very mobile in the water-filled cavity. The protein appears engineered so as to polarize water molecules inside the channel cavity. The resulting water induced dipole and the positively charged potassium ion within the cavity are the key ingredients for stabilizing the two K(+) ions in the binding sites. These two ions experience single file movements upon removal of the potassium in the cavity, confirming the role of the latter in ion transport through the channel.     

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.     

Ion hydration in nanopores and the molecular basis of selectivity

Using a simple model, it is shown that the cost of constraining a hydrated potassium ion inside a narrow pore is smaller than the cost of constraining hydrated sodium or lithium ions in pores of radius around 1.5 A. The opposite is true for pores of radius around 2.5 A. The reason for the selectivity in the first region is that the potassium ion allows for a greater distortion of its hydration shell and can therefore maintain a better coordination, and the reason for the reverse selectivity in the second region is that the smaller ions retain their hydration shells in these pores. This is relevant to the molecular basis of ion selective channels, and since this mechanism does not depend on the molecular details of the pore, it could also operate in all sorts of nanotubes.     

Influence of hydration on the dynamics of lysozyme

Quasielastic neutron and light-scattering techniques along with molecular dynamics simulations were employed to study the influence of hydration on the internal dynamics of lysozyme. We identified three major relaxation processes that contribute to the observed dynamics in the picosecond to nanosecond time range: 1), fluctuations of methyl groups; 2), fast picosecond relaxation; and 3), a slow relaxation process. A low-temperature onset of anharmonicity at T approximately 100 K is ascribed to methyl-group dynamics that is not sensitive to hydration level. The increase of hydration level seems to first increase the fast relaxation process and then activate the slow relaxation process at h approximately 0.2. The quasielastic scattering intensity associated with the slow process increases sharply with an increase of hydration to above h approximately 0.2. Activation of the slow process is responsible for the dynamical transition at T approximately 200 K. The dependence of the slow process on hydration correlates with the hydration dependence of the enzymatic activity of lysozyme, whereas the dependence of the fast process seems to correlate with the hydration dependence of hydrogen exchange of lysozyme.