Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K⁺ ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K⁺ channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K⁺ to Na⁺ substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na⁺ locally contracts the SF around Na⁺, releasing ligands that do not participate in Na⁺ coordination, and occluding the permeation pathway. K⁺ selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na⁺ ions due to basic structural differences between liquid water and the “quasi-liquid” SF matrix.     

Diverse gating in K channels: Differential role of the pore-helix glutamate in stabilizing the channel pore

The selectivity filter and adjacent regions in the bacterial KcsA and inwardly rectifying K⁺ (Kir) channels reveal significant conformational changes that cause the channel pore to transition from an activated to inactive state (C-type inactivation) once the channel is open. The meshwork of residues stabilizing the pore of KcsA involves Glu71–Asp80 carboxyl–carboxylate interaction ‘behind’ the selectivity filter. Interestingly, the Kir channels do not have this exact interaction, but instead have a Glu–Arg salt bridge where the Glu is in the same position but the Arg is one position N-terminal compared to the Asp in KcsA. Also, the Kir channels lack the Trp that hydrogen bonds to Asp80 in KcsA. Here, the sequence and structural information are combined to understand the dissimilarity in the role of the pore-helix Glu in stabilizing the pore structure in KcsA and Kir channels. This review illustrates that although Glu is quite conserved among both types of channels, the network of interactions is not translatable from one channel to the other; thereby suggesting a unique phenomenon of diverse gating patterns in K⁺ channels.     

K+ conduction in the selectivity filter of potassium channels is monitored by the charge distribution along their sequence

Potassium channels display a high conservation of sequence of the selectivity filter (SF), yet nature has designed a variety of channels that present a wide range of absolute rates of K(+) permeation. In KcsA, the structural archetype for K channels, under physiological concentrations, two K(+) ions reside in the SF in configurations 1,3 (up state) and 2,4 (down state) and ion conduction is believed to follow a throughput cycle involving a transition between these states. Using free-energy calculations of KcsA, Kv1.2, and mutant channels, we show that this transition is characterized by a channel-dependent energy barrier. This barrier is strongly influenced by the charges partitioned along the sequence of each channel. These results unveil therefore how, for similar structures of the SF, the rate of K(+) turnover may be fine-tuned within the family of potassium channels.     

Dynamics, Energetics, and Selectivity of the Low-K KcsA Channel Structure

Potassium channels are a diverse family of integral membrane proteins through which K⁺ can pass selectively. There is ongoing debate about the nature of conformational changes associated with the opening/closing and conductive/nonconductive states of potassium channels. The channels partly exert their function by varying their conductance through a mechanism known as C-type inactivation. Shortly after the activation of K⁺ channels, their selectivity filter stops conducting ions at a rate that depends on various stimuli. The molecular mechanism of C-type inactivation has not been fully understood yet. However, the X-ray structure of the KcsA channel obtained in the presence of low K⁺ concentration is thought to be representative of a K⁺ channel in the C-type inactivated state. Here, extensive, fully atomistic molecular dynamics and free-energy simulations of the low-K⁺ KcsA structure in an explicit lipid bilayer are performed to evaluate the stability of this structure and the selectivity of its binding sites. We find that the low-K⁺ KcsA structure is stable on the timescale of the molecular dynamics simulations performed, and that ions preferably remain in S1 and S4. In the absence of ions, the selectivity filter evolves toward an asymmetric architecture, as already observed in other computations of the high-K⁺ structure of KcsA and KirBac. The low-K⁺ KcsA structure is not permeable by Na⁺, K⁺, or Rb⁺, and the selectivity of its binding sites is different from that of the high-K⁺ structure.     

Comparative study of the energetics of ion permeation in Kv1.2 and KcsA potassium channels

Biological ion channels rely on a multi-ion transport mechanism for fast yet selective permeation of ions. The crystal structure of the KcsA potassium channel provided the first microscopic picture of this process. A similar mechanism is assumed to operate in all potassium channels, but the validity of this assumption has not been well investigated. Here, we examine the energetics of ion permeation in Shaker Kv1.2 and KcsA channels, which exemplify the six-transmembrane voltage-gated and two-transmembrane inward-rectifier channels. We study the feasibility of binding a third ion to the filter and the concerted motion of ions in the channel by constructing the potential of mean force for K⁺ ions in various configurations. For both channels, we find that a pair of K⁺ ions can move almost freely within the filter, but a relatively large free-energy barrier hinders the K⁺ ion from stepping outside the filter. We discuss the effect of the CMAP dihedral energy correction that was recently incorporated into the CHARMM force field on ion permeation dynamics.     

Pore-forming transmembrane domains control ion selectivity and selectivity filter conformation in the KirBac1.1 potassium channel

Potassium (K⁺) channels are membrane proteins with the remarkable ability to very selectively conduct K⁺ ions across the membrane. High-resolution structures have revealed that dehydrated K⁺ ions permeate through the narrowest region of the pore, formed by the backbone carbonyls of the signature selectivity filter (SF) sequence TxGYG. However, the existence of nonselective channels with similar SF sequences, as well as effects of mutations in other regions on selectivity, suggest that the SF is not the sole determinant of selectivity. We changed the selectivity of the KirBac1.1 channel by introducing mutations at residue I131 in transmembrane helix 2 (TM2). These mutations increase Na⁺ flux in the absence of K⁺ and introduce significant proton conductance. Consistent with K⁺ channel crystal structures, single-molecule FRET experiments show that the SF is conformationally constrained and stable in high-K⁺ conditions but undergoes transitions to dilated low-FRET states in high-Na⁺/low-K⁺ conditions. Relative to wild-type channels, I131M mutants exhibit marked shifts in the K⁺ and Na⁺ dependence of SF dynamics to higher K⁺ and lower Na⁺ concentrations. These results illuminate the role of I131, and potentially other structural elements outside the SF, in controlling ion selectivity, by suggesting that the physical interaction of these elements with the SF contributes to the relative stability of the constrained K⁺-induced SF configuration versus nonselective dilated conformations.     

Computational study of non-conductive selectivity filter conformations and C-type inactivation in a voltage-dependent potassium channel

C-type inactivation is a time-dependent process of great physiological significance that is observed in a large class of K⁺ channels. Experimental and computational studies of the pH-activated KcsA channel show that the functional C-type inactivated state, for this channel, is associated with a structural constriction of the selectivity filter at the level of the central glycine residue in the signature sequence, TTV(G)YGD. The structural constriction is allosterically promoted by the wide opening of the intracellular activation gate. However, whether this is a universal mechanism for C-type inactivation has not been established with certainty because similar constricted structures have not been observed for other K⁺ channels. Seeking to ascertain the general plausibility of the constricted filter conformation, molecular dynamics simulations of a homology model of the pore domain of the voltage-gated potassium channel Shaker were performed. Simulations performed with an open intracellular gate spontaneously resulted in a stable constricted-like filter conformation, providing a plausible nonconductive state responsible for C-type inactivation in the Shaker channel. While there are broad similarities with the constricted structure of KcsA, the hypothetical constricted-like conformation of Shaker also displays some subtle differences. Interestingly, those are recapitulated by the Shaker-like E71V KcsA mutant, suggesting that the residue at this position along the pore helix plays a pivotal role in determining the C-type inactivation behavior. Free energy landscape calculations show that the conductive-to-constricted transition in Shaker is allosterically controlled by the degree of opening of the intracellular activation gate, as observed with the KcsA channel. The behavior of the classic inactivating W434F Shaker mutant is also characterized from a 10-μs MD simulation, revealing that the selectivity filter spontaneously adopts a nonconductive conformation that is constricted at the level of the second glycine in the signature sequence, TTVGY(G)D.     

Non-Equilibrium Dynamics Contribute to Ion Selectivity in the KcsA Channel

The ability of biological ion channels to conduct selected ions across cell membranes is critical for the survival of both animal and bacterial cells. Numerous investigations of ion selectivity have been conducted over more than 50 years, yet the mechanisms whereby the channels select certain ions and reject others are not well understood. Here we report a new application of Jarzynski’s Equality to investigate the mechanism of ion selectivity using non-equilibrium molecular dynamics simulations of Na⁺ and K⁺ ions moving through the KcsA channel. The simulations show that the selectivity filter of KcsA adapts and responds to the presence of the ions with structural rearrangements that are different for Na⁺ and K⁺. These structural rearrangements facilitate entry of K⁺ ions into the selectivity filter and permeation through the channel, and rejection of Na⁺ ions. A mechanistic model of ion selectivity by this channel based on the results of the simulations relates the structural rearrangement of the selectivity filter to the differential dehydration of ions and multiple-ion occupancy and describes a mechanism to efficiently select and conduct K⁺. Estimates of the K⁺/Na⁺ selectivity ratio and steady state ion conductance for KcsA from the simulations are in good quantitative agreement with experimental measurements. This model also accurately describes experimental observations of channel block by cytoplasmic Na⁺ ions, the “punch through” relief of channel block by cytoplasmic positive voltages, and is consistent with the knock-on mechanism of ion permeation.     

Interaction of local anesthetics with the K+ channel pore domain

Local anesthetics and related drugs block ionic currents of Na⁺, K⁺ and Ca²⁺ conducted across the cell membrane by voltage-dependent ion channels. Many of these drugs bind in the permeation pathway, occlude the pore and stop ion movement. However channel-blocking drugs have also been associated with decreased membrane stability of certain tetrameric K⁺ channels, similar to the destabilization of channel function observed at low extracellular K⁺ concentration. Such drug-dependent stability may result from electrostatic repulsion of K⁺ from the selectivity filter by a cationic drug molecule bound in the central cavity of the channel. In this study we used the pore domain of the KcsA K⁺ channel protein to test this hypothesis experimentally with a biochemical assay of tetramer stability and theoretically by computational simulation of local anesthetic docking to the central cavity. We find that two common local anesthetics, lidocaine and tetracaine, promote thermal dissociation of the KcsA tetramer in a K⁺-dependent fashion. Docking simulations of these drugs with open, open-inactivated and closed crystal structures of KcsA yield many energetically favorable drug-channel complexes characterized by nonbonded attraction to pore-lining residues and electrostatic repulsion of K⁺. The results suggest that binding of cationic drugs to the inner cavity can reduce tetramer stability of K⁺ channels.     

Transferring knowledge towards understanding the pore stabilizing variations in K+ channels

Recent advances in structural biology underlying mechanisms of channel gating have strengthened our knowledge about how K⁺ channels can be inter-convertible between conductive and non-conductive states. We have reviewed and combined mutagenesis with biochemical, biophysical and structural information in order to understand the critical roles of the pore residues in stabilizing the pore structure and channel open state. We also discuss how the latest knowledge on the K⁺ channel KcsA may provide a step towards better understanding of distinct pore stabilizing differences among diversified K⁺ channels.     

The conserved potassium channel filter can have distinct ion binding profiles: Structural analysis of rubidium,cesium, and barium binding in NaK2K

Potassium channels are highly selective for K⁺ over the smaller Na⁺. Intriguingly, they are permeable to larger monovalent cations such as Rb⁺ and Cs⁺ but are specifically blocked by the similarly sized Ba²⁺. In this study, we used structural analysis to determine the binding profiles for these permeant and blocking ions in the selectivity filter of the potassium-selective NaK channel mutant NaK2K and also performed permeation experiments using single-channel recordings. Our data revealed that some ion binding properties of NaK2K are distinct from those of the canonical K⁺ channels KcsA and MthK. Rb⁺ bound at sites 1, 3, and 4 in NaK2K, as it does in KcsA. Cs⁺, however, bound predominantly at sites 1 and 3 in NaK2K, whereas it binds at sites 1, 3, and 4 in KcsA. Moreover, Ba²⁺ binding in NaK2K was distinct from that which has been observed in KcsA and MthK, even though all of these channels show similar Ba²⁺ block. In the presence of K⁺, Ba²⁺ bound to the NaK2K channel at site 3 in conjunction with a K⁺ at site 1; this led to a prolonged block of the channel (the external K⁺-dependent Ba²⁺ lock-in state). In the absence of K⁺, however, Ba²⁺ acts as a permeating blocker. We found that, under these conditions, Ba²⁺ bound at sites 1 or 0 as well as site 3, allowing it to enter the filter from the intracellular side and exit from the extracellular side. The difference in the Ba²⁺ binding profile in the presence and absence of K⁺ thus provides a structural explanation for the short and prolonged Ba²⁺ block observed in NaK2K.     

Sodium permeability and sensitivity induced by mutations in the selectivity filter of the KcsA channel towards Kir channels

The bacterial potassium (K⁺) channel KcsA provides an attractive model system to study ion permeation behavior in a selective K⁺-channel. We changed residue at the N-terminal end of the selectivity filter of KcsA (T74V) to its counterpart in inwardly rectifying K⁺-channels (Kir). The tetramer was found to be stable as unmodified KcsA. Under symmetrical and asymmetrical conditions, Na⁺ increased the inward current in the virtual absence of K⁺ however outward currents were nearly abolished which could be recovered upon internal K⁺ addition. Na⁺ also drastically increased the channel open time either in the presence or virtual absence of K⁺. Furthermore, the T74V mutation decreased the internal Ba²⁺ affinity of the channel possibly by binding to a K⁺ site in the pore. In additional experiments, another point mutation V76I in T74V mutant was carried out thus the selectivity filter resembled more the selectivity filter of Kir channels. The mutant tetramer was converted into monomers as determined by conventional gel electrophoresis. However, native like gel electrophoresis, Trp fluorescence and acrylamide quenching experiments indicated that this mutant still formed a tetramer and apparently adopted similar folding properties as unmodified KcsA. Single-channel experiments further demonstrated that the channel was selective for K⁺ over Na⁺ as Na⁺ blocked channel currents. These data suggest that single point mutation T74V alters the selectivity filter and allows simultaneous occupancy and conduction of K⁺ and Na⁺ probably via ion–ion interaction in the pore. In contrast, both mutations (T74V and V76I) in the same molecule seem to reorganize the pore conformation which controls the overall stability of a selective K⁺-channel.     

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.     

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.     

A multipoint hydrogen-bond network underlying KcsA C-type inactivation

In the prokaryotic potassium channel KcsA activation gating at the inner bundle gate is followed by C-type inactivation at the selectivity filter. Entry into the C-type inactivated state has been directly linked to the strength of the H-bond interaction between residues Glu-71 and Asp-80 behind the filter, and is allosterically triggered by the rearrangement of the inner bundle gate. Here, we show that H-bond pairing between residues Trp-67 and Asp-80, conserved in most K⁺ channels, constitutes another critical interaction that determines the rate and extent of KcsA C-type inactivation. Disruption of the equivalent interaction in Shaker (Trp-434-Asp-447) and Kv1.2 (Trp-366-Asp-379) leads also to modulation of the inactivation process, suggesting that these residues also play an analogous role in the inactivation gating of Kv channels. The present results show that in KcsA C-type inactivation gating is governed by a multipoint hydrogen-bond network formed by the triad Trp-67-Glu71-Asp-80. This triad exerts a critical role in the dynamics and conformational stability of the selectivity filter and might serve as a general modulator of selectivity filter gating in other members of the K⁺ channel family.     

Mutations in the K+-Channel KcsA Toward Kir Channels Alter Salt-Induced Clusterization and Blockade by Quaternary Alkylammonium Ions

Protein aggregation is a result of malfunction in protein folding, assembly, and transport, caused by protein mutation and/or changes in the cell environment, thus triggering many human diseases. We have shown that bacterial K⁺-channel KcsA, which acts as a representative model for ion channels, forms salt-induced large conductive complexes in a particular environment. In the present study, we investigated the effects of point mutations in the selectivity filter of KcsA on intrinsic stability, aggregation, and channel blocking behavior. First, we found that a low sodium chloride concentration in potassium-containing media induced fast transfer of single channels to a planar lipid bilayer. Second, increasing the sodium chloride concentration drastically increased the total channel current, indicating enhanced vesicle fusion and transfer of multiple channels to a planar lipid bilayer. However, such complexes exhibited high conductance as well as higher open probability compared to the unmodified KcsA behavior shown previously. Interestingly, the affinity of aggregated complexes for larger symmetric quaternary alkylammonium ions (QAs) was found to be much higher than that for tetraethylammonium, a classical blocker of the K⁺ channel. Based on these findings, we propose that mutant channel complexes exhibit larger pore dimensions, thus resembling more the topological properties of voltage-gated and inwardly rectifying K⁺ channels.     

A computational study of barium blockades in the KcsA potassium channel based on multi-ion potential of mean force calculations and free energy perturbation

Electrophysiological studies have established that the permeation of Ba²⁺ ions through the KcsA K⁺-channel is impeded by the presence of K⁺ ions in the external solution, while no effect is observed for external Na⁺ ions. This Ba²⁺ “lock-in” effect suggests that at least one of the external binding sites of the KcsA channel is thermodynamically selective for K⁺. We used molecular dynamics simulations to interpret these lock-in experiments in the context of the crystallographic structure of KcsA. Assuming that the Ba²⁺ is bound in site S₂ in the dominant blocked state, we examine the conditions that could impede its translocation and cause the observed “lock-in” effect. Although the binding of a K⁺ ion to site S₁ when site S₂ is occupied by Ba²⁺ is prohibitively high in energy (>10 kcal/mol), binding to site S₀ appears to be more plausible (ΔG > 4 kcal/mol). The 2D potential of mean force (PMF) for the simultaneous translocation of Ba²⁺ from site S₂ to site S₁ and of a K⁺ ion on the extracellular side shows a barrier that is consistent with the concept of external lock-in. The barrier opposing the movement of Ba²⁺ is very high when a cation is in site S₀, and considerably smaller when the site is unoccupied. Furthermore, free energy perturbation calculations show that site S₀ is selective for K⁺ by 1.8 kcal/mol when S₂ is occupied by Ba²⁺. However, the same site S₀ is nonselective when site S₂ is occupied by K⁺, which shows that the presence of Ba²⁺ affects the selectivity of the pore. A theoretical framework within classical rate theory is presented to incorporate the concentration dependence of the external ions on the lock-in effect.     

Conduction of Na+ and K+ through the NaK channel: molecular and Brownian dynamics studies

Conduction of ions through the NaK channel, with M0 helix removed, was studied using both Brownian dynamics and molecular dynamics. Brownian dynamics simulations predict that the truncated NaK has approximately a third of the conductance of the related KcsA K⁺ channel, is outwardly rectifying, and has a Michaelis-Menten current-concentration relationship. Current magnitude increases when the glutamine residue located near the intracellular gate is replaced with a glutamate residue. The channel is blocked by extracellular Ca²⁺. Molecular dynamics simulations show that, under the influence of a strong applied potential, both Na⁺ and K⁺ move across the selectivity filter, although conduction rates for Na⁺ ions are somewhat lower. The mechanism of conduction of Na⁺ differs significantly from that of K⁺ in that Na⁺ is preferentially coordinated by single planes of pore-lining carbonyl oxygens, instead of two planes as in the usual K⁺ binding sites. The water-containing filter pocket resulting from a single change in the selectivity filter sequence (compared to potassium channels) disrupts several of the planes of carbonyl oxygens, and thus reduces the filter's ability to discriminate against sodium.     

KcsA 通道对Na+、K+及Rb+离子选择性的统计热力学研究

钾离子的通透率至少比钠离子的通透率大10000倍,这个问题至今没有很好地解决.为了在分子水平阐释钾离子通道的选择性机制,以KcsA钾通道X射线衍射结构为基础,采用密度泛函理论计算了不同离子在离子通道中的位能.计算结果表明,Rb+离子具有与K+离子相类似的位能曲线,但是其在通透过程遇到的位垒要比K+离子的位垒高,因而所对应的通透率也就小于钾离子的通透率,而钠离子的的通透率仅仅是钾离子通透率的0.0067%.文中所涉及的系统仅仅包含269个原子,而用分子动力学虽然也可以得到相近的结果,但是它的系统大小为41 000个原子.     

The selectivity of K+ ion channels: testing the hypotheses

How K⁺ channels are able to conduct certain cations yet not others remains an important but unresolved question. The recent elucidation of the structure of NaK, an ion channel that conducts both Na⁺ and K⁺ ions, offers an opportunity to test the various hypotheses that have been put forward to explain the selectivity of K⁺ ion channels. We test the snug-fit, field-strength, and over-coordination hypotheses by comparing their predictions to the results of classical molecular dynamics simulations of the K⁺ selective channel KcsA and the less selective channel NaK embedded in lipid bilayers. Our results are incompatible with the so-called strong variant of the snug-fit hypothesis but are consistent with the over-coordination hypothesis and neither confirm nor refute the field-strength hypothesis. We also find that the ions and waters in the NaK selectivity filter unexpectedly move to a new conformation in seven K⁺ simulations: the two K⁺ ions rapidly move from site S4 to S2 and from the cavity to S4. At the same time, the selectivity filter narrows around sites S1 and S2 and the carbonyl oxygen atoms rotate 20°−40° inwards toward the ion. These motions diminish the large structural differences between the crystallographic structures of the selectivity filters of NaK and KcsA and appear to allow the binding of ions to S2 of NaK at physiological temperature.