Molecular dynamics study of the archaeal aquaporin AqpM
© Araya-Secchi et al; licensee BioMed Central Ltd. 2011
Published: 22 December 2011
Skip to main content
© Araya-Secchi et al; licensee BioMed Central Ltd. 2011
Published: 22 December 2011
Aquaporins are a large family of transmembrane channel proteins that are present throughout all domains of life and are implicated in human disorders. These channels, allow the passive but selective movement of water and other small neutral solutes across cell membranes. Aquaporins have been classified into two sub-families: i) strict aquaporins that only allow the passage of water and ii) the less selective aquaglyceroporins that transport water and other neutral solutes, such as glycerol, CO2 or urea. Recently, the identification and characterization of a number of archaeal and bacterial aquaporins suggested the existence of a third sub-family; one that is neither a strict aquaporin nor an aquaglyceroporin. The function and phylogeny of this third family is still a matter of debate.
Twenty nanosecond molecular dynamics (MD) simulation of a fully hydrated tetramer of AqpM embedded in a lipid bilayer permitted predictions to be made of key biophysical parameters including: single channel osmotic permeability constant (p f ), single channel diffusive permeability constant (p d ), channel radius, potential water occupancy of the channel and water orientation inside the pore. These properties were compared with those of well characterized representatives of the two main aquaporin sub-families. Results show that changes in the amino acid composition of the aromatic/arginine region affect the size and polarity of the selectivity filter (SF) and could help explain the difference in water permeability between aquaporins. In addition, MD simulation results suggest that AqpM combines characteristics of strict aquaporins, such as the narrow SF and channel radius, with those of aquaglyceroporins, such as a more hydrophobic and less polar SF.
MD simulations of AqpM extend previous evidence that this archaeal aquaporin exhibits hybrid features intermediate between the two known aquaporin sub-families, supporting the idea that it may constitute a member of a novel class of aquaporins.
Aquaporins are a large family of transmembrane channel proteins that are present throughout all domains of life and their malfunction has been implicated in several human disorders . Aquaporin channels allow the passive but selective movement of water and other small neutral solutes such as glycerol, CO2, or urea across cell membranes [1–4].
Structurally, Aquaporins present a homotetrameric organization in which each monomer forms an individual functional pore. The canonical fold of the aquaporin monomer is characterized by a right-handed helical bundle of six transmembrane α-helices (TM1 – TM6) connected by 5 loop regions (loops A to E), in which both amino and carboxyl termini face the cytoplasmic side of the membrane. Loops B and E are formed by a half-membrane spanning helical section (HB and HE respectively) and a non-helical section that contains the highly conserved Asn-Pro-Ala (NPA) motif, considered a signature of this protein family [5–10].
The aquaporin channel posses two major constriction zones, the NPA region, located at the centre of the pore and the selectivity filter (SF) containing the aromatic/arginine region, (ar/R), located ~8 Å above the NPA region extending towards the extracellular side of the channel. The ar/R region is formed by a residue from helix TM2 (H2 position), a residue from helix TM5 (H5 position) and two residues from loop E (LE1, LE2 position respectively) [11, 12]. This region plays a major role in solute selectivity and permeability .
Different permeability to water and other solutes and variations in the amino acid composition of the ar/R region has led to the classification of aquaporins into two main functional sub-families [14, 15]. The first are strict aquaporins, that only allow the passage of water molecules, in which the ar/R region contains a Phe residue from TM2 (H2 position), a His residue from TM5 (H5), an Arg residue from the loop E (LE2) and a fourth residue from loop E (LE1) that provides a backbone carbonyl oxygen, usually a Cys, Thr or Ala. In these aquaporins, the His and Arg residues of the SF are believed to provide donor hydrogen bonds for water molecules [14, 15]. Examples of strict aquaporins are: AqpZ from Escherichia coli [16, 17], and the human Aqp0 , Aqp1 and Aqp4. The second are aquaglyceroporins, less selective aquaporins that can transport water and other solutes, such as glycerol, urea and other uncharged small molecules [19, 20]. In these aquaporins, the Phe residue at position H2 is replaced by Trp residue, the His residue at position H5 is replaced by Gly residue and the small residue at position LE1 is replaced by a Phe, giving rise to a wider SF and a so called “hydrophobic corner” formed by Trp (H2) and Phe (LE1) residues [3, 20–23]. Similar to strict aquaporins, position LE2 is occupied by an Arg residue that is highly conserved in both aquaporin sub-families. Examples of aquaglyceroporins are: GlpF [19–21, 23] from E. coli, and human aquaporins Aqp7, Aqp3, Aqp9 and Aqp10 .
The recent identification of archaeal aquaporins that exhibit a different amino acid composition of their SF, and thus cannot be readily classified in either of the two aforementioned aquaporin sub-families, has led to the suggestion that a third functional sub-family of aquaporins may exist [24, 25]. It has been suggested that this family could contain special adaptations for the conduction of solutes specific for the life-style of the organisms or, alternatively, that they could be a primitive non-specialized aquaporin, an ancestor of the more specialized aquaporins [2, 24, 26]. The SF of these archaeal aquaporins exhibits the highly conserved Arg residue in position LE2, the Phe residue at position H2, and a Ser residue in position LE1 that provides its main chain carbonyl oxygen. However a key difference in the composition of the SF of these aquaporins lies in the presence of a medium size hydrophobic residue such as Ile, Leu or Val instead of the His residue highly conserved in strict aquaporins or the Gly residue present in aquaglyceroporins in position H5 [24, 25].
To date, the only representative of this type of aquaporins that has been experimentally studied is AqpM, found in the archaeon Methanothermobacter marburgensis [24, 25]. The osmotic water permeability constant, P f , obtained for AqpM, showed that its osmotic permeability was lower than the P f of the strict aquaporins AqpZ from E.coli, and slightly higher than the P f of the aquaglyceroporin from E. coli GlpF [24, 25]. (Pf (AqpM) = 57 μms-1 ; Pf (AqpZ) = 330 μms-1 ; Pf (GlpF) = 49 μms-1 ). Moreover, the transient glycerol permeability measured for AqpM resulted in values considerably lower than the glycerol permeability for GlpF . These previous results have led to the suggestion that AqpM could be a water transporter and also with the ability to transport other neutral solutes required by the archaeon such as CO2, the only carbon source available in the natural environment of M. marburgensis, or even H2S, the terminal electron acceptor in its energy production pathway .
The results of this study extend the information regarding the permeability and selectivity mechanisms of AqpM and provide further insight into the water permeability of aquaporins in general and the particular differences between the sub-families of aquaporins and the implications of these differences for their function.
Summary of results
Mean ± SD
p f / p d
19.9 ± 0.9
Comparison with other simulation results
p f /p d
Hashido et al. 2005
Jensen et al. 2006
Hashido et al.2007 
16 ± 3.0
Hashido et al. 2005
Zhu et al. 2002 
Jensen et al. 2006 
Hashido et al.2007 
Hashido et al.2007 
Hashido et al. 2005
Zhu et al. 2004 
Hashido et al.2007 
Garate, J.A. et al 
Kozono et al. 2003 
Lee et al. 2005 
Planar lipid bilayer
Pohl et al. 2001
Borgnia et al. 1999 
Borgnia et al. 2001 
Planar lipid bilayer
Pohl et al 2005 
Borgnia et al. 2001 
Zeidel et al. 1994
Walz et al. 1994 
Yang et al. 1997
Zeidel et al. 1992 
Tanimura et al. 2009
Yang et al. 1997 
3.5 - 9
Jung et al. 1994 
Tanimura et al. 2009
An examination of the published MD simulation results obtained for representative aquaporins from both sub-families, including AqpZ [28, 31, 32] and GlpF [28, 31–33] from E.coli, and mammalian Aqp1 [29, 31, 32] and Aqp4 [31, 34], revealed significant discrepancies on the values obtained for the osmotic and diffusive water permeability and the p f /p d ratio for the same aquaporin between different studies (see Table 2). These discrepancies may arise from the use of different force-fields, water models, simulation protocols, simulation times, and insufficient sampling due to the limited time scale of all-atom MD simulations . This issue together with the fact that only one representative from the aquaglyceroporins sub-family, GlpF from E. coli, has been subjected to MD-simulation studies, whereas four representatives from the strict aquaporins have been widely studied (Aqp0, Aqp1, Aqp4, AqpZ) exacerbates the problem of aquaporin sub-family comparisons in terms of water permeability and thus difficult the interpretation of the data obtained for newly described aquaporins such as AqpM in terms of water permeability. In an attempt to overcome these difficulties and to establish a clear point of reference, the present study of AqpM was performed using simulation protocols and analyses methods similar to those previously used in the study of AqpZ and GlpF [28, 35], including: i) the use of the same force-field (CHARMM22 for protein and CHARMM27 for lipids), and water model (TIP3P), ii) similar simulation protocol (NPT ensemble, Periodic Boundary Condition, no constraints applied in the production run), iii) similar simulation time (20ns) and temperature (310K) and iv) the same MD engine (NAMD) (see methods). The comparison results in the following trend: p f (AqpZ) < p f (AqpM) < p f (GlpF) (see table 2) indicating that the single channel osmotic water permeability of AqpM presents a value of p f that is intermediate between those obtained for AqpZ and GlpF. As a whole, this evidence suggests that AqpM exhibits water osmotic permeability that is intermediate between strict aquaporins and aquaglyceroporins and thus could belong to a third sub-family of aquaporins.
As noted for p f , different MD simulation studies of the same aquaporin show significant discrepancies between the values obtained for p d (see table 2). However, despite these discrepancies a trend is apparent, in which GlpF appears with the highest p d values with respect to all other aquaporins including AqpM. The observation that the higher transport rate appears to be in an opposite trend with respect to experimental results can be explained by the wider pore of GlpF that could allow more water-water interchanges and thus more permeation events per unit time, increasing the value of p d (see below). As mentioned above for p f , this fact can be mainly attributed to the nature of MD simulations, and specifically to the water model used, the van der Waals parameters of the forcefield, and non-physical breathing movements of the channel. Following the argument used for the comparison of p f values, we obtained for p d the same trend observed for p f (p d (AqpZ) < p d (AqpM) < p d (GlpF). This result indicates that the single channel diffusive water permeability of AqpM presents values that are intermediate between those obtained for AqpZ and GlpF, using equivalent simulation methods . This observation, together with the results observed for p f , provides the evidence to support the notion that AqpM could be a representative of a third aquaporin sub-family that exhibits properties that are intermediate between those of strict aquaporins and aquaglyceroporins.
According to the continuous time random walk (CTRW) model developed for single-file water transport , the p f / p d ratio is related to the number of steps that a water molecule must perform in order to traverse the channel and follows the relation: [37, 38] were is the average number of water molecules inside the channel. Thus the p f / p d ratio can be used as a measure of the “single-fileness” of water movement inside the channel, and for a channel in which the water molecules are perfectly aligned in single file [28, 37, 38]. Values of indicates that the single-fileness of the water molecules inside the channel is interrupted by occasional water-water interchanges that increase p d , i.e., more water molecules traverse the channel per unit time [28, 37, 38]. Water-water interchanges do not affect p f because the latter is related to the displacement of individual water molecules inside the channel and is not influenced by the occurrence of water-water interchanges. From the 20 ns simulation of the AqpM tetramer, we obtained 3.0 < p f / p d < 9.1 with an average value over the four monomers of p f / p d ≈ 6.1 and an average channel occupancy equal to (see below). These results suggest that water molecules move in single file in the AqpM channel with occasional water-water interchange events.
For AqpZ , Aqp1 and Aqp4, values of have been reported, whereas for the aquaglyceroporin GlpF [28, 31]. This indicates that strict-aquaporins, with their narrower pore posses a more idealized single-file water transport mechanism. Interestingly, this relation holds despite the aforementioned discrepancies between the values for p d and p f observed between different MD-simulation studies. Our results shows that, in terms of the p f / p d ratio, AqpM presents a similar level of single-fileness with respect to strict aquaporins (see table 2) but with an average value of suggesting that AqpM allows more water-water interchanges than strict aquaporins. This could be explained by the slightly wider ar/R region of AqpM due to the replacement of the His residue in position H5 conserved in strict aquaporins by the smaller and apolar Ile residue.
A bi-polar water orientation is clearly observed inside the AqpM constriction region in the course of our MD simulation (Fig.5b and Fig.6) and the hydrogen donor/hydrogen acceptor pattern appears to be formed by: i) the main chain carbonyl groups exposed to the channel lumen at both sides of the NPA region provided by residues G195(N182/G199), S196(T183/F200), S197(S184/A201) on the extra-cellular side, and residues H80(H61/H66), C79(G60/A65) on the intra-cellular side, that could act as hydrogen bond acceptors forming an hydrogen-bond ladder within the pore and ii) the Asn residues at the NPA region N199(N186/N203) and N82(N63/N68) that act as hydrogen bond donors. An increase in the values of P 1 (z) and negative values of P 2 (z) are observed in the SF of AqpM, indicating that water molecules in this region orient their dipoles perpendicular to the channel axis. This could reflect the interaction of water molecules with the polar hydrogen atoms Arg202(R189/R206):Hη and Arg202(R189/R206):Hε that causes a change in the dipole orientation of water molecules in this region. This interaction appears to be enhanced due to the narrowing of the channel in the SF of AqpM. Thus the unique amino acid composition of the SF of AqpM appears to influence the water orientation in this zone of the channel. Calculation of the same order parameters for other MD-simulated aquaporins indicate that, in the case of strict aquaporins like AqpZ and Aqp1, water molecules assume a more perpendicular orientation of their dipoles with respect to the channel axis around the SF due to the narrowing of the pore and the strong interaction between water molecules and the Arg (LE2 position) and His (H5 position) residues that comprise the SF of strict aquaporins [28, 31, 32]. On the other hand, the wider and hydrophobic SF of the aquaglyceroporin GlpF allows more water-water interactions and thus favours a water dipole orientation parallel to the channel axis [31, 32, 35].
As can be observed in Figure 7(a-d), there are significant differences between the compared aquaporins, regarding the residues that comprise their respective SFs and their spatial arrangement (residues highlighted in blue in Fig. 7d). Of particular interest, is the residue that occupies position H5 (TM5). In AqpZ, b-Aqp1 and h-Aqp4 this position is occupied by a His residue conserved among strict aquaporins (Fig.7d). In GlpF and other aquaglyceroporins position H5 is occupied by a Gly, a feature that allows the formation of the so-called “hydrophobic corner” by Trp48 (replaced by a Phe residue in AqpM and in strict aquaporins) in position H2 and Phe200 (a Ser residue in AqpM, a Tyr residue in AqpZ, a Cys residue in b-Aqp1 and an ala residue in h-Aqp4) in position LE1 that allows the accommodation of a glycerol molecule [19, 20, 49]. In AqpM, position H5 is occupied by Ile (I187), whereas in the other archaeal aquaporins it is occupied either by Ile, Val or Leu [24, 26] generating an SF that is narrower than the SF of aquaglyceroporins but wider and less polar than the SF of strict aquaporins. Thus, the SF of AqpM can be viewed as a slightly wider and more-hydrophobic version of the SF from strict aquaporins. Regarding the hydrophobic residues that line the pore below the SF (residues highlighted in red in red), it is interesting to note that among the compared aquaporins these residues correspond to Val, Ile or Leu, except for AqpZ that has two aromatic residues instead in two of these nine positions. These considerations underscore the notion that differences in water permeability between aquaporins are most likely related to the size and polarity of the SF resulting from differences in critically located amino acids. AqpM, with its unique SF, clearly differs from the two other sub-families, exhibiting characteristics that are hybrid between an aquaglyceroporin such as GlpF and a strict aquaporins like AqpZ. Whether this differences permits AqpM to carry out different functions or whether it is just a “molecular fossil” representing an early stage in the evolution of the other aquaporins, as previously suggested [24, 26, 50], is an issue that remains to be determined. However, the apparent absence of aquaporin genes in a large number of microbial genomes , together with the assertion that given the small size and resultant large surface-to-volume ratio of individual microbial cells, lipid membranes are sufficiently permeable to support growth without water channels that enhance water permeability [51, 52], have opened the debate about possible alternative roles for microbial aquaporins other than simple water transport. These potential roles include osmotic and turgor sensors , involvement in freeze tolerance mechanisms  and transport of other neutral solutes such as CO2, H2S, NH3 .
We have performed the first molecular dynamics simulation of the archaeal aquaporin AqpM, the most studied representative of a new group of aquaporins that are not easily classifiable within the current functional sub-families of aquaporins. From our MD simulation we have provided estimations for key biophysical parameters of AqpM including: single channel osmotic permeability constant (p f ), single channel diffusive permeability constant (p d ), channel radius, potential water occupancy of the channel and water orientation inside the pore. The results obtained from the MD simulations of AqpM were compared with those obtained by similar simulation methods for the well characterized microbial representatives of the two main aquaporin sub-families, AqpZ and GlpF, a strict aquaporin and an aquaglyceroporin respectively. Our results extend the evidence that the size and polarity of the residues that comprise the SF of aquaporins control much of the selectivity and water permeability of aquaporins and supports the notion that AqpM belongs to a third class of aquaporins that exhibits properties that are intermediate between those of the other two classes. However, it is clear that more experimental and theoretical evidence is needed to assess whether this third class of aquaporins differs from the other sub-families in other functional aspect besides water permeability rates i.e., if it exhibits permeability to other uncharged solutes, related to the specific life style of the organisms in which they are found or whether it could represent an earlier evolutionary form, a molecular fossil, which existed before the divergence of the other two sub-families.
The crystal structure of the AqpM monomer was obtained from the Protein Data Bank [PDB :2F2B] . Hydrogen atoms were added using the psfgen plug-in from VMD v1.86  (considering pH = 7.5 for protonation states). The tetrameric structure of AqpM was generated with VMD v1.86  using the coordinate transformation matrices provided in the PDB file. The simulated periodic cell was constructed using VMD v1.88  and comprised the AqpM tetramer embedded in a pre-equilibrated 1,2-dipalmitoylphosphatidylcholine (DPPC) membrane patch of dimensions 100x100x58 Å. The energy of the crystallographic waters of the monomer was evaluated with the program DOWSER , keeping the waters that were located inside the pore. The system was solvated with a water layer of 20Å above and below the membrane. Ions (Na+, Cl-) were randomly placed to neutralize the system, reaching a final concentration of 50mM. The final dimensions of the periodic cell were 100x100x96Å comprising a total of 85,019 atoms. The full system and the structural features of AqpM can be seen in Figure 1.
The system was minimized and subjected to MD for 0.5 ns with fixed protein. The protein was released keeping Cα atoms constrained with a force constant of 5 kcal/mol/Å2. The full system was minimized and a slow relaxation procedure was performed in which the constraint applied to Cα atoms of the protein was decreased at a rate of 0.5 kcal/mol/ps until no constraint was applied. Then 21 ns of NPT-MD simulations were performed with the first nanosecond considered equilibration and the last 20 ns used for analysis. The time evolution of Cα - RMSD with respect to the final conformation of the minimization/relaxation procedure is shown for each AqpM monomer in Figure S1(a) of Additional file 1. The fluctuations of Cα atoms during the production run are presented a Cα – RMSF plot for each AqpM monomer in Figure S1(b) of Additional file 1.
The program NAMD  with CHARMM27 parameter set [56–58] was used for the simulation. Periodic boundary conditions were imposed and the particle mesh Ewald method [59, 60] was used for electrostatic forces calculation. Constant temperature (310K) and pressure (1 atm) were maintained by using Langevin dynamics .
Where v w =V w /N A is the average volume of a water molecule (2.99x10-23 cm3).
where θ is the angle between the membrane normal, along the z-axis, and the dipole moment of water. z is the position of the oxygen atom of water along the channel axis.
The average structure considered for the structural alignment presented in section Structure / Sequence Comparisons (Results and discussion) and in Figure 7 was obtained by averaging the Cartesian coordinates of all the atoms of one monomer of AqpM from the last nanosecond of the simulation at a 2ps resolution (500 frames).
Residue in one letter code followed by the number in the AqpM sequence. Between parenthesis, the residue in one letter code and number of the equivalent residue of AqpZ and GlpF. Ej: AqpM(AqpZ/GlpF).
We wish to thank Leonardo Sepulveda (Center for Biophysics and Computational Biology University of Illinois, Urbana, IL 61801, USA) for his critical review of the paper.
RAS is a CONICYT fellow. This work was partially funded by CONICYT (PFB16 and PFB03), Fondecyt 1090451 and supported by the High Performance Computing infrastructure of the National Laboratory for High Performance Computing (NLHPC) “Powered@NLHPC”.
This article has been published as part of BMC Genomics Volume 12 Supplement 4, 2011: Proceedings of the 6th International Conference of the Brazilian Association for Bioinformatics and Computational Biology (X-meeting 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/12?issue=S4
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.