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.     

Role of the KcsA channel cytoplasmic domain in pH-dependent gating

The KcsA channel is a representative potassium channel that is activated by changes in pH. Previous studies suggested that the region that senses pH is entirely within its transmembrane segments. However, we recently revealed that the cytoplasmic domain also has an important role, because its conformation was observed to change dramatically in response to pH changes. Here, to investigate the effects of the cytoplasmic domain on pH-dependent gating, we made a chimera mutant channel consisting of the cytoplasmic domain of the KcsA channel and the transmembrane region of the MthK channel. The chimera showed a pH dependency similar to that of KcsA, indicating that the cytoplasmic domain can act as a pH sensor. To identify how this region detects pH, we substituted certain cytoplasmic domain amino acids that are normally negatively charged at pH 7 for neutral ones in the KcsA channels. These mutants opened independently of pH, suggesting that electrostatic charges have a major role in the cytoplasmic domain's ability to sense and respond to pH.     

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.     

Electrostatic state of the cytoplasmic domain influences inactivation at the selectivity filter of the KcsA potassium channel

KcsA is a proton-activated K⁺ channel that is regulated at two gates: an activation gate located in the inner entrance of the pore and an inactivation gate at the selectivity filter. Previously, we revealed that the cytoplasmic domain (CPD) of KcsA senses proton and that electrostatic changes of the CPD influences the opening and closing of the activation gate. However, our previous studies did not reveal the effect of CPD on the inactivation gate because we used a non-inactivating mutant (E71A). In the present study, we used mutants that did not harbor the E71A mutation, and showed that the electrostatic state of the CPD influences the inactivation gate. Three novel CPD mutants were generated in which some negatively charged amino acids were replaced with neutral amino acids. These CPD mutants conducted K⁺, but showed various inactivation properties. Mutants carrying the D149N mutation showed high open probability and slow inactivation, whereas those without the D149N mutation showed low open probability and fast inactivation, similar to wild-type KcsA. In addition, mutants with D149N showed poor K⁺ selectivity, and permitted Na⁺ to flow. These results indicated that electrostatic changes in the CPD by D149N mutation triggered the loss of fast inactivation and changes in the conformation of selectivity filter. Additionally, the loss of fast inactivation induced by D149N was reversed by R153A mutation, suggesting that not only the electrostatic state of D149, but also that of R153 affects inactivation.     

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.     

Characterization of the C-terminal domain of a potassium channel from Streptomyces lividans (KcsA)

KcsA, a potassium channel from Streptomyces lividans, is a good model for probing the general working mechanism of potassium channels. To date, the physiological activator of KcsA is still unknown, but in vitro studies showed that it could be opened by lowering the pH of the cytoplasmic compartment to 4. The C-terminal domain (CTD, residues 112-160) was proposed to be the modulator for this pH-responsive event. Here, we support this proposal by examining the pH profiles of: (a) thermal stability of KcsA with and without its CTD and (b) aggregation properties of a recombinant fragment of CTD. We found that the presence of the CTD weakened and enhanced the stability of KcsA at acidic and basic pH values, respectively. In addition, the CTD fragment oligomerized at basic pH values with a transition profile close to that of channel opening. Our results are consistent with the CTD being a pH modulator. We propose herein a mechanism on how this domain may contribute to the pH-dependent opening of KcsA.     

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.     

Molecular architecture of full-length KcsA: role of cytoplasmic domains in ion permeation and activation gating

The molecular architecture of the NH(2) and COOH termini of the prokaryotic potassium channel KcsA has been determined using site-directed spin-labeling methods and paramagnetic resonance EPR spectroscopy. Cysteine mutants were generated (residues 5-24 and 121-160) and spin labeled, and the X-band CW EPR spectra were obtained from liposome-reconstituted channels at room temperature. Data on probe mobility (DeltaHo(-1)), accessibility parameters (PiO(2) and PiNiEdda), and inter-subunit spin-spin interaction (Omega) were used as structural constraints to build a three-dimensional folding model of these cytoplasmic domains from a set of simulated annealing and restrained molecular dynamics runs. 32 backbone structures were generated and averaged using fourfold symmetry, and a final mean structure was obtained from the eight lowest energy runs. Based on the present data, together with information from the KcsA crystal structure, a model for the three-dimensional fold of full-length KcsA was constructed. In this model, the NH(2) terminus of KcsA forms an alpha-helix anchored at the membrane-water interface, while the COOH terminus forms a right-handed four-helix bundle that extend some 40-50 A towards the cytoplasm. Functional analysis of COOH-terminal deletion constructs suggest that, while the COOH terminus does not play a substantial role in determining ion permeation properties, it exerts a modulatory role in the pH-dependent gating mechanism.     

Molecular determinants of tetramerization in the KcsA cytoplasmic domain

The cytoplasmic C-terminal domain (CTD) of KcsA, a bacterial homotetrameric potassium channel, is an amphiphilic domain that forms a helical bundle with four-fold symmetry mediated by hydrophobic and electrostatic interactions. Previously we have established that a CTD-derived 34-residue peptide associates into a tetramer in a pH-dependent manner (Kamnesky et al., JMB 2012;418:237-247). Here we further investigate the molecular determinants of tetramer formation in the CTD by characterizing the kinetics of monomer-tetramer equilibrium for 10 alanine mutants using NMR, sedimentation equilibrium (SE) and molecular dynamics simulation. NMR and SE concur in finding single-residue contributions to tetramer stability to be in the 0.5 to 3.5 kcal/mol range. Hydrophobic interactions between residues lining the tetramer core generally contributed more to formation of tetramer than electrostatic interactions between residues R147, D149 and E152. In particular, alanine replacement of residue R147, a key contributor to inter-subunit salt bridges, resulted in only a minor effect on tetramer dissociation. Mutations outside of the inter-subunit interface also influenced tetramer stability by affecting the tetramerization on-rate, possibly by changing the inherent helical propensity of the peptide. These findings are interpreted in the context of established paradigms of protein-protein interactions and protein folding, and lay the groundwork for further studies of the CTD in full-length KcsA channels.     

The distal C-terminal region of the KcsA potassium channel is a pH-dependent tetramerization domain

The intracellular C-terminal domain (CTD) of KcsA, a bacterial homotetrameric potassium channel, is a 40-residue-long segment that natively adopts a helical bundle conformation with 4-fold symmetry. A hallmark of KcsA behavior is pH-induced conformational change, which leads to the opening of the channel at acidic pH. Previous studies have reached conflicting conclusions as to the role of the CTD in this transition. Here, we investigate the involvement of this domain in pH-mediated channel opening by NMR using a soluble peptide corresponding to residues 128-160 of the CTD (CTD34). At neutral pH, CTD34 exhibits concentration-dependent spectral changes consistent with oligomer formation. We prove this slowly tumbling species to be a tetramer with a dissociation constant of (2.0±0.5)×10(-)(11)?M(3) by NMR and sedimentation equilibrium experiments. Whereas monomeric CTD34 is only mildly helical, secondary chemical shifts prove that the tetrameric species adopts a tight native-like helical bundle conformation. The tetrameric species undergoes pH-dependent dissociation, and CTD34 is fully monomeric below pH?5.0. The structural basis for this phenomenon is the destabilization of the tetrameric CTD34 by protonation of residue H145 in the monomeric form of the peptide. We conclude that (i) the CTD34 peptide is independently capable of forming a tetrameric helical bundle, and (ii) this structurally significant conformational shift is modulated by the effects of solution pH on residue H145. Therefore, the involvement of this domain in the pH gating of the channel is strongly suggested.     

Binding of anionic lipids to at least three nonannular sites on the potassium channel KcsA is required for channel opening

In addition to the annular or boundary lipids that surround the transmembrane surface of the potassium channel KcsA from Streptomyces lividans, x-ray crystallographic studies have detected one anionic lipid molecule bound at each protein-protein interface in the homotetrameric structure, at sites referred to as nonannular sites. The binding constant for phosphatidylglycerol at the nonannular sites has been determined using fluorescence quenching methods with a mutant of KcsA lacking the normal three lipid-exposed Trp residues. Binding is weak, with a binding constant of 0.42 ± 0.06 in units of mol fraction, implying that the nonannular sites will only be ∼70% occupied in bilayers of 100% phosphatidylglycerol. However, the nonannular sites show high selectivity for anionic lipids over zwitterionic lipids, and it is suggested that a change in packing at the protein-protein interface leads to a closing of the nonannular binding site in the unbound state. Increasing the anionic lipid content of the membrane leads to a large increase in open channel probability, from ∼2.5% in the presence of 25 mol % phosphatidylglycerol to ∼62% in 100 mol % phosphatidylglycerol. The relationship between open channel probability and phosphatidylglycerol content shows cooperativity. The data are consistent with a model in which three or four of the four nonannular sites in the KcsA homotetramer have to be occupied by anionic lipid for the channel to open. The conductance of the open channel increases with increasing concentration of anionic lipid, an effect possibly due to effects of anionic lipid on the concentration of K⁺ close to the membrane surface.     

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.     

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 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.     

A comparison between two prokaryotic potassium channels (KirBac1.1 and KcsA) in a molecular dynamics (MD) simulation study

The two potassium ion channels KirBac1.1 and KcsA are compared in a Molecular Dynamics (MD) simulation study. The location and motion of the potassium ions observed in the simulations are compared to those in the X-ray structures and previous simulations. In our simulations several of the crystallography resolved ion sites in KirBac1.1 are occupied by ions. In addition to this, two in KirBac1.1 unresolved sites where occupied by ions at sites that are in close correspondence to sites found in KcsA. There is every reason to believe that the conserved alignment of the selectivity filter in the potassium ion channel family corresponds to a very similar mechanism for ion transport across the filter. The gate residues, Phe146 in KirBac1.1 and Ala111 in KcsA acted in the simulations as effective barriers which never were passed by ions nor water molecules.     

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.     

The protonation state of the Glu-71/Asp-80 residues in the KcsA potassium channel: a first-principles QM/MM molecular dynamics study

Although a few x-ray structures of the KcsA K(+) channel have been crystallized several issues concerning the mechanisms of the ionic permeation and the protonation state of the selectivity filter ionizable side chains are still open. Using a first-principles quantum mechanical/molecular mechanical simulation approach, we have investigated the protonation state of Glu-71 and Asp-80, two important residues located in the vicinity of the selectivity filter. Results from the dynamics show that a proton is shared between the two residues, with a slight preference for Glu-71. The proton is found to exchange on the picosecond timescale, an interesting phenomenon that cannot be observed in classical molecular dynamics. Simulations of different ionic loading states of the filter show that the probability for the proton transfer is correlated with the filter occupancy. In addition, the Glu-71/Asp-80 pair is able to modulate the potential energy profile experienced by a K(+) ion as it translates along the pore axis. These theoretical predictions, along with recent experimental results, suggest that changes of the filter structure could be associated with a shift in the Glu-Asp protonation state, which in turn would influence the ion translocation.     

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.     

Phosphatidic acid plays a special role in stabilizing and folding of the tetrameric potassium channel KcsA

In this study, we investigated how the presence of anionic lipids influenced the stability and folding properties of the potassium channel KcsA. By using a combination of gel electrophoresis, tryptophan fluorescence and acrylamide quenching experiments, we found that the presence of the anionic lipid phosphatidylglycerol (PG) in a phosphatidylcholine (PC) bilayer slightly stabilized the tetramer and protected it from trifluoroethanol-induced dissociation. Surprisingly, the presence of phosphatidic acid (PA) had a much larger effect on the stability of KcsA and this lipid, in addition, significantly influenced the folding properties of the protein. The data indicate that PA creates some specificity over PG, and that it most likely stabilizes the tetramer via both electrostatic and hydrogen bond interactions.     

Arranging the elements of the potassium channel: the T1 domain occludes the cytoplasmic face of the channel

The voltage-gated potassium channel is currently one of the few membrane proteins where functional roles have been mapped onto specific segments of sequence. Although high-resolution structures of the transmembrane portions of three bacterial potassium channels, the tetramerization domain and the cytoplasmic ball are available, their relative spatial arrangement in mammalian channels remains a matter of ongoing debate. Cryo-electron microscopic images of the six transmembrane voltage-gated Kv channel have been reconstructed at up to 18 Å resolution, revealing that the T1 domain tetramerizes and is suspended below the transmembrane segments. However, the resolution of these images is insufficient to reveal the location of the third piece of the puzzle, the inactivating ball domain. We have used the aberrant interactions observed in a series of chimæric channels to establish that an assembled T1 domain restricts access to the cytoplasmic face of the channel, suggesting that the N-terminal ball and chain may be confined in the space between the T1 domain and the transmembrane portion of the channel.