help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Originally published as Biophys J. BioFAST on March 31, 2006.
doi:10.1529/biophysj.106.084061
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.106.084061v1
90/10/L73    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chen, H.
Right arrow Articles by Voth, G. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chen, H.
Right arrow Articles by Voth, G. A.
Biophysical Journal 90:L73-L75 (2006)
© 2006 The Biophysical Society

Origins of Proton Transport Behavior from Selectivity Domain Mutations of the Aquaporin-1 Channel

Hanning Chen, Yujie Wu and Gregory A. Voth

Center for Biophysical Modeling and Simulation, Department of Chemistry, University of Utah, Salt Lake City, Utah, 84112-0850

Correspondence: Address reprint requests and inquiries to Gregory A. Voth, Tel.: 801-581-7272; Fax: 801-581-4353; E-mail: voth{at}chem.utah.edu.


   ABSTRACT
TOP
 ABSTRACT
ACKNOWLEDGEMENTS
 REFERENCES
 
The permeation free-energy profile and maximum ion conductance of proton transport along the channel of three aquaporin-1 (AQP1) mutants (H180A/R195V, H180A, and R195V) are calculated via molecular dynamics simulations and Poisson-Nernst-Planck theory. The proton dynamics was described by the multistate empirical valence bond (MS-EVB) model. The results reveal three major contributions to the overall free-energy barrier for proton transport in AQP1: 1), the bipolar field, 2), the electrostatic repulsion due to the Arg-195 residue, and 3), the dehydration penalty due to the narrow channel pore. The double mutation (H180A/R195V) drastically drops the overall free-energy barrier by roughly 20 kcal/mol via simultaneously relaxing the direct electrostatic interaction (by R195V) and dehydration effect (by H180A).

As the first identified and characterized integral transmembrane water channel (1Go), the aquaporin-1 (AQP1) channel has been studied extensively. An important property of this channel is its exceptionally high water permeability (~0.02 cm s–1) in liposomes (2Go,3Go), with efficient blockage of other neutral solutes and ions, including protons. The proton-blockage property of AQP1, and of the aquaporin superfamily in general, is especially interesting given the fact that the proton is the smallest and most mobile cation in aqueous solutions. Proton mobility is three times higher than that of water and is attributed mainly to its capability to shuttle rapidly along the hydrogen bond network of water molecules according to the Grotthuss mechanism (4Go). This interesting property of aquaporins has attracted considerable research effort in recent years, especially in the molecular simulation field (5Go–11Go). One focus of attention has been the electrostatic effects around the fingerprint Asn-Pro-Ala (NPA) motif where two half-membrane-spanning {alpha}-helices HB (at the cytoplasm side) and HE (at the extracellular side) (Fig. 1) meet with opposite orientations (12Go). It has been proposed that the opposite {alpha}-helical orientations might generate a bipolar field, preventing protons to move through the channel (5Go–9Go). This hypothesis has been partially supported by the finding that the main free-energy barrier is located at the NPA motif (5Go–9Go). A molecular simulation study has also revealed a clear secondary free-energy barrier located at the extracellular side roughly 8 Å from the main barrier (8Go). This barrier corresponds to a constriction region with a diameter of ~2 Å (12Go–13Go). The constriction region allows only one water molecule to pass through and is therefore believed to correspond to the selectivity filter (SF), which is formed by the following four residues: Phe-56, His-180, Cys-189, and Arg-195 (Fig. 1). The His-180 residue, whereas highly conserved in water-selective orthodox aquaporins, is usually replaced by Gly in aquaglyceroporins (14Go) to relax steric interaction with the alkyl chain of glycerol. Within the water permeation pathway, the Cys-189 residue is believed to be responsible for blocking Hg2+ by binding the ion to the backbone carbonyl oxygen that flanks the pore (15Go). The function of the aromatic hydrophobic Phe-56 is believed to help the water molecules to form stronger hydrogen bonds with the positively charged Arg-195 (16Go), which, if artificially made uncharged, would favor the entry of cations from the extracellular side (17Go).


Figure 1
View larger version (71K):
[in this window]
[in a new window]
  		FIGURE 1  A snapshot of the AQP1 double mutant (H180A/R195V). The Phe-56, Ala-180, Val-195, and Cys-189 are colored purple, red, yellow, and green, respectively. The HB and HE helices are modeled as highlighted with the indicated dipole field. The cytoplasm side is at the top.

 
Very recently, the conductance of water, neutral solute, and cations for several AQP1 mutants (H180A, R195V, H180A/R195V, F56A/H180A) was reported by Beitz et al. (18Go). Their results indicated that cations can pass through the channel, without affecting the water passage, simply by mutagenically removing the Arg-195 residue. Though the H180A mutation alone does not result in cation leakage, it does increase the cation conductance when combined with the R195V mutation (Fig. 1). These results indicate a very important role of the SF domain in cation blockage. Specifically, they suggest that the electrostatic repulsion due to the Arg-195 residue and the dehydration penalty (10Go,11Go) due to the narrow channel are also very important effects for cation blockage.

Inspired by this, we report here potential of mean force (PMF) results for the H180A/R195V, H180A, and R195V mutants (Fig. 2). The PMF is the free-energy profile for the cation (in this case excess proton) through the channel. To quantify the bipolarity contribution to the free-energy barrier, we also include the result for a wild-type (WT) channel that has the partial charges on the backbone of the HB and HE helices set to zero. This channel is hereinafter referred to as the no-backbone-charge (NBC) mutant. Such a computational mutant has also been studied for the GlpF aquaporin channel with a hybrid MD and continuum electrostatics approach (7Go). These calculations were molecular dynamics simulations as in our previous study (8Go,9Go,19Go). In particular, the dynamics of the Grotthuss shuttling excess proton was reliably described by the MS-EVB2 model (20Go,21Go), and the PMFs were calculated using the umbrella sampling and weighted histogram analysis methods. The initial structure of AQP1 was based on a crystal structure (Protein Data Bank entry code, 1J4N) (13Go), which had been equilibrated with a fully solvated lipid bilayer for 3 ns in a MD simulation.


Figure 2
View larger version (20K):
[in this window]
[in a new window]
  		FIGURE 2  The PMF of proton transport along the channel for the WT AQP1 and the H180A/R195V, R195V, H180A, and NBC mutants. (WT and NBC results from Chen et al. (9Go).) Average error bars for these curves are ±0.4 kcal/mol.

 
The PMF profile indicates that all of the four AQP1 variants have substantially lower free-energy barriers than that of the WT channel. The primary PMF peak at the NPA domain was lowered from 28 kcal/mol to 12 kcal/mol by solely turning off the dipoles of the HB and HE helices, and the bipolar orientation of the water file as observed in the WT and in the other mutants also disappeared. This indicates a clear contribution to the overall free-energy barrier from the NPA bipolarity whose origin can be attributed to the backbone dipoles of the two opposing helices. Compared to the other mutants, this computational mutation also results in a 4.0-Å shift of the primary barrier toward the center of the lipid bilayer, where the dehydration penalty is presumed to be maximal. For the H180A/R195V and R195V mutants, however, the removal of the Arg-195 residue at the narrow SF domain causes the secondary PMF peak of 17 kcal/mol to drop, which can further facilitate the proton to entry from the bulk water into the narrow AQP1 channel. By reducing the dehydration effect, the H180A mutation lowers the overall barrier by 13 kcal/mol, but it maintains an 11 kcal/mol barrier at the SF domain. Although the SF domain has been speculated for a long time to act as a solute filter because of its small pore radius, our study provides clear evidence that the SF domain is also essential in blocking proton transport mainly by direct electrostatic interaction between its positively charged Arg-195 and the excess proton and by forming a narrow pore.

To compare with the experimental proton conductance (18Go), the maximum ion conductance gmax for the H180A/R195V mutant was calculated based on the Poisson-Nernst-Planck theory (22Go). The gmax turns out to be 8 x 10–3 pS for the H180A/R195V mutant, in reasonable agreement with the experimental lower limit value ~10–4 pS (18Go).

Our results, taken together, suggest that the direct electrostatics, dehydration, and bipolar field are the three major contributing effects to the free-energy barrier to PT in AQP1.


   ACKNOWLEDGEMENTS
TOP
 ABSTRACT
 ACKNOWLEDGEMENTS
REFERENCES
 
We gratefully acknowledge the Center for High Performance Computing at the University of Utah, and the National Center for Supercomputing Applications for computational resources.

This research was supported by the National Institutes of Health (RO1-GM53148).

Submitted on February 25, 2006; accepted for publication March 13, 2006.


   REFERENCES
TOP
 ABSTRACT
 ACKNOWLEDGEMENTS
 REFERENCES
 
1. Denker, B. M., B. L. Smith, F. P. Kuhajda, and P. Agre. 1988. Identification, purification and partial characterization of a novel Mr 28,000 integral membrane protein from erythrocytes and renal tubules. J. Biol. Chem. 263:15634–15642.[Abstract/Free Full Text]

2. Zeidel, M. L., S. V. Ambudkar, B. L. Smith, and P. Agre. 1992. Reconstitution of functional water channels in liposomes containing purified red cell CHIP28 protein. Biochemistry. 31:7436–7440.[CrossRef][Medline]

3. Preston, G. M., T. P. Caroll, W. B. Guggino, and P. Agre. 1992. Appearance of water channels in Xenopus oocytes expressing red cell CHIP28 protein. Science. 256:385–387.[Abstract/Free Full Text]

4. Agmon, N. 1995. The Grotthuss mechanism. Chem. Phys. Lett. 244:456–462.[CrossRef]

5. de. Groot, B. L., T. Frigato, V. Helms, and H. Grubmüller. 2003. The mechanism of proton exclusion in the aquaporin-1 water channel. J. Mol. Biol. 333:279–293.[CrossRef][Medline]

6. Chakrabarti, N., E. Tajkhorshid, B. Roux, and R. Pomès. 2004. Molecular basis of proton blockage in aquaporins. Structure. 12:65–74.[Medline]

7. Chakrabarti, N., B. Roux, and R. Pomès. 2004. Structural determinants of proton blockage in aquaporins. J. Mol. Biol. 343:493–510.[CrossRef][Medline]

8. Ilan, B., E. Tajkhorshid, K. Schulten, and G. A. Voth. 2004. The mechanism of proton exclusion in aquaporin channel. Proteins. 55:223–228.[CrossRef][Medline]

9. Chen, H., B. Ilan, Y. Wu, C. J. Burnham, F. Zhu, E. Tajkhorshid, K. Schulten, and G. A. Voth. 2006. Electrostatics and charge delocalization in proton channels: the aquaporin channels and proton blockage. J. Mol. Biol. In press.

10. Burykin, A., and A. Warshel. 2003. What really prevents proton transport through Aquaporin? Charge self-energy versus proton wire proposals. Biophys. J. 85:3696–3706.[Abstract/Free Full Text]

11. Burykin, A., and A. Warshel. 2004. On the origin of the electrostatic barrier for proton transport in aquaporin. FEBS Lett. 570:41–46.[CrossRef][Medline]

12. Murata, K., K. Mitsuoka, T. Hirai, T. Walz, P. Agre, J. B. Heymann, A. Engel, and Y. Fujiyoshi. 2000. Structural determinants of water permeation through aquaporin-1. Nature. 407:599–605.[CrossRef][Medline]

13. Sui, H., B. G. Han, J. K. Lee, P. Walian, and B. K. Jap. 2001. Structural basis of water-specific transport through the AQP1 water channel. Nature. 414:872–878.[CrossRef][Medline]

14. Fu, D., A. Libson, L. J. W. Miercke, C. Weitzman, P. Nollert, J. Krucinski, and R. M. Stroud. 2000. Structure of a glycerol-conducting channel and the basis for its selectivity. Science. 290:525–530.

15. Preston, G. M., J. S. Jung, W. B. Guggino, and P. Agre. 1993. The mercury-sensitive residue at cysteine 189 in the CHIP28 water channel. J. Biol. Chem. 268:17–20.[Abstract/Free Full Text]

16. Liu, Z., J. M. Carbery, P. Agre, and B. P. Rosen. 2004. Arsenic trioxide uptake by human and rat aquaglyceroporins. Biochem. Biophys. Res. Commun. 316:1178–1185.[CrossRef][Medline]

17. Miloshevsky, G. V., and P. C. Jordan. 2004. Water and ion permeation in bAQP1 and GlpF channels: a kinetic Monte Carlo study. Biophys. J. 87:3690–3702.[Abstract/Free Full Text]

18. Beitz, E., B. Wu, L. M. Holm, J. E. Schultz, and T. Zeuthen. 2006. Point mutations in the aromatic/arginine region in aquaporin 1 allow passage of urea, glycerol, ammonia, and protons. Proc. Natl. Acad. Sci. USA. 103:269–274.[Abstract/Free Full Text]

19. Wu, Y., B. Ilan, and G. A. Voth. 2006. Electrostatics and proton delocalization in proton channels. II. The minimalist synthetic LS2 channel and proton selectivity. J. Mol. Biol. In press.

20. Day, T. J. F., A. V. Soudackov, M. Cuma, U. W. Schmitt, and G. A. Voth. 2002. A second generation multistate empirical valence bond model for proton transport in aqueous systems. J. Chem. Phys. 117:5839–5849.[CrossRef]

21. Schmitt, U. W., and G. A. Voth. 1999. The computer simulation of proton transport in water. J. Chem. Phys. 111:9361–9381.[CrossRef]

22. Levitt, D. 1991. General continuum theory for multi-ion channel. I. Theory. Biophys. J. 59:271–277.[Abstract/Free Full Text]




This article has been cited by other articles:


Home page
 Biophys. JHome page
H. Chen, Y. Wu, and G. A. Voth
Proton Transport Behavior through the Influenza A M2 Channel: Insights from Molecular Simulation
Biophys. J., November 15, 2007; 93(10): 3470 - 3479.
[Abstract] [Full Text] [PDF]

Home page
 Biophys. JHome page
M. Hashido, A. Kidera, and M. Ikeguchi
Water Transport in Aquaporins: Osmotic Permeability Matrix Analysis of Molecular Dynamics Simulations
Biophys. J., July 15, 2007; 93(2): 373 - 385.
[Abstract] [Full Text] [PDF]

Home page
 J. Gen. Physiol.Home page
A. B. Mamonov, R. D. Coalson, M. L. Zeidel, and J. C. Mathai
Water and Deuterium Oxide Permeability through Aquaporin 1: MD Predictions and Experimental Verification
J. Gen. Physiol., July 1, 2007; 130(1): 111 - 116.
[Abstract] [Full Text] [PDF]

Home page
 Biophys. JHome page
H. Chen, B. Ilan, Y. Wu, F. Zhu, K. Schulten, and G. A. Voth
Charge Delocalization in Proton Channels, I: The Aquaporin Channels and Proton Blockage
Biophys. J., January 1, 2007; 92(1): 46 - 60.
[Abstract] [Full Text] [PDF]

Home page
 Proc. Natl. Acad. Sci. USAHome page
L. Vrbka, J. Vondrasek, B. Jagoda-Cwiklik, R. Vacha, and P. Jungwirth
Quantification and rationalization of the higher affinity of sodium over potassium to protein surfaces
PNAS, October 17, 2006; 103(42): 15440 - 15444.
[Abstract] [Full Text] [PDF]

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.106.084061v1
90/10/L73    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Chen, H.
Right arrow Articles by Voth, G. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Chen, H.
Right arrow Articles by Voth, G. A.

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Copyright © 2006 by the Biophysical Society.