Water, Protons, and the Gating of Voltage-Gated Potassium Channels.

Ion channels are ubiquitous throughout all forms of life. Potassium channels are even found in viruses. Every cell must communicate with its surroundings, so all cells have them, and excitable cells, in particular, especially nerve cells, depend on the behavior of these channels. Every channel must be open at the appropriate time, and only then, so that each channel opens in response to the stimulus that tells that channel to open. One set of channels, including those in nerve cells, responds to voltage. There is a standard model for the gating of these channels that has a section of the protein moving in response to the voltage. However, there is evidence that protons are moving, rather than protein. Water is critical as part of the gating process, although it is hard to see how this works in the standard model. Here, we review the extensive evidence of the importance of the role of water and protons in gating these channels. Our principal example, but by no means the only example, will be the Kv1.2 channel. Evidence comes from the effects of D2O, from mutations in the voltage sensing domain, as well as in the linker between that domain and the gate, and at the gate itself. There is additional evidence from computations, especially quantum calculations. Structural evidence comes from X-ray studies. The hydration of ions is critical in the transfer of ions in constricted spaces, such as the gate region and the pore of a channel; we will see how the structure of the hydrated ion fits with the structure of the channel. In addition, there is macroscopic evidence from osmotic experiments and streaming current measurements. The combined evidence is discussed in the context of a model that emphasizes the role of protons and water in gating these channels.

Protons in Gating the Kv1.2 Channel: A Calculated Set of Protonation States in Response to Polarization/Depolarization of the Channel, with the Complete Proposed Proton Path from Voltage Sensing Domain to Gate.

We have in the past proposed that proton motion constitutes the gating current in the potassium channel Kv1.2 and is responsible for the gating mechanism. For this to happen, there must be a proton path between the voltage-sensing domain (VSD) and the channel gate, and here we present quantum calculations that lead to a specific pair of proton paths, defined at the molecular level, with well-defined water molecule linkages, and with hydrogen bonding between residues; there is also at least one interpath crossover, where protons can switch paths. Quantum calculations on the entire 563-atom system give the complete geometry, the energy, and atomic charges. Calculations show that three specific residues (in the pdb 3Lut numbering, H418, E327, R326), and the T1 intracellular moiety, all of which have been shown experimentally to be involved in gating, would necessarily be protonated or deprotonated in the path between the VSD and the gate. Hydroxyl reorientation of serine and threonine residues are shown to provide a means of adjusting proton directions of motion. In the deprotonated state for K312, a low energy state, our calculations come close to reproducing the X-ray structure. The demonstration of the existence of a double proton path between VSD and gate supports the proposed proton gating mechanism; when combined with our earlier demonstration of proton generation in the VSD, and comparison with other systems that are known to move protons, we are close to achieving the definition of a complete gating mechanism in molecular detail. The coupling of the paths to the VSD, and to the PVPV section that essentially forms the gate, can be easily seen from the results of the calculation. The gate itself remains for further computations.

Quantum Calculation of Proton and Other Charge Transfer Steps in Voltage Sensing in the Kv1.2 Channel.

Ion channels in cell membranes control entry and exit of ions; their gating (opening and closing) is key to their functioning. It is known that protons can pass through the voltage-sensing domain (VSD) of channels such as Kv1.2. Quantum calculations for a section of the VSD show the steps protons take in responding to voltage and show no major displacement of the protein backbone with voltage change; 70 amino acids are included, 42 with side chains (9 directly in the proton path), 28 as backbone only, and 24 water molecules. Protons provide much of the gating current, the capacitative current immediately preceding channel opening with significant additional contributions from charge transfer to other groups. Most gating models, in contrast, require major protein displacement during gating. Energy terms without classical analogues (exchange plus correlation energy, which are greater than thermal energy) show that quantum calculations are required. Energy as a function of voltage for a key proton transfer leads to, approximately, the correct voltage for channel opening. Calculated total charge transfer (not only protons) for gating is reasonable compared to experimental values. We are also able to account, at least qualitatively, for two mutations, one with the gating current curve left-shifted and one right-shifted, and show the alternate proton paths that are required to account for these.

The Role of Proton Transport in Gating Current in a Voltage Gated Ion Channel, as Shown by Quantum Calculations.

Over two-thirds of a century ago, Hodgkin and Huxley proposed the existence of voltage gated ion channels (VGICs) to carry Na⁺ and K⁺ ions across the cell membrane to create the nerve impulse, in response to depolarization of the membrane. The channels have multiple physiological roles, and play a central role in a wide variety of diseases when they malfunction. The first channel structure was found by MacKinnon and coworkers in 1998. Subsequently, the structure of a number of VGICs was determined in the open (ion conducting) state. This type of channel consists of four voltage sensing domains (VSDs), each formed from four transmembrane (TM) segments, plus a pore domain through which ions move. Understanding the gating mechanism (how the channel opens and closes) requires structures. One TM segment (S4) has an arginine in every third position, with one such segment per domain. It is usually assumed that these arginines are all ionized, and in the resting state are held toward the intracellular side of the membrane by voltage across the membrane. They are assumed to move outward (extracellular direction) when released by depolarization of this voltage, producing a capacitive gating current and opening the channel. We suggest alternate interpretations of the evidence that led to these models. Measured gating current is the total charge displacement of all atoms in the VSD; we propose that the prime, but not sole, contributor is proton motion, not displacement of the charges on the arginines of S4. It is known that the VSD can conduct protons. Quantum calculations on the Kv1.2 potassium channel VSD show how; the key is the amphoteric nature of the arginine side chain, which allows it to transfer a proton. This appears to be the first time the arginine side chain has had its amphoteric character considered. We have calculated one such proton transfer in detail: this proton starts from a tyrosine that can ionize, transferring to the NE of the third arginine on S4; that arginine's NH then transfers a proton to a glutamate. The backbone remains static. A mutation predicted to affect the proton transfer has been qualitatively confirmed experimentally, from the change in the gating current-voltage curve. The total charge displacement in going from a normal closed potential of -70 mV across the membrane to 0 mV (open), is calculated to be approximately consistent with measured values, although the error limits on the calculation require caution in interpretation.

Quantum effects in a simple ring with hydrogen bonds.

Complexes containing multiple arginines are common in proteins. The arginines are typically salt-bridged or hydrogen-bonded, so that their charges do not repel. Here we present a quantum calculation of a ring in which the components of a salt bridge composed of a guanidinium, the arginine side chain, and a carboxylic acid are separated by water molecules. When one water molecule is displaced from the ring, atomic charges of the other water molecule, as well as other properties, are significantly affected. The exchange and correlation energy differences between optimized and displaced rings are larger than thermal energy at room temperature, and larger than the sum of other energy differences. This suggests that calculations on proteins and other systems where such a ring may occur must take quantum effects into account; charges on certain atoms shift as substituents are added to the system: another water molecule, an -OH, or -CN bonded to either moiety. Also, charge shifts accompany proton shifts from the acid to guanidinium to ionize the salt bridge. The consequences of moving one water out of the ring give evidence for electron delocalization. Bond order and atomic charges are determined using natural bond orbital calculations. The geometry of the complex changes with ionization as well as the -OH and -CN additions but not in a simple manner. These results help in understanding the role of groups of arginines in salt-bridged clusters in proteins.

Caution is required in interpretation of mutations in the voltage sensing domain of voltage gated channels as evidence for gating mechanisms.

The gating mechanism of voltage sensitive ion channels is generally considered to be the motion of the S4 transmembrane segment of the voltage sensing domains (VSD). The primary supporting evidence came from R → C mutations on the S4 transmembrane segment of the VSD, followed by reaction with a methanethiosulfonate (MTS) reagent. The cys side chain is -SH (reactive form -S-); the arginine side chain is much larger, leaving space big enough to accommodate the MTS sulfonate head group. The cavity created by the mutation has space for up to seven more water molecules than were present in wild type, which could be displaced irreversibly by the MTS reagent. Our quantum calculations show there is major reorientation of three aromatic residues that face into the cavity in response to proton displacement within the VSD. Two phenylalanines reorient sufficiently to shield/unshield the cysteine from the intracellular and extracellular ends, depending on the proton positions, and a tyrosine forms a hydrogen bond to the cysteine sulfur with its side chain -OH. These could produce the results of the experiments that have been interpreted as evidence for physical motion of the S4 segment, without physical motion of the S4 backbone. The computations strongly suggest that the interpretation of cysteine substitution reaction experiments be re-examined in the light of these considerations.

The open gate of the K(V)1.2 channel: quantum calculations show the key role of hydration.

The open gate of the Kv1.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.

Voltage gated ion channel function: gating, conduction, and the role of water and protons.

Ion channels, which are found in every biological cell, regulate the concentration of electrolytes, and are responsible for multiple biological functions, including in particular the propagation of nerve impulses. The channels with the latter function are gated (opened) by a voltage signal, which allows Na(+) into the cell and K(+) out. These channels have several positively charged amino acids on a transmembrane domain of their voltage sensor, and it is generally considered, based primarily on two lines of experimental evidence, that these charges move with respect to the membrane to open the channel. At least three forms of motion, with greatly differing extents and mechanisms of motion, have been proposed. There is a "gating current", a capacitative current preceding the channel opening, that corresponds to several charges (for one class of channel typically 12-13) crossing the membrane field, which may not require protein physically crossing a large fraction of the membrane. The coupling to the opening of the channel would in these models depend on the motion. The conduction itself is usually assumed to require the "gate" of the channel to be pulled apart to allow ions to enter as a section of the protein partially crosses the membrane, and a selectivity filter at the opposite end of the channel determines the ion which is allowed to pass through. We will here primarily consider K(+) channels, although Na(+) channels are similar. We propose that the mechanism of gating differs from that which is generally accepted, in that the positively charged residues need not move (there may be some motion, but not as gating current). Instead, protons may constitute the gating current, causing the gate to open; opening consists of only increasing the diameter at the gate from approximately 6 Å to approximately 12 Å. We propose in addition that the gate oscillates rather than simply opens, and the ion experiences a barrier to its motion across the channel that is tuned by the water present within the channel. Our own quantum calculations as well as numerous experiments of others are interpreted in terms of this hypothesis. It is also shown that the evidence that supports the motion of the sensor as the gating current can also be consistent with the hypothesis we present.

Quantum mechanical calculations on selectivity in the KcsA channel: the role of the aqueous cavity.

We have carried out quantum calculations on selected residues at the intracellular side of the selectivity filter of the KcsA potassium channel, using the published X-ray coordinates as starting points. The calculations involved primarily the side chains of residues lining the aqueous cavity on the intracellular side of the selectivity filter, in addition to water molecules, plus a K+ or Na+ ion. The results showed unambiguously that Na+ significantly distorts the symmetry of the channel at the entrance to the selectivity filter (at the residue T75), while K+ does so to a much smaller extent. In all, three ion positions have been calculated: the S4 (lowest) position at the bottom of the selectivity filter, the top of the cavity, and the midpoint of the cavity; Na+ is trapped at the cavity top, while K+ is cosolvated by the selectivity filter carbonyl groups plus threonine hydroxyl groups so that it can traverse the filter. Only one water molecule remains in the K+ solvation shell at the upper position in the cavity; this solvation shell also contains four threonine (T75) hydroxyl oxygens and two backbone carbonyls, while Na+ is solvated by five molecules of water and one oxygen from threonine hydroxyls. T75 at the entrance to the selectivity filter has a key role in recognition of the alkali ion, and T74 has secondary importance. The energetic basis for the preferential bonding of potassium by these residues is briefly discussed, based on additional calculations. Taken together, the results suggest that Na+ would have difficulty entering the cavity, and if it did, it would not be able to enter the selectivity filter.

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