New subfamilies of major intrinsic proteins in fungi suggest novel transport properties in fungal channels: implications for the host-fungal interactions
© Verma et al.; BioMed Central Ltd 2014
Received: 9 May 2014
Accepted: 24 July 2014
Published: 12 August 2014
Aquaporins (AQPs) and aquaglyceroporins (AQGPs) belong to the superfamily of Major Intrinsic Proteins (MIPs) and are involved in the transport of water and neutral solutes across the membranes. MIP channels play significant role in plant-fungi symbiotic relationship and are believed to be important in host-pathogen interactions in human fungal diseases. In plants, at least five major MIP subfamilies have been identified. Fungal MIP subfamilies include orthodox aquaporins and five subgroups within aquaglyceroporins. XIP subfamily is common to both plants and fungi. In this study, we have investigated the extent of diversity in fungal MIPs and explored further evolutionary relationships with the plant MIP counterparts.
We have extensively analyzed the available fungal genomes and examined nearly 400 fungal MIPs. Phylogenetic analysis and homology modeling exhibit the existence of a new MIP cluster distinct from any of the known fungal MIP subfamilies. All members of this cluster are found in microsporidia which are unicellular fungal parasites. Members of this family are small in size, charged and have hydrophobic residues in the aromatic/arginine selectivity filter and these features are shared by small and basic intrinsic proteins (SIPs), one of the plant MIP subfamilies. We have also found two new subfamilies (δ and γ2) within the AQGP group. Fungal AQGPs are the most diverse and possess the largest number of subgroups. We have also identified distinguishing features in loops E and D in the newly identified subfamilies indicating their possible role in channel transport and gating.
Fungal SIP-like MIP family is distinct from any of the known fungal MIP families including orthodox aquaporins and aquaglyceroporins. After XIPs, this is the second MIP subfamily from fungi that may have possible evolutionary link with a plant MIP subfamily. AQGPs in fungi are more diverse and possess the largest number of subgroups. The aromatic/arginine selectivity filter of SIP-like fungal MIPs and the δ AQGPs are unique, hydrophobic in nature and are likely to transport novel hydrophobic solutes. They can be attractive targets for developing anti-fungal drugs. The evolutionary pattern shared with their plant counterparts indicates possible involvement of new fungal MIPs in plant-fungi symbiosis and host-pathogen interactions.
KeywordsAquaglyceroporin Aquaporin Aromatic/arginine selectivity filter Host-pathogen interactions Intra-helical salt-bridge Membrane transport Plant-fungi symbiosis Selective transport
The superfamily of Major Intrinsic Proteins (MIPs) contains channel proteins that transport water and other neutral solutes –. These integral membrane proteins are found from bacteria to humans and are abundantly present in plants ,. Aquaporin (AQP) and aquaglyceroporin (AQGP) are the prototype members of the MIP superfamily ,. While in mammals three major subfamilies are detected , plant MIPs are found to have at least five subfamilies –. In humans, MIP members play significant role in kidney nephron, epithelial fluid secretion, maintaining brain water balance, cell migration, skin hydration, adipocyte metabolism and neuroexcitation ,–. They are implicated in various human diseases such as glaucoma, epilepsy, cancer, obesity, nephrogenic diabetes insipidus and neuromyelitis optica –. In plants, multiple roles have been recognized in terms of plant development, growth and physiology ,,. They have been shown to be important in stress tolerance in plants and they play vital role in their interactions with soil microorganisms, fungi and pathogens –. In addition to water and glycerol, expression studies in Xenopus oocytes showed that MIPs transport diverse neutral solutes and gases such as ammonia and CO2,,. MIP members also serve as an important component in host-parasite interactions and several protozoan parasite aquaporins have been identified including those from Plasmodium, Trypanosoma and Leishmania species –. Functional studies suggest that protozoan aquaporins are likely to transport both water and glycerol efficiently . These parasite channel proteins could either serve as potential drug targets or vehicles for transporting cytotoxic compounds ,,. We have recently identified more than 1000 MIP genes from more than 340 organisms and details about the genes, protein products and structural models are available in MIPModDB database .
Symbiotic relationship between plants and fungi occurs through mycorrhiza and these interactions help in translocating soil nutrients to the host plants ,. This mutualistic association also aids in the transfer of organic carbon from plants to fungal partners . In the mycorrhized plants, movement of water and other nutrients are greatly influenced by the expression of both plant and fungal aquaporins and this plays a major role in the drought resistance of plants ,. The expression of plant aquaporins in mycorrhized and non-mycorrhized forms have been investigated in several plants under different stress conditions ,,–. On the contrary, very few studies have been carried out on the role of fungal aquaporins in mycorrhizal symbiosis ,. In addition to their role in root water transport in plants through symplastic pathway, fungal aquaporins present in pathogenic fungal species may act as attractive targets for antifungal drugs ,,. With increasing drug resistance reported in human fungal pathogens ,, it has become even more important to identify new drug targets in fungal organisms.
Fungal aquaporins have been identified in individual species such as Laccaria bicolor (an ectomycorrhizal fungus) Glomus intraradices (an arbuscular mycorrhizal fungus)  and Encephalitozoon cuniculi (a microsporidia pathogenic to humans) . Their transport properties, gene expression profiles and role in mycorrhizia or pathogenesis have been investigated. Available literature shows that only a few fungal MIPs have been studied in detail. Fungal species are diverse and form a large group of eukaryotic organisms and hence it is anticipated that the MIP channels present in fungi also could be diverse with different transport specificities and unique regulatory mechanisms and/or physiological functions. In an earlier study, Pettersson et al.  analyzed genomes of yeast and filamentous fungi. Their phylogenetic analysis with a small number of 55 MIPs revealed three different subgroups within fungal AQGPs in addition to members belonging to orthodox AQP family. A new subfamily of fungal aquaporins, called X-intrinsic proteins (XIPs), has been discovered in our laboratory which is found to be common between fungi and plants . More recently, Dietz et al.  and Zwiazek and coworkers  have identified even a larger number of 135 and 229 fungal MIPs respectively from database searches. Both their phylogenetic studies indicated at least four clusters including the orthodox AQPs, XIPs and two to four distinct subgroups of AQGPs. With the improved technologies and falling costs of sequencing, fungal genomes are rapidly sequenced and they are available in public databases. In this study, we have systematically searched completely and partially sequenced fungal genomes and identified a large number of fungal MIPs. With a comprehensive analysis of these MIP sequences using different bioinformatics tools, we have identified a new group of fungal MIPs which forms a distinct cluster separately from the orthodox AQPs and AQGPs. Members of this new family share several features with plant SIP (small and basic intrinsic proteins) subfamily. Our analysis also revealed two new subfamilies of aquaglyceroporins which are separate from the other subgroups of AQGPs. Distinct features of these three new families of fungal MIPs, their possible evolutionary relationship with plants and lower organisms and the nature of solutes that are likely to be transported across these channels are presented in detail.
Identification of fungal MIPs and their taxonomic distributions
Taxonomical distribution of fungal MIPs
Taxonomical distribution of different clusters and the subgroups are also indicated in Figure 2. The orthodox fungal AQPs are present mainly in the species of Ascomycota (yeast and filamentous) and the rest are found in filamentous Basidiomycota. Only two examples are found in Glomeromycota and a lone representative is from Mucoromycotina. All XIPs are found in Ascomycota filamentous fungi. Among the fungal AQGPs, we have identified two additional subgroups in this study apart from the known subgroups. They have different taxonomic distributions. The phylogenetic tree of all fungal AQGPs is presented in Figure 2b. The Fps1-like AQGPs are found invariably in yeast Ascomycota while the Yfl054-like AQGPs, also known as facultative aquaporins , are predominantly identified in yeast and filamentous Ascomycota. The other major subgroups classified as α and β are mainly found in Ascomycota (filamentous) and Basidomycota (filamentous) respectively. There are two small subgroups which we have designated as γ1 and γ2. The γ1 cluster has been previously recognized  and is mainly found in the species of Mucoromycotina. The second cluster γ2 is newly recognized in this study and is found in filamentous Ascomycota. The major fungal AQGP subgroup newly identified from this analysis consists of 49 members and the members of this group are found in filamentous Ascomycota and Basidiomycota (yeast & filamentous). We have named this cluster as δ. The new fungal MIP group, which is distinct from orthodox AQPs and AQGPs, are found exclusively in Microsporidia and these parasitic species are known to infect their animal hosts. In this context, it is even more important to understand and characterize this MIP group. We have designated this group as “SIP-like” fungal MIPs since they possess several characteristics of plant SIP subgroup  (see below). From this analysis, it appears that some fungal MIP subfamilies such as XIPs and the newly identified SIP-like members are exclusive to a single fungal phylum (Figure 2). However, one has to keep in mind that the present analysis was carried out with limited number of available fungal genome sequences. When more fungal genome sequences are available, it is possible that MIP members from these subfamilies may be found in other fungal phyla also.
Characterization of new fungal MIP cluster that is distinct from AQPs, AQGPs and XIPs
Sixteen fungal MIPs form a cluster separate from AQPs, AQGPs and XIPs. Average pairwise sequence identity of all these 16 sequences is about 54% and similarity is close to 70%. We also calculated the average pairwise sequence identity between the sequences from the newly identified cluster and all the sequences belonging to the clusters of orthodox AQPs and XIPs. A similar exercise was carried out for all the subgroups of AQGPs also. The results of pairwise sequence analysis are given in Additional file 2: Table S2. Inter-group average sequence identity varies from 16 to 27% and similarity varies from 32 to 44% indicating that this group of sequences exhibits a large degree of divergence from all other fungal MIP sequences.
Analysis of Ar/R selectivity filter residues
Aromatic/arginine selectivity filter residues and substitutions in the conserved NPA motifs in different fungal MIP subgroups
Aromatic/arginine selectivity filterc
Fungal MIP subgroupa
Number of sequences/Total number of sequencesb
NPA substitution (%)d
FPS-1 like AQGPs
Aromatic/arginine selectivity filter residues and substitutions in the conserved NPA motifs in plant subfamilies TIPs and SIPs
Aromatic/arginine selectivity filterc
Plant MIP subfamilya
Number of sequences/Total number of sequencesb
NPA substitutions (%)d
Substitutions in the highly conserved NPA boxes
Characterization of δ subfamily of fungal AQGPs
The newly identified SIP-like fungal MIP family is distinct from AQPs, AQGPs and XIPs within fungal MIP family. Within the fungal AQGP group, five different subgroups have already been identified in a recently published study . In the current study, we have also identified two additional subgroups distinct from the previously identified fungal AQGP clusters. The biggest cluster is δ with 49 members (Figure 2b). The average pairwise sequence identity and similarity within this group are about 51% and 67% respectively. This indicates that the sequences within the group are closely related. However, the same cannot be said for AQGP members from other subgroups. The inter-group average sequence identity between δ cluster and other AQGP groups varies from 24% to 30% indicating that δ subgroup members have diverged significantly and are distantly related to other fungal AQGP groups (Additional file 2: Table S2). Hence, our study mainly focused on features that are unique to δ cluster AQGPs and hence, we present our analysis by comparing δ subgroup with other fungal AQGP groups.
γ2 cluster of AQGPs
This small group contains only 5 members forming a separate cluster within the fungal AQGPs (Figure 2b) and all of them belong to filamentous Ascomycota. Their selectivity filter contains a Tyr residue at LE1 position which seems to be a feature shared by many members of α and β clusters (Table 1). The TM2 position is preferred by a bulky residue (Trp/Met/Ile) and LE2 position is occupied by Arg. There is no strikingly different feature observed in this small group. Perhaps, a larger dataset would help to identify unique characteristics of small clusters like γ1 and γ2.
Fungal MIPs are as diverse as plant MIPs
The diversity of plant MIPs has been established in many earlier studies –. Only few studies identified and addressed the extent of diversity in fungal MIPs ,,. In the present study, by identifying two major clusters and one minor cluster within the fungal MIP superfamily, we have determined that fungal MIPs are as diverse as their plant counterparts. Till date, the only MIP subgroup common to both plant and fungal species has been found to be XIPs . However, by systematically searching the genome sequences of fungal species, we have recognized another group in fungal MIPs which has common features with plant SIP subfamily. Although phylogenetic analysis cluster the newly identified “SIP-like” fungal MIPs separately from the plant SIPs, several features are shared by the two MIP families from the two different species group. The small size of the channel protein, the nature of selectivity filter residues and substitutions in the conserved NPA boxes indicate that the new fungal MIP cluster shares some characteristic features with plant SIPs.
Unique selectivity filter suggests that SIP-like MIPs are likely to transport larger hydrophobic solutes
Very few functional studies have been carried on fungal MIPs ,,,. Experimental, structural and simulation studies have emphasized the role of the aromatic/arginine selectivity filter in MIP channel’s transport and selectivity –,–. In this context, it would be interesting to at least make a speculation regarding the transport properties of SIP-like fungal MIPs based on the selectivity filter residues. Even for plant SIPs, not many studies have been reported regarding their function. When Arabidopsis SIPs (AtSIPs) were expressed in yeast, at least one of the two members displayed water channel activity . This water transport has been reported in spite of the fact that the aromatic/arginine selectivity filter of AtSIPs lacks the arginine residue. In a recent study , an aquaporin AQP2 from the protozoan parasite Trypanosoma brucei gambiense or Trypanosoma brucei rhodesiense has been shown to have an unusual selectivity filter comprising of large hydrophobic residues IVLL. T. brucei causes sleeping sickness in humans and with all hydrophobic residues in the selectivity filter without any aromatic and arginine residue, this aquaporin has been shown to play a vital role in increased cross resistance to the drugs melarsoprol and pentamidine . Based on their studies, Baker et al.  have suggested that T. brucei aquaporin AQP2 may function as a transporter of pentamidine and melarsoprol drugs. Both drugs with molecular mass above 330 Da are larger molecules compared to glycerol (92 Da). It is suggested that the unique selectivity filter of AQP2 from T. Brucei with bulky hydrophobic residues may facilitate the transport of large molecules.
We have also examined all non-plant MIPs whose selectivity filter lacks arginine residue. In MIPModDB  which has more than 1000 MIPs, we could find only 52 examples in which arginine is absent in the selectivity filter. Among them, only 17 of them including mammalian AQP12 have features similar to that found in SIP-like fungal MIPs (two bulky hydrophobic residues and two small residues in the selectivity filter). We found that in mammalian AQP12, both TM2 and LE2 positions are occupied by Leu residues and Ala is found in TM5 and LE1 positions. Although no direct functional studies are reported for AQP12, it has been speculated that AQP12 may regulate proper secretion of pancreatic fluid or digestive enzymes ,. Thus, plant SIPs, AQP12 and AQP2 of T. Brucei also have selectivity filters similar to fungal SIP-like MIPs with all four positions occupied by hydrophobic residues and arginine is absent. Based on the available literature on other MIPs with similar features, we conclude that the SIP-like fungal MIPs are likely to transport solutes which are novel, hydrophobic and larger in size. The same argument may hold good for the plant SIPs also.
Loop D in SIP-like members may play a role in regulating the channel function
The SIP-like fungal MIPs exhibit some characteristic features in loop D which connects the TM4 and TM5 helices. The average length of loop D is the longest in SIP-like cluster compared to all other fungal MIP subgroups (Figure 5). Moreover, large number of charged residues is found in loop D in SIP-like members. This observation indicates that loop D in SIP-like MIPs can play a special role in the function of these channels. Loop D in some mammalian and plant aquaporin channels have been shown to be important in channel gating and regulation of channel function. Experimental and computational studies in AQP1, spinach aquaporin SoPIP2;1 and AQP4 have suggested a role for loop D in channel gating , or binding site for channel agonists  or antagonists . Specific residues in loop D which respond to pH have been shown to be responsible for the closed state in AQP1 channel . The arginine-rich loop D in AQP1 has also been shown to bind cGMP which seems to initiate channel opening . Loop D in AQP4 serves as a possible region for metal-binding activity and has also been implicated in regulating the channel activity ,. Structural studies on spinach aquaporin characterized the open and closed states and exhibited two different conformations of loop D revealing that the gating mechanism at molecular level involves a large movement of this loop . The available experimental and simulation studies provide a clue for the role of loop D in fungal SIP-like MIPs. The charged residues in this loop could be involved in binding ions or molecules that could trigger the channel opening and closing. The longer loop in this subfamily can easily adopt different conformations that can regulate the channel transport. Both possibilities can also simultaneously exist implying that the loop D in SIP-like channels should be given specific attention while investigating the structure-function relationship of these subgroup members.
SIP-like channels could be attractive targets for some anti-fungal diseases
It is interesting to note that all the SIP-like members are from the lower fungi Microsporidia. We did not find even a single example from higher fungi. Although Micorsporidia were earlier classified as part of the kingdom protozoa, these spore-forming unicellular parasites are now known to be fungi. They infect both humans and animals and cause microsporidiosis in humans . As human pathogens they cause diarrhea and infections in immune-compromised individuals (Additional file 6: Table S4). Hence, understanding the host-pathogen interactions and the infection mechanism is extremely important that will aid in developing anti-fungal drugs. SIP-like members in Microsporidia are unique in several respects. They are small in size and the channels formed by them have hydrophobic selectivity filters. The loop D in these channels is longer and is highly charged. The nature of molecules that are transported through these channels and the role of loop D in regulating the channel function have to be first experimentally investigated. Although human AQP12 has similar selectivity filter, phylogenetic analysis clearly clusters them into separate clades. Hence, SIP-like MIP channels in disease causing Microsporidial fungi can be considered as attractive drug targets.
Fungal AQGPs have the largest number of subfamilies
With the identification of δ and γ2 AQGP subfamilies in this study, fungal AQGPs can be classified into seven different subgroups (Fps1-like, facultative, α, β, γ1, γ2 and δ). To our knowledge, the kingdom of fungi possesses the largest subfamilies within the major group of AQGPs. This makes the fungal AQGPs as the most diverse among all known AQGPs compared to any major organism groups. In majority of fungal AQGPs, Trp is found to be preferred in the TM2 position of the selectivity filter and Arg is absolutely conserved in LE2 position of all fungal AQGPs (Table 1). Most of the variations are seen in TM5 and LE1 positions. Members of γ1 AQGPs and the newly identified δ subgroup of AQGPs defy the trend observed in TM2 position. While Trp is absent in members of both the clusters, a small polar residue (Thr, Ser or Asn) is found in γ1 at TM2 position. However in δ subgroup, bulky residues of hydrophobic nature (Phe, Tyr, Met, Val, Ile and Leu) are found in TM2, TM5 and LE1 positions. As in SIP-like MIPs, the selectivity filter of δ AQGPs appears to be unique. We searched MIPModDB database with more than 1000 MIP members looking for examples with similar features in their selectivity filters. We could find only four MIPs in which the channel selectivity filter is formed by three bulky hydrophobic residues and the fourth position is occupied by an arginine residue.
As mentioned earlier, δ subgroup of AQGPs has been found in filamentous Ascomycota and yeast and filamentous Basidiomycota. Many examples from this species group are known to be plant pathogens  (Additional file 7: Table S5). Members belonging to this group have been identified to cause serious diseases in wide range of plants including agriculturally important cereals such as rice, wheat, rye and barley. They also affect economically important crops like cotton. These pathogens display resistance to fungicides and δ subgroup of AQGPs can be considered as an important target to contain these plant pathogens.
Loop E as a marker to study the evolution of MIPs
Loop E provides two residues for the selectivity filter. It also possesses one of the two conserved NPA motifs. While the other conserved NPA motif resides in loop B, majority of the positions in this loop show variations (Additional file 5: Figure S2). In the case of loop E, this region seems to show some distinctive features depending upon the subgroups. Six out of seven AQGPs exhibit a conserved motif RDxGPR next to the NPA motif in loop E (Figure 4). A very high conservation of acidic and basic residues separated by four positions will enable the formation of salt-bridge interaction in the half-helix. However, as mentioned earlier, the presence of two helix destabilizing residues Gly and Pro will weaken the stability of the same helix. In δ AQGPs, the motif present is RDxGxR and the Pro residue next to Gly is mostly substituted by other residues. Hence, we can think of the half-helix in δ AQGPs as relatively more stable compared to all other fungal AQGP subgroups. In XIPs, the acidic and basic residues are replaced by a strictly conserved cysteine and mostly bulky residues (Tyr, Phe or Met) respectively. The motif in XIPs is found to be RCx[G/A]xx. SIP-like members exhibit no such conserved motif in this region. In orthodox AQPs, the acidic and basic residues are replaced by small residues and the conservation pattern in the same region is R[S/A]xG[P/A][A/C/D/S]. It should be noted that the absence of intra-helical salt-bridge and the conservation of Gly and to some extent Pro in AQPs will render this half-helix relatively less stable compared to the same region in AQGPs in general and δ-AQGPs in particular. An extensive bioinformatics analysis of MIPs from microbial organisms to mammals in loop E is necessary and this may spring some surprising insights into the evolution of MIPs.
The role of different fungal proteins in the symbiotic relationship with plants needs to be clearly established. Similarly, host-pathogen interactions at the time of fungi-induced infection require understanding at the molecular level. In this context, understanding the evolution, function and diversity of fungal MIP channels has become very significant. Experimental studies clearly reveal that MIP channels are important players in plant-fungi interactions. Identification of new targets for anti-fungal drugs is a major goal for several human diseases caused by fungi and MIP channels could be considered as attractive anti-fungal targets. In the present study, analysis of fungal genome sequences has identified additional MIP channels in different species groups within the kingdom of fungi. Phylogenetic analysis of nearly 400 fungal MIP channels has revealed the existence of a new MIP cluster completely distinct from the orthodox AQP and AQGP channels. Further sequence analysis and homology modeling studies indicate features that are shared between the new fungal MIP cluster and the plant SIP subfamily. The size of the protein, chemical nature of the residues that form the narrow aromatic/arginine selectivity filter and the substitutions found in the conserved NPA motifs are the common characteristics between the new family and the plant SIP channels. Hence in addition to the XIPs, this “SIP-like” fungal MIP channels can possibly be another evolutionary link between the plants and fungi. Since SIP-like channels are observed only in Microsporidia which are the unicellular parasites, they can also be considered as drug targets for developing anti-fungal drugs in human infections caused by fungi.
We have also identified one major subgroup and another minor cluster within the fungal AQGP family. The δ-subgroup with 49 members have unique selectivity filter very rarely found in other MIP channels. With the discovery of these two new AQGP subfamilies, it appears that fungal AQGPs are the most diverse among all known AQGPs. Many fungal species possessing δ-AQGPs also are known to act as plant pathogens. The newly identified δ-group MIP channels can be exploited to contain some of the serious infections affecting agriculturally and economically important crops.
With very few functional studies available in fungal MIP channels, the present study has given a picture of the diverse fungal organisms evolved with two new MIP subfamilies with unique selectivity filter residues. The possible solutes that are transported through these narrow regions within the channel interior, the role of loop D in SIP-like channels and the importance of these channels in the fungal life cycles are some of the immediate questions that have to be addressed by researchers in this field.
Identification of fungal MIP sequences
The non-redundant database of NCBI (http://www.ncbi.nlm.nih.gov) and the UniProtKB (http://www.uniprot.org) database  were searched using known AQP and AQGP sequences as query. The UniProt accession IDs of these sequences are P0AER0, P47865, P60844, Q6J8I9, Q9C4Z5, Q41372, Q8WPZ6, P55064, P55087 and F2QVG4. These sequences belong to different species and organism groups and they come from bacteria (Escherichia coli), mammals (Bos taurus, Ovis aries and Homo sapiens), archaea (Methanothermobacter marburgensis), plant (Spinacia oleracea), protozoan parasite (Plasmodium falciparum) and yeast (Komagataella pastoris). The three-dimensional structures of all these channels have been determined. In addition to these sequences, we have also considered the sequences belonging to the XIP subfamily from fungal species as query sequences previously identified in our laboratory . We used PSI-BLAST  on NCBI non-redundant database to identify fungal MIP sequences. For each sequence, three iterations were performed with threshold e-value 0.001. When going from previous iteration to the next iteration, we included only those hits whose expect value was 0.001 or less to generate the position-specific scoring matrix that can be used in the next iteration. This is to make sure that the PSSM is not corrupted and at the same time diverse MIP sequences can be obtained as hits. BLASTp was used with UniProt database with the same query sequences. Additionally, we also used tBLASTn against NCBI GenBank database to retrieve additional fungal MIP sequences. The newly identified sequences from PSI-BLAST and tBLASTn searches were again used as query sequences and second round of search was carried out using the same tools. This exercise was repeated until no new MIP sequences were found.
We then used the program CD-HIT , on the set all fungal MIP sequences thus obtained to remove any redundancy. It is possible that some of the MIP sequences in this non-redundant set may be partial or may not have all the features associated with a MIP channel protein. In order to find out the partial sequences or sequences that don’t have all MIP features, we adopted the following strategy. We used a multiple sequence alignment program PRALINE (http://www.ibi.vu.nl/programs/pralinewww/) to align all the fungal MIP sequences . There are advantages using PRALINE. It uses profile-based approach to align sequences that have low sequence identity. It also guides the sequence alignment using secondary structure information by employing different scoring matrices for helix, strand or coil. For each MIP sequence, we used TMHMM  and PSIPRED  which are available as part of the PRALINE toolkit to predict transmembrane regions and secondary structures respectively. The multiple sequence alignment thus created was examined to find out the following features specific to MIP family sequences. These are (a) presence of two NPA or NPA-like motifs, (b) presence of six transmembrane segments and two functionally important loops possessing the features characteristically present in MIP channels and (c) group-based conservation of small and weakly polar residues in most of the 17 positions that typically occur in the helix-helix interface of MIP hour-glass helical fold . Only those sequences which satisfy all the three conditions were considered further for analysis.
We identified 326 fungal MIP sequences from the database search described above. We also included additional 69 fungal MIP sequences from the recently published work of Xu et al.  which fulfilled all our criteria. In total, 395 fungal MIP sequences were considered for phylogenetic analysis.
To understand the fungal MIP diversity, the evolution of MIP family and to enable comparison with other groups of species, we performed phylogenetic analysis on all fungal MIP sequences using MEGA 5.0 suite of software . Two different clustering algorithms, namely neighbor-joining and maximum parsimony methods, were used to derive a phylogenetic tree. Reliability of individual branches of the tree was estimated by performing bootstrapping with 1000 replicates. We applied 50% majority rule so that branches with less than 50% confidence level were collapsed to get the final tree topology. The fungal MIP subgroups obtained from this tree was considered for further analysis. The groupings obtained from this tree were validated by comparing the trees obtained from the two different clustering methods. We also used different multiple sequence alignments as input for MEGA 5.0 obtained from various MSA tools like Clustal-W , Clustal-Ω  and MUSCLE  and the resulting trees generated from the newly obtained multiple sequence alignments were compared.
Homology modeling of fungal MIP sequences
Three-dimensional structures of all fungal MIP sequences were modeled using the protocol developed in our laboratory earlier ,. Briefly, the software MODELLER version 9.10 , was used to build the structure and the experimentally determined high-resolution MIP structures from mammalian AQP1 (PDB ID:1J4N) , bacterial GlpF (PDB ID: 1FX8)  and archael AQPM (PDB ID: 2F2B)  were used as templates. The target-template sequence alignment was manually examined to find out whether there are gaps in the transmembrane helical regions or in the functionally important loop regions which possess the highly conserved NPA motifs. We also looked at the high conservation of at least one residue in each transmembrane helical segment (see above). Among the ten models built for each fungal MIP sequence, the one with the optimum MODELLER objective function was selected. Loops and side-chain conformations of non-conserved residues were further refined using the MODELLER’s loop optimization protocol and the program SCWRL3  respectively. The resultant model was minimized using GROMACS ver. 4.5  and the quality of the model was examined using PROCHECK .
No human or animal experiments were carried out in this study.
Major intrinsic proteins
Small and basic intrinsic proteins
Tonoplast intrinsic proteins
We gratefully acknowledge the support received from the Centre of Excellence for Chemical Biology at IIT-Kanpur funded by MHRD, Government of India. RS is holding USV Chair Professorship at IIT-Kanpur. RKV thanks CSIR, Government of India for a Senior Research Fellowship. NDP thanks IIT-Kanpur for a fellowship. We acknowledge Dr. Anjali Bansal Gupta and Mr. Manu Vajpai for their help. We thank our lab members for useful discussions.
- Gomes D, Agasse A, Thiebaud P, Delrot S, Geros H, Chaumont F: Aquaporins are multifunctional water and solute transporters highly divergent in living organisms. Biochim Biophys Acta. 2009, 1788: 1213-1228. 10.1016/j.bbamem.2009.03.009.PubMedView ArticleGoogle Scholar
- Maurel C: Plant aquaporins: Novel functions and regulation properties. FEBS Lett. 2007, 581: 2227-2236. 10.1016/j.febslet.2007.03.021.PubMedView ArticleGoogle Scholar
- Wu B, Beitz E: Aquaporins with selectivity for unconventional permeants. Cell Mol Life Sci. 2007, 64: 2413-2421. 10.1007/s00018-007-7163-2.PubMedView ArticleGoogle Scholar
- Bienert GP, Schussler MD, Jahn TP: Metalloids: essential, beneficial or toxic? Major intrinsic proteins sort it out. Trends Biochem Sci. 2008, 33: 20-26. 10.1016/j.tibs.2007.10.004.PubMedView ArticleGoogle Scholar
- Abascal F, Irisarri I, Zardoya R: Diversity and evolution of membrane intrinsic proteins. Biochim Biophys Acta. 2014, 1840: 1468-1481. 10.1016/j.bbagen.2013.12.001.PubMedView ArticleGoogle Scholar
- Gupta AB, Verma RK, Agarwal V, Vajpai M, Bansal V, Sankararamakrishnan R: MIPModDB: A central resource for the superfamily of major intrinsic proteins. Nucleic Acids Res. 2012, 40: D362-D369. 10.1093/nar/gkr914.PubMedPubMed CentralView ArticleGoogle Scholar
- Borgnia MJ, Agre P: Reconstitution and functional comparison of purified GlpF and AqpZ, the glycerol and water channels from Escherichia coli. Proc Natl Acad Sci U S A. 2001, 98: 2888-2893. 10.1073/pnas.051628098.PubMedPubMed CentralView ArticleGoogle Scholar
- Rojek A, Praetorius J, Frokiaer J, Nielsen S, Fenton RA: A current view of the mammalian aquaglyceroporins. Annu Rev Physiol. 2008, 70: 301-327. 10.1146/annurev.physiol.70.113006.100452.PubMedView ArticleGoogle Scholar
- Ishibashi K, Hara S, Kondo S: Aquaporin water channels in mammals. Clin Exp Nephrol. 2009, 13: 107-117. 10.1007/s10157-008-0118-6.PubMedView ArticleGoogle Scholar
- Johanson U, Karlsson M, Johansson I, Gustavsson S, Sjovall S, Fraysse L, Weig AR, Kjellbom P: The complete set of genes encoding major intrinsic proteins in Arabidopsis provides a framework for a new nomenclature for major intrinsic proteins in plants. Plant Physiol. 2001, 126: 1358-1369. 10.1104/pp.126.4.1358.PubMedPubMed CentralView ArticleGoogle Scholar
- Johanson U, Gustavsson S: A new subfamily of major intrinsic proteins in plants. Mol Biol Evol. 2002, 19: 456-461. 10.1093/oxfordjournals.molbev.a004101.PubMedView ArticleGoogle Scholar
- Gupta AB, Sankararamakrishnan R: Genome-wide analysis of major intrinsic proteins in the tree plant Populus trichocarpa: Characterization of XIP subfamily of aquaporins from evolutionary perspective.BMC Plant Biol 2009, 9:Art. no. 134.,
- Nielsen S, Frokiaer J, Marples D, Kwon TH, Agre P, Knepper MA: Aquaporins in kidney: From molecules to medicine. Physiol Rev. 2002, 82: 205-244.PubMedView ArticleGoogle Scholar
- Papadopoulos MC, Verkman AS: Aquaporin water channels in the nervous system. Nature Rev Neurosci. 2013, 14: 265-277. 10.1038/nrn3468.View ArticleGoogle Scholar
- Verkman AS: Aquaporins in clinical medicine. Annu Rev Med. 2012, 63: 303-316. 10.1146/annurev-med-043010-193843.PubMedPubMed CentralView ArticleGoogle Scholar
- Agre P, King LS, Yasui M, Guggino WB, Otterson OP, Fujiyoshi Y, Engel A, Nielsen S: Aquaporin water channels - from atomic structure to clinical medicine. J Physiol London. 2002, 542: 3-16. 10.1113/jphysiol.2002.020818.PubMedPubMed CentralView ArticleGoogle Scholar
- Lennon VA, Kryzer TJ, Pittock SJ, Verkman AS, Hinson SR: IgG marker of optic-spinal multiple sclerosis binds to the aquaporin-4 water channel. J Exp Med. 2005, 202: 473-477. 10.1084/jem.20050304.PubMedPubMed CentralView ArticleGoogle Scholar
- Verkman AS, Hara-Chikuma M, Papadopoulos MC: Aquaporins - new players in cancer biology. J Mol Med. 2008, 86: 523-529. 10.1007/s00109-008-0303-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Verkman AS: Role of aquaporin water channels in eye function. Exp Eye Res. 2003, 76: 137-143. 10.1016/S0014-4835(02)00303-2.PubMedView ArticleGoogle Scholar
- Binder DK, Nagelhus EA, Ottersen OP: Aquaporin-4 and epilepsy. Glia. 2012, 60: 1203-1214. 10.1002/glia.22317.PubMedView ArticleGoogle Scholar
- Hibuse T, Maeda N, Funahashi T, Yamamoto K, Nagasawa A, Mizunoya W, Kishida K, Inoue K, Kuriyama H, Nakamura T, Fushiki T, Kihara S, Shimomura L: Aquaporin 7 deficiency is associated with development of obesity through activation of adipose glycerol kinase. Proc Natl Acad Sci U S A. 2005, 102: 10993-10998. 10.1073/pnas.0503291102.PubMedPubMed CentralView ArticleGoogle Scholar
- Loonen AJM, Knoers NVAM, van Os CH, Deen PMT: Aquaporin 2 mutations in nephrogenic diabetes insipidus. Seminars Nephrol. 2008, 26: 252-265. 10.1016/j.semnephrol.2008.03.006.View ArticleGoogle Scholar
- Maurel C, Verdoucq L, Luu DT, Santoni V: Plant aquaporins: Membrane channels with multiple integrated functions. Annu Rev Plant Biol. 2008, 59: 595-624. 10.1146/annurev.arplant.59.032607.092734.PubMedView ArticleGoogle Scholar
- Johansson I, Karlsson M, Johanson U, Larsson C, Kjellbom P: The role of aquaporins in cellular and whole plant water balance. Biochim Biophys Acta. 2000, 1465: 324-342. 10.1016/S0005-2736(00)00147-4.PubMedView ArticleGoogle Scholar
- Tyerman SD, Niemietz CM, Bramley H: Plant aquaporins: multifunctional water and solute channels with expanding roles. Plant Cell Environ. 2002, 25: 173-194. 10.1046/j.0016-8025.2001.00791.x.PubMedView ArticleGoogle Scholar
- Sade N, Vinocur BJ, Diber A, Shatil A, Ronen G, Nissan H, Wallach R, Karchi H, Moshelion M: Improving plant stress tolerance and yield production: Is the tonoplat aquaporin SITIP2;2 a key to isohydric to anisohydric conversion?. New Phytol. 2009, 181: 651-661. 10.1111/j.1469-8137.2008.02689.x.PubMedView ArticleGoogle Scholar
- Marjanovic Z, Uehlein N, Kaldenhoff R, Zwiazek JJ, Weiss M, Hampp R, Nehls U: Aquaporins in poplar: What a difference a symbiont makes!. Planta. 2005, 222: 258-268. 10.1007/s00425-005-1539-z.PubMedView ArticleGoogle Scholar
- Uehlein N, Fileschi K, Eckert M, Bienert GP, Bertl A, Kaldenhoff R: Arbuscular mycorrhizal symbiosis and plant aquaporin expression. Phytochemistry. 2007, 68: 122-129. 10.1016/j.phytochem.2006.09.033.PubMedView ArticleGoogle Scholar
- Dowd C, Wilson LW, McFadden H: Gene expression profile changes in cotton root and hypocotyl tissues in response to infection with Fusarium oxysporum f.sp vasinfectum. Mol Plant Microbe Interact. 2004, 17: 654-667.PubMedView ArticleGoogle Scholar
- Fadiel A, Isokpehi RD, Stambouli N, Hamza A, Benammar-Elgaaied A, Scalise TJ: Protozoan parasite aquaporins. Expert Rev Proteomics. 2009, 6: 199-211. 10.1586/epr.09.10.PubMedView ArticleGoogle Scholar
- Beitz E: Aquaporins from pathogenic protozoan parasites: structure, function and potential for chemotherapy. Biol Cell. 2005, 97: 373-383. 10.1042/BC20040095.PubMedView ArticleGoogle Scholar
- Hansen M, Kun JFJ, Schultz JE, Beitz E: A single, bi-functional aquaglyceroporin in blood-stage Plasmodium falciparum malaria parasites. J Biol Chem. 2002, 277: 4874-4882. 10.1074/jbc.M110683200.PubMedView ArticleGoogle Scholar
- Montalvetti A, Rohloff P, Docampo R: A functional aquaporin co-localizes with the vacuolar proton pyrophosphatase to acidocalcisomes and the contractile vacuole complex of Trypanosoma cruzi. J Biol Chem. 2004, 279: 38673-38682. 10.1074/jbc.M406304200.PubMedView ArticleGoogle Scholar
- Baker N, Glover L, Munday JC, Andres DA, Barrett MP, de Koning HP, Horn D: Aquaglyceroporin 2 controls susceptibility to melarsoprol and pentamidine in African tyrpanosomes. Proc Natl Acad Sci U S A. 2012, 109: 10996-11001. 10.1073/pnas.1202885109.PubMedPubMed CentralView ArticleGoogle Scholar
- Kun JF, de Carvalho EG: Novel therapeutic targets in Plasmodium falciparum: aquaglyceroporins. Expert Opin Ther Targets. 2009, 13: 385-394. 10.1517/14728220902817839.PubMedView ArticleGoogle Scholar
- Bonfante P, Genre A: Mechanisms underlying beneficial plant-fungus interactions in mycorrhizal symbiosis.Nature Commn 2010, 1:Art. No. 48.,
- Casieri L, Lahmidi NA, Doidy J, Veneault-Fourrey C, Migeon A, Bonneau L, Courty PE, Garcia K, Charbonnier M, Delteil A, Brun A, Zimmermann S, Plassard C, Wipf D: Biotrophic transportome in mutualistic plant-fungal interactions. Mycorrhiza. 2013, 23: 597-625. 10.1007/s00572-013-0496-9.PubMedView ArticleGoogle Scholar
- Parniske M: Arbuscular mycorrhiza: the mother of plant root endosymbioses. Nature Rev Microbiol. 2008, 6: 763-775. 10.1038/nrmicro1987.View ArticleGoogle Scholar
- Li T, Hu YJ, Hao ZP, Li H, Wang YS, Chen BD: First cloning and characterization of two functional aquaporin genes from an arbuscular mycorrhizal fungus Glomus intraradices. New Phytol. 2013, 197: 617-630. 10.1111/nph.12011.PubMedView ArticleGoogle Scholar
- Giovannetti M, Balestrini R, Volpe V, Guether M, Straub D, Costa A, Ludewig U, Bonfante P: Two putative aquaporin genes are differently expressed during arbuscular mycorrhizal symbiosis inLotus japonicus.BMC Plant Biol 2012, 12:Art. No: 186.,
- Navarro-Rodenas A, Ruiz-Lozano JM, Kaldenhoff R, Morte A: The aquaporin TcAQP1 of the desert truffle Terfezia claveryi is a membrane pore for water and CO2 transport. Mol Plant-Microbe Interact. 2012, 25: 259-266. 10.1094/MPMI-07-11-0190.PubMedView ArticleGoogle Scholar
- Lee SH, Calvo-Polanco M, Chung GC, Zwiazek JJ: Role of aquaporins in root water transport of ectomycorrhizal jack pine (Pinus banksiana) seedlings exposed to NaCl and fluoride. Plant Cell Environ. 2010, 33: 769-780.PubMedGoogle Scholar
- Aroca R, Porcel R, Ruiz-Lozano JM: How does arbuscular mycorrhizal symbiosis regulate root hydraulic properties and plasma membrane aquaporins in Phaseolus vulgaris under drought, cold or salinity stresses?. New Phytol. 2007, 173: 808-816. 10.1111/j.1469-8137.2006.01961.x.PubMedView ArticleGoogle Scholar
- Dietz S, von Bulow J, Baker N, Nehls U: The aquaporin gene family of the ectomycorrhizal fungus Laccaria bicolor: lessons for symbiotic functions. New Phytol. 2011, 190: 927-940. 10.1111/j.1469-8137.2011.03651.x.PubMedView ArticleGoogle Scholar
- Aroca R, Bago A, Sutka M, Paz JA, Cano C, Amodeo G, Ruiz-Lozano JM: Expression analysis of the first arbuscular mycorrhizal fungi aquaporin described reveals concerted gene expression between salt-stressed and nonstressed mycelium. Mol Plant Microbe Interact. 2009, 22: 1169-1178.PubMedView ArticleGoogle Scholar
- Castle NA: Aquaporins as targets for drug discovery. Drug Disc Today. 2005, 10: 485-493. 10.1016/S1359-6446(05)03390-8.View ArticleGoogle Scholar
- Kontoyiannis DP, Lewis RE: Antifungal drug resistance of pathogenic fungi. Lancet. 2002, 359: 1135-1144. 10.1016/S0140-6736(02)08162-X.PubMedView ArticleGoogle Scholar
- Sanglard D: Resistance of human fungal pathogens to antifungal drugs. Curr Opin Microbiol. 2002, 5: 379-385. 10.1016/S1369-5274(02)00344-2.PubMedView ArticleGoogle Scholar
- Ghosh K, Cappiello CD, McBride SM, Occi JL, Cali A, Takvorian PM, McDonald TV, Weiss LM: Functional characterization of a putative aquaporin from Encephalitozoon cuniculi, a microsporidia pathogenic to humans. Int J Parasitol. 2006, 36: 57-62. 10.1016/j.ijpara.2005.08.013.PubMedPubMed CentralView ArticleGoogle Scholar
- Pettersson N, Filipsson C, Becit E, Brive L, Hohmann S: Aquaporins in yeasts and filamentous fungi. Biol Cell. 2005, 97: 487-500. 10.1042/BC20040144.PubMedView ArticleGoogle Scholar
- Xu H, Cooke JEK, Zwiazek JJ: Phylogenetic analysis of fungal aquaporins provide insight into their possible role in water transport of mycorrhizal associations. Botany-Botanique. 2013, 91: 495-504. 10.1139/cjb-2013-0041.View ArticleGoogle Scholar
- Consortium TU: Activities at the Universal Protein Resource (UniProt). Nucleic Acids Res. 2014, 42: D191-D198. 10.1093/nar/gkt1140.View ArticleGoogle Scholar
- Kirk PM, Cannon PF, Minter DW, Stalpers JA: Dictionary of the fungi. 2008, CAB International, Wallingford, U. K.Google Scholar
- Jones MDM, Forn I, Gadelha C, Egan MJ, Bass D, Massana R, Richards TA: Discovery of novel intermediate forms redefines the fungal tree of life. Nature. 2011, 474: 200-203. 10.1038/nature09984.PubMedView ArticleGoogle Scholar
- Azad AK, Yoshikawa N, Ishikawa T, Sawa Y, Shibata H: Substitution of a single amino acid residue in the aromatic/arginine selectivity filter alters the transport profiles of tonoplast aquaporin homologs. Biochim Biophys Acta. 1818, 2012: 1-11.Google Scholar
- Mitani-Ueno N, Yamaji N, Zhao FJ, Ma JF: The aromatic/arginine selectivity filter of NIP aquaporins plays a critical role in substrate selectivity for silicon, boron, and arsenic. J Exp Botany. 2011, 62: 4391-4398. 10.1093/jxb/err158.View ArticleGoogle Scholar
- Beitz E, Wu B, Holm LM, Schultz JE, Zeuthen T: Point mutations in the aromatic/arginine region in aquaporin 1 allow passage of urea, glycerol, ammonia, and protons. Proc Natl Acad Sci U S A. 2006, 103: 269-274. 10.1073/pnas.0507225103.PubMedPubMed CentralView ArticleGoogle Scholar
- Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: A sequence logo generator. Genome Res. 2004, 14: 1186-1190. 10.1101/gr.849004.View ArticleGoogle Scholar
- Bansal A, Sankararamakrishnan R: Homology modeling of major intrinsic proteins in rice, maize and Arabidopsis: comparative analysis of transmembrane helix association and aromatic/arginine selectivity filters.BMC Struct Biol 2007, 7:Art. No. 27.,
- Verma RK, Jain A, Sankararamakrishnan R: Distinguishing features of aquaglyceroporin in Plasmodium falciparum: Comparative molecular dynamics simulations of three aquaporins.Biophys J 2012, 102:452A.,
- Li H, Chen HN, Steinbronn C, Wu BH, Beitz E, Zeuthen T, Voth GA: Enhancement of proton conductance by mutations of the selectivity filter of aquaporin-1. J Mol Biol. 2011, 407: 607-620. 10.1016/j.jmb.2011.01.036.PubMedView ArticleGoogle Scholar
- Savage DF, O’Connell JD, Miercke LJW, Finer-Moore J, Stroud RM: Structural context shapes the aquaporin selectivity filter. Proc Natl Acad Sci U S A. 2010, 107: 17164-17169. 10.1073/pnas.1009864107.PubMedPubMed CentralView ArticleGoogle Scholar
- Bienert GP, Bienert MD, Jahn TP, Boutry M, Chaumont F: Solanaceae XIPs are plasma membrane aquaporins that facilitate the transport of many uncharged substrates. Plant J. 2011, 66: 306-317. 10.1111/j.1365-313X.2011.04496.x.PubMedView ArticleGoogle Scholar
- Hub JS, de Groot BL: Does CO2 permeate through aquaporin-1?. Biophys J. 2006, 91: 842-848. 10.1529/biophysj.106.081406.PubMedPubMed CentralView ArticleGoogle Scholar
- Ishikawa F, Suga S, Uemura T, Sato MH, Maeshima M: Novel type aquaporin SIPs are mainly localized to the ER membrane and show cell-specific expression in Arabidopsis thaliana. FEBS Lett. 2005, 579: 5814-5820. 10.1016/j.febslet.2005.09.076.PubMedView ArticleGoogle Scholar
- Graf FE, Ludin P, Wenzler T, Kaiser M, Brun R, Pyana PP, Buscher P, de Koning HP, Horn D, Maser P: Aquaporin 2 mutations in Trypanosoma brucei gambiense field isolates correlate with decreased susceptibility to pentamidine and melarsoprol.PLoS Neg Trop Dis 2013, 7:Art. No. e2475.,
- Ohta E, Itoh T, Nemoto T, Kumagai J, Ko SBH, Ishibashi K, Ohno M, Uchida K, Ohta A, Sohara E, Uchida S, Sasaki S, Rai T: Pancreas-specific aquaporin 12 null mice showed increased susceptibility to caerulein-induced acute pancreatitis. Am J Physiol Cell Physiol. 2009, 297: C1368-C1378.PubMedView ArticleGoogle Scholar
- Itoh T, Rai T, Kuwahara M, Ko SBH, Uchida S, Sasaki S, Ishibashi K: Identification of a novel aquaporin, AQP12, expressed in pancreatic acinar cells. Biochem Biophys Res Commn. 2005, 330: 832-838. 10.1016/j.bbrc.2005.03.046.View ArticleGoogle Scholar
- Nyblom M, Frick A, Wang Y, Ekvall M, Hallgren K, Hedfalk K, Neutze R, Tajkhorshid E, Tornroth-Horsefield S: Structural and functional analysis of SoPIP2;1 mutants adds insight into plant aquaporin gating. J Mol Biol. 2009, 387: 653-668. 10.1016/j.jmb.2009.01.065.PubMedView ArticleGoogle Scholar
- Yukutake Y, Tsuji S, Hirano Y, Adachi T, Takahashi T, Fujihara K, Agre P, Yasui M, Suematsu M: Mercury chloride decreases the water permeability of aquaporin-4-reconstituted proteoliposomes. Biol Cell. 2008, 100: 355-363. 10.1042/BC20070132.PubMedView ArticleGoogle Scholar
- Yool AJ, Morelle J, Cnops Y, Verbavatz JM, Campbell EM, Beckett EAH, Booker GW, Flynn G, Devuyst O: AqF206 is a pharmacologic agonist of the water channel aquaporin-1. J Am Soc Nephrol. 2013, 24: 1045-1052. 10.1681/ASN.2012080869.PubMedPubMed CentralView ArticleGoogle Scholar
- Yool AJ, Brown EA, Flynn GA: Roles for novel pharmacological blockers of aquaporins in the treatment of brain oedema and cancer. Clin Exp Pharmacol Physiol. 2010, 37: 403-409. 10.1111/j.1440-1681.2009.05244.x.PubMedView ArticleGoogle Scholar
- Frick A, Jarva M, Tornroth-Horsefield S: Structural basis for pH gating of plant aquaporins. FEBS Lett. 2013, 587: 989-993. 10.1016/j.febslet.2013.02.038.PubMedView ArticleGoogle Scholar
- Yu J, Yool AJ, Schulten K, Tajkhorshid E: Mechanism of gating and ion conductivity of a possible tetrameric pore in aquaporin-1. Structure. 2006, 14: 1411-1423. 10.1016/j.str.2006.07.006.PubMedView ArticleGoogle Scholar
- Yukutake Y, Yasui M: Regulation of water permeability through aquaporin-4. Neuroscience. 2010, 168: 885-891. 10.1016/j.neuroscience.2009.10.029.PubMedView ArticleGoogle Scholar
- Yukutake Y, Hirano Y, Suematsu M, Yasui M: Rapid and reversible inhibition of aquaporin-4 by zinc. Biochemistry. 2009, 48: 12059-12061. 10.1021/bi901762y.PubMedView ArticleGoogle Scholar
- Tornroth-Horsefield S, Wang Y, Hedfalk K, Johanson U, Karlsson M, Tajkhorshid E, Neutze R, Kjellbom P: Structural mechanism of plant aquaporin gating. Nature. 2006, 439: 688-694. 10.1038/nature04316.PubMedView ArticleGoogle Scholar
- Didier ES: Microsporidiosis: An emerging and opportunistic infection in humans and animals. Acta Trop. 2005, 94: 61-76. 10.1016/j.actatropica.2005.01.010.PubMedView ArticleGoogle Scholar
- Oide S, Moeder W, Krasnoff S, Gibson D, Haas H, Yoshioka K, Turgeon BG: NPS6, encoding a nonribosomal pepide synthetase involved in siderophore-mediated iron metabolism, is a conserved virulence determinant of plant pathogenic ascomycetes. Plant Cell. 2006, 18: 2836-2853. 10.1105/tpc.106.045633.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Li WZ, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22: 1658-1659. 10.1093/bioinformatics/btl158.PubMedView ArticleGoogle Scholar
- Huang Y, Niu BF, Gao Y, Fu LM, Li WZ: CD-HIT suite: a web server for clustering and comparing biological sequences. Bioinformatics. 2010, 26: 680-682. 10.1093/bioinformatics/btq003.PubMedPubMed CentralView ArticleGoogle Scholar
- Simossis VA, Heringa J: PRALINE: a multiple sequence alignment toolbox that integrates homology-extended and secondary structure information. Nucleic Acids Res. 2005, 33: W289-W294. 10.1093/nar/gki390.PubMedPubMed CentralView ArticleGoogle Scholar
- Krogh A, Larsson B, von Heijne G, Sonnhammer ELL: Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J Mol Biol. 2001, 305: 567-580. 10.1006/jmbi.2000.4315.PubMedView ArticleGoogle Scholar
- McGuffin LJ, Bryson K, Jones DT: The PSIPRED protein structure prediction server. Bioinformatics. 2000, 16: 404-405. 10.1093/bioinformatics/16.4.404.PubMedView ArticleGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nel M, Kumar S: MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMedPubMed CentralView ArticleGoogle Scholar
- Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31: 3497-3500. 10.1093/nar/gkg500.PubMedPubMed CentralView ArticleGoogle Scholar
- Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Soding J, Thompson JD, Higgins DG: Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega.Mol Systems Biol 2011, 7:Art. No. 539.,
- Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity.BMC Bioinformatics 2004, 5:Art. No. 113.,
- Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993, 234: 779-815. 10.1006/jmbi.1993.1626.PubMedView ArticleGoogle Scholar
- Eswar N, Eramian D, Webb B, Shen M-Y, Sali A: Protein structure modeling with MODELLER. Methods Mol Biol. 2008, 426: 145-159. 10.1007/978-1-60327-058-8_8.PubMedView ArticleGoogle Scholar
- Sui H, Han BG, Lee JK, Walian P, Jap BK: Structural basis of water-specific transport through the AQP1 water channel. Nature. 2001, 414: 872-878. 10.1038/414872a.PubMedView ArticleGoogle Scholar
- Fu D, Libson A, Miercke LJ, Weitzman C, Nollert P, Krucinski J, Stroud RM: Structure of a glycerol-conducting channel and the basis for its selectivity. Science. 2000, 290: 481-486. 10.1126/science.290.5491.481.PubMedView ArticleGoogle Scholar
- Lee JK, Kozono D, Remis J, Kitagawa Y, Agre P, Stroud RM: Structural basis for conductance by the archaeal aquaporin AqpM at 1.68 A. Proc Natl Acad Sci U S A. 2005, 102: 18932-18937. 10.1073/pnas.0509469102.PubMedPubMed CentralView ArticleGoogle Scholar
- Canutescu AA, Shelenkov AA, Dunbrack RL: A graph-theory algorithm for rapid protein side-chain prediction. Proc Natl Acad Sci U S A. 2003, 12: 2001-2014.Google Scholar
- Hess B, Kutzner C, van der Spoel D, Lindahl E: GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation. J Chem Theoy Comput. 2008, 4: 435-447.View ArticleGoogle Scholar
- Laskowski RA, MacArthur MW, Moss DS, Thornton JM: PROCHECK: A program to check the stereochemical quality of protein structures. J Appl Cryst. 1993, 26: 283-291. 10.1107/S0021889892009944.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.