|
|
||||||||
1 Department of Physiology and Cellular Biophysics, College of Physicians and Surgeons, Columbia University, New York, New York 10032, USA
2 Department of Opthalmology, College of Physicians and Surgeons, Columbia University, New York, New York 10032, USA
Reprint requests to: Dr. Jorge Fischbarg, Department of Physiology and Cellular Biophysics, College of P. & S., Columbia University, 630 West 168th Street, New York, NY 10032, USA; e-mail: jf20{at}columbia.edu; fax: 212-305-2461.
(RECEIVED February 7, 2001; FINAL REVISION April 16, 2001; ACCEPTED May 16, 2001)
3 Present address: Department of Ophthalmology, Weifang Medical College, Shandong, China. ![]()
4 Present address: Department of Radiation Oncology, 2 Donner, Hospital of the University of Pennsylvania, 3400 Spruce Street, Philadelphia, PA 19104, USA. ![]()
Article and publication are at www.proteinscience.org/cgi/doi/10.1110/ps.5901
| Abstract |
|---|
|
|
|---|
6Å wide), suggest that residues 7173 at the pore are accessible to extracellular mercurials. A 30-ps molecular dynamics simulation (at 300 K) starting with crystallographic coordinates of AQP1 showed that the width of the pore bottleneck (between Connolly surfaces) can vary (wavg = 3.9 Å,
= 0.75; hydrated AQP1). Thus, although the pore width would be
6 Å only for 0.0026 of the time, this might suffice for pCMBS to reach residues 7173. Alternative explanations such as passage of pCMBS across the AQP1 tetramer center or other unspecified transmembrane pathways cannot be excluded. Keywords: Water channel; pore; hourglass model; osmotic permeability
| Introduction |
|---|
|
|
|---|
The wild-type AQP1 water channel is characteristically sensitive to mercurials placed outside of cells. It has been established that one single cysteine residue in its chain (Cys 189) confers that property (Preston et al. 1993; Zhang et al. 1993), presumably by steric blockage of the aqueous pore through the protein. Replacing C189 by a non-sulfhydryl containing amino-acid-like serine (C189S) makes AQP1 mercurial insensitive (Preston et al. 1993; Zhang et al. 1993). Hence, by use of C189S as a template, further substitution of Cys residues at strategic locations can be used to test the accessibility to mercurials for the residues replaced. For this work, we substituted cysteine one at a time (Akabas et al. 1992; Xu and Akabas 1993) for three consecutive residues (7173) near the first NPA repeat (7678). The cysteine AQP1 mutants were expressed in Xenopus laevis oocytes and the osmotic permeability (Pf) of oocytes expressed was determined under control conditions and in the presence of the mercurial reagents HgCl2 and p-chloromercuribenzene-sulfonic acid (pCMBS). Unexpectedly, we found that all three mutants were sensitive to both mercurials. To study such behavior, we ran molecular dynamics (MD)simulations of AQP1 at 300 K, and found that its pore bottleneck can vary in width over a surprisingly wide range (avg. width = 3.9 Å,
= 0.75; hydrated AQP1). This seems consistent with pCMBS (width
6 Å) possibly traversing the pore if the incubation time is long enough. We discuss how residues 7173 could be accessible to external mercurials by this or alternative routes, and the possible implications for MIP protein selectivities.
| Results |
|---|
|
|
|---|
|
|
Molecular dynamics
We utilized the AQP1 coordinates (at 3.8 Å resolution) communicated by Murata el al. (2000). Prior to the simulations, a recommended step of energy minimization was performed as described in Materials and Methods. Figures 57![]()
![]()
, below, show the results of the MD simulations. The simulations ran for 30 ps, during which time the positions of the atoms oscillate within a steady-state range throughout the simulation. The events during the first
20 ps depicted here are adequate to describe the process.
|
|
|
Figure 6
shows the progression of pore widths during 20 ps (sampled at 100-fs intervals) for both simulations (AQP1 in vacuum, and hydrated). To define pore widths, three residues of interest (Phe 24, Leu 149, and Ile 191) were shown with Connolly (1983) surfaces (with InsightII; atom radius scale: 0.7 of VDW; probe radius 1.4 Å). The surfaces were at 0.8 Å of the H centers, from which pore widths were estimated accordingly from the interatomic distances. As can be seen, the estimated widths vary considerably during the short time interval examined.
Figure 7
shows the distribution of pore widths for the two simulations with hydrated AQP1 given in Figure 6
. As can be calculated from those data, although the average pore width was 3.9 ± 0.75 Å, there was a 0.003 probability that the pore would be
6 Å.
| Discussion |
|---|
|
|
|---|
The Pf values we determined for our mutants G72C/C189S and A73C/C189S were relatively high, but somewhat lower than those obtained with wild-type AQP1. Hence, these mutations did not appear to affect the overall protein fold. The
30% decrease in Pf we observed for them compared with the wild type (Fig. 2
, second and third bar groups) might be due to increased steric hindrance as C is bulkier than G or A. In support of this hypothesis, mutants A73F and A73M yielded much-reduced Pf in a prior study (Jung et al. 1994). However, different expression levels of the mutants cannot be excluded as an explanation.
It seems also interesting that the mutation S71C resulted in a decrease in Pf to one-third of that in wild-type or C189S mutant-expressing oocytes (Fig. 2
, fourth bar group ). S71 seems to be important to stabilize the molecular structure (Y. Fujiyoshi, pers. comm.), therefore, it is conceivable that replacement of this residue might lead to complex structural changes and lowered Pf.
Rather surprisingly, our main observation is that the three C mutants (from 71to 73) recovered the sensitivity to both HgCl2 and pCMBS added to the outside (see Fig. 2
). In support of our interpretations, we note that Jung et al. (1994) preceded us in showing that A73C/C189S exhibited mercurial sensitivity (
60% inhibition for 1 mM HgCl2, compared with
90% inhibition in the present data). From this result, they concluded that residue 73 could be placed near the cytoplasmic face of the single aqueous pathway across the protein, which appears remarkably prescient in view of recent structural information.
In our results, the membrane-impermeant sulfhydryl reagent pCMBS (see Fig. 2
) blocks the Pf corresponding to mutant expression by some 70% (just as it does in the case of wild-type AQP1). How the mercurials, especially the larger one, pCMBS, would reach these residues is less clear. One possibility is that the mutations in residues 7173 alter the stability of the molecule and produce slight conformation changes either enhancing the diameter of the pore or exposing cysteines that before were inaccessible to mercurials. For example, this might happen with Cys 102, which is close to the narrowest region of the pore (cf. Fig. 4
). However, another explanation that deserves attention is that the pore dimensions might vary. In the MD simulations at 300 K we report here, the width of the pore went from the 3 Å described by Murata et al (2000) to, if briefly, as much as
6.1 Å (Fig. 6
), a size large enough to allow pCMBS (
6 Å wide, cf. Fig. 3
) to traverse the pore. There is a range of values for the diameter of the AQP1 pore in different reports [
3.0 Å (Murata et al. 2000);
4.5 Å (Mitsuoka et al. 1999);
6.5 Å (Ren et al. 2000);
4.0 ± 0.5 Å (Ren et al. 2001)]. This variation might be due mostly to the limited resolutions achieved this far, but perhaps might also reflect the existence of different conformations. In this connection, it seems interesting that there is evidence that in lactose permease, as temperature increases, additional regions are exposed to pCMBS inhibition (Venkatesan et al. 2000). As a note of caution, given that the starting point of our simulations is a relatively low-resolution model, there may be an inherent inaccuracy, especially for the positions of the side chains. On the other hand, we observed comparatively large mobility for the helical segments with both the esff and the cvff force fields (not shown), so this may correspond to the general behavior of the protein. In a related subject, the time sequence for the hydrated AQP1 in Figure 6
shows what are roughly two overall pore width levels. Interestingly, this might indicate that the pore could be gated. More work is required to refine these aspects.
|
|
6 Å appears to be only 0.003. Hence, for pCMBS to migrate all of the length of the pore to reach the intracellular end, comparatively long incubation times would be necessary. Yet, this is what happens in practice as oocytes are preincubated with pCMBS typically for 2030 min. If pore widths are that variable, one would then need to contend with the possibility that substrates normally excluded by AQP1 such as urea and glycerol might occasionally traverse the pore. On one hand, such events are rather rare, and, hence, almost undetectable. On the other hand, an increase in pore width might not necessarily increase the selectivity for such substrates. Such selectivity might rather stem crucially from favorable binding sites rather than steric factors. That this may be so can be gleaned from the cogent arguments given to explain glycerol selectivity by the glycerol facilitator GlpF (a glyceroaquaporin), for which crystallographic coordinates at 2.2 Å resolution have been recently communicated (Fu et al. 2000). In that case, several characteristic arrangements along the pore would add up to contribute to enhance glycerol selectivity. In addition, we note that another glyceroaquaporin, MIP26, has a water conductance of only
15% that of AQP1 (Zampighi et al. 1995), whereas its pore diameter has to be more than adequate for water passage. Therefore, conversely, it could be that in some glyceroaquaporins, the motifs that enhance water selectivity in aquaporins would be either absent or have less weight.
Another question that this hypothesis poses is whether given substrates could interact with the protein and contribute to modify its dynamic behavior, enlarging the pore width as they move. Passage of pCMBS across AQP1 might be possible this way. Even if one contemplates the possibility of passage of pCMBS across another pathway such as the center of the AQP1 tetramer, from recent evidence, the corresponding width is only
3 Å (Ren et al. 2001), which is inadequate unless dynamic changes are once more postulated.
Lastly, it is conceivable that pCMBS would traverse the membrane through pathways other than AQP1. In this view, a carrier permeable to pCMBS may be coexpressed with AQP1, or the expression of AQP1 might lead to unspecified local changes in the lipid allowing the mercurial(s) across. These possibilities appear somewhat unlikely to us, but have not been ruled out by the present work.
| Materials and methods |
|---|
|
|
|---|
2600 bases (which included the full-length AQP1 cDNA) was ligated into pAlter-1 cut with KpnI and BamHI. This vector had its tetracycline site active, which was used to select and amplify the ligation product (AQP-pAlter) by use of Escherichia coli cells. Mutation oligonucleotides (27 to 33 bases) were designed (by use of the Gene Runner program, Hastings Software) so as to generate the appropriate cysteine replacement or substitution while simultaneously introducing a silent restriction site; this usually requires a 34 base mismatch (Table 1
|
Preparation of oocytes and measurement of Pf
Oocytes were removed from the ovaries of Xenopus laevis (NASCO), defolliculated, and separated (Echevarria et al. 1993). The largest undamaged oocytes (stages 5 and 6) were transferred to Barth's medium. Oocytes were injected with 50 nL of either water or cRNA (25 ng/µl) as described (Echevarria et al. 1993) and were incubated at 18°C for 72 h. Oocytes were then transferred to a glass-bottom chamber containing Barth's medium at room temperature. The oocyte equatorial cross-section was viewed with a Nikon TMs inverted microscope equipped with a video camera (model NC-65, Dage-MTI) connected to a monitor screen. Oocytes were superfused at first in isotonic Barth's solution (178 mOsm) for a period of 60 s, and then in hypotonic solution (by reducing NaCl, 15 mOsm) for another 100 s. A frame grabber recorded an image every 10 s, and a computer calculated oocyte area and volume from such image. Pf values were calculated from the induced change in oocyte volume,
![]() |
in which A was the area (assumed spherical) of the oocyte at zero time, dV/dt was the rate of volume change at zero time, and
C was the osmolarity gradient at zero time.
Molecular dynamics
MD were performed using as the starting position the AQP1 coordinates reported by Murata et al. (2000). We used the Discover 3 module of InsightII (Molecular Simulations Inc.; version 97.2); the force field was ESFF (Extensible Systemic Force Field). We ran simulations as follows: (1) with AQP1 in vacuum, and (2) with the exofacial and endofacial vestibules plus the pore of AQP1 hydrated (using the assembly/soak command of InsightII, which added 254 water molecules). Initially, before running dynamics, the energy was optimized with the steepest descent and conjugate algorithms (1000 steps each). Parameters for MD were as follows: equilibration 5000 steps, run 30,000 steps, temperature 300 K, time step 1 femtosecond, and no constraints. Otherwise, settings were default. History files for the 30-ps runs recorded the position of the atoms every 50 steps.
| Acknowledgments |
|---|
The publication costs of this article were defrayed in part by payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 USC section 1734 solely to indicate this fact.
| References |
|---|
|
|
|---|
Cheng, A., van Hoekm A.N., Yeager, M., Verkman, A.S., and Mitra, A.K. 1997. Three-dimensional organization of a human water channel. Nature 387: 627630.[CrossRef][Medline]
Chepelinsky, A.B. 1994. The MIP transmembrane channel family. In Handbook of membrane channels: Molecular and cellular physiology, (ed. C. Peracchia), pp. 413432, Academic Press, San Diego, CA.
Connolly, M.L. 1983. Solvent-accessible surfaces of proteins and nucleic acids. Science 221: 709713.
Echevarria, M., Kuangm, K., Iserovich, P., Li, J., Preston, G.M., Agre, P., and Fischbarg, J. 1993. Cultured bovine corneal endothelial cells express CHIP28 water channels. Am. J. Physiol. 265: C1349C1355.
Fu, D., Libson, A., Miercke, L.J., Weitzman, C., Nollert, P., Krucinski, J., and Stroud, R.M. 2000. Structure of a glycerol-conducting channel and the basis for its selectivity. Science 290: 481486.
Heymann, J.B. and Engel, A. 2000. Structural clues in the sequences of the aquaporins. J. Mol. Biol. 295: 10391053.[CrossRef][Medline]
Jung, J.S., Preston, G.M., Smith, B.L., Guggino, W.B., and Agre, P. 1994. Molecular structure of the water channel through aquaporin CHIP. The hourglass model. J. Biol. Chem. 269: 1464814654.
Mitsuoka, K., Murata, K., Walz, T., Hirai, T., Agre, P., Heymann, J.B., Engel, A., and Fujiyoshi, Y. 1999. The structure of aquaporin-1 at 4.5-A resolution reveals short alpha- helices in the center of the monomer. J. Struct. Biol. 128: 3443.[CrossRef][Medline]
Murata, K., Mitsuoka, K., Hirai, T., Walz, T., Agre, P., Heymann, J.B., Engel, A., and Fujiyoshi, Y. 2000. Structural determinants of water permeation through aquaporin-1. Nature 407: 599605.[CrossRef][Medline]
Preston, G.M. and Agre, P. 1991. Isolation of the cDNA for erythrocyte membrane protein of 28 kilodaltons: Member of an ancient channel familiy. Proc. Natl. Acad. Sci. 88: 1111011114.
Preston, G.M., Jung, J.S., Guggino, W.B., and Agre, P. 1993. The mercury-sensitive residue at cysteine 189 in the CHIP28 water channel. J. Biol. Chem. 268: 1720.
Ren, G., Cheng, A., Reddy, V., Melnyk, P., and Mitra, A.K. 2000. Three-dimensional fold of the human AQP1 water channel determined at 4 A resolution by electron crystallography of two-dimensional crystals embedded in ice. J. Mol. Biol. 301: 369387.[CrossRef][Medline]
Ren, G., Reddy, S., Cheng, A., Melnyk, P., and Mitra, A.K. 2001. Visualization of a water-selective pore by electron crystallography in vitreous ice. Proc. Natl. Acad. Sci. 98: 13981403.
Venkatesan, P., Hu, Y., and Kaback, H.R. 2000. Site-directed sulfhydryl labeling of the lactose permease of Escherichia coli: Helix X. Biochemistry 39: 1065610661.[CrossRef][Medline]
Xu, M. and Akabas, M.H. 1993. Amino acids lining the channel of the gamma-aminobutyric acid type A receptor identified by cysteine substitution. J. Biol. Chem. 268: 2150521508.
Zampighi, G.A., Kreman, M., Boorer, K.J., Loo, D.D.F., Bezanilla, F., Chandy, G., Hall, J.E., and Wright, E.M. 1995. A method for determining the unitary functional capacity of cloned channels and transporters expressed in Xenopus laevis oocytes. J. Membr. Biol. 148: 6578.[Medline]
Zhang, R., van-Hoek, A.N., Biwersi, J., and Verkman, A.S. 1993. A point mutation at cysteine 189 blocks the water permeability of rat kidney water channel CHIP28k. Biochemistry 32: 29382941.[CrossRef][Medline]
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
This article has been cited by other articles:
![]() |
G. J. Weber, A. P. Mehr, J. C. Sirota, S. G. Aller, S. E. Decker, D. C. Dawson, and J. N. Forrest Jr. Mercury and zinc differentially inhibit shark and human CFTR orthologues: involvement of shark cysteine 102 Am J Physiol Cell Physiol, March 1, 2006; 290(3): C793 - C801. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Kong and J. Ma Dynamic mechanisms of the membrane water channel aquaporin-1 (AQP1) PNAS, November 15, 2001; (2001) 251507998. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Kong and J. Ma Dynamic mechanisms of the membrane water channel aquaporin-1 (AQP1) PNAS, December 4, 2001; 98(25): 14345 - 14349. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |