Switch from an aquaporin to a glycerol channel by two amino acids substitution

The MIP (major intrinsic protein) proteins constitute a channel family of currently 150 members that have been identified in cell membranes of organisms ranging from bacteria to man. Among these proteins, two functionally distinct subgroups are characterized: aquaporins that allow specific water transfer and glycerol channels that are involved in glycerol and small neutral solutes transport. Since the flow of small molecules across cell membranes is vital for every living organism, the study of such proteins is of particular interest. For instance, aquaporins located in kidney cell membranes are responsible for reabsorption of 150 liters of water/day in adult human. To understand the molecular mechanisms of solute transport specificity, we analyzed mutant aquaporins in which highly conserved residues have been substituted by amino acids located at the same positions in glycerol channels. Here, we show that substitution of a tyrosine and a tryptophan by a proline and a leucine, respectively, in the sixth transmembrane helix of an aquaporin leads to a switch in the selectivity of the channel, from water to glycerol.     

The structure of aquaporins

The ubiquitous members of the aquaporin (AQP) family form transmembrane pores that are either exclusive for water (aquaporins) or are also permeable for other small neutral solutes such as glycerol (aquaglyceroporins). The purpose of this review is to provide an overview of our current knowledge of AQP structures and to describe the structural features that define the function of these membrane pores. The review will discuss the mechanisms governing water conduction, proton exclusion and substrate specificity, and how the pore permeability is regulated in different members of the AQP family.     

Prediction of functional residues in water channels and related proteins.

In this paper, we present an updated classification of the ubiquitous MIP (Major Intrinsic Protein) family proteins, including 153 fully or partially sequenced members available in public databases. Presently, about 30 of these proteins have been functionally characterized, exhibiting essentially two distinct types of channel properties: (1) specific water transport by the aquaporins, and (2) small neutral solutes transport, such as glycerol by the glycerol facilitators. Sequence alignments were used to predict amino acids and motifs discriminant in channel specificity. The protein sequences were also analyzed using statistical tools (comparisons of means and correspondence analysis). Five key positions were clearly identified where the residues are specific for each functional subgroup and exhibit high dissimilar physico-chemical properties. Moreover, we have found that the putative channels for small neutral solutes clearly differ from the aquaporins by the amino acid content and the length of predicted loop regions, suggesting a substrate filter function for these loops. From these results, we propose a signature pattern for water transport.     

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".     

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.     

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.     

Architecture and selectivity in aquaporins: 2.5 a X-ray structure of aquaporin Z

Aquaporins are a family of water and small molecule channels found in organisms ranging from bacteria to animals. One of these channels, the E. coli protein aquaporin Z (AqpZ), has been shown to selectively conduct only water at high rates. We have expressed, purified, crystallized, and solved the X-ray structure of AqpZ. The 2.5 A resolution structure of AqpZ suggests aquaporin selectivity results both from a steric mechanism due to pore size and from specific amino acid substitutions that regulate the preference for a hydrophobic or hydrophilic substrate. This structure provides direct evidence on the molecular mechanisms of specificity between water and glycerol in this family of channels from a single species. It is to our knowledge the first atomic resolution structure of a recombinant aquaporin and so provides a platform for combined genetic, mutational, functional, and structural determinations of the mechanisms of aquaporins and, more generally, the assembly of multimeric membrane proteins.     

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

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

Accurate Simulation and Detection of Coevolution Signals in Multiple Sequence Alignments

Background

While the conserved positions of a multiple sequence alignment (MSA) are clearly of interest, non-conserved positions can also be important because, for example, destabilizing effects at one position can be compensated by stabilizing effects at another position. Different methods have been developed to recognize the evolutionary relationship between amino acid sites, and to disentangle functional/structural dependencies from historical/phylogenetic ones.

Methodology/Principal Findings

We have used two complementary approaches to test the efficacy of these methods. In the first approach, we have used a new program, MSAvolve, for the in silico evolution of MSAs, which records a detailed history of all covarying positions, and builds a global coevolution matrix as the accumulated sum of individual matrices for the positions forced to co-vary, the recombinant coevolution, and the stochastic coevolution. We have simulated over 1600 MSAs for 8 protein families, which reflect sequences of different sizes and proteins with widely different functions. The calculated coevolution matrices were compared with the coevolution matrices obtained for the same evolved MSAs with different coevolution detection methods. In a second approach we have evaluated the capacity of the different methods to predict close contacts in the representative X-ray structures of an additional 150 protein families using only experimental MSAs.

Conclusions/Significance

Methods based on the identification of global correlations between pairs were found to be generally superior to methods based only on local correlations in their capacity to identify coevolving residues using either simulated or experimental MSAs. However, the significant variability in the performance of different methods with different proteins suggests that the simulation of MSAs that replicate the statistical properties of the experimental MSA can be a valuable tool to identify the coevolution detection method that is most effective in each case.     

Structure, function and translational relevance of aquaporin dual water and ion channels

Aquaporins have been assumed to be selective for water alone, and aquaglyceroporins are accepted as carrying water and small uncharged solutes including glycerol. This review presents an expanded view of aquaporins as channels with more complex mechanisms of regulation and diverse repertoires of substrate permeabilities than were originally appreciated in the early establishment of the field. The role of aquaporins as dual water and gated ion channels is likely to have physiological and potentially translational relevance, and can be evaluated with newly developed molecular and pharmacological tools. Ion channel activity has been shown for Aquaporins -0, -1, and -6, Drosphila Big Brain, and plant Nodulin-26. Although the concept of ion channel function in aquaporins remains controversial, research advances are beginning to define not only the ion channel function but also the detailed molecular mechanisms that govern and mediate the multifunctional capabilities. With regard to physiological relevance, the adaptive benefit of expression of ion channel activity in aquaporins, implied by amino acid sequence conservation of the ion channel gating domains, suggests they provide more than water or glycerol and solute transport. Dual ion and water channels are of interest for understanding the modulation of transmembrane fluid gradients, volume regulation, and possible signal transduction in tissues expressing classes of aquaporins that have the dual function capability. Other aquaporin classes might be found in future work to have ion channel activities, pending identification of the possible signaling pathways that could govern activation.     

Homology modeling of representative subfamilies of Arabidopsis major intrinsic proteins. Classification based on the aromatic/arginine selectivity filter

Major intrinsic proteins (MIPs) are a family of membrane channels that facilitate the bidirectional transport of water and small uncharged solutes such as glycerol. The 35 full-length members of the MIP family in Arabidopsis are segregated into four structurally homologous subfamilies: plasma membrane intrinsic proteins (PIPs), tonoplast intrinsic proteins (TIPs), nodulin 26-like intrinsic membrane proteins (NIPs), and small basic intrinsic proteins (SIPs). Computational methods were used to construct structural models of the putative pore regions of various plant MIPs based on homology modeling with the atomic resolution crystal structures of mammalian aquaporin 1 and the bacterial glycerol permease GlpF. Based on comparisons of the narrow selectivity filter regions (the aromatic/Arg [ar/R] filter), the members of the four phylogenetic subfamilies of Arabidopsis MIPs can be classified into eight groups. PIPs possess a uniform ar/R signature characteristic of high water transport aquaporins, whereas TIPs are highly diverse with three separate conserved ar/R regions. NIPs possess two separate conserved ar/R regions, one that is similar to the archetype, soybean (Glycine max) nodulin 26, and another that is characteristic of Arabidopsis NIP6;1. The SIP subfamily possesses two ar/R subgroups, characteristic of either SIP1 or SIP2. Both SIP ar/R residues are divergent from all other MIPs in plants and other kingdoms. Overall, these findings suggest that higher plant MIPs have a common fold but show distinct differences in proposed pore apertures, potential to form hydrogen bonds with transported molecules, and amphiphilicity that likely results in divergent transport selectivities.     

Architecture and Selectivity in Aquaporins: 2.5 Å X-Ray Structure of Aquaporin Z

Aquaporins are a family of water and small molecule channels found in organisms ranging from bacteria to animals. One of these channels, the E. coli protein aquaporin Z (AqpZ), has been shown to selectively conduct only water at high rates. We have expressed, purified, crystallized, and solved the X-ray structure of AqpZ. The 2.5 Å resolution structure of AqpZ suggests aquaporin selectivity results both from a steric mechanism due to pore size and from specific amino acid substitutions that regulate the preference for a hydrophobic or hydrophilic substrate. This structure provides direct evidence on the molecular mechanisms of specificity between water and glycerol in this family of channels from a single species. It is to our knowledge the first atomic resolution structure of a recombinant aquaporin and so provides a platform for combined genetic, mutational, functional, and structural determinations of the mechanisms of aquaporins and, more generally, the assembly of multimeric membrane proteins.     

Structural insights into eukaryotic aquaporin regulation

Aquaporin-mediated water transport across cellular membranes is an ancient, ubiquitous mechanism within cell biology. This family of integral membrane proteins includes both water selective pores (aquaporins) and transport facilitators of other small molecules such as glycerol and urea (aquaglyceroporins). Eukaryotic aquaporins are frequently regulated post-translationally by gating, whereby the rate of flux through the channel is controlled, or by trafficking, whereby aquaporins are shuttled from intracellular storage sites to the plasma membrane. A number of high-resolution X-ray structures of eukaryotic aquaporins have recently been reported and the new structural insights into gating and trafficking that emerged from these studies are described. Basic structural themes reoccur, illustrating how the problem of regulation in diverse biological contexts builds upon a limited set of possible solutions.     

Architecture and Selectivity in Aquaporins: 2.5 Å X-Ray Structure of Aquaporin Z

Aquaporins are a family of water and small molecule channels found in organisms ranging from bacteria to animals. One of these channels, the E. coli protein aquaporin Z (AqpZ), has been shown to selectively conduct only water at high rates. We have expressed, purified, crystallized, and solved the X-ray structure of AqpZ. The 2.5 Å resolution structure of AqpZ suggests aquaporin selectivity results both from a steric mechanism due to pore size and from specific amino acid substitutions that regulate the preference for a hydrophobic or hydrophilic substrate. This structure provides direct evidence on the molecular mechanisms of specificity between water and glycerol in this family of channels from a single species. It is to our knowledge the first atomic resolution structure of a recombinant aquaporin and so provides a platform for combined genetic, mutational, functional, and structural determinations of the mechanisms of aquaporins and, more generally, the assembly of multimeric membrane proteins.     

Folding and stability of the aquaglyceroporin GlpF: Implications for human aqua(glycero)porin diseases

Aquaporins are highly selective polytopic transmembrane channel proteins that facilitate the permeation of water across cellular membranes in a large diversity of organisms. Defects in aquaporin function are associated with common diseases, such as nephrogenic diabetes insipidus, congenital cataract and certain types of cancer. In general, aquaporins have a highly conserved structure; from prokaryotes to humans. The conserved structure, together with structural dynamics and the structural framework for substrate selectivity is discussed. The folding pathway of aquaporins has been a topic of several studies in recent years. These studies revealed that a conserved protein structure can be reached by following different folding pathways. Based on the available data, we suggest a complex folding pathway for aquaporins, starting from the insertion of individual helices up to the formation of the tetrameric aquaporin structure. The consequences of some known mutations in human aquaporin-encoding genes, which most likely affect the folding and stability of human aquaporins, are discussed.     

Identification of functional motions in the adenylate kinase (ADK) protein family by computational hybrid approaches

Based on integrative computational hybrid approaches that combined statistical coupling analysis (SCA), molecular dynamics (MD), and normal mode analysis (NMA), evolutionarily coupled residues involved in functionally relevant motion in the adenylate kinase protein family were identified. The hybrids identified four top-ranking site pairs that belong to a conserved hydrogen bond network that is involved in the enzyme's flexibility. A second group of top-ranking site pairs was identified in critical regions for functional dynamics, such as those related to enzymatic turnover. The high consistency of the results obtained by SCA with NMA (SCA.NMA) and by SCA.MD hybrid analyses suggests that suitable replacement of the matrix of cross-correlation analysis of atomic fluctuations (derived by using NMA) with those based on MD contributes to the identification of such sites by means of a fast computational calculation. The analysis presented here strongly supports the hypothesis that evolutionary forces, such as coevolution at the sequence level, have promoted functional dynamic properties of the adenylate kinase protein family. Finally, these hybrid approaches can be used to identify, at the residue level, protein motion coordination patterns not previously observed, such as in hinge regions.     

Genome-Wide Identification and Expression Analysis of Aquaporins in Tomato

The family of aquaporins, also called water channels or major intrinsic proteins, is characterized by six transmembrane domains that together facilitate the transport of water and a variety of low molecular weight solutes. They are found in all domains of life, but show their highest diversity in plants. Numerous studies identified aquaporins as important targets for improving plant performance under drought stress. The phylogeny of aquaporins is well established based on model species like Arabidopsis thaliana, which can be used as a template to investigate aquaporins in other species. In this study we comprehensively identified aquaporin encoding genes in tomato (Solanum lycopersicum), which is an important vegetable crop and also serves as a model for fleshy fruit development. We found 47 aquaporin genes in the tomato genome and analyzed their structural features. Based on a phylogenetic analysis of the deduced amino acid sequences the aquaporin genes were assigned to five subfamilies (PIPs, TIPs, NIPs, SIPs and XIPs) and their substrate specificity was assessed on the basis of key amino acid residues. As ESTs were available for 32 genes, expression of these genes was analyzed in 13 different tissues and developmental stages of tomato. We detected tissue-specific and development-specific expression of tomato aquaporin genes, which is a first step towards revealing the contribution of aquaporins to water and solute transport in leaves and during fruit development.     

Using information theory to search for co-evolving residues in proteins

MOTIVATION: Some functionally important protein residues are easily detected since they correspond to conserved columns in a multiple sequence alignment (MSA). However important residues may also mutate, with compensatory mutations occurring elsewhere in the protein, which serve to preserve or restore functionality. It is difficult to distinguish these co-evolving sites from other non-conserved sites. RESULTS: We used Mutual Information (MI) to identify co-evolving positions. Using in silico evolved MSAs, we examined the effects of the number of sequences, the size of amino acid alphabet and the mutation rate on two sources of background MI: finite sample size effects and phylogenetic influence. We then assessed the performance of various normalizations of MI in enhancing detection of co-evolving positions and found that normalization by the pair entropy was optimal. Real protein alignments were analyzed and co-evolving isolated pairs were often found to be in contact with each other. AVAILABILITY: All data and program files can be found at http://www.biochem.uwo.ca/cgi-bin/CDD/index.cgi     

Methylation of aquaporins in plant plasma membrane

A thorough analysis, using MS, of aquaporins expressed in plant root PM (plasma membrane) was performed, with the objective of revealing novel post-translational regulations. Here we show that the N-terminal tail of PIP (PM intrinsic protein) aquaporins can exhibit multiple modifications and is differentially processed between members of the PIP1 and PIP2 subclasses. Thus the initiating methionine was acetylated or cleaved in native PIP1 and PIP2 isoforms respectively. In addition, several residues were detected to be methylated in PIP2 aquaporins. Lys3 and Glu6 of PIP2;1, one of the most abundant aquaporins in the PM, occurred as di- and mono-methylated residues respectively. Ectopic expression in Arabidopsis suspension cells of PIP2;1, either wild-type or with altered methylation sites, revealed an interplay between methylation at the two sites. Measurements of water transport in PM vesicles purified from these cells suggested that PIP2;1 methylation does not interfere with the aquaporin intrinsic water permeability. In conclusion, the present study identifies methylation as a novel post-translational modification of aquaporins, and even plant membrane proteins, and may represent a critical advance towards the identification of new regulatory mechanisms of membrane transport.     

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.