Water wires in atomistic models of the Hv1 proton channel

Omega currents in voltage-gated ion channels: what can we learn from uncovering the voltage-sensing mechanism using MD simulations?

Ion channels conduct charged species through otherwise impermeable biological membranes. Their activity supports a number of physiological processes, and genetic mutations can disrupt their function dramatically. Among these channels, voltage gated cation channels (VGCCs) are ubiquitous transmembrane proteins involved in electrical signaling. In addition to their selectivity for ions, their function requires membrane-polarization-dependent gating. Triggered by changes in the transmembrane voltage, the activation and deactivation of VGCCs proceed through a sensing mechanism that prompts motion of conserved positively charged (basic) residues within the S4 helix of a four-helix bundle, the voltage sensor domain (VSD). Decades of experimental investigations, using electrophysiology, molecular biology, pharmacology, and spectroscopy, have revealed details about the function of VGCCs. However, in 2005, the resolution of the crystal structure of the activated state of one member of the mammalian voltage gated potassium (Kv) channels family (the Kv1.2) enabled researchers to make significant progress in understanding the structure-function relationship in these proteins on a molecular level. In this Account, we review the use of a complementary technique, molecular dynamics (MD) simulations, that has offered new insights on this timely issue. Starting from the "open-activated state" crystal structure, we have carried out large-scale all atom MD simulations of the Kv1.2 channel embedded in its lipidic environment and submitted to a hyperpolarizing (negative) transmembrane potential. We then used steered MD simulations to complete the full transition to the resting-closed state. Using these procedures, we have followed the operation of the VSDs and uncovered three intermediate states between their activated and deactivated conformations. Each conformational state is characterized by its network of salt bridges and by the occupation of the gating charge transfer center by a specific S4 basic residue. Overall, the global deactivation mechanism that we have uncovered agrees with proposed kinetic models and recent experimental results that point towards the presence of several intermediate states. The understanding of these conformations has allowed us to examine how mutations of the S4 basic residues analogous to those involved in genetic diseases affect the function of VGCCs. In agreement with electrophysiology experiments, mutations perturb the VSD structure and trigger the appearance of state-dependent "leak" currents. The simulation results unveil the key elementary molecular processes involved in these so-called "omega" currents. We generalize these observations to other members of the VGCC family, indicating which type of residues may generate such currents and which conditions might cause leaks that prevent proper function of the channel. Today, the understanding of the intermediate state conformations enables researchers to confidently tackle other key questions such as the mode of action of toxins or modulation of channel function by lipids.

Peregrination of the selectivity filter delineates the pore of the human voltage-gated proton channel hHV1

Extraordinary selectivity is crucial to all proton-conducting molecules, including the human voltage-gated proton channel (hH_V1), because the proton concentration is >10⁶ times lower than that of other cations. Here we use “selectivity filter scanning” to elucidate the molecular requirements for proton-specific conduction in hH_V1. Asp¹¹², in the middle of the S1 transmembrane helix, is an essential part of the selectivity filter in wild-type (WT) channels. After neutralizing Asp¹¹² by mutating it to Ala (D112A), we introduced Asp at each position along S1 from 108 to 118, searching for “second site suppressor” activity. Surprisingly, most mutants lacked even the anion conduction exhibited by D112A. Proton-specific conduction was restored only with Asp or Glu at position 116. The D112V/V116D channel strikingly resembled WT in selectivity, kinetics, and ΔpH-dependent gating. The S4 segment of this mutant has similar accessibility to WT in open channels, because R211H/D112V/V116D was inhibited by internally applied Zn²⁺. Asp at position 109 allowed anion permeation in combination with D112A but did not rescue function in the nonconducting D112V mutant, indicating that selectivity is established externally to the constriction at F150. The three positions that permitted conduction all line the pore in our homology model, clearly delineating the conduction pathway. Evidently, a carboxyl group must face the pore directly to enable conduction. Molecular dynamics simulations indicate reorganization of hydrogen bond networks in the external vestibule in D112V/V116D. At both positions where it produces proton selectivity, Asp frequently engages in salt linkage with one or more Arg residues from S4. Surprisingly, mean hydration profiles were similar in proton-selective, anion-permeable, and nonconducting constructs. That the selectivity filter functions in a new location helps to define local environmental features required to produce proton-selective conduction.

Mapping of sites facing aqueous environment of voltage-gated proton channel at resting state: a study with PEGylation protection

Hv1 (also named, voltage-sensor only protein, VSOP) lacks an authentic pore domain, and its voltage sensor domain plays both roles in voltage sensing and proton permeation. The activities of a proton channel are intrinsic to protomers of Hv1, while Hv1 is dimeric in biological membranes; cooperative gating is exerted by interaction between two protomers. As the signature pattern conserved among voltage-gated channels and voltage-sensing phosphatase, Hv1 has multiple arginines intervened by two hydrophobic residues on the fourth transmembrane segment, S4. S4 moves upward relative to other helices upon depolarization, causing conformational change possibly leading to the formation of a proton-selective conduction pathway. However, detailed mechanisms of proton-selectivity and gating of Hv1 are unknown. Here we took an approach of PEGylation protection assay to define residues facing the aqueous environment of mouse Hv1 (mHv1). Accessibilities of two maleimide molecules, N-ethylmaleimide (NEM) and 4-acetamido-4'-maleimidylstilbene-2,2'-disulfonic acid (AMS), were examined on cysteine introduced into individual sites. Only the first arginine on S4 (R1: R201) was inaccessible by NEM and AMS in mHv1. This is consistent with previous results of electrophysiology on the resting state channel, suggesting that the accessibility profile represents the resting state of mHv1. D108, critical for proton selectivity, was accessible by AMS and NEM, suggesting that D108 faces the vestibule. F146, a site critical for blocking by a guanidinium-reagent, was accessible by NEM, suggesting that F146 also faces the inner vestibule. These findings suggest an inner vestibule lined by several residues on S2 including F146, D108 on S1, and the C-terminal half of S4.

Construction and validation of a homology model of the human voltage-gated proton channel hHV1

The topological similarity of voltage-gated proton channels (H_V1s) to the voltage-sensing domain (VSD) of other voltage-gated ion channels raises the central question of whether H_V1s have a similar structure. We present the construction and validation of a homology model of the human H_V1 (hH_V1). Multiple structural alignment was used to construct structural models of the open (proton-conducting) state of hH_V1 by exploiting the homology of hH_V1 with VSDs of K⁺ and Na⁺ channels of known three-dimensional structure. The comparative assessment of structural stability of the homology models and their VSD templates was performed using massively repeated molecular dynamics simulations in which the proteins were allowed to relax from their initial conformation in an explicit membrane mimetic. The analysis of structural deviations from the initial conformation based on up to 125 repeats of 100-ns simulations for each system reveals structural features consistently retained in the homology models and leads to a consensus structural model for hH_V1 in which well-defined external and internal salt-bridge networks stabilize the open state. The structural and electrostatic properties of this open-state model are compatible with proton translocation and offer an explanation for the reversal of charge selectivity in neutral mutants of Asp¹¹². Furthermore, these structural properties are consistent with experimental accessibility data, providing a valuable basis for further structural and functional studies of hH_V1. Each Arg residue in the S4 helix of hH_V1 was replaced by His to test accessibility using Zn²⁺ as a probe. The two outermost Arg residues in S4 were accessible to external solution, whereas the innermost one was accessible only to the internal solution. Both modeling and experimental data indicate that in the open state, Arg²¹¹, the third Arg residue in the S4 helix in hH_V1, remains accessible to the internal solution and is located near the charge transfer center, Phe¹⁵⁰.

Voltage-sensing domain of voltage-gated proton channel Hv1 shares mechanism of block with pore domains

Voltage-gated sodium, potassium, and calcium channels are made of a pore domain (PD) controlled by four voltage-sensing domains (VSDs). The PD contains the ion permeation pathway and the activation gate located on the intracellular side of the membrane. A large number of small molecules are known to inhibit the PD by acting as open channel blockers. The voltage-gated proton channel Hv1 is made of two VSDs and lacks the PD. The location of the activation gate in the VSD is unknown and open channel blockers for VSDs have not yet been identified. Here, we describe a class of small molecules which act as open channel blockers on the Hv1 VSD and find that a highly conserved phenylalanine in the charge transfer center of the VSD plays a key role in blocker binding. We then use one of the blockers to show that Hv1 contains two intracellular and allosterically coupled gates.

Molecular dynamics simulations reveal proton transfer pathways in cytochrome C-dependent nitric oxide reductase

Nitric oxide reductases (NORs) are membrane proteins that catalyze the reduction of nitric oxide (NO) to nitrous oxide (N₂O), which is a critical step of the nitrate respiration process in denitrifying bacteria. Using the recently determined first crystal structure of the cytochrome c-dependent NOR (cNOR) [Hino T, Matsumoto Y, Nagano S, Sugimoto H, Fukumori Y, et al. (2010) Structural basis of biological N2O generation by bacterial nitric oxide reductase. Science 330: 1666–70.], we performed extensive all-atom molecular dynamics (MD) simulations of cNOR within an explicit membrane/solvent environment to fully characterize water distribution and dynamics as well as hydrogen-bonded networks inside the protein, yielding the atomic details of functionally important proton channels. Simulations reveal two possible proton transfer pathways leading from the periplasm to the active site, while no pathways from the cytoplasmic side were found, consistently with the experimental observations that cNOR is not a proton pump. One of the pathways, which was newly identified in the MD simulation, is blocked in the crystal structure and requires small structural rearrangements to allow for water channel formation. That pathway is equivalent to the functional periplasmic cavity postulated in cbb₃ oxidase, which illustrates that the two enzymes share some elements of the proton transfer mechanisms and confirms a close evolutionary relation between NORs and C-type oxidases. Several mechanisms of the critical proton transfer steps near the catalytic center are proposed.

Denitrification is an anaerobic process performed by several bacteria as an alternative to aerobic respiration. A key intermediate step is catalyzed by the nitric oxide reductase (NOR) enzyme, which is situated in the cytoplasmic membrane. Proton delivery to the catalytic site inside NOR is an important part of its functioning. In this work we use molecular dynamics simulations to describe water distribution and to identify proton transfer pathways in cNOR. Our results reveal two channels from the periplasmic side of the membrane and none from the cytoplasmic side, indicating that cNOR is not a proton pump. It is our hope that these results will provide a basis for further experimental and computational studies aimed to understand details of the NOR mechanism. Furthermore, this work sheds light on the molecular evolution of respiratory enzymes.

Molecular dynamics simulations of voltage-gated cation channels: insights on voltage-sensor domain function and modulation

Since their discovery in the 1950s, the structure and function of voltage-gated cation channels (VGCC) has been largely understood thanks to results stemming from electrophysiology, pharmacology, spectroscopy, and structural biology. Over the past decade, computational methods such as molecular dynamics (MD) simulations have also contributed, providing molecular level information that can be tested against experimental results, thereby allowing the validation of the models and protocols. Importantly, MD can shed light on elements of VGCC function that cannot be easily accessed through “classical” experiments. Here, we review the results of recent MD simulations addressing key questions that pertain to the function and modulation of the VGCC’s voltage-sensor domain (VSD) highlighting: (1) the movement of the S4-helix basic residues during channel activation, articulating how the electrical driving force acts upon them; (2) the nature of the VSD intermediate states on transitioning between open and closed states of the VGCC; and (3) the molecular level effects on the VSD arising from mutations of specific S4 positively charged residues involved in certain genetic diseases.

Voltage-gated proton channels

Voltage-gated proton channels, HV1, have vaulted from the realm of the esoteric into the forefront of a central question facing ion channel biophysicists, namely, the mechanism by which voltage-dependent gating occurs. This transformation is the result of several factors. Identification of the gene in 2006 revealed that proton channels are homologues of the voltage-sensing domain of most other voltage-gated ion channels. Unique, or at least eccentric, properties of proton channels include dimeric architecture with dual conduction pathways, perfect proton selectivity, a single-channel conductance approximately 10(3) times smaller than most ion channels, voltage-dependent gating that is strongly modulated by the pH gradient, ΔpH, and potent inhibition by Zn(2+) (in many species) but an absence of other potent inhibitors. The recent identification of HV1 in three unicellular marine plankton species has dramatically expanded the phylogenetic family tree. Interest in proton channels in their own right has increased as important physiological roles have been identified in many cells. Proton channels trigger the bioluminescent flash of dinoflagellates, facilitate calcification by coccolithophores, regulate pH-dependent processes in eggs and sperm during fertilization, secrete acid to control the pH of airway fluids, facilitate histamine secretion by basophils, and play a signaling role in facilitating B-cell receptor mediated responses in B-lymphocytes. The most elaborate and best-established functions occur in phagocytes, where proton channels optimize the activity of NADPH oxidase, an important producer of reactive oxygen species. Proton efflux mediated by HV1 balances the charge translocated across the membrane by electrons through NADPH oxidase, minimizes changes in cytoplasmic and phagosomal pH, limits osmotic swelling of the phagosome, and provides substrate H(+) for the production of H2O2 and HOCl, reactive oxygen species crucial to killing pathogens.

Evolution of the voltage sensor domain of the voltage-sensitive phosphoinositide phosphatase VSP/TPTE suggests a role as a proton channel in eutherian mammals

The voltage-sensitive phosphoinositide phosphatases provide a mechanism to couple changes in the transmembrane electrical potential to intracellular signal transduction pathways. These proteins share a domain architecture that is conserved in deuterostomes. However, gene duplication events in primates, including humans, give rise to the paralogs TPTE and TPTE2 that retain protein domain organization but, in the case of TPTE, have lost catalytic activity. Here, we present evidence that these human proteins contain a functional voltage sensor, similar to that in nonmammalian orthologs. However, domains of these human proteins can also generate a noninactivating outward current that is not observed in zebra fish or tunicate orthologs. This outward current has the anticipated characteristics of a voltage-sensitive proton current and is due to the appearance of a single histidine residue in the S4 transmembrane segment of the voltage sensor. Histidine is observed at this position only during the eutherian radiation. Domains from both human paralogs generate proton currents. This apparent gain of proton channel function during the evolution of the TPTE protein family may account for the conservation of voltage sensor domains despite the loss of phosphatase activity in some human paralogs.

Voltage sensor of ion channels and enzymes

Placed in the cell membrane (a two-dimensional environment), ion channels and enzymes are able to sense voltage. How these proteins are able to detect the difference in the voltage across membranes has attracted much attention, and at times, heated debate during the last few years. Sodium, Ca2+ and K+ voltage-dependent channels have a conserved positively charged transmembrane (S4) segment that moves in response to changes in membrane voltage. In voltage-dependent channels, S4 forms part of a domain that crystallizes as a well-defined structure consisting of the first four transmembrane (S1-S4) segments of the channel-forming protein, which is defined as the voltage sensor domain (VSD). The VSD is tied to a pore domain and VSD movements are allosterically coupled to the pore opening to various degrees, depending on the type of channel. How many charges are moved during channel activation, how much they move, and which are the molecular determinants that mediate the electromechanical coupling between the VSD and the pore domains are some of the questions that we discuss here. The VSD can function, however, as a bona fide proton channel itself, and, furthermore, the VSD can also be a functional part of a voltage-dependent phosphatase.

The pore of the voltage-gated proton channel

In classical tetrameric voltage-gated ion channels four voltage-sensing domains (VSDs), one from each subunit, control one ion permeation pathway formed by four pore domains. The human Hv1 proton channel has a different architecture, containing a VSD, but lacking a pore domain. Since its location is not known, we searched for the Hv permeation pathway. We find that mutation of the S4 segment's third arginine R211 (R3) compromises proton selectivity, enabling conduction of a metal cation and even of the large organic cation guanidinium, reminiscent of Shaker's omega pore. In the open state, R3 appears to interact with an aspartate (D112) that is situated in the middle of S1 and is unique to Hv channels. The double mutation of both residues further compromises cation selectivity. We propose that membrane depolarization reversibly positions R3 next to D112 in the transmembrane VSD to form the ion selectivity filter in the channel's open conformation.