Electrostatic tuning of permeation and selectivity in aquaporin water channels

Water permeation and electrostatic interactions between water and channel are investigated in the Escherichia coli glycerol uptake facilitator GlpF, a member of the aquaporin water channel family, by molecular dynamics simulations. A tetrameric model of the channel embedded in a 16:0/18:1c9-palmitoyloleylphosphatidylethanolamine membrane was used for the simulations. During the simulations, water molecules pass through the channel in single file. The movement of the single file water molecules through the channel is concerted, and we show that it can be described by a continuous-time random-walk model. The integrity of the single file remains intact during the permeation, indicating that a disrupted water chain is unlikely to be the mechanism of proton exclusion in aquaporins. Specific hydrogen bonds between permeating water and protein at the channel center (at two conserved Asp-Pro-Ala "NPA" motifs), together with the protein electrostatic fields enforce a bipolar water configuration inside the channel with dipole inversion at the NPA motifs. At the NPA motifs water-protein electrostatic interactions facilitate this inversion. Furthermore, water-water electrostatic interactions are in all regions inside the channel stronger than water-protein interactions, except near a conserved, positively charged Arg residue. We find that variations of the protein electrostatic field through the channel, owing to preserved structural features, completely explain the bipolar orientation of water. This orientation persists despite water translocation in single file and blocks proton transport. Furthermore, we find that for permeation of a cation, ion-protein electrostatic interactions are more unfavorable at the conserved NPA motifs than at the conserved Arg, suggesting that the major barrier against proton transport in aquaporins is faced at the NPA motifs.     

Enhancement of proton conductance by mutations of the selectivity filter of aquaporin-1

Prevention of cation permeation in wild-type aquaporin-1 (AQP1) is believed to be associated with the Asn-Pro-Ala (NPA) region and the aromatic/arginine selectivity filter (SF) domain. Previous work has suggested that the NPA region helps to impede proton permeation due to the protein backbone collective macrodipoles that create an environment favoring a directionally discontinuous channel hydrogen-bonded water chain and a large electrostatic barrier. The SF domain contributes to the proton permeation barrier by a spatial restriction mechanism and direct electrostatic interactions. To further explore these various effects, the free-energy barriers and the maximum cation conductance for the permeation of various cations through the AQP1-R195V and AQP1-R195S mutants are predicted computationally. The cations studied included the hydrated excess proton that utilizes the Grotthuss shuttling mechanism, a model “classical” charge localized hydronium cation that exhibits no Grotthuss shuttling, and a sodium cation. The hydrated excess proton was simulated using a specialized multi-state molecular dynamics method including a proper physical treatment of the proton shuttling and charge defect delocalization. Both AQP1 mutants exhibit a surprising cooperative effect leading to a reduction in the free-energy barrier for proton permeation around the NPA region due to altered water configurations in the SF region, with AQP1-R195S having a higher conductance than AQP1-R195V. The theoretical predictions are experimentally confirmed in wild-type AQP1 and the mutants expressed in Xenopus oocytes. The combined results suggest that the SF domain is a specialized structure that has evolved to impede proton permeation in aquaporins.     

The barrier for proton transport in aquaporins as a challenge for electrostatic models: the role of protein relaxation in mutational calculations

The origin of the barrier for proton transport through the aquaporin channel is a problem of general interest. It is becoming increasingly clear that this barrier is not attributable to the orientation of the water molecules across the channel but rather to the electrostatic penalty for moving the proton charge to the center of the channel. However, the reason for the high electrostatic barrier is still rather controversial. It has been argued by some workers that the barrier is due to the so-called NPA motif and/or to the helix macrodipole or to other specific elements. However, our works indicated that the main reason for the high barrier is the loss of the generalized solvation upon moving the proton charge from the bulk to the center of the channel and that this does not reflect a specific repulsive electrostatic interaction but the absence of sufficient electrostatic stabilization. At this stage it seems that the elucidation and clarification of the origin of the electrostatic barrier can serve as an instructive test case for electrostatic models. Thus, we reexamine the free-energy surface for proton transport in aquaporins using the microscopic free-energy perturbation/umbrella sampling (FEP/US) and the empirical valence bond/umbrella sampling (EVB/US) methods as well as the semimacroscopic protein dipole Langevin dipole model in its linear response approximation version (the PDLD/S-LRA). These extensive studies help to clarify the nature of the barrier and to establish the "reduced solvation effect" as the primary source of this barrier. That is, it is found that the barrier is associated with the loss of the generalized solvation energy (which includes of course all electrostatic effects) upon moving the proton charge from the bulk solvent to the center of the channel. It is also demonstrated that the residues in the NPA region and the helix dipole cannot be considered as the main reasons for the electrostatic barrier. Furthermore, our microscopic and semimacroscopic studies clarify the problems with incomplete alternative calculations, illustrating that the effects of various electrostatic elements are drastically overestimated by macroscopic calculations that use a low dielectric constant and do not consider the protein reorganization. Similarly, it is pointed out that microscopic potential of mean force calculations that do not evaluate the electrostatic barrier relative to the bulk water cannot be used to establish the origin of the electrostatic barrier. The relationship between the present study and calculations of pK(a)s in protein interiors is clarified, pointing out that approaches that are applied to study the aquaporin barrier should be validated by pK(a)s calculations. Such calculations also help to clarify the crucial role of solvation energies in establishing the barrier in aquaporins.     

Structural determinants of proton blockage in aquaporins

Aquaporins are an important class of membrane channels selective for water and linear polyols but impermeable to ions, including protons. Recent computational studies have revealed that the relay of protons through the water-conduction pathway of aquaporin channels is opposed by a substantial free energy barrier peaking at the signature NPA motifs. Here, free-energy simulations and continuum electrostatic calculations are combined to examine the nature and the magnitude of the contribution of specific structural elements to proton blockage in the bacterial glycerol uptake facilitator, GlpF. Potential of mean-force profiles for both hop and turn steps of structural diffusion in the narrow pore are obtained for artificial variants of the GlpF channel in which coulombic interactions between the pore contents and conserved residues Asn68 and Asn203 at the NPA signature motifs, Arg206 at the selectivity filter, and the peptidic backbone of the two half-helices M3 and M7, which are arranged in head-to-head fashion around the NPA motifs, are turned off selectively. A comparison of these results with electrostatic energy profiles for the translocation of a probe cation throughout the water permeation pathway indicates that the free-energy profile for proton movement inside the narrow pore is dominated by static effects arising from the distribution of charged and polar groups of the channel, whereas dielectric effects contribute primarily to opposing the access of H+ to the pore mouths (desolvation penalty). The single most effective way to abolish the free-energy gradients opposing the movement of H+ around the NPA motif is to turn off the dipole moments of helices M3 and M7. Mutation of either of the two NPA Asn residues to Asp compensates for charge-dipole and dipole-dipole effects opposing the hop and turn steps of structural diffusion, respectively, and dramatically reduces the free energy barrier of proton translocation, suggesting that these single mutants could leak protons.     

Proton exclusion by an aquaglyceroprotein: a voltage clamp study

BACKGROUND INFORMATION: In silico both orthodox aquaporins and aquaglyceroporins are shown to exclude protons. Supporting experimental evidence is available only for orthodox aquaporins. In contrast, the subset of the aquaporin water channel family that is permeable to glycerol and certain small, uncharged solutes has not yet been shown to exclude protons. Moreover, different aquaglyceroporins have been reported to conduct ions when reconstituted in planar bilayers. RESULTS: To clarify these discrepancies, we have measured proton permeability through the purified Escherichia coli glycerol facilitator (GlpF). Functional reconstitution into planar lipid bilayers was demonstrated by imposing an osmotic gradient across the membrane and detecting the resulting small changes in ionic concentration close to the membrane surface. The osmotic water flow corresponds to a GlpF single channel water permeability of 0.7x10(-14) cm(3).subunit(-1).s(-1). Proton conductivity measurements carried out in the presence of a pH gradient (1 unit) revealed an upper limit of the H(+) (OH(-)) to H(2)O molecules transport stoichiometry of 2x10(-9). A significant GlpF-mediated ion conductivity was also not detectable. CONCLUSIONS: The lack of a physiologically relevant GlpF-mediated proton conductivity agrees well with predictions made by molecular dynamics simulations.     

On the origin of the electrostatic barrier for proton transport in aquaporin

The nature of the electrostatic barrier for proton transport in aquaporins is analyzed by semimacroscopic and microscopic models. It is found that the barrier is associated with the loss of the generalized solvation energy upon moving from the bulk solvent to the center of the channel. It is clarified that our solvation concept includes the effect of the protein polar groups and ionized residues. The nature of the contributions to the solvation barrier is examined by using the linear response approximation. It is found that the residues in the NPA region contribute much less than what would be deduced from calculations that do not consider the protein reorganization. It is clarified that the contributions of different structural or electrostatic elements to the solvation barrier can be established by removing these elements and examining the corresponding effect on the barrier height. Using this definition and “mutating” the NPA residues to their non-polar analogues establishes that these residues do not provide the major contribution to the solvation barrier.     

Oxygen and proton pathways in cytochrome c oxidase

Cytochrome c oxidase is a redox-driven proton pump, which couples the reduction of oxygen to water to the translocation of protons across the membrane. The recently solved x-ray structures of cytochrome c oxidase permit molecular dynamics simulations of the underlying transport processes. To eventually establish the proton pump mechanism, we investigate the transport of the substrates, oxygen and protons, through the enzyme. Molecular dynamics simulations of oxygen diffusion through the protein reveal a well-defined pathway to the oxygen-binding site starting at a hydrophobic cavity near the membrane-exposed surface of subunit I, close to the interface to subunit III. A large number of water sites are predicted within the protein, which could play an essential role for the transfer of protons in cytochrome c oxidase. The water molecules form two channels along which protons can enter from the cytoplasmic (matrix) side of the protein and reach the binuclear center. A possible pumping mechanism is proposed that involves a shuttling motion of a glutamic acid side chain, which could then transfer a proton to a propionate group of heme α₃. Proteins 30:100–107, 1998. © 1998 Wiley-Liss, Inc.     

Concerted action of two cation filters in the aquaporin water channel

Aquaporin (AQP) facilitated water transport is common to virtually all cell membranes and is marked by almost perfect specificity and high flux rates. Simultaneously, protons and cations are strictly excluded to maintain ionic transmembrane gradients. Yet, the AQP cation filters have not been identified experimentally. We report that three point mutations turned the water-specific AQP1 into a proton/alkali cation channel with reduced water permeability and the permeability sequence: H⁺ ≫K⁺ >Rb⁺ >Na⁺ >Cs⁺ >Li⁺. Contrary to theoretical models, we found that electrostatic repulsion at the central asn-pro-ala (NPA) region does not suffice to exclude protons. Full proton exclusion is reached only in conjunction with the aromatic/arginine (ar/R) constriction at the pore mouth. In contrast, alkali cations are blocked by the NPA region but leak through the ar/R constriction. Expression of alkali-leaking AQPs depolarized membrane potentials and compromised cell survival. Our results hint at the alkali-tight but solute-unselective NPA region as a feature of primordial channels and the proton-tight and solute-selective ar/R constriction variants as later adaptations within the AQP superfamily.     

Does CO2 permeate through aquaporin-1?

Aquaporins facilitate water permeation across biological membranes. Additionally, glycerol and other small neutral solutes are permeated by related aquaglyceroporins. The role of aquaporins in gas permeation has been a long-standing and controversially discussed issue. We present an extensive set of atomistic molecular dynamics simulations that address the question of CO(2) permeation through human aquaporin-1. Free energy profiles derived from the simulations display a barrier of approximately 23 kJ/mol in the aromatic/arginine constriction region of the water pore, whereas a barrier of approximately 4 kJ/mol was observed for a palmitoyloleoylphosphatidylethanolamine lipid bilayer membrane. The results indicate that significant aquaporin-1-mediated CO(2) permeation is to be expected only in membranes with a low intrinsic CO(2) permeability.     

Water and ion permeation in bAQP1 and GlpF channels: a kinetic Monte Carlo study

The kinetic Monte Carlo reaction-path-following technique is applied to determine the lowest-energy water pathway and the coordinating amino acids in bAQP1 and GlpF channels, both treated as rigid. In bAQP1, water molecules pass through the pore between the asparagine-proline-alanine (NPA) and selectivity filter (SF) sites one at a time. The water chain is interrupted at the SF where one water forms three stable hydrogen bonds with protein atoms. In this SF, water's conformation depends on the protonation locus of H182. In GlpF, two water molecules bond simultaneously to the NPA asparagines and pass through the SF in zigzag fashion. No water single-file forms in rigid GlpF. To accommodate a single file of waters requires narrowing the GlpF pore. Our results reveal that in both proteins a proposed bipolar water arrangement is thermally disrupted in the NPA region, especially in the cytoplasmic part of the pore. The equilibrium hydrogen-bonded chain is occasionally interrupted in the hydrophobic zones adjacent to the NPA motifs. The permeation of alkali cations through bAQP1 and GlpF is barred due to a large free-energy barrier in the NPA region as well as a large energy barrier blocking entry from the cytoplasm. Permeation of halides is prevented due to two large energy barriers in the cytoplasmic and periplasmic pores as well as a large free-energy barrier barring entry from the periplasm. Our results, based on modeling charge permeation, support an electrostatic rather than orientational basis for proton exclusion. Binding within the aquaporin pore cannot compensate sufficiently for dehydration of the protonic charge; there is also an electrostatic barrier in the NPA region blocking proton transport. The highly ordered single file of waters, which is drastically interrupted at the SF of bAQP1, may also contribute to proton block.     

Dynamics and energetics of water permeation through the aquaporin channel

Structural properties of water inside bovine aquaporin-1 are investigated by molecular simulation. The calculations, which are based on the recently determined X-ray structure at 2.2 A resolution (Sui et al., Nature 2001;414:872-878), are carried out on one monomeric subunit immersed in a water-n-octane-water bilayer. Molecular dynamics (MD) simulations suggest that His182, a fully conserved residue in the channel pore, is protonated in the delta position. Furthermore, they reveal a highly ordered water structure in the channel, induced by the electrostatic properties of the protein. Multiple-steering MD simulations are used to calculate the free-energy of water diffusion. To the best of our knowledge, this represents the first free-energy calculation based on the new, high-resolution structure of the pore. The calculated barrier is 2.5 kcal/mol, and it is associated to water permeation through the Asn-Pro-Ala (NPA) region of the pore, where water molecules are only hydrogen-bonded with themselves. These findings are fully consistent with those based on the previous MD studies on the human protein (de Groot and Grubmüller, Science 2001;294:2353-2357).     

Dynamic water networks in cytochrome cbb(3) oxidase

Heme-copper oxidases (HCOs) are terminal electron acceptors in aerobic respiration. They catalyze the reduction of molecular oxygen to water with concurrent pumping of protons across the mitochondrial and bacterial membranes. Protons required for oxygen reduction chemistry and pumping are transferred through proton uptake channels. Recently, the crystal structure of the first C-type member of the HCO superfamily was resolved [Buschmann et al. Science 329 (2010) 327-330], but crystallographic water molecules could not be identified. Here we have used molecular dynamics (MD) simulations, continuum electrostatic approaches, and quantum chemical cluster calculations to identify proton transfer pathways in cytochrome cbb(3). In MD simulations we observe formation of stable water chains that connect the highly conserved Glu323 residue on the proximal side of heme b(3) both with the N- and the P-sides of the membrane. We propose that such pathways could be utilized for redox-coupled proton pumping in the C-type oxidases. Electrostatics and quantum chemical calculations suggest an increased proton affinity of Glu323 upon reduction of high-spin heme b(3). Protonation of Glu323 provides a mechanism to tune the redox potential of heme b(3) with possible implications for proton pumping.     

Molecular basis of proton blockage in aquaporins

Water transport channels in membrane proteins of the aquaporin superfamily are impermeable to ions, including H+ and OH-. We examine the molecular basis for the blockage of proton translocation through the single-file water chain in the pore of a bacterial aquaporin, GlpF. We compute the reversible thermodynamic work for the two complementary steps of the Grotthuss "hop-and-turn" relay mechanism: consecutive transfers of H+ along the hydrogen-bonded chain (hop) and conformational reorganization of the chain (turn). In the absence of H+, the strong preference for the bipolar orientation of water around the two Asn-Pro-Ala (NPA) motifs lining the pore over both unidirectional polarization states of the chain precludes the reorganization of the hydrogen-bonded network. Inversely, translocation of an excess proton in either direction is opposed by a free-energy barrier centered at the NPA region. Both hop and turn steps of proton translocation are opposed by the electrostatic field of the channel.     

Ion exclusion mechanism in aquaporin at an atomistic level

Although the mechanism of proton exclusion in aquaporin is investigated by many researchers, the detailed molecular mechanism for ion exclusion in aquaporin is still not completely understood. In the present work, a detailed mechanism for ion exclusion in aquaporin-1 (AQP1) at an atomistic level is investigated by calculating the free energy for transport of ions in AQP1 using an atomistic molecular dynamics simulation. For this purpose, sodium and chloride ions are chosen as representatives for nonprotonic ions. The simulation shows that the free energy barrier showing its maximum is located at the NPA region for sodium ion while it is located at both the front and the rear for chloride ion and that the barrier height is 18 and 9 kcal/mol, respectively, indicating that the ions are not able to pass through aquaporin. Analysis of the pair interaction energy between the permeating ion and its environment reveals that sodium ion is excluded by the positive charge generated by two alpha-helical macro-dipoles, while chloride ion is expelled by carbonyl oxygen atoms protruding from pore-making residues before it reaches the NPA motif. It is also found that the number of water molecules hydrating the ions is reduced as the ions enter the pore, implying that the energetic cost for detaching water molecules from a permeating ion also contributes to the free energy barriers of ion transport in AQP1.     

Molecular dynamics of a protein surface: ion-residues interactions

Time-resolved measurements indicated that protons could propagate on the surface of a protein or a membrane by a special mechanism that enhanced the shuttle of the proton toward a specific site. It was proposed that a suitable location of residues on the surface contributes to the proton shuttling function. In this study, this notion was further investigated by the use of molecular dynamics simulations, where Na(+) and Cl(-) are the ions under study, thus avoiding the necessity for quantum mechanical calculations. Molecular dynamics simulations were carried out using as a model a few Na(+) and Cl(-) ions enclosed in a fully hydrated simulation box with a small globular protein (the S6 of the bacterial ribosome). Three independent 10-ns-long simulations indicated that the ions and the protein's surface were in equilibrium, with rapid passage of the ions between the protein's surface and the bulk. However, it was noted that close to some domains the ions extended their duration near the surface, thus suggesting that the local electrostatic potential hindered their diffusion to the bulk. During the time frame in which the ions were detained next to the surface, they could rapidly shuttle between various attractor sites located under the electrostatic umbrella. Statistical analysis of the molecular dynamics and electrostatic potential/entropy consideration indicated that the detainment state is an energetic compromise between attractive forces and entropy of dilution. The similarity between the motion of free ions next to a protein and the proton transfer on the protein's surface are discussed.     

A molecular modeling approach defines a new group of Nodulin 26-like aquaporins in plants

The three-dimensional models built for the Nod26-like aquaporins all exhibit the typical α-helical fold of other aquaporins containing the two ar/R and NPA constriction filters along the central water channel. Besides these structural homologies, they readily differ with respect to the amino acid residues forming the ar/R selective filter. According to these discrepancies in both the hydrophilicity and pore size of the ar/R filter, Nod26-like aquaporins can be distributed in three subgroups corresponding to NIP-1, NIP-II and a third subgroup of Nod26-like aquaporins exhibiting a highly hydrophilic and widely open filter. However, all Nod26-like aquaporins display a bipartite distribution of electrostatic charges along the water channel with an electropositive extracellular vestibular portion followed by an electronegative cytosolic vestibular portion. The specific transport of water, non-ionic solutes (glycerol, urea, ammoniac), ions and gas (NH₃) across the Nod26-like obviously depends on the electrostatic and conformational properties of their central water channel.     

Protons @ interfaces: implications for biological energy conversion

The review focuses on the anisotropy of proton transfer at the surface of biological membranes. We consider (i) the data from "pulsed" experiments, where light-triggered enzymes capture or eject protons at the membrane surface, (ii) the electrostatic properties of water at charged interfaces, and (iii) the specific structural attributes of proton-translocating enzymes. The pulsed experiments revealed that proton exchange between the membrane surface and the bulk aqueous phase takes as much as about 1 ms, but could be accelerated by added mobile pH-buffers. Since the accelerating capacity of the latter decreased with the increase in their electric charge, it was concluded that the membrane surface is separated from the bulk aqueous phase by a barrier of electrostatic nature. The barrier could arise owing to the water polarization at the negatively charged membrane surface. The barrier height depends linearly on the charge of penetrating ions; for protons, it has been estimated as about 0.12 eV. While the proton exchange between the surface and the bulk aqueous phase is retarded by the interfacial barrier, the proton diffusion along the membrane, between neighboring enzymes, takes only microseconds. The proton spreading over the membrane is facilitated by the hydrogen-bonded networks at the surface. The membrane-buried layers of these networks can eventually serve as a storage/buffer for protons (proton sponges). As the proton equilibration between the surface and the bulk aqueous phase is slower than the lateral proton diffusion between the "sources" and "sinks", the proton activity at the membrane surface, as sensed by the energy transducing enzymes at steady state, might deviate from that measured in the adjoining water phase. This trait should increase the driving force for ATP synthesis, especially in the case of alkaliphilic bacteria.     

On the definition, nomenclature and classification of water channel proteins (aquaporins and relatives)

A water channel protein (WCP) or a water channel can be defined as a transmembrane protein that has a specific three-dimensional structure with a pore that provides a pathway for water permeation across biological membranes. The pore is formed by two highly conserved regions in the amino acid sequence, called NPA boxes (or motifs) with three amino acid residues (asparagine-proline-alanine, NPA) and several surrounding amino acids. The NPA boxes have been called the "signature" sequence of WCPs. WCPs are a family of proteins belonging to the Membrane Intrinsic Proteins (MIPs) superfamily. In addition, in the MIP superfamily (with more than 1000 members) there are also proteins with no channel activity. The WCP family include three subfamilies: aquaporins, aquaglyceroporins and S-aquaporins. (1) The aquaporins (AQPs) are water selective or specific water channels, also named by various authors as "orthodox", "ordinary", "conventional", "classical", "pure", "normal", or "sensu strictu" aquaporins); (2) The aquaglyceroporins are permeable to water, but also to other small uncharged molecules, in particular glycerol; this family includes the glycerol facilitators, abbreviated as GlpFs, from glycerol permease facilitators. The "signature" sequence for aquaglyceroporins is the aspartic acid residue (D) in the second NPA box. (3) The third subfamily of WCPs have little conserved amino acid sequences around the NPA boxes, unclassifiable to the first two subfamilies. I recommend to use always for this subfamily the name S-aquaporins. They are also named "superaquaporins", "aquaporins with unusual (or deviated) NPA boxes", "subcellular aquaporins", or "sip-like aquaporins". I also recommend to use always the spelling aquaporin (not aquaporine), and, for various AQPs, the abbreviation AQP followed immediately by the number, (e.g. AQP1), with no space or - which might create confusions with "minus".     

The mechanism of glycerol conduction in aquaglyceroporins.

BACKGROUND: The E. coli glycerol facilitator, GlpF, selectively conducts glycerol and water, excluding ions and charged solutes. The detailed mechanism of the glycerol conduction and its relationship to the characteristic secondary structure of aquaporins and to the NPA motifs in the center of the channel are unknown. RESULTS: Molecular dynamics simulations of GlpF reveal spontaneous glycerol and water conduction driven, on a nanosecond timescale, by thermal fluctuations. The bidirectional conduction, guided and facilitated by the secondary structure, is characterized by breakage and formation of hydrogen bonds for which water and glycerol compete. The conduction involves only very minor changes in the protein structure, and cooperativity between the GlpF monomers is not evident. The two conserved NPA motifs are strictly linked together by several stable hydrogen bonds and their asparagine side chains form hydrogen bonds with the substrates passing the channel in single file. CONCLUSIONS: A complete conduction of glycerol through the GlpF was deduced from molecular dynamics simulations, and key residues facilitating the conduction were identified. The nonhelical parts of the two half-membrane-spanning segments expose carbonyl groups towards the channel interior, establishing a curve-linear pathway. The conformational stability of the NPA motifs is important in the conduction and critical for selectivity. Water and glycerol compete in a random manner for hydrogen bonding sites in the protein, and their translocations in single file are correlated. The suggested conduction mechanism should apply to the whole family.     

Computer simulation of ion channel gating: the M(2) channel of influenza A virus in a lipid bilayer

The transmembrane fragment of the influenza virus M(2) protein forms a homotetrameric channel that transports protons. In this paper, we use molecular dynamics simulations to help elucidate the mechanism of channel gating by four histidines that occlude the channel lumen in the closed state. We test two competing hypotheses. In the "shuttle" mechanism, the delta nitrogen atom on the extracellular side of one histidine is protonated by the incoming proton, and, subsequently, the proton on the epsilon nitrogen atom is released on the opposite side. In the "water-wire" mechanism, the gate opens because of electrostatic repulsion between four simultaneously biprotonated histidines. This allows for proton transport along the water wire that penetrates the gate. For each system, composed of the channel embedded in a hydrated phospholipid bilayer, a 1.3-ns trajectory was obtained. It is found that the states involved in the shuttle mechanism, which contain either single-protonated histidines or a mixture of single-protonated histidines plus one biprotonated residue, are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfer. In contrast, the fully biprotonated state is not stable. Additional simulations show that if only two histidines are biprotonated, the channel deforms but the gate remains closed. These results support the shuttle mechanism but not the gate-opening mechanism of proton gating in M(2).