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.     

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.     

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.     

Role of C-terminal domain and transmembrane helices 5 and 6 in function and quaternary structure of major intrinsic proteins: analysis of aquaporin/glycerol facilitator chimeric proteins

We previously observed that aquaporins and glycerol facilitators exhibit different oligomeric states when studied by sedimentation on density gradients following nondenaturing detergent solubilization. To determine the domains of major intrinsic protein (MIP) family proteins involved in oligomerization, we constructed protein chimeras corresponding to the aquaporin AQPcic substituted in the loop E (including the proximal part of transmembrane domain (TM) 5) and/or the C-terminal part (including the distal part of TM 6) by the equivalent domain of the glycerol channel aquaglyceroporin (GlpF) (chimeras called AGA, AAG, and AGG). The analogous chimeras of GlpF were also constructed (chimeras GAG, GGA, and GAA). cRNA corresponding to all constructs were injected into Xenopus oocytes. AQPcic, GlpF, AAG, AGG, and GAG were targeted to plasma membranes. Water or glycerol membrane permeability measurements demonstrated that only the AAG chimera exhibited a channel function corresponding to water transport. Analysis of all proteins expressed either in oocytes or in yeast by velocity sedimentation on sucrose gradients following solubilization by 2% n-octyl glucoside indicated that only AQPcic and AAG exist in tetrameric forms. GlpF, GAG, and GAA sediment in a monomeric form, whereas GGA and AGG were found mono/dimeric. These data bring new evidence that, within the MIP family, aquaporins and GlpFs behave differently toward nondenaturing detergents. We demonstrate that the C-terminal part of AQPcic, including the distal half of TM 6, can be substituted by the equivalent domain of GlpF (AAG chimera) without modifying the transport specificity. Our results also suggest that interactions of TM 5 of one monomer with TM 1 of the adjacent monomer are crucial for aquaporin tetramer stability.     

Multiple pore lining residues modulate water permeability of GlpF

The water permeability of aquaporins (AQPs) varies by more than an order of magnitude even though the pore structure, geometry, as well as the channel lining residues are highly conserved. However, channel gating by pH, divalent ions or phosphorylation was only shown for a minority of AQPs. Structural and in silico indications of water flux modulation by flexible side chains of channel lining residues have not been experimentally confirmed yet. Hence, the aquaporin “open state” is still considered to be a continuously open pore with water molecules permeating in a single‐file fashion. Using protein mutations outside the selectivity filter in the aqua(glycerol)facilitator GlpF of Escherichia coli we, to the best of our knowledge, for the first time, modulate the position of the highly conserved Arg in the selectivity filter. This in turn enhances or reduces the unitary water permeability of GlpF as shown in silico by molecular dynamics (MD) simulations and in vitro with purified and reconstituted GlpF. This finding suggests that AQP water permeability can indeed be regulated by lipid bilayer asymmetry and the transmembrane potential. Strikingly, our long‐term MD simulations reveal that not only the conserved Arg in the selectivity filter, but the position and dynamics of multiple other pore lining residues modulate water passage through GlpF. This finding is expected to trigger a wealth of future investigations on permeability and regulation of AQPs among others with the aim to tune water permeability for biotechnological applications.     

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.     

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.     

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.     

Diffusion of glycerol through Escherichia coli aquaglyceroporin GlpF

The glycerol uptake facilitator, GlpF, a major intrinsic protein found in Escherichia coli, selectively conducts water and glycerol across the inner membrane. The free energy landscape characterizing the assisted transport of glycerol by this homotetrameric aquaglyceroporin has been explored by means of equilibrium molecular dynamics over a timescale spanning 0.12 μs. To overcome the free energy barriers of the conduction pathway, an adaptive biasing force is applied to the glycerol molecule confined in each of the four channels. The results illuminate the critical role played by intramolecular relaxation on the diffusion properties of the permeant. These free energy calculations reveal that glycerol tumbles and isomerizes on a timescale comparable to that spanned by its adaptive-biasing-force-assisted conduction in GlpF. As a result, reorientation and conformational equilibrium of glycerol in GlpF constitute a bottleneck in the molecular simulations of the permeation event. A profile characterizing the position-dependent diffusion of the permeant has been determined, allowing reaction rate theory to be applied for investigating conduction kinetics based on the measured free energy landscape.     

Crystal structure of SecB from Escherichia coli

The chaperone SecB from Escherichia coli is primarily involved in passing precursor proteins into the Sec system via specific interactions with SecA. The crystal structure of SecB from E. coli has been solved to 2.35 A resolution. The structure shows flexibility in the crossover loop and the helix-connecting loop, regions that have been implicated to be part of the SecB substrate-binding site. Moreover conformational variability of Trp36 is observed as well as different loop conformations for the different monomers. Based on this, we speculate that SecB can regulate the access or extent of its hydrophobic substrate-binding site, by modulating the conformation of the crossover loop and the helix-connecting loop. The structure also clearly explains why the tetrameric equilibrium is shifted towards the dimeric state in the mutant SecBCys76Tyr. The buried cysteine residue is crucial for tight packing, and mutations are likely to disrupt the tetramer formation but not the dimer formation.     

A conserved Trp residue in HwBR contributes to its unique tolerance toward acidic environments

Bacteriorhodopsin (BR) is a light-driven outward proton pump found mainly in halophilic archaea. A BR from an archaeon Haloquadratum walsbyi (HwBR) was found to pump protons under more acidic conditions compared with most known BR proteins. The atomic structural study on HwBR unveiled that a pair of hydrogen bonds between the BC and FG loop in its periplasmic region may be a factor in such improved pumping capability. Here, we further investigated the retinal-binding pocket of HwBR and found that Trp94 contributes to the higher acid tolerance. Through single mutations in a BR from Halobacterium salinarum and HwBR, we examined the conserved tryptophan residues in the retinal-binding pocket. Among these residues of HwBR, mutagenesis at Trp94 facing the periplasmic region caused the most significant disruption to optical stability and proton-pumping capability under acidic conditions. The other tryptophan residues of HwBR exerted little impact on both maximum absorption wavelength and pH-dependent proton pumping. Our findings suggest that the residues from Trp94 to the hydrogen bonds at the BC loop confer both optical stability and functionality on the overall protein in low-pH environments.     

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.     

Homology modelling and molecular dynamics simulations: comparative studies of human aquaporin-1

The structures of the mammalian water transport protein Aqp1 and of its bacterial homologue GlpF enables us to test whether homology models can be used to explore relationships between structure, dynamics and function in mammalian transport proteins. Molecular dynamics simulations (totalling almost 40 ns) were performed starting from: the X-ray structure of Aqp1; a homology model of Aqp1 based on the GlpF structure; and intermediate resolution structures of Aqp1 derived from electron microscopy. Comparisons of protein RMSDs vs. time suggest that the homology models are of comparable conformational stability to the X-ray structure, whereas the intermediate resolution structures exhibit significant conformation drift. For simulations based on the X-ray structure and on homology models, the flexibility profile vs. residue number correlates well with the crystallographic B-values for each residue. In the simulations based on intermediate resolution structures, mobility of the highly conserved NPA loops is substantially higher than in the simulations based on the X-ray structure or the homology models. Pore radius profiles remained relatively constant in the X-ray and homology model simulations but showed substantial fluctuations (reflecting the higher NPA loop mobility) in the intermediate resolution simulations. The orientation of the dipoles of water molecules within the pore is of key importance in maintaining low proton permeability through Aqp1. This property seems to be quite robust to the starting model used in the simulation. These simulations suggest that homology models based on bacterial homologues may be used to derive functionally relevant information on the structural dynamics of mammalian transport proteins.     

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.     

Role of Tryptophan Residues in Toxicity of Cry1Ab Toxin from Bacillus thuringiensis

Bacillus thuringiensis produces insecticidal proteins (Cry protoxins) during the sporulation phase as parasporal crystals. During intoxication, the Cry protoxins must change from insoluble crystals into membrane-inserted toxins which form ionic pores. The structural changes of Cry toxins during oligomerization and insertion into the membrane are still unknown. The Cry1Ab toxin has nine tryptophan residues; seven are located in domain I, the pore-forming domain, and two are located in domain II, which is involved in receptor recognition. Eight Trp residues are highly conserved within the whole family of three-domain Cry proteins, suggesting an essential role for these residues in the structural folding and function of the toxin. In this work, we analyzed the role of Trp residues in the structure and function of Cry1Ab toxin. We replaced the Trp residues with phenylalanine or cysteine using site-directed mutagenesis. Our results show that W65 and W316 are important for insecticidal activity of the toxin since their replacement by Phe reduced the toxicity against Manduca sexta. The presence of hydrophobic residue is important at positions 117, 219, 226, and 455 since replacement by Cys affected either the crystal formation or the insecticidal activity of the toxin in contrast to replacement by Phe in these positions. Additionally, some mutants in positions 219, 316, and 455 were also affected in binding to brush border membrane vesicles (BBMV). This is the first report that studies the role of Trp residues in the activity of Cry toxins.     

Antimonite is accumulated by the glycerol facilitator GlpF in Escherichia coli.

In a search for genes responsible for the accumulation of antimonite in Escherichia coli, TnphoA was used to create a pool of random insertional mutants, from which one antimonite-resistant mutant was isolated. Sequence analysis showed that the TnphoA insertion was located in the glpF gene, coding for the glycerol facilitator GlpF. The mutant was shown to be defective in polyol transport by GlpF. These results suggest that in solution Sb(III) is recognized as a polyol by the glycerol facilitator.     

A rule governing the FtsH-mediated proteolysis of the MgtC virulence protein from <Emphasis Type="Italic">Salmonella enterica</Emphasis> serovar Typhimurium

A tightly controlled turnover of membrane proteins is required for lipid bilayer stability, cell metabolism, and cell viability. Among the energy-dependent AAA⁺ proteases in Salmonella, FtsH is the only membrane-bound protease that contributes to the quality control of membrane proteins. FtsH preferentially degrades the C-terminus or N-terminus of misfolded, misassembled, or damaged proteins to maintain physiological functions. We found that FtsH hydrolyzes the Salmonella MgtC virulence protein when we substitute the MgtC 226^th Trp, which is well conserved in other intracellular pathogens and normally protects MgtC from the FtsH-mediated proteolysis. Here we investigate a rule determining the FtsH-mediated proteolysis of the MgtC protein at Trp226 residue. Substitution of MgtC tryptophan 226^th residue to alanine, glycine, or tyrosine leads to MgtC proteolysis in a manner dependent on the FtsH protease whereas substitution to phenylalanine, methionine, isoleucine, leucine, or valine resists MgtC degradation by FtsH. These data indicate that a large and hydrophobic side chain at 226^th residue is required for protection from the FtsH-mediated MgtC proteolysis.     

Role of conserved residues in structure and stability: tryptophans of human serum retinol-binding protein, a model for the lipocalin superfamily

Serum retinol binding protein (RBP) is a member of the lipocalin family, proteins with up-and-down beta-barrel folds, low levels of sequence identity, and diverse functions. Although tryptophan 24 of RBP is highly conserved among lipocalins, it does not play a direct role in activity. To determine if Trp24 and other conserved residues have roles in stability and/or folding, we investigated the effects of conservative substitutions for the four tryptophans and some adjacent residues on the structure, stability, and spectroscopic properties of apo-RBP. Crystal structures of recombinant human apo-RBP and of a mutant with substitutions for tryptophans 67 and 91 at 1.7 A and 2.0 A resolution, respectively, as well as stability measurements, indicate that these relatively exposed tryptophans have little influence on structure or stability. Although Trp105 is largely buried in the wall of the beta-barrel, it can be replaced with minor effects on stability to thermal and chemical unfolding. In contrast, substitutions of three different amino acids for Trp24 or replacement of Arg139, a conserved residue that interacts with Trp24, lead to similar large losses in stability and lower yields of native protein generated by in vitro folding. The results and the coordinated nature of natural substitutions at these sites support the idea that conserved residues in functionally divergent homologs have roles in stabilizing the native relative to misfolded structures. They also establish conditions for studies of the kinetics of folding and unfolding by identifying spectroscopic signals for monitoring the formation of different substructures.