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.     

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.     

K(+)/Na(+) selectivity of the KcsA potassium channel from microscopic free energy perturbation calculations

Microscopic molecular dynamics free energy perturbation calculations of the K(+)/Na(+) selectivity in the KcsA potassium channel, based on its experimental three-dimensional structure, are reported. The relative binding free energies for K(+) and Na(+) in the most relevant ion occupancy states of the four-site selectivity filter are calculated. The previously proposed mechanism for ion permeation through the KcsA channel is predicted, in agreement with available experimental data, to have a significant selectivity for K(+) over Na(+). The calculations also show that the individual 'binding site' selectivities are generally not additive and the doubly loaded states of the filter thus display cooperative effects. The only site that is not K(+) selective is that which is located at the entrance to the internal water cavity, suggesting the possibility that internal Na(+) could block outward currents.     

Brownian dynamics study of an open-state KcsA potassium channel.

Three-dimensional Brownian dynamics simulations are used to study conductance of the KcsA potassium channel using the known crystallographic structure. Employing an open-state channel created by molecular dynamics simulations, current-voltage and current-concentration curves broadly consistent with experimental measurements are obtained. In the absence of an applied potential, the channel houses three potassium ions at positions that are in close agreement with X-ray diffraction maps.     

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.     

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.     

Modelling the effect of a GHz electric field on the dynamics of K+ ions in KcsA potassium channel

The dynamics of potassium ions in a KcsA channel, located within a stochastically fluctuating medium, is modelled via the application of the molecular dynamics simulation method. We investigate the effect of presence and absence of an applied electric field, either constant or periodic, on the dynamics of the channel. It is found that the ions undergo a hopping motion when the channel is exposed to a constant electric field of strength 0.03 V/nm. Furthermore, an alternating electric field in the GHz range, normally present in the daily environment and encountered by most biological systems, is applied to the channel, showing that in this frequency range, the rigidity of the atomic bonds of the filter is increased, which in turn disturbs the ionic passage rate through the filter. Consequently, in this frequency range, the application of electric fields may affect the function of such channels.     

Conducting-state properties of the KcsA potassium channel from molecular and Brownian dynamics simulations.

The mechanisms underlying transport of ions across the potassium channel are examined using electrostatic calculations and three-dimensional Brownian dynamics simulations. We first build open-state configurations of the channel with molecular dynamics simulations, by pulling the transmembrane helices outward until the channel attains the desired interior radius. To gain insights into ion permeation, we construct potential energy profiles experienced by an ion traversing the channel in the presence of other resident ions. These profiles reveal that in the absence of an applied field the channel accommodates three potassium ions in a stable equilibrium, two in the selectivity filter and one in the central cavity. In the presence of a driving potential, this three-ion state becomes unstable, and ion permeation across the channel is observed. These qualitative explanations are confirmed by the results of three-dimensional Brownian dynamics simulations. We find that the channel conducts when the ionizable residues near the extracellular entrance are fully charged and those near the intracellular side are partially charged. The conductance increases steeply as the radius of the intracellular mouth of the channel is increased from 2 A to 5 A. Our simulation results reproduce several experimental observations, including the current-voltage curves, conductance-concentration relationships, and outward rectification of currents.     

The ionization state and the conformation of Glu-71 in the KcsA K(+) channel.

The side chain of Glu-71 of the KcsA K(+) channel, an important residue in the vicinity of the selectivity filter, was not resolved in the crystallographic structure of Doyle et al. (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). Its atomic coordinates are undetermined and its ionization state is unknown. For meaningful theoretical and computational studies of the KcsA K(+) channel, it is essential to address questions about the conformation and the ionization state of this residue in detail. In previous MD simulations in which the side chain of Glu-71 is protonated and forming a strong hydrogen bond with Asp-80 it was observed that the channel did not deviate significantly from the crystallographic structure (Bernèche, S., and B. Roux. 2000. Biophys. J. 78:2900-2917). In contrast, we show here that the structure of the selectivity filter of the KcsA channel is significantly disrupted when these side chains are fully ionized on each of the four monomers. To further resolve questions about the ionization state of Glu-71 we calculated the pK(a) value of this residue using molecular dynamics free energy simulations (MD/FES) with a fully flexible system including explicit solvent and membrane and finite-difference Poisson-Boltzmann (PB) continuum electrostatics. It is found that the pK(a) of Glu-71 is shifted by approximately +10 pK(a) units. These results strongly suggest that Glu-71 is protonated under normal conditions.     

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.     

Interactions of phospholipids with the potassium channel KcsA

The potassium channel KcsA from Streptomyces lividans has been reconstituted into bilayers of phosphatidylcholines and fluorescence spectroscopy has been used to characterize the response of KcsA to changes in bilayer thickness. The Trp residues in KcsA form two bands, one on each side of the membrane. Trp fluorescence emission spectra and the proportion of the Trp fluorescence intensity quenchable by I(-) hardly vary in the lipid chain length range C10 to C24, suggesting efficient hydrophobic matching between KcsA and the lipid bilayer over this range. Measurements of fluorescence quenching for KcsA reconstituted into mixtures of brominated and nonbrominated phospholipids have been analyzed to give binding constants of lipids for KcsA, relative to that for dioleoylphosphatidylcholine (di(C18:1)PC). Relative lipid binding constants increase by only a factor of three with increasing chain length from C10 to C22 with a decrease from C22 to C24. Strongest binding to di(C22:1)PC corresponds to a state in which the side chains of the lipid-exposed Trp residues are likely to be located within the hydrocarbon core of the lipid bilayer. It is suggested that matching of KcsA to thinner bilayers than di(C24:1)PC is achieved by tilting of the transmembrane alpha-helices in KcsA. Measurements of fluorescence quenching of KcsA in bilayers of brominated phospholipids as a function of phospholipid chain length suggest that in the chain length range C14 to C18 the Trp residues move further away from the center of the lipid bilayer with increasing chain length, which can be partly explained by a decrease in helix tilt angle with increasing bilayer thickness. In the chain length range C18 to C24 it is suggested that the Trp residues become more buried within the hydrocarbon core of the bilayer.     

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.     

A mutant KcsA K(+) channel with altered conduction properties and selectivity filter ion distribution

The selectivity filter of K(+) channels is comprised of a linear queue of four equal-spaced ion-binding sites spanning a distance of 12A. Each site is formed of eight oxygen atoms from the protein. The first three sites, numbered 1-3 from the extracellular side, are made of exclusively main-chain carbonyl oxygen atoms. The fourth site, closest to the intracellular side, is made of four main-chain carbonyl oxygen atoms and four threonine side-chain hydroxyl oxygen atoms. Here we characterize the effects of mutating the threonine to cysteine on the distribution of ions in the selectivity filter and on the conduction of ions through the filter. The mutation influences the occupancy of K(+) at sites 2 and 4 and it reduces the maximum rate of conduction in the limit of high K(+) concentration. The mutation does not affect the conduction of Rb(+). These results can be understood in the context of a conduction mechanism in which a pair of K ions switch between energetically balanced 1,3 and 2,4 configurations.     

Rapid intracellular TEA block of the KcsA potassium channel

Intracellular tetraethylammonium (TEA) inhibition was studied at the single-channel level in the KcsA potassium channel reconstituted in planar lipid bilayers. TEA acts as a fast blocker (resulting in decreased current amplitude) with an affinity in the 75 mM range even at high bandwidth. Studies over a wide voltage range reveal that TEA block has a complex voltage-dependence that also depends on the ionic conditions. These observations are examined in the context of permeation models to extend our understanding of the coupling between permeant ions and TEA blockade.     

Intermolecular functional coupling between phosphoinositides and the potassium channel KcsA

Biological membranes are composed of a wide variety of lipids. Phosphoinositides (PIPns) in the membrane inner leaflet only account for a small percentage of the total membrane lipids but modulate the functions of various membrane proteins, including ion channels, which play important roles in cell signaling. KcsA, a prototypical K⁺ channel that is small, simple, and easy to handle, has been broadly examined regarding its crystallography, in silico molecular analysis, and electrophysiology. It has been reported that KcsA activity is regulated by membrane phospholipids, such as phosphatidylglycerol. However, there has been no quantitative analysis of the correlation between direct lipid binding and the functional modification of KcsA, and it is unknown whether PIPns modulate KcsA function. Here, using contact bubble bilayer recording, we observed that the open probability of KcsA increased significantly (from about 10% to 90%) when the membrane inner leaflet contained only a small percentage of PIPns. In addition, we found an increase in the electrophysiological activity of KcsA correlated with a larger number of negative charges on PIPns. We further analyzed the affinity of the direct interaction between PIPns and KcsA using microscale thermophoresis and observed a strong correlation between direct lipid binding and the functional modification of KcsA. In conclusion, our approach was able to reconstruct the direct modification of KcsA by PIPns, and we propose that it can also be applied to elucidate the mechanism of modification of other ion channels by PIPns.     

Brownian dynamics theory for predicting internal and external blockages of tetraethylammonium in the KcsA potassium channel

The theory of Brownian dynamics is used to model permeation and the blocking of KcsA potassium channels by tetraethylammonium (TEA). A novel Brownian dynamics simulation algorithm is implemented that comprises two free energy profiles; one profile is seen by the potassium ions and the other by the TEA molecules whose shape is approximated by a sphere. Our simulations reveal that internally applied TEA blocks the passage of K⁺ ions by physically occluding the pore. A TEA molecule in the external reservoir encounters an attractive energy-well created by four tyrosine residues at position 82, in addition to all other attractive and repulsive forces impinging on it. Using Brownian dynamics, we investigate how deep the energy-well needs to be to reproduce the experimentally determined inhibitory constant k_i for the TEA blockade of KcsA or the mutant Shaker T449Y. The one-dimensional free energy profile obtained from molecular dynamics is first converted into a one-dimensional potential energy profile, and is then transformed into a three-dimensional free energy profile in Brownian dynamics by adding the short-range potential from the channel walls. When converted, the free energy profile calculated from molecular dynamics gives a well-depth of ∼10 kT. We systematically alter the depths of the profiles, and then use Brownian dynamics simulations to numerically determine the current versus TEA-concentration curves. We show that the sequence of binding and unbinding events of the TEA molecule to the binding pocket can be modeled by a first-order Markov process. The Brownian dynamics simulations also reveal that the probability of a TEA molecule binding to the binding pocket in KcsA potassium channels increases exponentially with TEA concentration and depends also on the applied potential and the K⁺ concentration in the simulation assembly.     

Efficacy of external tetraethylammonium block of the KcsA potassium channel: Molecular and Brownian dynamics studies

Blockade of the KcsA potassium channel by externally applied tetraethylammonium is investigated using molecular dynamics calculations and Brownian dynamics simulations. In KcsA, the aromatic rings of four tyrosine residues located just external to the selectivity filter create an attractive energy well or a binding cage for a tetraethylammonium molecule. We first investigate the effects of re-orienting the four tyrosine residues such that the centers of the aromatic rings face the tetraethylammonium molecule directly. Then, we systematically move the residues inward in both orientations so that the radius of the binding cage formed by them becomes smaller. For each configuration, we construct a one-dimensional free energy profile by bringing in a tetraethylammonium molecule from the external reservoir toward the selectivity filter. The free energy profile is then converted to a one-dimensional potential energy profile, taking the available space between the tyrosine residues and the tetraethylammonium molecule into account. Incorporating this potential energy profile into the Brownian dynamics algorithm, we determine the conductance properties of the channel under various conditions, construct the current-tetraethylammonium-concentration curve and compare it with the experimentally determined inhibitory constant k_i for externally applied tetraethylammonium. We show that the experimentally determined binding affinity for externally applied tetraethylammonium can be replicated when each of the four tyrosine residues is moved inward by about 0.7 Å, irrespective of orientation of their aromatic rings.     

Inactivation of the KcsA potassium channel explored with heterotetramers

The tetrameric prokaryotic potassium channel KcsA is activated by protons acting on the intracellular aspect of the protein and inactivated through conformational changes in the selectivity filter. Inactivation is modulated by a network of interactions within each protomer between the pore helix and residues at the external entrance of the channel. Inactivation is suppressed by the E71A mutation, which perturbs the stability of this network. Here, cell-free protein synthesis followed by protein purification by sodium dodecyl sulfate–polyacrylamide gel electrophoresis was used to produce heterotetramers of KcsA that contain different combinations of wild-type and E71A subunits. Single-channel recordings from these heterotetramers reveal how the network of interactions in individual protomers affects ionic conduction and channel inactivation, suggesting that the latter is a cooperative process.     

Mixed modes in opening of KcsA potassium channel from a targeted molecular dynamics simulation

Potassium channels conduct K⁺ flow selectively across the membrane through a central pore. During a process called gating, the potassium channels undergo a conformational change that opens or closes the ion-conducting pore. The potassium channel KcsA has been structurally determined in its closed state. However, the dynamic mechanism of the gating transition of the KcsA channel is still being investigated. Here, a targeted molecular dynamics simulation up to 150 ns is performed to investigate the detailed opening process of the KcsA channel with an open Kv1.2 structure serving as the target. The channel arrived at a self-determined quasi-stable state within 60 ns. The rigid-body and hinge-bending modes are observed mixed together in the remaining 90 ns long quasi-stable state. The mixed-mode movement seems come from the competition between the helix rigidity and the biased-applied gating force.