|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Department of Chemistry, Swarthmore College, Swarthmore, Pennsylvania 19081, USA2 Department of Biochemistry and Cell Biology, Center for Structural Biology, and 3 Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794-5115, USA4 Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
Reprint requests to: Kathleen P. Howard, Department of Chemistry, Swarthmore College, Swarthmore, PA 19081, USA; e-mail: khoward1{at}swarthmore.edu; fax: (610) 328-7355; or Steven O. Smith, Department of Biochemistry and Cell Biology, Center for Structural Biology, Stony Brook University, Stony Brook, NY 11794-5115, USA; e-mail: steven. o.smith{at}sunysb.edu; fax: (631) 632-8575.
(RECEIVED September 10, 2004; FINAL REVISION November 30, 2004; ACCEPTED December 3, 2004)
| Abstract |
|---|
|
|
|---|
Keywords: transmembrane peptides; helix association; 2H MAS NMR; ATR-FTIR
Article published online ahead of print. Article and publication date are at http://www.proteinscience.org/cgi/doi/10.1110/ps.041110605.
| Introduction |
|---|
|
|
|---|
Scanning mutagenesis along the length of the hydrophobic sequence, followed by gel electrophoresis, has been used to locate interfacial regions in TM helical bundles (Simmerman et al. 1996; Ruan et al. 2004). An alternative method that avoids the need for extensive peptide production is the use of 2H NMR spectroscopy (Ying et al. 2000). In this case, the motion of deuterated amino acid side-chains is measured around at least one turn of a TM helix. The 2H NMR line shape provides direct information on side-chain dynamics in hydrated lipid bilayers above the lipid phase transition temperature. The side-chains packed within a helixhelix interface exhibit restricted motion relative to side-chains oriented toward the surrounding lipid bilayer (Ying et al. 2000).
In this study, 2H NMR is used to establish the rotational orientation of a TM peptide designed to form a homo-oligomeric helical bundle. The peptide belongs to a family of model peptides derived from a membrane-solubilized version of the water-soluble coiled-coil GCN4-P1 (Zhou et al. 2001; Lear et al. 2003). The variant studied here (GCN4MS1-N7N14) forms a stable trimer, as demonstrated by analytical ultracentrifugation and SDSpolyacrylamide gel electrophoresis (Lear et al. 2003). The specificity and strength of association of this particular peptide is attributed to two asparagine residues located at positions 7 and 14 that are located in "a" positions of the helical repeat and are thought to engage in interhelical hydrogen bonding (Fig. 1
). For the 2H NMR studies, three different peptides were synthesized, each with a single deuterium-labeled leucine at either position 8, 9, or 10 in the peptide sequence. These leucines correspond to positions "b," "c," and "d" of the helical repeat and are located between the two critical asparagines. L8 and L9 are in "b" and "c" positions and are expected to face the lipid bilayer, while L10 is in the "d" position and is predicted to be in the helical interface.
|
| Results and Discussion |
|---|
|
|
|---|
-helical secondary structure (Tamm and Tatulian 1997). Although both the amide I and II normal modes are sensitive to protein secondary structure, the amide I vibration is the most widely used to determine protein conformation. The dichroic ratio of the amide I band is sensitive to the orientation of the GCN4MS1-N7N14 helix relative to the plane of membrane. The observed dichroic ratio of 3.4 corresponds to a helix tilt of ~13° relative to the membrane normal (see Materials and Methods). These data indicate that the GCN4MS1-N7N14 peptides have a helical conformation with a TM orientation in DMPC bilayers.
|
The strategy in these experiments is to deuterate consecutive leucines in the sequence of the GCN4MS1-N7N14 peptide. A single deuterated methyl group at the end of the highly flexible leucine side-chain provides a sensitive probe of molecular motion. We have previously shown that the motion of amino acids located in TM helix interfaces is restricted relative to amino acids oriented toward the lipid membrane (Ying et al. 2000). By characterizing side-chain motion as a function of the TM sequence, one can map out the interface of interacting TM helices.
2H MAS spectra of the three leucine-labeled peptides incorporated into DMPC bilayers are shown in Figure 3
. The MAS side bands are spaced at 3 kHz (the spinning frequency), and their intensities trace out the deuterium line shape. The intense resonance in the center of each spectrum is due to residual HDO, while the rest of the signal arises from a
-methyl labeled (CD3) leucine. Because of rapid methyl group rotation, the three deuterium sites in CD3 are equivalent. Rapid methyl group rotation reduces the quadrupole coupling constant to a maximum value of ~40 kHz compared with the maximum theoretical value of ~127 kHz in the rigid limit (Keniry 1989). Additional rotations and librations about the C
C
and C
C
bonds, as well as rotational diffusion about the axis of the trimer, further reduce the quadrupolar interaction and produce the observed line shapes (Batchelder et al. 1982). These motions result in narrowing of the spinning sideband manifold and a line shape with an asymmetry parameter (
) close to 1.0.
|
|
-helix is minor for the length of peptides studied here.
Figure 5
shows a cross section of a molecular model of GCN4MS1-N7N14 at the level of L8 to L10 that is consistent with the 2H NMR data and the designed structure as pictured in Figure 1
. L8 and L9 at the "b" and "c" positions of the helical wheel face outward from the bundle. L10 at the "d" position is packed in the interior of the helical bundle. N7 and N14 at the "a" positions above and below the targeted leucines form interhelical hydrogen bonded rings in the low energy structure.
|
Conclusions
In summary, we have used 2H MAS NMR spectroscopy to determine the rotational orientation of the TM helices within the trimer structure formed by the designed peptide GCN4MS1-N7N14. The 2H NMR line shapes of the three consecutive leucines are consistent with a trimer structure that is stabilized by interhelical hydrogen bonding of as-paragines at positions 7 and 14. More generally, the study describes a simple approach for mapping the interface of interacting TM helices. The approach complements gel electrophoresis, analytical ultracentrifugation, and fluorescence resonant energy transfer, which provide information on the oligomerization states of membrane proteins and TM peptides. Helix interactions are important in the folding of polytopic membrane proteins and in the oligomerization of membrane proteins and peptides having a single TM helix. De novo design in combination with 2H MAS NMR provides a way to understand and test the interactions that govern the association of the TM helices in these proteins.
| Materials and methods |
|---|
|
|
|---|
Reconstitution into lipid bilayers
HPLC purified peptides were cosolubilized with DMPC (Avanti Polar Lipids) in ethanol at a peptidelipid molar ratio of 1:100. Ethanol was removed under a stream of nitrogen and the peptide/ lipid films placed under high vacuum overnight to remove traces of solvent. The samples were then hydrated with deuterium-depleted water, vortexed, and incubated above the gelliquid crystalline phase temperature of DMPC.
Polarized infrared spectroscopy
Polarized ATR FTIR spectra were obtained on a Bruker IFS 66V/S spectrometer. Multilamellar vesicles at a concentration of 10 mg lipid/mL were layered on a germanium internal reflection element using a slow flow of air directed at an oblique angle to the ATR plate to form an oriented multilamellar lipid-peptide film. ATR FTIR spectra were acquired at a resolution of 4 cm1 using parallel and perpendicularly polarized light. Each spectrum represents the average of 1000 transients, and a background spectrum of the germanium plate without sample was subtracted in each case.
The dichroic ratio (RATR) is defined as the ratio between absorption of parallel (A||) and perpendicular (A
) polarized light (Smith et al. 2002) and was used to calculate the order parameter Smeas using the equation,
![]() | (1) |
where
is the angle between the helix director and the normal of the internal-reflection element, and
is the angle between the helix director and the transition-dipole moment of the amide I vibrational mode. The order parameter depends on the electric field amplitudes, Ex = 1.399, Ey = 1.514, and Ez = 1.621, in thick film limit, which is applicable for the concentration of lipid used here.
For our calculations, the value of
(42.4°) was derived from parallel experiments on bacteriorhodopsin, where the average helix tilt (14°) is known independently from the crystal structure (Goormaghtigh et al. 1999). Smem is an order parameter that accounts for disorder in the orientation of the reconstituted peptide. The value of
used for our calculations implies a value of Smem between 0.8 and 0.9, corresponding the order parameter calculated for bacteriorhodopsin from measurements of mosaic spread (see Smith et al. 2002). We assume that the amount of order cannot be more than that for the highly ordered two-dimensional purple membranes containing bacteriorhodopsin. As a result, the value we estimate for the helix tilt of GCN4MS1-N7N14 (~13°) is an upper limit for the angle between the helix axis and the bilayer normal. If the peptide orientation is more disordered than bacteriorhodopsin, then the true helix tilt for GCN4MS1-N7N14 will be <13°. The dichroic ratios of the asymmetric and symmetric lipid CH2 stretching modes at 2924 cm1 and 2852 cm1 were 1.19 and 1.09, respectively, which correspond to order parameters of 0.64 and 0.68, and are consistent with well-oriented lipids (Smith et al. 1994).
Deuterium NMR spectroscopy
Deuterium NMR spectra were obtained at a 2H frequency of 55.2 MHz on a Bruker Avance NMR spectrometer using a 4-mm MAS probe. A MAS frequency of 3 kHz was used to increase the number of spinning side bands. Single pulse excitation was employed using a 4.2-µsec 90° pulse, followed by a 30-µsec delay before data acquisition. The repetition delay was 0.25-sec. A total of 600,000800,000 transients were averaged for each spectrum and processed using a 200-Hz exponential line broadening function. Spectra were obtained at 30°C.
MAS deuterium spectra were simulated using the program SIMPSON version 1.1.0 (Bak et al. 2000) with a MAS frequency of 3 kHz. The asymmetry parameter was set to 1.0 to reflect extensive motional averaging. Only the quadrupole coupling constant was varied to quantify the breadth of the line shape.
Molecular modeling
Computational searches were performed to obtain the lowest energy structure of the GCN4MS1-N7N14 trimer consistent with the NMR data. A library of possible trimer structures was generated from bundles of three ideal
-helices (3.5 residues/turn) that differed in the radius of the bundle and the rotational orientation of the individual helices. A supercoil pitch of 180 Å was subsequently applied to each bundle (Betz and DeGrado 1996; Dieckmann and DeGrado 1997; Pinto et al. 1997). Two sets of searches were performed: The first sampled over radii between 5.0 and 8.0 Å in 0.25 Å increments and over 360° of helix rotation in 10° increments, and the second search expanded on the global minimum of the first, sampling between 6.3 and 6.7 Å in 0.1 Å increments and over 20° of rotational space in 2° increments. For each backbone sampled, side-chain conformations were selected using the Goldstein Dead End Elimination criterion (Desmet et al. 1992; Goldstein 1994) followed by exhaustive enumeration of the remaining conformations. A backbone-dependent rotamer library was used to limit the conformational space to
-helix-compatible conformations (Bower et al. 1997). Only on-rotamer conformations were sampled, and the same library was used in both coarse-and fine-grained conformational searches. Energies were calculated using a van der Waals (vdW) 126 potential and a hydrogen bonding 1210 potential function and a cosine angular dependence term based on donor and acceptor hybridization state (Gordon et al. 1999). Atom radii and well-depth parameters for vdW energies were taken from the united-atom set of AMBER parameters (Weiner et al. 1984). Hydrogen bonds were uniformly modeled with an ideal distance of 2.8 Å and a well depth of 2.0 kcal/mol. No electrostatics terms were included. All backbone generation and side-chain optimization was performed using the protCAD library of protein manipulation tools (Summa 2002). The lowest energy structure from the second search was subject to further minimization using DISCOVER97.0 (MSI/Biosym) and the AMBER force field (Weiner et al. 1984). The final structure had a radius of 6.4 Å and the helix rotation shown in Figure 5
.
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
Bak, M., Rasmussen, J.T., and Nielsen, N.C. 2000. SIMPSON: A general simulation program for solid-state NMR spectroscopy. J. Magn. Res. 147: 296330.[CrossRef][Medline]
Batchelder, L.S., Sullivan, C.E., Jelinski, L.W., and Torchia, D.A. 1982. Characterization of leucine side-chain reorientation in collagen-fibrils by solid state 2H NMR. Proc. Natl. Acad. Sci. 79: 386389.
Betz, S.F. and DeGrado, W.F. 1996. Controlling topology and native-like behavior of de novo designed transmembrane peptides: Design and characterization of antiparallel four-stranded coiled coils. Biochemistry 35: 69556962.[CrossRef][Medline]
Bloom, M. and Smith, I.C.P. 1985. Manifestations of lipidprotein interactions in deuterium NMR. In Progress in proteinlipid interactions, Vol. 1 (eds. A. Watts and J.J.H.H.M. De Pont), pp. 6188. Elsevier, Amsterdam.
Bower, M.J., Cohen, F.E., and Dunbrack, R.L. 1997. Prediction of protein side-chain rotamers from a backbone-dependent rotamer library: A new homology modeling tool. J. Mol. Biol. 267: 12681282.[CrossRef][Medline]
Chothia, C., Levitt, M., and Richardson, D. 1981. Helix to helix packing in proteins. J. Mol. Biol. 145: 215250.[CrossRef][Medline]
DeGrado, W.F., Gratkowski, H., and Lear, J.D. 2003. How do helixhelix interactions help determine the folds of membrane proteins?: Perspectives from the study of homo-oligomeric helical bundles. Protein Sci. 12: 647665.
Desmet, J., Demaeyer, M., Hazes, B., and Lasters, I. 1992. The dead-end elimination theorem and its use in protein side-chain positioning. Nature 356: 539542.[CrossRef]
Dieckmann, G.R. and DeGrado, W.F. 1997. Modeling transmembrane helical oligomers. Curr. Opin. Struc. Biol. 7: 486494.[CrossRef][Medline]
Fringeli, U.P. 1993. In situ infrared attenuated total reflection membrane spectroscopy. In Internal reflection spectroscopy: Theory and applications (ed. F.M. Mirabella), pp. 255324. Marcel Dekker Inc., New York.
Goldstein, R.F. 1994. Efficient rotamer elimination applied to protein side-chains and related spin-glasses. Biophys. J. 66: 13351340.
Goormaghtigh, E., Raussens, V., and Ruysschaert, J. 1999. Attenuated total reflection infrared spectroscopy of proteins and lipids in biological membranes. Biochim. Biophys. Acta 1422: 105185.[Medline]
Gordon, D.B., Marshall, S.A., and Mayo, S.L. 1999. Energy functions for protein design. Curr. Opin. Struc. Biol. 9: 509513.[CrossRef][Medline]
Gratkowski, H., Dai, Q., Wand, J., DeGrado, W.F., and Lear, J.D. 2002. Cooperativity and specificity of association of a designed transmembrane peptide. Biophys. J. 83: 16131619.
Keniry, M.A. 1989. Solid-state deuterium nuclear magnetic resonance spectroscopy of proteins. Meth. Enzymol. 176: 376386.[Medline]
Lear, J.D., Gratkowski, H., Adamian, L., Liang, J., and DeGrado, W.F. 2003. Position-dependence of stabilizing polar interactions of asparagine in transmembrane helical bundles. Biochemistry 42: 64006407.[CrossRef][Medline]
Pinto, L.H., Dieckmann, G.R., Gandhi, C.S., Papworth, C.G., Braman, J., Shaughnessy, M.A., Lear, J.D., Lamb, R.A., and DeGrado, W.F. 1997. A functionally defined model for the M2 proton channel of influenza A virus suggests a mechanism for its ion selectivity. Proc. Natl. Acad. Sci. 94: 1130111306.
Ruan, W., Linder, E., and Langosch, D. 2004. The interface of a membrane-spanning leucine zipper mapped by asparagine-scanning mutagenesis. Protein Sci. 13: 555559.
Simmerman, H.K.E., Kobayashi, Y.M., Autry, J.M., and Jones, L.R. 1996. A leucine zipper stabilizes the pentameric membrane domain of phospholamban and forms a coiled-coil pore structure. J. Biol. Chem. 271: 59415946.
Smith, S.O.S., Jonas, R., Braiman, M., and Bormann, B.J. 1994. Structure and orientation of the transmembrane domain of glycophorin A in lipid bilayers. Biochemistry 33: 63346341.[CrossRef][Medline]
Smith, S.O.S., Eilers, M., Song, D., Crocker, E., Ying, W., Groesbeck, M., Metz, G., Ziliox, M., and Aimoto, S. 2002. Implications of threonine hydrogen bonding in the glycophorin A transmembrane helix dimer. Biophys. J. 82: 24762486.
Summa, C.M. 2002. "Computational methods and their applications for de novo functional protein design and membrane protein solubilization." Ph.D. thesis, University of Pennsylvania School of Medicine, Philadelphia.
Tamm, L.K. and Tatulian, S.A. 1997. Infrared spectroscopy of proteins and peptides in lipid bilayers. Q. Rev. Biophys. 30: 365429.[CrossRef][Medline]
Weiner, S.J., Kollman, P.A., Case, D.A., Singh, U.C., Ghio, C., Alagona, G., Profeta, S., and Weiner, P. 1984. A new force field for molecular mechanical simulation of nucleic acids and proteins. J. Am. Chem. Soc. 106: 765784.[CrossRef]
Ying, W., Irvine, S.E., Beekman, R.A., Siminovitch, D.J., and Smith, S.O. 2000. Deuterium NMR reveals helix packing interactions in phospholamban. J. Am. Chem. Soc. 122: 1112511128.[CrossRef]
Zhou, F.X., Cocco, M.J., Russ, W.P., Brunger, A.T., and Engelman, D.M. 2001. Polar residues drive association of polyleucine transmembrane helices. Proc. Natl. Acad. Sci. 98: 22502255.
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |