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.     

Widening the scope of constriction

JGP modeling study suggests that selectivity filter constriction is a plausible mechanism for C-type inactivation of the Shaker voltage-gated potassium channel.

In response to prolonged activation, many K⁺ channels spontaneously reduce the membrane conductance by undergoing C-type inactivation, a kinetic process crucial for the pacing of cardiac action potentials and the modulation of neuronal firing patterns. In the pH-activated bacterial channel KcsA, C-type inactivation appears to involve constriction of the channel’s selectivity filer that prohibits ion conduction, but whether voltage-gated channels like Drosophila Shaker use a similar mechanism is controversial (1). In this issue of JGP, a computational study by Li et al. suggests that filter constriction is indeed a plausible mechanism for the C-type inactivation of Shaker (2).(Left to right) Jing Li, Benoît Roux, and colleagues use computational modeling to show that selectivity filter constriction, allosterically promoted by opening of the intracellular activation gate, is a plausible mechanism for the C-type inactivation of voltage-gated K⁺ channels such as Drosophila Shaker. The selectivity filter is conductive (left) when the intracellular gate is partially open, but adopts a constricted conformation (right) when the gate is open wide.Various structural approaches have shown that C-type inactivation of KcsA channels is associated with the symmetrical constriction of all four channel subunits at the level of the central glycine residue in the selectivity filter. Benoît Roux and colleagues at The University of Chicago used MD simulations to show that the KcsA pore can transition from the conductive to the constricted conformation on an appropriate timescale, and that this transition is allosterically promoted by the wide opening of the pore’s intracellular gate (3). Modeling by Roux and colleagues suggests that C-type inactivation of cardiac hERG channels could also involve selectivity filter constriction, though in this case it appears to be an asymmetric process in which only two of the channel’s subunits move closer together (4).“In view of the high similarity between the pore domains of Shaker and KcsA (almost 40% sequence identity), we wanted to examine if it’s possible for the Shaker selectivity filter to constrict and, if so, how similar it is to KcsA,” Roux explains. Led by first author Jing Li—now an assistant professor at the University of Mississippi—Roux and colleagues developed several homology models of the Shaker pore domain with the intracellular gate open to various degrees (2).MD simulations and free energy calculations revealed that the Shaker selectivity filter can dynamically transition from a conductive to a constricted conformation, and that this transition is allosterically coupled to the intracellular gate; the constricted conformation is stable when the gate is wide open. “Our computations strongly suggest that constriction is a plausible mechanism for the C-type inactivation of Shaker,” Roux says. “There’s no reason based on the currently available information to reject the existence of a constricted state in Shaker channels.”As with KcsA, Shaker channels appear to constrict symmetrically at the level of the selectivity filter’s central glycine. But Li et al.’s simulations revealed some small variations between the two channels, including differences in the number of water molecules bound to each channel subunit and the arrangement of the hydrogen-bond network they form to stabilize the constricted state.Li et al. also modeled the pore domain of the Shaker W434F mutant, which is widely assumed to be trapped in a C-type inactivated state. The simulation suggests that the mutant channel’s filter adopts a stable constricted conformation even when the intracellular gate is only partially open, although the constriction is asymmetric and occurs at the level of a different filter residue (2).Constriction may therefore be a universal mechanism of C-type inactivation, even if the exact conformation varies from channel to channel. But, says Roux, confirming this will require more experimental work using the right conditions and mutations to capture the structure of inactivated channels.     

Potassium and sodium ions in a potassium channel studied by molecular dynamics simulations

We have performed simulations of both a single potassium ion and a single sodium ion within the pore of the bacterial potassium channel KcsA. For both ions there is a dehydration energy barrier at the cytoplasmic mouth suggesting that the crystal structure is a closed conformation of the channel. There is a potential energy barrier for a sodium ion in the selectivity filter that is not seen for potassium. Radial distribution functions for both ions with the carbonyl oxygens of the selectivity filter indicate that sodium may interact more tightly with the filter than does potassium. This suggests that the key to the ion selectivity of KcsA is the greater dehydration energy of Na⁺ ions, and helps to explain the block of KcsA by internal Na⁺ ions.     

Voltage-insensitive Gating after Charge-neutralizing Mutations in the S4 Segment of Shaker Channels

Shaker channel mutants, in which the first (R362), second (R365), and fourth (R371) basic residues in the S4 segment have been neutralized, are found to pass potassium currents with voltage-insensitive kinetics when expressed in Xenopus oocytes. Single channel recordings clarify that these channels continue to open and close from −160 to +80 mV with a constant opening probability (P _o). Although P _o is low (∼0.15) in these mutants, mean open time is voltage independent and similar to that of control Shaker channels. Additionally, these mutant channels retain characteristic Shaker channel selectivity, sensitivity to block by 4-aminopyridine, and are partially blocked by external Ca²⁺ ions at very negative potentials. Furthermore, mean open time is approximately doubled, in both mutant channels and control Shaker channels, when Rb⁺ is substituted for K⁺ as the permeant ion species. Such strong similarities between mutant channels and control Shaker channels suggests that the pore region has not been substantially altered by the S4 charge neutralizations. We conclude that single channel kinetics in these mutants may indicate how Shaker channels would behave in the absence of voltage sensor input. Thus, mean open times appear primarily determined by voltage-insensitive transitions close to the open state rather than by voltage sensor movement, even in control, voltage-sensitive Shaker channels. By contrast, the low and voltage-insensitive P _o seen in these mutant channels suggests that important determinants of normal channel opening derive from electrostatic coupling between S4 charges and the pore domain.     

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.     

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.     

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.     

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.     

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.     

The Molecular Basis of Polyunsaturated Fatty Acid Interactions with the Shaker Voltage-Gated Potassium Channel

Voltage-gated potassium (K_V) channels are membrane proteins that respond to changes in membrane potential by enabling K⁺ ion flux across the membrane. Polyunsaturated fatty acids (PUFAs) induce channel opening by modulating the voltage-sensitivity, which can provide effective treatment against refractory epilepsy by means of a ketogenic diet. While PUFAs have been reported to influence the gating mechanism by electrostatic interactions to the voltage-sensor domain (VSD), the exact PUFA-protein interactions are still elusive. In this study, we report on the interactions between the Shaker K_V channel in open and closed states and a PUFA-enriched lipid bilayer using microsecond molecular dynamics simulations. We determined a putative PUFA binding site in the open state of the channel located at the protein-lipid interface in the vicinity of the extracellular halves of the S3 and S4 helices of the VSD. In particular, the lipophilic PUFA tail covered a wide range of non-specific hydrophobic interactions in the hydrophobic central core of the protein-lipid interface, while the carboxylic head group displayed more specific interactions to polar/charged residues at the extracellular regions of the S3 and S4 helices, encompassing the S3-S4 linker. Moreover, by studying the interactions between saturated fatty acids (SFA) and the Shaker K_V channel, our study confirmed an increased conformational flexibility in the polyunsaturated carbon tails compared to saturated carbon chains, which may explain the specificity of PUFA action on channel proteins.     

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.     

Absence of ion-binding affinity in the putatively inactivated low-[K+] structure of the KcsA potassium channel

Potassium channels are membrane proteins that selectively conduct K(+) across cellular membranes. The narrowest part of their pore, the selectivity filter, is responsible for distinguishing K(+) from Na(+), and can also act as a gate through a mechanism known as C-type inactivation. It has been proposed that a conformation of the KcsA channel obtained by crystallization in presence of low concentration of K(+) (PDB 1K4D) could correspond to the C-type inactivated state. Here, we show using molecular mechanics simulations that such conformation has little ion-binding affinity and that ions do not contribute to its stability. The simulations suggest that, in this conformation, the selectivity filter is mostly occupied by water molecules. Whether such ion-free state of the KcsA channel is physiologically accessible and representative of the inactivated state of eukaryotic channels remains unclear.     

Functional Equilibrium of the KcsA Structure Revealed by NMR

KcsA is a tetrameric K⁺ channel that is activated by acidic pH. Under open conditions of the helix bundle crossing, the selectivity filter undergoes an equilibrium between permeable and impermeable conformations. Here we report that the population of the permeable conformation (p_perm) positively correlates with the tetrameric stability and that the population in reconstituted high density lipoprotein, where KcsA is surrounded by the lipid bilayer, is lower than that in detergent micelles, indicating that dynamic properties of KcsA are different in these two media. Perturbation of the membrane environment by the addition of 1–3% 2,2,2-trifluoroethanol increases p_perm and the open probability, revealed by NMR and single-channel recording analyses. These results demonstrate that KcsA inactivation is determined not only by the protein itself but also by the surrounding membrane environments.     

Bupivacaine blocks N-type inactivating Kv channels in the open state: no allosteric effect on inactivation kinetics

Local anesthetics bind to ion channels in a state-dependent manner. For noninactivating voltage-gated K channels the binding mainly occurs in the open state, while for voltage-gated inactivating Na channels it is assumed to occur mainly in inactivated states, leading to an allosterically caused increase in the inactivation probability, reflected in a negative shift of the steady-state inactivation curve, prolonged recovery from inactivation, and a frequency-dependent block. How local anesthetics bind to N-type inactivating K channels is less explored. In this study, we have compared bupivacaine effects on inactivating (Shaker and K_v3.4) and noninactivating (Shaker-IR and K_v3.2) channels, expressed in Xenopus oocytes. Bupivacaine was found to block these channels time-dependently without shifting the steady-state inactivation curve markedly, without a prolonged recovery from inactivation, and without a frequency-dependent block. An analysis, including computational testing of kinetic models, suggests binding to the channel mainly in the open state, with affinities close to those estimated for corresponding noninactivating channels (300 and 280 μM for Shaker and Shaker-IR, and 60 and 90 μM for K_v3.4 and K_v3.2). The similar magnitudes of K_d, as well as of blocking and unblocking rate constants for inactivating and noninactivating Shaker channels, most likely exclude allosteric interactions between the inactivation mechanism and the binding site. The relevance of these results for understanding the action of local anesthetics on Na channels is discussed.     

Role of polyphosphate in regulation of the Streptomyces lividans KcsA channel

We examine the hypotheses that the Streptomyces lividans potassium channel KcsA is gated at neutral pH by the electrochemical potential, and that its selectivity and conductance are governed at the cytoplasmic face by interactions between the KcsA polypeptides and a core molecule of inorganic polyphosphate (polyP). The four polypeptides of KcsA are postulated to surround the end unit of the polyP molecule with a collar of eight arginines, thereby modulating the negative charge of the polyP end unit and increasing its preference for binding monovalent cations. Here we show that KcsA channels can be activated in planar lipid bilayers at pH 7.4 by the chemical potential alone. Moreover, one or both of the C-terminal arginines are replaced with residues of progressively lower basicity-lysine, histidine, valine, asparagine-and the effects of these mutations on conductance and selectivity for K⁺ over Mg²⁺ is tested in planar bilayers as a function of Mg²⁺ concentration and pH. As the basicity of the C-terminal residues decreases, Mg²⁺ block increases, and Mg²⁺ becomes permeant when medium pH is greater than the pI of the C-terminal residues. The results uphold the premise that polyP and the C-terminal arginines are decisive elements in KcsA channel regulation.     

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.     

Intermediate Conductances during Deactivation of Heteromultimeric Shaker Potassium Channels

A previous study of the T442S mutant Shaker channel revealed activation-coupled subconductance levels that apparently represent kinetic intermediates in channel activation (Zheng, J., and F.J. Sigworth. 1997. J. Gen. Physiol. 110:101–117). We have now extended the study to heteromultimeric channels consisting of various numbers of mutant subunits as well as channels without mutant subunits, all in the background of a chimeric Shaker channel having increased conductance. It has been found that activation-coupled sublevels exist in all these channel types, and are traversed in at least 80% of all deactivation time courses. In symmetric K⁺ solutions, the currents in the two sublevels have a linear voltage dependence, being 23–44% and 54–70% of the fully open conductance. Sublevels in different channel types share similar voltage dependence of the mean lifetime and similar ion selectivity properties. However, the mean lifetime of each current level depends approximately geometrically on the number of mutant subunits in the channel, becoming shorter in channels having fewer mutant subunits. Each mutant subunit appears to stabilize all of the conducting states by ∼0.5 kcal/mol. Consistent with previous results in the mutant channel, sublevels in channels with two or no mutant subunits also showed ion selectivities that differ from that of the fully open level, having relatively higher K⁺ than Rb⁺ conductances. A model is presented in which Shaker channels have two coupled activation gates, one associated with the selectivity filter and a second associated with the S6 helix bundle.     

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.     

Homology modeling and molecular dynamics simulation studies of an inward rectifier potassium channel

A homology model has been generated for the pore-forming domain of Kir6.2, a component of an ATP-sensitive K channel, based on the x-ray structure of the bacterial channel KcsA. Analysis of the lipid-exposed and pore-lining surfaces of the model reveals them to be compatible with the known features of membrane proteins and Kir channels, respectively. The Kir6.2 homology model was used as the starting point for nanosecond-duration molecular dynamics simulations in a solvated phospholipid bilayer. The overall drift from the model structure was comparable to that seen for KcsA in previous similar simulations. Preliminary analysis of the interactions of the Kir6.2 channel model with K(+) ions and water molecules during these simulations suggests that concerted single-file motion of K(+) ions and water through the selectivity filter occurs. This is similar to such motion observed in simulations of KcsA. This suggests that a single-filing mechanism is conserved between different K channel structures and may be robust to changes in simulation details. Comparison of Kir6.2 and KcsA suggests some degree of flexibility in the filter, thus complicating models of ion selectivity based upon a rigid filter.     

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.