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.     

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.     

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.     

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.     

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.     

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

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

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.     

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.     

Interactions of anionic phospholipids and phosphatidylethanolamine with the potassium channel KcsA

Fluorescence quenching methods have been used to study interactions of anionic phospholipids with the potassium channel KcsA from Streptomyces lividans. Quenching of the Trp fluorescence of KcsA reconstituted into mixtures of dioleoylphosphatidylcholine (DOPC) and an anionic phospholipid with dibromostearoyl chains is more marked at low mole fractions of the brominated anionic phospholipid than is quenching in mixtures of dibromostearoylphosphatidylcholine and nonbrominated anionic lipid. The quenching data are consistent with two classes of binding site for lipid on KcsA, one set corresponding to annular binding sites around KcsA to which DOPC and two-chain anionic phospholipids bind with similar affinities, the other set (non-annular sites) corresponding to sites at which anionic phospholipids can bind but from which DOPC is either excluded or binds with very low affinity. The binding constant for tetraoleoylcardiolipin at the annular sites is significantly less than that for DOPC, being comparable to that for dioleoylphosphatidylethanolamine. Tetraoleoylcardiolipin binds with highest affinity to the non-annular sites, the affinity for dioleoylphosphatidylglycerol being the lowest. The affinity for dioleoylphosphatidylserine decreases at high ionic strength, suggesting that electrostatic interactions between the anionic phospholipid headgroup and positively charged residues on KcsA are important for binding at the non-annular site. The effect of ionic strength on the binding of phosphatidic acid is less marked than on phosphatidylserine. The value of the binding constant for the non-annular site depends on the extent of Trp fluorescence quenching following from binding at the non-annular site. It is suggested that the non-annular site to which binding is detected in the fluorescence quenching experiments corresponds to the binding site for phosphatidylglycerol detected at monomer-monomer interfaces in x-ray diffraction studies.     

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.     

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.     

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.     

Influence of C-terminal protein domains and protein-lipid interactions on tetramerization and stability of the potassium channel KcsA

KcsA is a prokaryotic potassium channel formed by the assembly of four identical subunits around a central aqueous pore. Although the high-resolution X-ray structure of the transmembrane portion of KcsA is known [Doyle, D. A., Morais, C. J., Pfuetzner, R. A., Kuo, A., Gulbis, J. M., Cohen, S. L., Chait, B. T., and MacKinnon, R. (1998) Science 280, 69-77], the identification of the molecular determinant(s) involved in promoting subunit tetramerization remains to be determined. Here, C-terminal deletion channel mutants, KcsA Delta125-160 and Delta120-160, as well as 1-125 KcsA obtained from chymotrypsin cleavage of full-length 1-160 KcsA, have been used to evaluate the role of the C-terminal segment on the stability and tetrameric assembly of the channel protein. We found that the lack of the cytoplasmic C-terminal domain of KcsA, and most critically the 120-124 sequence stretch, impairs tetrameric assembly of channel subunits in a heterologous E. coli expression system. Molecular modeling of KcsA predicts that, indeed, such sequence stretch provides intersubunit interaction sites by hydrogen bonding to amino acid residues in N- and C-terminal segments of adjacent subunits. However, once the KcsA tetramer is assembled, its remarkable in vitro stability to detergent or to heat-induced dissociation into subunits is not greatly influenced by whether the entire C-terminal domain continues being part of the protein. Finally and most interestingly, it is observed that, even in the absence of the C-terminal domain involved in tetramerization, reconstitution into membrane lipids promotes in vitro KcsA tetramerization very efficiently, an event which is likely mediated by allowing proper hydrophobic interactions involving intramembrane protein domains.     

Ion binding affinity in the cavity of the KcsA potassium channel

The hydrophobic cell membrane interior presents a large energy barrier for ions to permeate. Potassium channels reduce this barrier by creating a water-filled cavity at the middle of their ion conduction pore to allow ion hydration and by directing the C-terminal "end charge" of four alpha-helices toward the water-filled cavity. Here we have studied the interaction of monovalent cations with the cavity of the KcsA K(+) channel using X-ray crystallography. In these studies, Tl(+) was used as an analogue for K(+) and the total ion-stabilization energy for Tl(+) in the cavity was estimated by measuring its binding affinity. Binding affinity for the Na(+) ion was also measured, revealing a weak selectivity ( approximately 7-fold) favoring Tl(+) over Na(+). The structures of the cavity containing Na(+), K(+), Tl(+), Rb(+), and Cs(+) are compared. These results are consistent with a fairly large (more negative than -100 mV) electrostatic potential inside the cavity, and they also imply the presence of a weak nonelectrostatic component to a cation's interaction with the cavity.     

NMR study of the tetrameric KcsA potassium channel in detergent micelles

Nuclear magnetic resonance (NMR) studies of large membrane-associated proteins are limited by the difficulties in preparation of stable protein-detergent mixed micelles and by line broadening, which is typical of these macroassemblies. We have used the 68-kDa homotetrameric KcsA, a thermostable N-terminal deletion mutant of a bacterial potassium channel from Streptomyces lividans, as a model system for applying NMR methods to membrane proteins. Optimization of measurement conditions enabled us to perform the backbone assignment of KcsA in SDS micelles and establish its secondary structure, which was found to closely agree with the KcsA crystal structure. The C-terminal cytoplasmic domain, absent in the original structure, contains a 14-residue helix that could participate in tetramerization by forming an intersubunit four-helix bundle. A quantitative estimate of cross- relaxation between detergent and KcsA backbone amide protons, together with relaxation and light scattering data, suggests SDS-KcsA mixed micelles form an oblate spheroid with approximately 180 SDS molecules per channel. K(+) ions bind to the micelle-solubilized channel with a K(D) of 3 +/- 0.5 mM, resulting in chemical shift changes in the selectivity filter. Related pH-induced changes in chemical shift along the "outer" transmembrane helix and the cytoplasmic membrane interface hint at a possible structural explanation for the observed pH-gating of the potassium channel.     

KcsA closed and open: modelling and simulation studies

Bacterial homologues of mammalian potassium channels provide structures of two states of a gated K channel. Thus, the crystal structure of KcsA represents a closed state whilst that of MthK represents an open state. Using homology modelling and molecular dynamics simulations we have built a model of the transmembrane domain of KcsA in an open state and have compared its conformational stability with that of the same domain of KcsA in a closed state. Approximate Born energy calculations of monovalent cations within the two KcsA channel states suggest that the intracellular hydrophobic gate in the closed state provides a barrier of height ~5 kT to ion permeation, whilst in the open state the barrier is absent. Simulations (10 ns duration) in an octane slab (a simple membrane mimetic) suggest that closed- and open-state models are of comparable conformational stability, both exhibiting conformational drifts of ~3.3 Å C RMSD relative to the respective starting models. Substantial conformational fluctuations are observed in the intracellular gate region during both simulations (closed state and open state). In the simulation of open-state KcsA, rapid (<5 ns) exit of all three K⁺ ions occurs through the intracellular mouth of the channel. Helix kink and swivel motion is observed at the molecular hinge formed by residue G99 of the M2 helix. This motion is more substantial for the open- than for the closed-state model of the channel.     

Multiple binding sites for fatty acids on the potassium channel KcsA

Interactions of fatty acids with the potassium channel KcsA were studied using Trp fluorescence quenching and electron paramagnetic resonance (EPR) techniques. The brominated analogue of oleic acid was shown to bind to annular sites on KcsA and to the nonannular sites at each protein-protein interface in the homotetrameric structure with binding constants relative to dioleoylphosphatidylcholine of 0.67 ± 0.04 and 0.87 ± 0.08, respectively. Mutation of the two Arg residues close to the nonannular binding sites had no effect on fatty acid binding. EPR studies with a spin-labeled analogue of stearic acid detected a high-affinity binding site for the fatty acid with strong immobilization. Fluorescence quenching studies with the spin-labeled analogue showed that the binding site detected in the EPR experiments could not be one of the annular or nonannular binding sites. Instead, it is proposed that the EPR studies detect binding to the central hydrophobic cavity of the channel, with a binding constant in the range of ~0.1-1 μM.     

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.