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.     

Anionic Lipids Modulate the Activity of the Aquaglyceroporin GlpF

The structure and composition of a biological membrane can severely influence the activity of membrane-embedded proteins. Here, we show that the E. coli aquaglyceroporin GlpF has only little activity in lipid bilayers formed from native E. coli lipids. Thus, at first glance, GlpF appears to not be optimized for its natural membrane environment. In fact, we found that GlpF activity was severely affected by negatively charged lipids regardless of the exact chemical nature of the lipid headgroup, whereas GlpF was not sensitive to changes in the lateral membrane pressure. These observations illustrate a potential mechanism by which the activity of an α-helical membrane protein is modulated by the negative charge density around the protein.     

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.     

Molecular determinants of ammonia and urea conductance in plant aquaporin homologs

Aquaporins and/or aquaglyceroporins regulate the permeability of plant membranes to water and small, uncharged molecules. Using molecular simulations with a plant plasma membrane aquaporin tetramer, the residues in the channel constriction region were identified as the crucial determinants of ammonia and urea conductance. The impact of these residues was experimentally verified using AtPIP2;1 pore mutants. Several, but not all, mutants with a NIP-like selectivity filter promoted yeast growth on urea or ammonia as sole sources of nitrogen. TIP-like mutants conducted urea but not NH(3), and a residue without direct contact to the pore lumen was critical for conduction in the mutants.     

The 6.9-A structure of GlpF: a basis for homology modeling of the glycerol channel from Escherichia coli

The three-dimensional structure of GlpF, the glycerol facilitator of Escherichia coli, was determined by cryo-electron microscopy. The 6.9-A density map calculated from images of two-dimensional crystals shows the GlpF helices to be similar to those of AQP1, the erythrocyte water channel. While the helix arrangement of GlpF does not reflect the larger pore diameter as seen in the projection map, additional peripheral densities observed in GlpF are compatible with the 31 additional residues in loops C and E, which accordingly do not interfere with the inner channel construction. Therefore, the atomic structure of AQP1 was used as a basis for homology modeling of the GlpF channel, which is predicted to be free of bends, wider, and more vertically oriented than the AQP1 channel. Furthermore, the residues facing the GlpF channel exhibit an amphiphilic nature, being hydrophobic on one side and hydrophilic on the other side. This property may partially explain the contradiction of glycerol diffusion but limited water permeation capacity.     

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.     

Single-channel water permeabilities of Escherichia coli aquaporins AqpZ and GlpF

From equilibrium molecular dynamics simulations we have determined single-channel water permeabilities for Escherichia coli aquaporin Z (AqpZ) and aquaglyceroporin GlpF with the channels embedded in lipid bilayers. GlpF's osmotic water permeability constant pf exceeds by 2-3 times that of AqpZ and the diffusive permeability constant (pd) of GlpF is found to exceed that of AqpZ 2-9-fold. Achieving complete water selectivity in AqpZ consequently implies lower transport rates overall relative to the less selective, wider channel of GlpF. For AqpZ, the ratio pf/pd congruent with 12 is close to the average number of water molecules in the channel lumen, whereas for GlpF, pf/pd congruent with 4. This implies that single-file structure of the luminal water is more pronounced for AqpZ, the narrower channel of the two. Electrostatics profiles across the pore lumens reveal that AqpZ significantly reinforces water-channel interactions, and weaker water-water interactions in turn suppress water-water correlations relative to GlpF. Consequently, suppressed water-water correlations across the narrow selectivity filter become a key structural determinant for water permeation causing luminal water to permeate slower across AqpZ.     

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.     

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.     

Permeability properties of ENaC selectivity filter mutants.

The epithelial Na(+) channel (ENaC), located in the apical membrane of tight epithelia, allows vectorial Na(+) absorption. The amiloride-sensitive ENaC is highly selective for Na(+) and Li(+) ions. There is growing evidence that the short stretch of amino acid residues (preM2) preceding the putative second transmembrane domain M2 forms the outer channel pore with the amiloride binding site and the narrow ion-selective region of the pore. We have shown previously that mutations of the alphaS589 residue in the preM2 segment change the ion selectivity, making the channel permeant to K(+) ions. To understand the molecular basis of this important change in ionic selectivity, we have substituted alphaS589 with amino acids of different sizes and physicochemical properties. Here, we show that the molecular cutoff of the channel pore for inorganic and organic cations increases with the size of the amino acid residue at position alpha589, indicating that alphaS589 mutations enlarge the pore at the selectivity filter. Mutants with an increased permeability to large cations show a decrease in the ENaC unitary conductance of small cations such as Na(+) and Li(+). These findings demonstrate the critical role of the pore size at the alphaS589 residue for the selectivity properties of ENaC. Our data are consistent with the main chain carbonyl oxygens of the alphaS589 residues lining the channel pore at the selectivity filter with their side chain pointing away from the pore lumen. We propose that the alphaS589 side chain is oriented toward the subunit-subunit interface and that substitution of alphaS589 by larger residues increases the pore diameter by adding extra volume at the subunit-subunit interface.     

Topology of the P segments in the sodium channel pore revealed by cysteine mutagenesis.

The P segments of the voltage-dependent Na+ channel line the outer mouth and selectivity filter of the pore. The residues that form the cytoplasmic mouth of the pore of the channel have not been identified. To study the structure of the inner pore mouth, the presumed selectivity filter residues (D400, E755, K1237, and A1529), and three amino acids just amino-terminal to each of these residues in the rat skeletal muscle Na+ channel, were mutated to cysteine and expressed in tsA 201 cells. These amino acids are predicted (by analogy to K+ channels) to be on the cytoplasmic side of the putative selectivity filter residues. Inward and outward Na+ currents were measured with the whole-cell configuration of the patch-clamp technique. Cysteinyl side-chain accessibility was gauged by sensitivity to Cd2+ block and by reactivity with methanethiosulfonate (MTS) reagents applied to both the inside and the outside of the cell. Outward currents through the wild-type and all of the mutant channels were unaffected by internal Cd2+ (100 microM). Similarly, 1 mM methanethiosulfonate ethylammonium (MTSEA) applied to the inside of the membrane did not affect wild-type or mutant outward currents. However, two mutants amino-terminal to the selectivity position in domain III (F1236C and T1235C) and one in domain IV (S1528C) were blocked with high affinity by external Cd2+. The Na+ current through F1236C and S1528C channels was inhibited by MTSEA applied to the outside of the cell. The accessibility of these mutants to externally applied cysteinyl ligands indicates that the side chains of the mutated residues face outward rather than inward. The K+ channel model of the P segments as protein loops that span the selectivity region is not applicable to the Na+ channel.     

水甘油通道蛋白结构和功能与疾病研究进展

水甘油通道蛋白(aquaglyceroporins)属于内在膜蛋白(major intrinsic protein,MIP)家族成员,是既可以运输水分子,又可以运输甘油等小分子物质穿透细胞膜等生理屏障的双功能水通道蛋白(aquaporin,AQP)。研究发现其在人和其他生物的许多脏器中广泛分布,研究表明,其生理状态的改变与许多疾病的发生、发展及预后密切相关。本文从水甘油通道蛋白的基本结构、功能和相关疾病的关系,以及针对甘油水通道蛋白靶向治疗等方面的研究进展进行综述。     

Analysis of side-chain rotamers in transmembrane proteins

We measured the frequency of side-chain rotamers in 14 alpha-helical and 16 beta-barrel membrane protein structures and found that the membrane environment considerably perturbs the rotamer frequencies compared to soluble proteins. Although there are limited experimental data, we found statistically significant changes in rotamer preferences depending on the residue environment. Rotamer distributions were influenced by whether the residues were lipid or protein facing, and whether the residues were found near the N- or C-terminus. Hydrogen-bonding interactions with the helical backbone perturbs the rotamer populations of Ser and His. Trp and Tyr favor side-chain conformations that allow their side chains to extend their polar atoms out of the membrane core, thereby aligning the side-chain polarity gradient with the polarity gradient of the membrane. Our results demonstrate how the membrane environment influences protein structures, providing information that will be useful in the structure prediction and design of transmembrane proteins.     

Insights into the mechanisms of the selectivity filter of Escherichia coli aquaporin Z

Aquaporin Z (AQPZ) is a tetrameric protein that forms water channels in the cell membrane of Escherichia coli. The histidine residue (residue 174) in the selectivity filter (SF) region plays an important role in the transport of water across the membrane. In this work, we perform equilibrium molecular dynamics (MD) simulations to illustrate the gating mechanism of the SF and the influences of residue 174 in two different protonation states: Hsd174 with the proton at Nδ, and Hse174 with the proton at Nε. We calculate the pore radii in the SF region versus the simulation time. We perform steered MD to compute the free-energy profile, i.e., the potential of mean force (PMF) of a water molecule through the SF region. We conduct a quantum mechanics calculation of the binding energy of one water molecule with the residues in the SF region. The hydrogen bonds formed between the side chain of Hsd174 and the side chain of residue 189 (Arg189) play important roles in the selectivity mechanism of AQPZ. The radii of the pores, the hydrogen-bond analysis, and the free energies show that it is easier for water molecules to permeate through the SF region of AQPZ with residue 174 in the Hse state than in the Hsd state.     

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.     

A new subfamily of major intrinsic proteins in plants

The major intrinsic proteins (MIPs) form a large protein family of ancient origin and are found in bacteria, fungi, animals, and plants. MIPs act as channels in membranes to facilitate passive transport across the membrane. Some MIPs allow small polar molecules like glycerol or urea to pass through the membrane. However, the majority of MIPs are thought to be aquaporins (AQPs), i.e., they are specific for water transport. Plant MIPs can be subdivided into the plasma membrane intrinsic protein, tonoplast intrinsic protein, and NOD26-like intrinsic protein subfamilies. By database mining and phylogenetic analyses, we have identified a new subfamily in plants, the Small basic Intrinsic Proteins (SIPs). Comparisons of sequences from the new subfamily with conserved amino acid residues in other MIPs reveal characteristic features of SIPs. Possible functional consequences of these features are discussed in relation to the recently solved structures of AQP1 and GlpF. We suggest that substitutions at conserved and structurally important positions imply a different substrate specificity for the new subfamily.     

Molecular Dynamics Simulation of WSK-3, a Computationally Designed, Water-Soluble Variant of the Integral Membrane Protein KcsA

Poor solubility and low expression levels often make membrane proteins difficult to study. An alternative to the use of detergents to solubilize these aggregation-prone proteins is the partial redesign of the sequence so as to confer water solubility. Recently, computationally assisted membrane protein solubilization (CAMPS) has been reported, where exposed hydrophobic residues on a protein's surface are computationally redesigned. Herein, the structure and fluctuations of a designed, water-soluble variant of KcsA (WSK-3) were studied using molecular dynamics simulations. The root mean square deviation of the protein from its starting structure, where the backbone coordinates are those of KcsA, was 1.8 angstroms. The structure of salt bridges involved in structural specificity and solubility were examined. The preferred configuration of ions and water in the selectivity filter of WSK-3 was consistent with the reported preferences for KcsA. The structure of the selectivity filter was maintained, which is consistent with WSK-3 having an affinity for agitoxin2 comparable to that of wild-type KcsA. In contrast to KcsA, the central cavity's side chains were observed to reorient, allowing water diffusion through the side of the cavity wall. These simulations provide an atomistic analysis of the CAMPS strategy and its implications for further investigations of membrane proteins.     

Extent of the selectivity filter conferred by the sixth transmembrane region in the CFTR chloride channel pore

Point mutations within the pore region of the cystic fibrosis transmembrane conductance regulator (CFTR) Cl^? channel have previously been shown to alter the selectivity of the channel between different anions, suggesting that part of the pore may form an anion 'selectivity filter'. However, the full extent of this selectivity filter region and the location of anion binding sites in the pore are currently unclear. As a result, comparisons between CFTR and other classes of Cl^? channel of known structure are difficult. We compare here the effects of point mutations at each of eight consecutive amino acid residues (arginine 334-serine 341) in the crucial sixth transmembrane region (TM6) of CFTR. Anion selectivity was determined using patch-clamp recording from inside-out membrane patches excised from transiently transfected mammalian cell lines. The results suggest that selectivity is predominantly controlled by a single site involving adjacent residues phenylalanine 337 and threonine 338, and that the selectivity conferred by this 'filter' region is modified by anion binding to flanking sites involving the more extracellular arginine 334 and the more intracellular serine 341. Other residues within this part of the pore play only minor roles in controlling anion permeability and conductance. Our results support a model in which specific TM6 residues make important contributions to a single, localized anion selectivity filter in the CFTR pore, and also contribute to multiple anion binding sites both within and on either side of the filter region.     

Normal-mode refinement of anisotropic thermal parameters for potassium channel KcsA at 3.2 A crystallographic resolution

We report a normal-mode method for anisotropic refinement of membrane-protein structures, based on a hypothesis that the global near-native-state disordering of membrane proteins in crystals follows low-frequency normal modes. Thus, a small set of modes is sufficient to represent the anisotropic thermal motions in X-ray crystallographic refinement. By applying the method to potassium channel KcsA at 3.2 A, we obtained a structural model with an improved fit with the diffraction data. Moreover, the improved electron density maps allowed for large structural adjustments for 12 residues in each subunit, including the rebuilding of 3 missing side chains. Overall, the anisotropic KcsA structure at 3.2 A was systematically closer to a 2.0 A KcsA structure, especially in the selectivity filter. Furthermore, the anisotropic thermal ellipsoids from the refinement revealed functionally relevant structural flexibility. We expect this method to be a valuable tool for structural refinement of many membrane proteins with moderate-resolution diffraction data.     

Extent of the selectivity filter conferred by the sixth transmembrane region in the CFTR chloride channel pore

Point mutations within the pore region of the cystic fibrosis transmembrane conductance regulator (CFTR) Cl(-) channel have previously been shown to alter the selectivity of the channel between different anions, suggesting that part of the pore may form an anion 'selectivity filter'. However, the full extent of this selectivity filter region and the location of anion binding sites in the pore are currently unclear. As a result, comparisons between CFTR and other classes of Cl(-) channel of known structure are difficult. We compare here the effects of point mutations at each of eight consecutive amino acid residues (arginine 334-serine 341) in the crucial sixth transmembrane region (TM6) of CFTR. Anion selectivity was determined using patch-clamp recording from inside-out membrane patches excised from transiently transfected mammalian cell lines. The results suggest that selectivity is predominantly controlled by a single site involving adjacent residues phenylalanine 337 and threonine 338, and that the selectivity conferred by this 'filter' region is modified by anion binding to flanking sites involving the more extracellular arginine 334 and the more intracellular serine 341. Other residues within this part of the pore play only minor roles in controlling anion permeability and conductance. Our results support a model in which specific TM6 residues make important contributions to a single, localized anion selectivity filter in the CFTR pore, and also contribute to multiple anion binding sites both within and on either side of the filter region.