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.     

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.     

Simulations of ion current in realistic models of ion channels: the KcsA potassium channel

Realistic studies of ion current in biologic channels present a major challenge for computer simulation approaches. All-atom molecular dynamics simulations involve serious time limitations that prevent their use in direct evaluation of ion current in channels with significant barriers. The alternative use of Brownian dynamics (BD) simulations can provide the current for simplified macroscopic models. However, the time needed for accurate calculations of electrostatic energies can make BD simulations of ion current expensive. The present work develops an approach that overcomes some of the above challenges and allows one to simulate ion currents in models of biologic channels. Our method provides a fast and reliable estimate of the energetics of the system by combining semimacroscopic calculations of the self-energy of each ion and an implicit treatment of the interactions between the ions, as well as the interactions between the ions and the protein-ionizable groups. This treatment involves the use of the semimacroscopic version of the protein dipole Langevin dipole (PDLD/S) model in its linear response approximation (LRA) implementation, which reduces the uncertainties about the value of the protein "dielectric constant." The resulting free energy surface is used to generate the forces for on-the-fly BD simulations of the corresponding ion currents. Our model is examined in a preliminary simulation of the ion current in the KcsA potassium channel. The complete free energy profile for a single ion transport reflects reasonable energetics and captures the effect of the protein-ionized groups. This calculated profile indicates that we are dealing with the channel in its closed state. Reducing the barrier at the gate region allows us to simulate the ion current in a reasonable computational time. Several limiting cases are examined, including those that reproduce the observed current, and the nature of the productive trajectories is considered. The ability to simulate the current in realistic models of ion channels should provide a powerful tool for studies of the biologic function of such systems, including the analysis of the effect of mutations, pH, and electric potentials.     

Measurement of 15N relaxation in the detergent-solubilized tetrameric KcsA potassium channel

A set of TROSY-HNCO (tHNCO)-based 3D experiments is presented for measuring ¹⁵N relaxation parameters in large, membrane-associated proteins, characterized by slow tumbling times and significant spectral overlap. Measurement of backbone ¹⁵N R ₁, R _1ρ, ¹⁵N–{¹H} NOE, and ¹⁵N CSA/dipolar cross correlation is demonstrated and applied to study the dynamic behavior of the homotetrameric KcsA potassium channel in SDS micelles under conditions where this channel is in the closed state. The micelle-encapsulated transmembrane domain, KcsA^TM, exhibits a high degree of order, tumbling as an oblate ellipsoid with a global rotational correlation time, τ_c = 38 ± 2.5 ns, at 50 °C and a diffusion anisotropy, , corresponding to an aspect ratio a/b ≥ 1.4. The N- and C-terminal intracellular segments of KcsA exhibit considerable internal dynamics (S ² values in the 0.2–0.45 range), but are distinctly more ordered than what has been observed for unstructured random coils. Relaxation behavior in these domains confirms the position of the C-terminal helix, and indicates that in SDS micelles, this amphiphilic helix does not associate into a stable homotetrameric helical bundle. The relaxation data indicate the absence of elevated backbone dynamics on the ps–ns time scale for the 5-residue selectivity filter, which selects K⁺ ions to enter the channel. Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at . An erratum to this article can be found at     

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.     

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.     

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.     

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.     

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.     

Anionic phospholipid interactions with the potassium channel KcsA: simulation studies

Molecular dynamics (MD) simulations have been used to unmask details of specific interactions of anionic phospholipids with intersubunit binding sites on the surface of the bacterial potassium channel KcsA. Crystallographic data on a diacyl glycerol fragment at this site were used to model phosphatidylethanolamine (PE), or phosphatidylglycerol (PG), or phosphatidic acid (PA) at the intersubunit binding sites. Each of these models of a KcsA-lipid complex was embedded in phosphatidyl choline bilayer and explored in a 20 ns MD simulation. H-bond analysis revealed that in terms of lipid-protein interactions PA > PG > PE and revealed how anionic lipids (PG and PA) bind to a site provided by two key arginine residues (R(64) and R(89)) at the interface between adjacent subunits. A 27 ns simulation was performed in which KcsA (without any lipids initially modeled at the R(64)/R(89) sites) was embedded in a PE/PG bilayer. There was a progressive specific increase over the course of the simulation in the number of H-bonds of PG with KcsA. Furthermore, two specific PG binding events at R(64)/R(89) sites were observed. The phosphate oxygen atoms of bound PG formed H-bonds to the guanidinium group of R(89), whereas the terminal glycerol H-bonded to R(64). Overall, this study suggests that simulations can help identify and characterize sites for specific lipid interactions on a membrane protein surface.     

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.     

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.     

Energy distribution and ion selectivity of the bacterial potassium channel

The ion selectivity of the bacterial potassium channel K_CSA is explained upon comparing the energy characteristics of the interaction of cations (Li⁺, Na⁺, K⁺) with atoms of the selectivity filter of the protein pore. Quantum-chemical calculations reveal a deeper potential well for potassium ions, which accounts for preferred K⁺ permeation. It is shown that the conventional methods with AMBER, CHARMM, OPLS force fields in standard parametrization as well as partial re-parametrization give incorrect estimates of ion energy distribution in the channel.     

Permeation mechanism of a two-state potassium channel

A two-state hopping model was proposed to study the permeation of ion channel.The Nemst equation in equilibrium and the Michaelis-Menten relation in steady state were derived from the two-state kinetic model.The currentvoltage relationship obtained in the symmetrical solutions case was linear when the applied potential was less than 100 mV,which met Ohm's law.The conductance-concentration relationship exhibited the saturation property.Moreover,the characteristic time reaching the steady state of the KcsA channel was also discussed.     

Permeation mechanism of a two-state potassium channel

A two-state hopping model was proposed to study the permeation of ion channel. The Nernst equation in equilibrium and the Michaelis-Menten relation in steady state were derived from the two-state kinetic model. The current-voltage relationship obtained in the symmetrical solutions case was linear when the applied potential was less than 100 mV, which met Ohm’s law. The conductance-concentration relationship exhibited the saturation property. Moreover, the characteristic time reaching the steady state of the KcsA channel was also discussed. Translated from Acta Biophysica Sinica, 2005, 21(4): 289–294 [译自: 生物物理学报]     

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.     

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.     

Mechanisms of tetraethylammonium ion block in the KcsA potassium channel

We report results from automated docking and microscopic molecular dynamics simulations of the tetraethylammonium (TEA) complexes with KcsA. Binding modes and energies for TEA binding at the external and internal sides of the channel pore are examined utilising the linear interaction energy method. Effects of the channel ion occupancy (based on our previous results for the ion permeation mechanisms) on the binding energies are considered. Calculations show that TEA forms stable complexes at both the external and internal entrances of the selectivity filter. Furthermore, the effects of the Y82V mutation are evaluated and the results show, in agreement with experimental data, that the mutant has a significantly reduced binding affinity for TEA at the external binding site, which is attributed to stabilising hydrophobic interactions between the ligand and the tyrosines.     

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.     

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.