Aquaglyceroporins,one channel for two molecules

In the light of the recently published structure of GlpF and AQP1, we have analysed the nature of the residues which could be involved in the formation of the selectivity filter of aquaporins, glycerol facilitators and aquaglyceroporins. We demonstrate that the functional specificity for major intrinsic protein (MIP) channels can be explained on one side by analysing the polar environment of the residues that form the selective filter. On the other side, we show that the channel selectivity could be associated with the oligomeric state of the membrane protein. We conclude that a non-polar environment in the vicinity of the top of helix 5 could allow aquaglyceroporins and GlpF to exist as monomers within the hydrophobic environment of the membrane.     

The 3.7 A projection map of the glycerol facilitator GlpF: a variant of the aquaporin tetramer

GlpF, the glycerol facilitator protein of Escherichia coli, is an archetypal member of the aquaporin superfamily. To assess its structure, recombinant histidine-tagged protein was overexpressed, solubilized in octylglucoside and purified to homogeneity. Negative stain electron microscopy of solubilized GlpF protein revealed a tetrameric structure of ~80 Å side length. Scanning transmission electron microscopy yielded a mass of 170 kDa, corroborating the tetrameric nature of GlpF. Reconstitution of GlpF in the presence of lipids produced highly ordered two-dimensional crystals, which diffracted electrons to 3.6 Å resolution. Cryoelectron microscopy provided a 3.7 Å projection map exhibiting a unit cell comprised of two tetramers. In projection, GlpF is similar to AQP1, the erythrocyte water channel. However, the major density minimum within each monomer is distinctly larger in GlpF than in AQP1.     

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.     

Selectivity and conductance among the glycerol and water conducting aquaporin family of channels

The atomic structures of a transmembrane water plus glycerol conducting channel (GlpF), and now of aquaporin Z (AqpZ) from the same species, Escherichia coli, bring the total to three atomic resolution structures in the aquaporin (AQP) family. Members of the AQP family each assemble as tetramers of four channels. Common helical axes support a wider channel in the glycerol plus water channel paradigm, GlpF. Water molecules form a single hydrogen bonded file throughout the 28 A long channel in AqpZ. The basis for absolute exclusion of proton or hydronium ion conductance through the line of water is explored using simulations.     

Arsenic trioxide uptake by human and rat aquaglyceroporins

Aquaglyceroporins are channels that allow downhill movement of uncharged solutes such as glycerol and urea. Arsenic trioxide has recently been shown to be translocated by mouse mAQP7 and rat rAQP9. In this study we examined the ability of the four known human members of the aquaglyceroporin family, hAQP3, hAQP7, hAQP9, and hAQP10, to facilitate As(OH)(3) movement in Xenopus oocytes. The order of effectiveness as an As(III) transporter was found to be hAQP9 > hAQP7, with little or no transport by hAQP3 or hAQP10. From comparison with the crystal structure of the bacterial homologue GlpF and the bovine erythrocyte water channel bAQP1, AQP9 residues Phe-64 and Arg-219 are predicted to serve as part of the selectivity filter. The requirement for Phe-64 and Arg-219 in arsenic trioxide translocation was examined by site-directed mutagenesis of rAQP9, taking advantage of the fact that rat AQP9 catalyzes (73)As(OH)(3) uptake in Saccharomyces cerevisiae and in oocytes. R219A, R219K, F64A, F64T, and F64W were expressed in both yeast and oocytes, and permeability of arsenic trioxide and glycerol was measured. A lysine but not an alanine residue could substitute for the highly conserved Arg-219, indicating that a positive charge is required at the entry to the channel. In contrast, the phenylalanine residue, which is believed to position substrates near the conserved arginine, was not required for either arsenic trioxide or glycerol uptake. The results support the hypothesis that arsenic trioxide and glycerol use the same translocation pathway in AQP9.     

A water channel of the nematode C.elegans and its implications for channel selectivity of MIP proteins

A genome project focusingon the nematode Caenorhabditis elegans has demonstrated thepresence of eight cDNAs belonging to the major intrinsic proteinsuperfamily. We functionally characterized one of these cDNAs namedC01G6.1. Injection of C01G6.1 cRNA increased the osmotic waterpermeability (P_f) of Xenopusoocytes 11-fold and the urea permeability 4.5-fold but failed toincrease the glycerol permeability. It has been speculated that the MIPfamily may be separated into two large subfamilies based on thepresence or absence of two segments of extra amino acid residues (~15amino acids) at the second and third extracellular loops. BecauseC01G6.1 (designated AQP-CE1), AQP3, and glycerol facilitator (GlpF) all have these two segments, we replaced the segments of AQP-CE1 with thoseof AQP3 and GlpF to identify their roles. The functional characteristics of these mutants were principally similar to that ofwild-type AQP-CE1, although the values of P_f andurea permeability were decreased by 39-74% and 28-65%,respectively. These results suggest that the two segments of extraamino acid residues may not contribute to channel selectivity orformation of the route for small solutes.

     

Glycerol conductance and physical asymmetry of the Escherichia coli glycerol facilitator GlpF

The aquaglyceroporin GlpF is a transmembrane channel of Escherichia coli that facilitates the uptake of glycerol by the cell. Its high glycerol uptake rate is crucial for the cell to survive in very low glycerol concentrations. Although GlpF allows both influx and outflux of glycerol, its structure, similar to the structure of maltoporin, exhibits a significant degree of asymmetry. The potential of mean force characterizing glycerol in the channel shows a corresponding asymmetry with an attractive vestibule only at the periplasmic side. In this study, we analyze the potential of mean force, showing that a simplified six-step model captures the kinetics and yields a glycerol conduction rate that agrees well with observation. The vestibule improves the conduction rate by 40% and 75% at 10- micro M and 10-mM periplasmic glycerol concentrations, respectively. In addition, neither the conduction rate nor the conduction probability for a single glycerol (efficiency) depends on the orientation of GlpF. GlpF appears to conduct equally well in both directions under physiological conditions.     

Comparative simulations of aquaporin family: AQP1, AQPZ, AQP0 and GlpF

Molecular dynamics simulations were performed for four members of the aquaporin family (AQP1, AQPZ, AQP0, and GlpF) in the explicit membrane environment. The single-channel water permeability, pf, was evaluated to be GlpF approximately AQPZ > AQP1 > AQP0, while their relative pore sizes were GlpF > AQP1 > AQPZ > AQP0. This relation between pf and pore size indicates that water permeability was determined not only by the channel radius, but also another competing factor. Analysis of water dynamics revealed that this factor was the single-file nature of water transport.     

Projection map of aquaporin-9 at 7 A resolution

Aquaporin-9, an aquaglyceroporin present in diverse tissues, is unique among aquaporins because it is not only permeable to water, urea and glycerol, but also allows passage of larger uncharged solutes. Single particle analysis of negatively stained recombinant rat aquaporin-9 revealed a particle size characteristic of the tetrameric organization of all members of the aquaporin family. Reconstitution of aquaporin-9 into two-dimensional crystals enabled us to calculate a projection map at 7 A resolution. The projection structure indicates a tetrameric structure, similar to GlpF, with each square-like monomer forming a pore. A comparison of the pore-lining residues between the crystal structure of GlpF and a homology model of aquaporin-9 locates substitutions in these residues predominantly to the hydrophobic edge of the tripathic pore of GlpF, providing first insights into the structural basis for the broader substrate specificity of aquaporin-9.     

The structure of the aquaporin-1 water channel: a comparison between cryo-electron microscopy and X-ray crystallography

Three different medium-resolution structures of the human water channel aquaporin-1 (AQP1) have been solved by cryo-electron microscopy (cryo-EM) during the last two years. Recently, the structure of the strongly related bovine AQP1 was solved by X-ray crystallography at higher resolution, allowing a validation of the original medium-resolution structures, and providing a good indication for the strengths and limitations of state of the art cryo-EM methods. We present a detailed comparison between the different models, which shows that overall, the structures are highly similar, deviating less than 2.5 A from each other in the helical backbone regions. The two original cryo-EM structures, however, also show a number of significant deviations from the X-ray structure, both in the backbone positions of the transmembrane helices and in the location of the amino acid side-chains facing the pore. In contrast, the third cryo-EM structure that included information from the X-ray structure of the homologous bacterial glycerol facilitator GlpF and that was subsequently refined against cryo-EM AQP1 data, shows a root mean square deviation of 0.9A from the X-ray structure in the helical backbone regions.     

What makes an aquaporin a glycerol channel? A comparative study of AqpZ and GlpF

The recent availability of high-resolution structures of two structurally highly homologous, but functionally distinct aquaporins from the same species, namely Escherichia coli AqpZ, a pure water channel, and GlpF, a glycerol channel, presents a unique opportunity to understand the mechanism of substrate selectivity in these channels. Comparison of the free energy profile of glycerol conduction through AqpZ and GlpF reveals a much larger barrier in AqpZ (22.8 kcal/mol) than in GlpF (7.3 kcal/mol). In either channel, the highest barrier is located at the selectivity filter. Analysis of substrate-protein interactions suggests that steric restriction of AqpZ is the main contribution to this large barrier. Another important difference is the presence of a deep energy well at the periplasmic vestibule of GlpF, which was not found in AqpZ. The latter difference can be attributed to the more pronounced structural asymmetry of GlpF, which may play a role in attracting glycerol.     

The GlpF residue Trp219 is part of an amino-acid cluster crucial for aquaglyceroporin oligomerization and function

The vestibule loop regions of aquaglyceroporins are involved in accumulation of glycerol inside the channel pore. Even though most loop regions do not show high sequence similarity among aquaglyceroporins, loop E is highly conserved in aquaglyceroporins, but not in members of the homologous aquaporins. Specifically, a tryptophan residue is extremely conserved within this loop. We have investigated the role of this residue (Trp219) that deeply protrudes into the protein and potentially interacts with adjacent loops, using the E. coli aqualgyeroporin GlpF as a model. Replacement of Trp219 affects the activity of GlpF and impairs the stability of the tetrameric protein. Furthermore, we have identified an amino acid cluster involving Trp219 that stabilizes the GlpF tetramer. Based on our results we propose that Trp219 is key for formation of a defined vestibule structure, which is crucial for glycerol accumulation as well as for the stability of the active GlpF tetramer.     

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.     

Transport of water and glycerol in aquaporin 3 is gated by H(+).

Aquaporins (AQPs) were expressed in Xenopus laevis oocytes in order to study the effects of external pH and solute structure on permeabilities. For AQP3 the osmotic water permeability, L(p), was abolished at acid pH values with a pK of 6.4 and a Hill coefficient of 3. The L(p) values of AQP0, AQP1, AQP2, AQP4, and AQP5 were independent of pH. For AQP3 the glycerol permeability P(Gl), obtained from [(14)C]glycerol uptake, was abolished at acid pH values with a pK of 6.1 and a Hill coefficient of 6. Consequently, AQP3 acts as a glycerol and water channel at physiological pH, but predominantly as a glycerol channel at pH values around 6.1. The pH effects were reversible. The interactions between fluxes of water and straight chain polyols were inferred from reflection coefficients (sigma). For AQP3, water and glycerol interacted by competing for titratable site(s): sigma(Gl) was 0.15 at neutral pH but doubled at pH 6.4. The sigma values were smaller for polyols in which the -OH groups were free to form hydrogen bonds. The activation energy for the transport processes was around 5 kcal mol(-1). We suggest that water and polyols permeate AQP3 by forming successive hydrogen bonds with titratable sites.     

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.     

A Phylogenetic Framework for the Aquaporin Family in Eukaryotes

A comprehensive evolutionary analysis of aquaporins, a family of intrinsic membrane proteins that function as water channels, was conducted to establish groups of homology (i.e., to identify orthologues and paralogues) within the family and to gain insights into the functional constraints acting on the structure of the aquaporin molecule structure. Aquaporins are present in all living organisms, and therefore, they provide an excellent opportunity to further our understanding of the broader biological significance of molecular evolution by gene duplication followed by functional and structural specialization. Based on the resulting phylogeny, the 153 channel proteins analyzed were classified into six major paralogous groups: (1) GLPs, or glycerol-transporting channel proteins, which include mammalian AQP3, AQP7, and AQP9, several nematode paralogues, a yeast paralogue, and Escherichia coli GLP; (2) AQPs, or aquaporins, which include metazoan AQP0, AQP1, AQP2, AQP4, AQP5, and AQP6; (3) PIPs, or plasma membrane intrinsic proteins of plants, which include PIP1 and PIP2; (4) TIPs, or tonoplast intrinsic proteins of plants, which include alphaTIP, gammaTIP, and deltaTIP; (5) NODs, or nodulins of plants; and (6) AQP8s, or metazoan aquaporin 8 proteins. Of these groups, AQPs, PIPs, and TIPs cluster together. According to the results, the capacity to transport glycerol shown by several members of the family was acquired only early in the history of the family. The new phylogeny reveals that several water channel proteins are misclassified and require reassignment, whereas several previously undetermined ones can now be classified with confidence. The deduced phylogenetic framework was used to characterize the molecular features of water channel proteins. Three motifs are common to all family members: AEF (Ala-Glu-Phe), which is located in the N-terminal domain; and two NPA (Asp-Pro-Ala) boxes, which are located in the center and C-terminal domains, respectively. Other residues are found to be conserved within the major groups but not among them. Overall, the PIP subfamily showed the least variation. In general, no radical amino acid replacements affecting tertiary structure were identified, with the exception of Ala-->Ser in the TIP subfamily. Constancy of rates of evolution was demonstrated within the different paralogues but rejected among several of them (GLP and NOD).     

A refined structure of human aquaporin-1

A refined structure of the human water channel aquaporin-1 is presented. The model rests on the high resolution X-ray structure of the homologous bacterial glycerol transporter GlpF, electron crystallographic data at 3.8 A resolution and a multiple sequence alignment of the aquaporin superfamily. The crystallographic R and free R values (36.7% and 37.8%) for the refined structure are significantly lower than for previous models. Improved geometry and enhanced stability in molecular dynamics simulations demonstrate a significant improvement of the aquaporin-1 structure. Comparison with previous aquaporin-1 models shows significant differences, not only in the loop regions, but also in the core of the water channel.     

Biophysical Assessment of Human Aquaporin-7 as a Water and Glycerol Channel in 3T3-L1 Adipocytes

The plasma membrane aquaporin-7 (AQP7) has been shown to be expressed in adipose tissue and its role in glycerol release/uptake in adipocytes has been postulated and correlated with obesity onset. However, some studies have contradicted this view. Based on this situation, we have re-assessed the precise localization of AQP7 in adipose tissue and analyzed its function as a water and/or glycerol channel in adipose cells. Fractionation of mice adipose tissue revealed that AQP7 is located in both adipose and stromal vascular fractions. Moreover, AQP7 was the only aquaglyceroporin expressed in adipose tissue and in 3T3-L1 adipocytes. By overexpressing the human AQP7 in 3T3-L1 adipocytes it was possible to ascertain its role as a water and glycerol channel in a gain-of-function scenario. AQP7 expression had no effect in equilibrium cell volume but AQP7 loss of function correlated with higher triglyceride content. Furthermore it is also reported for the first time a negative correlation between water permeability and the cell non-osmotic volume supporting the observation that AQP7 depleted cells are more prone to lipid accumulation. Additionally, the strong positive correlation between the rates of water and glycerol transport highlights the role of AQP7 as both a water and a glycerol channel and reflects its expression levels in cells. In all, our results clearly document a direct involvement of AQP7 in water and glycerol transport, as well as in triglyceride content in adipocytes.     

Water transport in aquaporins: osmotic permeability matrix analysis of molecular dynamics simulations

Single-channel osmotic water permeability (p(f)) is a key quantity for investigating the transport capability of the water channel protein, aquaporin. However, the direct connection between the single scalar quantity p(f) and the channel structure remains unclear. In this study, based on molecular dynamics simulations, we propose a p(f)-matrix method, in which p(f) is decomposed into contributions from each local region of the channel. Diagonal elements of the p(f) matrix are equivalent to the local permeability at each region of the channel, and off-diagonal elements represent correlated motions of water molecules in different regions. Averaging both diagonal and off-diagonal elements of the p(f) matrix recovers p(f) for the entire channel; this implies that correlated motions between distantly-separated water molecules, as well as adjacent water molecules, influence the osmotic permeability. The p(f) matrices from molecular dynamics simulations of five aquaporins (AQP0, AQP1, AQP4, AqpZ, and GlpF) indicated that the reduction in the water correlation across the Asn-Pro-Ala region, and the small local permeability around the ar/R region, characterize the transport efficiency of water. These structural determinants in water permeation were confirmed in molecular dynamics simulations of three mutants of AqpZ, which mimic AQP1.