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.     

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.     

Requirement for asparagine in the aquaporin NPA sequence signature motifs for cation exclusion

Two highly conserved NPA motifs are a hallmark of the aquaporin (AQP) family. The NPA triplets form N-terminal helix capping structures with the Asn side chains located in the centre of the water or solute-conducting channel, and are considered to play an important role in AQP selectivity. Although another AQP selectivity filter site, the aromatic/Arg (ar/R) constriction, has been well characterized by mutational analysis, experimental data concerning the NPA region--in particular, the Asn position--is missing. Here, we report on the cloning and mutational analysis of a novel aquaglyceroporin carrying one SPA motif instead of the NPA motif from Burkholderia cenocepacia, an epidemic pathogen of cystic fibrosis patients. Of 1357 AQP sequences deposited in RefSeq, we identified only 15 with an Asn exchange. Using direct and phenotypic permeability assays, we found that Asn and Ser are freely interchangeable at both NPA sites without affecting protein expression or water, glycerol and methylamine permeability. However, other mutations in the NPA region led to reduced permeability (S186C and S186D), to nonfunctional channels (N64D), or even to lack of protein expression (S186A and S186T). Using electrophysiology, we found that an analogous mammalian AQP1 N76S mutant excluded protons and potassium ions, but leaked sodium ions, providing an argument for the overwhelming prevalence of Asn over other amino acids. We conclude that, at the first position in the NPA motifs, only Asn provides efficient helix cap stabilization and cation exclusion, whereas other small residues compromise structural stability or cation exclusion but not necessarily water and solute permeability.     

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.     

Aquaporin-11 containing a divergent NPA motif has normal water channel activity

Recently, two novel mammalian aquaporins (AQPs), AQPs 11 and 12, have been identified and classified as members of a new AQP subfamily, the "subcellular AQPs". In members of this subfamily one of the two asparagine-proline-alanine (NPA) motifs, which play a crucial role in selective water conduction, are not completely conserved. Mouse AQP11 (mAQP11) was expressed in Sf9 cells and purified using the detergent Fos-choline 10. The protein was reconstituted into liposomes, which were used for water conduction studies with a stopped-flow device. Single water permeability (pf) of AQP11 was measured to be 1.72+/-0.03x10(-13) cm(3)/s, suggesting that other members of the subfamily with incompletely conserved NPA motifs may also function as water channels.     

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.     

Aquaporin-11 containing a divergent NPA motif has normal water channel activity

Recently, two novel mammalian aquaporins (AQPs), AQPs 11 and 12, have been identified and classified as members of a new AQP subfamily, the “subcellular AQPs”. In members of this subfamily one of the two asparagine-proline-alanine (NPA) motifs, which play a crucial role in selective water conduction, are not completely conserved. Mouse AQP11 (mAQP11) was expressed in Sf9 cells and purified using the detergent Fos-choline 10. The protein was reconstituted into liposomes, which were used for water conduction studies with a stopped-flow device. Single water permeability (pf) of AQP11 was measured to be 1.72 ± 0.03 × 10^− 13 cm³/s, suggesting that other members of the subfamily with incompletely conserved NPA motifs may also function as water channels.     

Water transport in human aquaporin-4: molecular dynamics (MD) simulations

Aquaporin-4 (AQP4) is the predominant water channel in the central nervous system, where it has been reported to be involved in many pathophysiological roles including water transport. In this paper, the AQP4 tetramer was modeled from its PDB structure file, embedded in a palmitoyl-oleoyl-phosphatidyl-choline (POPC) lipid bilayer, solvated in water, then minimized and equilibrated by means of molecular dynamics simulations. Analysis of the equilibrated structure showed that the central pore along the fourfold axis of the tetramers is formed with hydrophobic amino acid residues. In particular, Phe-195, Leu-191 and Leu-75, form the narrowest part of the pore. Therefore water molecules are not expected to transport through the central pore, which was confirmed by MD simulations. Each monomer of the AQP4 tetramers forms a channel whose walls consist mostly of hydrophilic residues. There are eight water molecules in single file observed in each of the four channels, transporting through the selectivity filter containing Arg-216, His-201, Phe-77, Ala-210, and the two conserved Asn-Pro-Ala (NPA) motifs containing Asn-213 and Asn-97. By using Brownian dynamics fluctuation–dissipation-theorem (BD-FDT), the overall free-energy profile was obtained for water transporting through AQP4 for the first time, which gives a complete map of the entire channel of water permeation.     

Blockade of sodium and potassium channels in the node of Ranvier by ajmaline and N-propyl ajmaline

The inhibition of sodium and potassium currents in frog myelinated fibres by ajmaline (AM) and its quaternary derivative, N-propyl ajmaline (NPA), depends on voltage-clamp pulses and the state of channel gating mechanisms. The permanently charged NPA and protonated AM interact only (or mainly) with open channels, while unprotonated AM affects preferently inactivated Na channels. Inhibition of Na currents by NPA and AM does not depend on the current direction and Na ion concentration in external or internal media. In contrast only the outward potassium currents can be blocked by NPA and AM; the inward potassium currents in high K+ ions external media are resistant to the blocking action of these drugs. The voltage dependence of ionic current inhibition by charged drugs suggests the location of their binding sites in the inner mouths of Na and K channels. Judging by the kinetics of current restoration after cessation of pulsing, the drug-binding site complex is much more stable in Na than in potassium channels. Batrachotoxin and aconitine, unlike veratridine and sea anemone toxin, decrease greatly the affinity of Na channel binding sites to NPA and AM. The effects of NPA and AM are compared with those of local anesthetics and other amine blocking drugs.     

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.     

Analysis of the pore of the unusual major intrinsic protein channel, yeast Fps1p.

Fps1p is a glycerol efflux channel from Saccharomyces cerevisiae. In this atypical major intrinsic protein neither of the signature NPA motifs of the family, which are part of the pore, is preserved. To understand the functional consequences of this feature, we analyzed the pseudo-NPA motifs of Fps1p by site-directed mutagenesis and assayed the resultant mutant proteins in vivo. In addition, we took advantage of the fact that the closest bacterial homolog of Fps1p, Escherichia coli GlpF, can be functionally expressed in yeast, thus enabling the analysis in yeast cells of mutations that make this typical major intrinsic protein more similar to Fps1p. We observed that mutations made in Fps1p to "restore" the signature NPA motifs did not substantially affect channel function. In contrast, when GlpF was mutated to resemble Fps1p, all mutants had reduced activity compared with wild type. We rationalized these data by constructing models of one GlpF mutant and of the transmembrane core of Fps1p. Our model predicts that the pore of Fps1p is more flexible than that of GlpF. We discuss the fact that this may accommodate the divergent NPA motifs of Fps1p and that the different pore structures of Fps1p and GlpF may reflect the physiological roles of the two glycerol facilitators.     

Presence of aquaporin and V-ATPase on the contractile vacuole of Amoeba proteus.

BACKGROUND INFORMATION: The results of water permeability measurements suggest the presence of an AQP (aquaporin) in the membrane of the CV (contractile vacuole) in Amoeba proteus [Nishihara, Shimmen and Sonobe (2004) Cell Struct. Funct. 29, 85-90]. RESULTS: In the present study, we cloned an AQP gene from A. proteus [ApAQP (A. proteus AQP)] that encodes a 295-amino-acid protein. The protein has six putative TMs (transmembrane domains) and two NPA (Asn-Pro-Ala) motifs, which are conserved among various AQPs and are thought to be involved in the formation of water channels that span the lipid bilayer. Using Xenopus oocytes, we have demonstrated that the ApAQP protein product can function as a water channel. Immunofluorescence microscopy with anti-ApAQP antibody revealed that ApAQP is detected on the CV membrane and on the vesicles around the CV. The presence of V-ATPase (vacuolar H+-ATPase) on the vesicle membrane around the CV was also detected. CONCLUSIONS: Our data on ApAQP allow us to provide the first informed explanation of the high water permeability of the CV membrane in amoeba. Moreover, the results suggest that vesicles possessing V-ATPase are involved in generating an osmotic gradient. Based on our findings, we propose a new hypothesis for the mechanism of CV function.     

Cloning, heterologous expression, and characterization of three aquaglyceroporins from Trypanosoma brucei

Trypanosoma brucei, causative for African sleeping sickness, relies exclusively on glycolysis for ATP production. Under anaerobic conditions, glucose is converted to equimolar amounts of glycerol and pyruvate, which are both secreted from the parasite. As we have shown previously, glycerol transport in T. brucei occurs via specific membrane proteins (Wille, U., Schade, B., and Duszenko, M. (1998) Eur. J. Biochem. 256, 245-250). Here, we describe cloning and biochemical characterization of the three trypanosomal aquaglyceroporins (AQP; TbAQP1-3), which show a 40-45% identity to mammalian AQP3 and -9. AQPs belong to the major intrinsic protein family and represent channels for small non-ionic molecules. Both TbAQP1 and TbAQP3 contain two highly conserved NPA motifs within the pore-forming region, whereas TbAQP2 contains NSA and NPS motifs instead, which are only occasionally found in AQPs. For functional characterization, all three proteins were heterologously expressed in yeast and Xenopus oocytes. In the yeast fps1Delta mutant, TbAQPs suppressed hypoosmosensitivity and rendered cells to a hyper-osmosensitive phenotype, as expected for unregulated glycerol channels. Under iso- and hyperosmotic conditions, these cells constitutively released glycerol, consistent with a glycerol efflux function of TbAQP proteins. TbAQP expression in Xenopus oocytes increased permeability for water, glycerol and, interestingly, dihydroxyacetone. Except for urea, TbAQPs were virtually impermeable for other polyols; only TbAQP3 transported erythritol and ribitol. Thus, TbAQPs represent mainly water/glycerol/dihydroxyacetone channels involved in osmoregulation and glycerol metabolism in T. brucei. This function and especially the so far not investigated transport of dihydroxyacetone may be pivotal for the survival of the parasite survival under non-aerobic or osmotic stress conditions.     

AQP1 expression analysis in human diseases: implications for proteomic characterization

Aquaporin (AQP)1 belongs to a ubiquitous family of water channel proteins characterized by sequence similarity and the presence of two NPA (Asp-Pro-Ala) motifs existing in almost all organs and tissues. Currently, 13 human AQPs are known and they are divided into two subgroups according to their ability to transport only water molecules, such as AQP1, or also glycerol and other small solutes. The genomic, structural and functional aspects of AQP1 are briefly described. An in-depth discussion is devoted to proteomic approaches that are useful for identifying and characterizing AQP1, mainly through electrophoretic techniques combined with different extraction procedures followed by mass spectrometry analysis. Moreover, the relevance of AQP1 in human diseases is also explained. Its role in human tumors and, in particular, those of the kidney (e.g., clear cell renal carcinoma) is discussed.     

Functional characterization of a putative aquaporin from Encephalitozoon cuniculi, a microsporidia pathogenic to humans

The microsporidia are a group of obligate intracellular parasitic protists that have been implicated as both human and veterinary pathogens. The infectious process of these organisms is believed to be dependent upon the rapid influx of water into spores, presumably via aquaporins (AQPs), transmembrane channels that facilitate osmosis. An AQP-like sequence of the microsporidium Encephalitozoon cuniculi (EcAQP), when cloned and expressed in oocytes of Xenopus laevis, rendered these oocytes highly permeable to water. No permeability to the solutes glycerol or urea was observed. Pre-treatment of EcAQP-expressing oocytes with HgCl(2) failed to inhibit their osmotic permeability, as predicted from EcAQP's lack of mercury-sensitive cysteine residues near the NPA motifs which line the AQP aqueous pore. EcAQP exhibits sequence identity to AQP A of Dictyostelium discoideum (26%) and human AQP 2 (24%). Further study of AQPs in microsporidia and their potential inhibitors may yield novel therapeutic agents for microsporidian infections.     

The mechanism of proton exclusion in the aquaporin-1 water channel

Aquaporins are efficient, yet strictly selective water channels. Remarkably, proton permeation is fully blocked, in contrast to most other water-filled pores which are known to conduct protons well. Blocking of protons by aquaporins is essential to maintain the electrochemical gradient across cellular and subcellular membranes. We studied the mechanism of proton exclusion in aquaporin-1 by multiple non-equilibrium molecular dynamics simulations that also allow proton transfer reactions. From the simulations, an effective free energy profile for the proton motion along the channel was determined with a maximum-likelihood approach. The results indicate that the main barrier is not, as had previously been speculated, caused by the interruption of the hydrogen-bonded water chain, but rather by an electrostatic field centered around the fingerprint Asn-Pro-Ala (NPA) motif. Hydrogen bond interruption only forms a secondary barrier located at the ar/R constriction region. The calculated main barrier height of 25-30 kJ mol(-1) matches the barrier height for the passage of protons across pure lipid bilayers and, therefore, suffices to prevent major leakage of protons through aquaporins. Conventional molecular dynamics simulations additionally showed that negatively charged hydroxide ions are prevented from being trapped within the NPA region by two adjacent electrostatic barriers of opposite polarity.     

Phosphorylation-dependent regulation of ryanodine receptors: a novel role for leucine/isoleucine zippers

Ryanodine receptors (RyRs), intracellular calcium release channels required for cardiac and skeletal muscle contraction, are macromolecular complexes that include kinases and phosphatases. Phosphorylation/dephosphorylation plays a key role in regulating the function of many ion channels, including RyRs. However, the mechanism by which kinases and phosphatases are targeted to ion channels is not well understood. We have identified a novel mechanism involved in the formation of ion channel macromolecular complexes: kinase and phosphatase targeting proteins binding to ion channels via leucine/isoleucine zipper (LZ) motifs. Activation of kinases and phosphatases bound to RyR2 via LZs regulates phosphorylation of the channel, and disruption of kinase binding via LZ motifs prevents phosphorylation of RyR2. Elucidation of this new role for LZs in ion channel macromolecular complexes now permits: (a) rapid mapping of kinase and phosphatase targeting protein binding sites on ion channels; (b) predicting which kinases and phosphatases are likely to regulate a given ion channel; (c) rapid identification of novel kinase and phosphatase targeting proteins; and (d) tools for dissecting the role of kinases and phosphatases as modulators of ion channel function.     

Intra-helical salt-bridge and helix destabilizing residues within the same helical turn: Role of functionally important loop E half-helix in channel regulation of major intrinsic proteins

The superfamily of major intrinsic proteins (MIPs) includes aquaporin (AQP) and aquaglyceroporin (AQGP) and it is involved in the transport of water and neutral solutes across the membrane. Diverse MIP sequences adopt a unique hour-glass fold with six transmembrane helices (TM1 to TM6) and two half-helices (LB and LE). Loop E contains one of the two conserved NPA motifs and contributes two residues to the aromatic/arginine selectivity filter. Function and regulation of majority of MIP channels are not yet characterized. We have analyzed the loop E region of 1468 MIP sequences and their structural models from six different organism groups. They can be phylogenetically clustered into AQGPs, AQPs, plant MIPs and other MIPs. The LE half-helix in all AQGPs contains an intra-helical salt-bridge and helix-breaking residues Gly/Pro within the same helical turn. All non-AQGPs lack this salt-bridge but have the helix destabilizing Gly and/or Pro in the same positions. However, the segment connecting LE half-helix and TM6 is longer by 10–15 residues in AQGPs compared to all non-AQGPs. We speculate that this longer loop in AQGPs and the LE half-helix of non-AQGPs will be relatively more flexible and this could be functionally important. Molecular dynamics simulations on glycerol-specific GlpF, water-transporting AQP1, its mutant and a fungal AQP channel confirm these predictions. Thus two distinct regions of loop E, one in AQGPs and the other in non-AQGPs, seem to be capable of modulating the transport. These regions can also act in conjunction with other extracellular residues/segments to regulate MIP channel transport.