The Open Gate of the KV1.2 Channel: Quantum Calculations Show the Key Role of Hydration

The open gate of the K_v1.2 voltage-gated potassium channel can just hold a hydrated K⁺ ion. Quantum calculations starting from the x-ray coordinates of the channel confirm this, showing little change from the x-ray coordinates for the protein. Water molecules not in the x-ray coordinates, and the ion itself, are placed by the calculation. The water molecules, including their orientation and hydrogen bonding, with and without an ion, are critical for the path of the ion, from the solution to the gate. A sequence of steps is postulated in which the potential experienced by the ion in the pore is influenced by the position of the ion. The gate structure, with and without the ion, has been optimized. The charges on the atoms and bond lengths have been calculated using natural bond orbital calculations, giving K⁺ ∼0.77 charges, rather than 1.0. The PVPV hinge sequence has been mutated in silico to PVVV (P407V in the 2A79 numbering). The water structure around the ion becomes discontinuous, separated into two sections, above and below the ion. PVPV conservation closely relates to maintaining the water structure. Finally, these results have implications concerning gating.     

Quantum mechanical calculations of charge effects on gating the KcsA channel

A series of ab initio (density functional) calculations were carried out on side chains of a set of amino acids, plus water, from the (intracellular) gating region of the KcsA K⁺ channel. Their atomic coordinates, except hydrogen, are known from X-ray structures [D.A. Doyle, J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, R. MacKinnon, The structure of the potassium channel: molecular basis of K⁺ conduction and selectivity, Science 280 (1998) 69-77; R. MacKinnon, S.L. Cohen, A. Kuo, A. Lee, B.T. Chait, Structural conservation in prokaryotic and eukaryotic potassium channels, Science 280 (1998) 106-109; Y. Jiang, A. Lee, J. Chen, M. Cadene, B.T. Chait, R. MacKinnon, The open pore conformation of potassium channels. Nature 417 (2001) 523-526], as are the coordinates of some water oxygen atoms. The 1k4c structure is used for the starting coordinates. Quantum mechanical optimization, in spite of the starting configuration, places the atoms in positions much closer to the 1j95, more tightly closed, configuration. This state shows four water molecules forming a “basket” under the Q119 side chains, blocking the channel. When a hydrated K⁺ approaches this “basket”, the optimized system shows a strong set of hydrogen bonds with the K⁺ at defined positions, preventing further approach of the K⁺ to the basket. This optimized structure with hydrated K⁺ added shows an ice-like 12 molecule nanocrystal of water. If the water molecules exchange, unless they do it as a group, the channel will remain blocked. The “basket” itself appears to be very stable, although it is possible that the K⁺ with its hydrating water molecules may be more mobile, capable of withdrawing from the gate. It is also not surprising that water essentially freezes, or forms a kind of glue, in a nanometer space; this agrees with experimental results on a rather different, but similarly sized (nm dimensions) system [K.B. Jinesh, J.W.M. Frenken, Capillary condensation in atomic scale friction: how water acts like a glue, Phys. Rev. Lett. 96 (2006) 166103/1-4]. It also agrees qualitatively with simulations on channels [A. Anishkin, S. Sukharev, Water dynamics and dewetting transitions in the small mechanosensitive channel MscS, Biophys. J. 86 (2004) 2883-2895; O. Beckstein, M.S.P. Sansom, Liquid-vapor oscillations of water in hydrophobic nanopores, Proc. Natl Acad. Sci. U. S. A. 100 (2003) 7063-7068] and on featureless channel-like systems [J. Lu, M.E. Green, Simulation of water in a pore with charges: application to a gating mechanism for ion channels, Prog. Colloid Polym. Sci. 103 (1997) 121-129], in that it forms a boundary on water that is not obvious from the liquid state. The idea that a structure is stable, even if individual molecules exchange, is well known, for example from the hydration shell of ions. We show that when charges are added in the form of protons to the domains (one proton per domain), the optimized structure is open. No stable water hydrogen bonds hold it together; an opening of 11.0 Å appears, measured diagonally between non-neighboring domains as glutamine 119 carbonyl O-O distance. This is comparable to the opening in the MthK potassium channel structure that is generally agreed to be open. The appearance of the opening is in rather good agreement with that found by Perozo and coworkers. In contrast, in the uncharged structure this diagonal distance is 6.5 Å, and the water “basket” constricts the uncharged opening still further, with the ice-like structure that couples the K⁺ ion to the gating region freezing the entrance to the channel. Comparison with our earlier model for voltage gated channels suggests that a similar mechanism may apply in those channels.     

Molecular dynamics study of the KcsA channel at 2.0-Å resolution: stability and concerted motions within the pore

The stability of the KcsA channel accommodating more than one ion in the pore has been studied with molecular dynamics. We have used the very last X-ray structure of the KcsA channel at 2.0-Å resolution determined by Zhou et al. [Nature 414 (2001) 43]. In this channel, six of the seven experimentally evidenced sites have been considered. We show that the protein remains very stable in the presence of four K⁺ ions (three in the selectivity filter and one in the cavity). The locations and the respective distances of the different K⁺ ions and water molecules (W), calculated within our KWKWKK sequence, also fits well with the experimental observations. The analysis of the K⁺ ions and water molecules displacements shows concerted file motions on the simulated time scale (≈1 ns), which could act as precursor to the diffusion of K⁺ ions inside the channel. A simple one-dimensional dynamical model is used to interpret the concerted motions of the ions and water molecules in the pore leading ultimately to ion transfer.     

Quantum calculations on water in the KcsA channel cavity with permeant and non-permeant ions

Different ions in the pore of the KcsA channel behave differently, and we relate this to their solvation. We show that the selectivity is dependent, in part, on the solvation in the cavity (sometimes referred to as the vestibule, it is the region containing water molecules between the intracellular gate and the selectivity filter at the extracellular end of the pore). We have shown earlier that potassium is more dependent at the upper end of the cavity region on solvation by the threonines there, while sodium ion has more water molecules as ligands. In addition, sodium ion is placed asymmetrically, while potassium is nearly exactly symmetric with respect to the four-fold symmetry of the channel. We have now extended these calculations to rubidium and cesium ions, and find that rubidium solvation resembles that of potassium (and both are permeant ions), while cesium resembles sodium (and both are non-permeant), in terms of the geometry of up to eight hydrating, and four non-hydrating, water molecules. In each case, a maximum of 12 water molecules are relevant to the calculation. The placement of the water molecules in the two cases is essentially the same as found from the electron density in the X-ray structure of Zhou and MacKinnon. For Na⁺ and K⁺, we show that energy decreases from bulk to the cavity to the lowest position in the selectivity filter (accurate energy could not be calculated for the heavier ions). A separate calculation shows that fixing the Na⁺ ion at the position of the K⁺ minimum, followed by re-optimization produced a significantly modified system, not something that could be produced by thermal fluctuations. Moving the K⁺ into the Na⁺ position in the upper cavity led to a small increase in energy, ≈ 3 k_BT, but was accompanied by large shifts in the positions of hydrating waters, which would create a major kinetic barrier. Therefore, thermal fluctuations could not invalidate the conclusions of the main calculations.     

Coupling of activation and inactivation gate in a K+‐channel: potassium and ligand sensitivity

Potassium (K⁺)‐channel gating is choreographed by a complex interplay between external stimuli, K⁺ concentration and lipidic environment. We combined solid‐state NMR and electrophysiological experiments on a chimeric KcsA–Kv1.3 channel to delineate K⁺, pH and blocker effects on channel structure and function in a membrane setting. Our data show that pH‐induced activation is correlated with protonation of glutamate residues at or near the activation gate. Moreover, K⁺ and channel blockers distinctly affect the open probability of both the inactivation gate comprising the selectivity filter of the channel and the activation gate. The results indicate that the two gates are coupled and that effects of the permeant K⁺ ion on the inactivation gate modulate activation‐gate opening. Our data suggest a mechanism for controlling coordinated and sequential opening and closing of activation and inactivation gates in the K⁺‐channel pore.     

Gating at the selectivity filter of ion channels that conduct Na+ and K+ ions

The NaK channel is a cation selective channel with similar permeability for K⁺ and Na⁺. The available crystallographic structure of wild-type (WT) NaK is usually associated with a conductive state of the channel. Here, potential of mean force for complete conduction events of Na⁺ and K⁺ ions through NaK show that: i), large energy barriers prevent the passage of ions through the WT NaK structure, ii), the barriers are correlated to the presence of a hydrogen bond between Asp-66 and Asn-68, and iii), the structure of NaK mutated to mimic cyclic nucleotide-gated channels conducts Na⁺ and K⁺. These results support the hypothesis that the filter of cation selective channels can adopt at least two different structures: a conductive one, represented by the x-ray structures of the NaK-CNG chimeras, and a closed one, represented by the x-ray structures of the WT NaK.     

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.     

Determinants of cation transport selectivity: Equilibrium binding and transport kinetics

ConclusionThe equilibrium ion-binding properties of ion channels and transporters can be difficult to discern from crystal structures alone, as proteins often adopt different lowest energy states depending on the ions bound. In cases where transport is slow, their inherent ion-binding preferences can be used to infer their transport preferences. However, in cases where transport is fast, the transport selectivity can hide their equilibrium preferences by accentuating the kinetics of ions hopping through a channel over its inherent ion-binding preferences. Thus, depending on the arrangement of ion-binding sites in a channel’s selectivity filter, one can achieve either selective or nonselective ion transport.The equilibrium K⁺ selectivity of some nonselective channels suggests a potential mechanism whereby they could evolve into a fast K⁺-selective channel. K⁺ channels and nonselective channels like CNG and HCN are related to one another in both sequence and structure, suggesting an evolutionary link between them. Swap experiments show that only a few mutations separate a nonselective channel from a K⁺-selective channel. One might imagine an evolutionary path between these channels in which the equilibrium preference for a K⁺ ion in a nonselective channel evolves into a K⁺-selective channel through these few mutations to create the selective ion queue. Alternatively, a slow single-ion channel with an equilibrium and transport preference for K⁺ ions could be transformed into a fast multi-ion channel through mutations that create a queue of K⁺-selective ion-binding sites, as is seen in most K⁺ channels studied to date.In the case of multi-ion selectivity filters, such as those found in K⁺ channels, the selectivity filter can be viewed as the active site that interacts with different queues of ions and water molecules. At least three properties emerge from multi-ion queues: (1) high conductance by reducing the affinity of multiple bound ions versus single ions; (2) high selectivity by allowing disfavored ions time to dissociate back into solution; and, consequently, (3) robust selectivity in an environment where ion concentrations can change. For transporters and carriers, the equilibrium preference and slow transport naturally create robust selectivity. In all these cases, equilibrium-based ion selectivity is achieved by slowing transport enough so that the disfavored ion is able to dissociate back into solution before transport takes place.     

An ion gating mechanism of gastric H,K-ATPase based on molecular dynamics simulations

Gastric H,K-ATPase is an electroneutral transmembrane pump that moves protons from the cytoplasm of the parietal cell into the gastric lumen in exchange for potassium ions. The mechanism of transport against the established electrochemical gradients includes intermediate conformations in which the transferred ions are trapped (occluded) within the membrane domain of the pump. The pump cycle involves switching between the E1 and E2P states. Molecular dynamics simulations on homology models of the E2P and E1 states were performed to investigate the mechanism of K⁺ movement in this enzyme. We performed separate E2P simulations with one K⁺ in the luminal channel, one K⁺ ion in the occlusion site, two K⁺ ions in the occlusion site, and targeted molecular dynamics from E2P to E1 with two K⁺ ions in the occlusion site. The models were inserted into a lipid bilayer system and were stable over the time course of the simulations, and K⁺ ions in the channel moved to a consistent location near the center of the membrane domain, thus defining the occlusion site. The backbone carbonyl oxygen from residues 337 through 342 on the nonhelical turn of M4, as well as side-chain oxygen from E343, E795, and E820, participated in the ion occlusion. A single water molecule was stably bound between the two K⁺ ions in the occlusion site, providing an additional ligand and partial shielding the positive charges from one another. Targeted molecular dynamics was used to transform the protein from the E2P to the E1 state (two K⁺ ions to the cytoplasm). This simulation identified the separation of the water column in the entry channel as the likely gating mechanism on the luminal side. A hydrated exit channel also formed on the cytoplasmic side of the occlusion site during this simulation. Hence, water molecules became available to hydrate the ions. The movement of the M1M2 transmembrane segments, and the displacement of residues Q159, E160, Q110, and T152 during the conformational change, as well as the motions of E343 and L346, acted as the cytoplasmic-side gate.     

Ionic interactions of Ba2+ blockades in the MthK K+ channel

The movement and interaction of multiple ions passing through in single file underlie various fundamental K⁺ channel properties, from the effective conduction of K⁺ ions to channel blockade by Ba²⁺ ions. In this study, we used single-channel electrophysiology and x-ray crystallography to probe the interactions of Ba²⁺ with permeant ions within the ion conduction pathway of the MthK K⁺ channel. We found that, as typical of K⁺ channels, the MthK channel was blocked by Ba²⁺ at the internal side, and the Ba²⁺-blocking effect was enhanced by external K⁺. We also obtained crystal structures of the MthK K⁺ channel pore in both Ba²⁺–Na⁺ and Ba²⁺–K⁺ environments. In the Ba²⁺–Na⁺ environment, we found that a single Ba²⁺ ion remained bound in the selectivity filter, preferably at site 2, whereas in the Ba²⁺–K⁺ environment, Ba²⁺ ions were predominantly distributed between sites 3 and 4. These ionic configurations are remarkably consistent with the functional studies and identify a molecular basis for Ba²⁺ blockade of K⁺ channels.     

Structural themes in ion channels

The recent crystal structure of the prokaryotic inwardly rectifying potassium channel, KirBac1.1, revealed for the first time the structure of a K⁺ channel in the closed state plus the location of the activation gate. Comparison of the KirBac1.1 structure with other known ion channels reveals a number of common structural features. These common characteristics include the formation of the ion conduction pathway at the interface between adjacent subunits, non-fixed charges forming part of the ion pathway, electrostatic sinks drawing ions into the channel, helix dipoles, and hydrophobic gates that ultimately prevent ion movement. This review describes in detail common structural themes present in ion channels.Presented at the Biophysical Society Meeting on Ion channels – from structure to disease held in May 2003, Rennes, France     

Hydration valve controlled non-selective conduction of Na and K in the NaK channel

The Na⁺ and K⁺ channels are essential to neural signaling, but our current knowledge at the atomic level is mainly limited to the conducting mechanism of K⁺. Unlike a K⁺ channel having four equivalent K⁺-binding sites in its selectivity filter, a NaK channel has a vestibule in the middle part of its selectivity filter, and can conduct both Na⁺ and K⁺ ions. However, the underlying mechanism for non-selective ion conduction in NaK remains elusive. Here we find four small grottos connecting with the vestibule of the NaK selectivity filter, which form a vestibule-grotto complex perpendicular to the filter pore with a few water molecules within it. It is shown that two or more of the water molecules coming to the vestibule to coordinate the cation are necessary for conducting both Na⁺ and K⁺ ions, while only one water molecule in the vestibule will obstruct ion permeation. Thus, the complex with the aid of interior water movement forms a dynamic hydration valve which is flexible in conveying different cations through the vestibule. Similar exquisite hydration valve mechanisms are expected to be utilized by other non-selective cation channels, and the results should shed new light on the importance of water in neural signaling.     

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.     

A Nanosecond Molecular Dynamics Study of Antiparallel d(G)7 Quadruplex Structures: Effect of the Coordinated Cations

Abstract

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G₇) quadruplex structures with different coordinated ions, namely Na⁺ and K⁺ ion, water and Na⁺ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of ～ 1.5 Å from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K⁺ ions are quite immobile. MD studies at 400K indicate that K⁺ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K⁺ ion coordinated structure is more favourable than that for the Na⁺ ion coordinated antiparallel quadruplex structure.     

More gating charges are needed to open a Shaker K+ channel than are needed to open an rBIIA Na+ channel

This study presents what is, to our knowledge, a novel technique by means of which the ratio of the single gating charges of voltage-gated rat brain IIA (rBIIA) sodium and Shaker potassium ion channels was estimated. In the experiment, multiple tandems of enhanced green fluorescent protein were constructed and inserted into the C-terminals of Na⁺ and K⁺ ion channels. cRNA of Na⁺ and K⁺ ion channels was injected and expressed in Xenopus laevis oocytes. The two electrode voltage-clamp technique allowed us to determine the total gating charge of sodium and potassium ion channels, while a relative measure of the amount of expressed channels could be established on the basis of the quantification of the fluorescence intensity of membrane-bound channels marked by enhanced green fluorescent proteins. As a result, gating charge and fluorescence intensity were found to be positively correlated. A relative comparison of the single gating charges of voltage-gated sodium and potassium ion channels could thus be established: the ratio of the single gating charges of the Shaker potassium channel and the rBIIA sodium channel was found to be 2.5 ± 0.4. Assuming the single channel gating charge of the Shaker K⁺ channel to be ∼13 elementary charges (well supported by other studies), this leads to approximately six elementary charges for the rBIIA sodium channel, which includes a fraction of gating charge that is missed during inactivation.     

On the Theory of Ion Transport across the Nerve Membrane: V. Two Models for the Cole-Moore K+ Hyperpolarization Delay

Two illustrative molecular models, designed to explain the Cole-Moore K⁺ hyperpolarization delay, are proposed and analyzed. Both introduce a process supplementary to the usual Hodgkin-Huxley (HH) one for a K⁺ channel. In both cases the new process becomes involved as a consequence of the conditioning hyperpolarization of the membrane and would account for the observed delay time in the K⁺ current after depolarization to near ENa. The first model uses adsorption or desorption of phospholipid molecules on the surface of the assumed protein K⁺ channel or gate. The second model involves the translocation of the charged subunits of the channel in the hyperpolarizing electric field.     

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.     

An overview of the potassium channel family

Potassium channels, tetrameric integral membrane proteins that form aqueous pores through which K⁺ can flow, are found in virtually all organisms; the genomes of humans, Drosophila, and Caenorhabditis elegans contain 30-100 K⁺ channel genes each. The structure of a bacterial K⁺ channel, sequence comparisons with other channels and electrophysiological measurements have enabled conclusions about the mechanism of gating and ion flow to be drawn for many other channels.     

Energy distribution and ion selectivity of the bacterial potassium channel

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

Selectivity of the Ca2+-activated and Light-dependent K+ Channels for Monovalent Cations

The ionic selectivity of the Ca²⁺-activated K⁺ channel of Aplysia neurons and of the light-dependent K⁺ channel of Pecten photoreceptors to metal and organic cations was studied. The selectivity sequence determined from reversal potential measurements is T1⁺ K⁺ > Rb⁺ > NH⁺₄ > Cs⁺ > Na⁺, Li⁺ and is identical to the sequence determined previously for voltage-dependent K⁺ channels in a variety of tissues. Our results suggest that some physical aspect of the K⁺ channel is conserved in phyllogenetically different tissues and cells.