help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Originally published as Biophys J. BioFAST on May 27, 2005.
doi:10.1529/biophysj.105.059857
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.105.059857v1
89/2/1120    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 Google Scholar
Google Scholar
Right arrow Articles by Soubias, O.
Right arrow Articles by Réat, V.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Soubias, O.
Right arrow Articles by Réat, V.
Biophysical Journal 89:1120-1131 (2005)
© 2005 The Biophysical Society

Determination of the Orientation and Dynamics of Ergosterol in Model Membranes Using Uniform 13C Labeling and Dynamically Averaged 13C Chemical Shift Anisotropies as Experimental Restraints

O. Soubias *, F. Jolibois {dagger}, S. Massou {ddagger}, A. Milon * and V. Réat *

* Institut de Pharmacologie et de Biologie Structurale, UMR 5089, Toulouse, France; {dagger} Laboratoire de Physique Quantique, UMR 5626, Institut de Recherché sur les Systèmes Atomiques et Moléculaires Complexes, Université P. Sabatier, Toulouse, France; and {ddagger} Service de Résonance Magnétique Nucléaire, Fédération de Recherche Toulousaine en Chimie Moléculaire, Université P. Sabatier, Toulouse, France

Correspondence: Address reprint requests to Valerie Réat, Tel.: 33-5-6117-5418; Fax: 33-5-6117-5424; E-mail: valerie.reat{at}ipbs.fr.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 MATERIAL AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
A new strategy was established to determine the average orientation and dynamics of ergosterol in dimyristoylphosphatidylcholine model membranes. It is based on the analysis of chemical shift anisotropies (CSAs) averaged by the molecular dynamics. Static 13C CSA tensors were computed by quantum chemistry, using the gauge-including atomic-orbital approach within Hartree-Fock theory. Uniformly 13C-labeled ergosterol was purified from Pichia pastoris cells grown on labeled methanol. After reconstitution into dimyristoylphosphatidylcholine lipids, the complete 1H and 13C assignment of ergosterol's resonances was performed using a combination of magic-angle spinning two-dimensional experiments. Dynamically averaged CSAs were determined by standard side-band intensity analysis for isolated 13C resonances (C3 and ethylenic carbons) and by off-magic-angle spinning experiments for other carbons. A set of 18 constraints was thus obtained, from which the sterol's molecular order parameter and average orientation could be precisely defined. The validity of using computed CSAs in this strategy was verified on cholesterol model systems. This new method allowed us to quantify ergosterol's dynamics at three molar ratios: 16 mol % (Ld phase), 30 mol % (Lo phase), and 23 mol % (mixed phases). Contrary to cholesterol, ergosterol's molecular diffusion axis makes an important angle (14°) with the inertial axis of the rigid four-ring system.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
MATERIAL AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Cholesterol and ergosterol are well-known sterols, which play an essential role in the cell membrane of eukaryotic organisms. Cholesterol is the major sterol found in animal cells membrane and ergosterol is the predominant sterol in fungi and yeast; and they are absolutely required for viability and cell proliferation (1Go–4Go). One of the primary roles of sterols in eukaryotic cells is to modulate the physical properties of the plasma membrane phospholipid bilayer (2Go). For example, they reduce molecular surface area (5Go–9Go) and membrane permeability (2Go,10Go–12Go), alter lateral diffusion rates of proteins and lipids (13Go–17Go), broaden gel to liquid-crystalline phase transition (18Go–21Go) and modulate acyl-chain order in both the gel and liquid-crystalline phases (22Go–25Go). Many studies have been carried out to gain understanding of the molecular basis of cholesterol's biological properties (for excellent recent reviews on physical studies of cholesterol-phospholipid interactions, see Refs. 26Go and 27Go). In contrast, there are very few detailed studies on effects of ergosterol on physical properties of phospholipid bilayer membranes. Early monolayer studies (28Go) noted that ergosterol was much less effective than cholesterol in condensing egg-phosphatidylcholine (PC) monolayers, whereas Demel and De Kruyff (10Go) found that ergosterol was less effective than cholesterol in reducing the permeability of egg-PC liposomes to glycerol, glucose, and rubidium ions. Semer and Gelerinter (24Go), using electron-paramagnetic-resonance of doxyl-fatty acid and cholestane spin-label probes in egg-PC membrane, reported the observation that ergosterol, in contrast with cholesterol, only ordered the acyl chain up to 15–20 mol %; above this percentage, ergosterol induced disorder of the acyl-chains. Ergosterol has a very low solubility (<3%) into soybean PC (containing 63% of di-unsaturated fatty acid) compared to cholesterol, which can be incorporated up to 23% (29Go). Urbina et al. (25Go) found that the ordering effects are slightly higher for ergosterol than for cholesterol, in sterol/dimyristoylphosphatidylcholine (DMPC) membrane at 30 mol % and at 25°C. More recently, Hsueh et al. (30Go) have determined the complete phase diagram of ergosterol/dipalmitoylphosphatidylcholine mixtures by deuterium nuclear-magnetic-resonance (NMR) and calorimetry.

Three main features of sterols have been linked to their characteristic effects on lipid bilayer membrane: a planar four-fused-ring motif, a 3ß-OH group, and a hydrophobic side chain linked to C17. However, the mechanisms by which the slight differences in structures between cholesterol and ergosterol may be responsible for modification in the physical properties of the bilayer and in the sterol's dynamics are still largely unknown. Although the average orientation and dynamics of cholesterol have been well established by deuterium NMR (see Ref. 22Go and references cited herein), very little is known concerning ergosterol on this aspect.

Solid-state NMR is uniquely suited to the determination of orientation, dynamics, and molecular structure of membrane compounds. Combination of solid NMR methods and stable isotope 15N labeling of the amide nitrogens (31Go–35Go) or 13C labeling of the amide carbonyls (36Go,37Go), are currently used to obtain orientational and structural constraints for membrane peptides and proteins. In the case of sterol, it has been shown that cholesterol order parameters can be extracted from the quadrupole splittings of specifically deuterated positions (22Go). This strategy has proven to be very powerful but requires critical steps such as chemical synthesis of specifically deuterated positions and assignment of the obtained quadrupolar splittings to be performed. 13C NMR could be an interesting alternative to study such compounds due to the large carbon chemical shift range and the sensitivity enhancement brought about by cross-polarization (CP) from neighboring protons. Moreover, modulation of biosynthetic labeling strategies can be used to increase both sensitivity and selectivity. In particular, methylotrophic yeast has been shown to be an excellent tool for the cost-efficient production of uniformly labeled biomass and membrane proteins (38Go–40Go). Another advantage for NMR is that this yeast can be easily grown on a minimal medium with methanol, i.e., a C1 compound as a sole carbon source. Randomly distributed 13C labeling to any level is therefore possible, and in particular low level 13C labeling (11%) is suitable to measure carbon chemical shifts or H–C dipolar couplings on isolated spin pairs with a reasonable sensitivity.

The 13C chemical shift anisotropy (CSA) is dependent on both molecular structure and dynamics and therefore represents a powerful probe of these two parameters. Furthermore, recent developments and implementation in quantum chemistry allow the accurate determination of theoretical chemical shift tensors (for reviews, see Refs. 41Go and 42Go).

In this article, by combining quantum chemistry calculations with various magic-angle spinning (MAS) solid-state NMR experiments, carbon chemical shift anisotropies were used as restraints to probe orientation and dynamics of uniformly 13C-labeled ergosterol in DMPC membrane. The complete proton and carbon assignment of ergosterol in membrane has been achieved, based on a combination of one-dimensional and two-dimensional 13C–1H heteronuclear (J- and dipolar-heteronuclear correlation, i.e., HETCOR) and 13C–13C homonuclear (incredible-natural-abundance-double-quantum-transfer experiment; i.e., INADEQUATE) MAS experiments. Using this strategy we could determine ergosterol's order parameter and diffusion axis in the two major coexisting phases in ergosterol/DMPC mixtures, i.e., Ld and Lo phases. We analyzed at 313 K three different ergosterol/DMPC molar ratio; 16 mol % ergosterol (pure Ld phase), 30 mol % (pure Lo phase), and 23 mol % (two-phase coexistence) according to Hsueh et al. (30Go).


   MATERIAL AND METHODS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Growth condition and culture media
Fermentation of Pichia pastoris methylotrophic yeast was used to produce ergosterol. The inoculum for fermentation was obtained by growing this microorganism in a liquid medium containing 2% (v/v) methanol, 3 g/l KH2PO4, 0.2 g/l MgSO4 7H2O, 1.5 g/l (NH4)2SO4, 2 ml/l of a trace element solution (50 g/l EDTA, 6 g/l ZnSO4 7H2O, 5.4 g/l CuCl2 2H2O, 0.39 g/l MnSO4 H2O, 4 g/l FeSO4 7H2O, 0.06 g/l (NH4)2 Mo2O7, 1.2 g/l CuSO4 5H2O, 0.8 g/l CoCl2 6H2O) and 1 ml/l vitamin solution (5 mg/l biotin and 300 mg/l thiamine chlorohydrate). The methanol and the vitamin solutions were sterilized by filtration. The pH was set to 6 with KOH 3N. Shake flask culture (200 ml) was performed at 37°C until OD640 reached 1. The production of large amount of 11% randomly labeled ergosterol was carried out in a 2.5-l fermentor. The initial fermentation culture media contained 10 g/l KH2PO4, 1 g/l MgSO4, 10 g/l (NH4)2SO4, 12.5 ml/l of trace element solution, 1 ml/l of vitamin solution, and 0.5% (v/v) of methanol solution (10% of 13C methanol and 90% of natural abundance methanol). The pH was regulated at 6.0 by addition of KOH 3N. The temperature was set to 37°C. During the fed batch phases, the flow rate of methanol was increased from 3 g/h to 30 g/h. At the end of fermentation, the methanol total consumption was 300 g. When the OD640 was 50, 50g of (NH4)2SO4 were added. The O2 pressure was regulated by increasing airflow and stirring speed. To maintain the dissolved O2 level above 25%, 40% oxygen-enriched air was used when the OD640 reached 70. At the end of fermentation, the cells were harvested by centrifugation and lyophilized.

Purification from Pichia pastoris yeast cells, and chemical characterization
Pichia pastoris yeast cells (10 g) were broken by freeze/thaw cycles and delipidated in methanol/chloroform/acetone solvent mixture (1:1:2 v/v/v) during 4 h under reflux and argon atmosphere. Polar and neutral lipids were separated by precipitation in acetone (2 days, 4°C). The mixture was then centrifuged at 4000g for 5 min and the supernatant saponified in a 1:1:4 potassium hydroxide/water/ethanol mixture. Sterols were extracted with hexane and further purified using silicic column chromatography (20 g silica-gel 60; Merck, Whitehouse Station, NJ).

Before the gas chromatography/mass spectrometry (GC/MS) experiment, ergosterol was acetylated in a 1:1 anhydrous pyridine/acetic anhydride mixture. GC/MS spectra were acquired on an HP 5989X glass column and interfaced to a Hewlett-Packard 5989A mass spectrometer (Palo Alto, CA). Acetylated ergosterol derivatives were analyzed by electron-impact/mass spectroscopy (ion source temperature 250°C, ionization voltage 70 eV). Labeling extent (10.7%) was calculated from GC/MS spectra on molecular ion (m/z = 441) and deacetylated fragment (m/z = 381). From GC/MS spectra, another minor sterol was identified (ergosta-5-7-22-24-tetraen-3ß-ol, m/z = 379), which represented 4.8% of the total sterol content.

Sample preparation
Samples were prepared by drying down under nitrogen chloroform solution containing perdeuterated acyl chains dimyristoylphosphatidylcholine (DMPC-d54; Avanti Polar Lipids, Alabaster, AL) and purified ergosterol in the desired proportions (8.4:1.6, 7.7:2.3, and 7.0:3.0 mol/mol), then eliminating residual solvent under high vacuum overnight. The lipids were resuspended using deuterium-depleted water at 33 wt % solids, then subjected to five mixing and freeze (253 K) thaw (323 K) cycles, to ensure homogeneity of the multilamellar vesicles as described in Urbina et al. (25Go). Before NMR spectroscopy, sample hydration was adjusted at 50% weight by adding the appropriate amount of deuterium-depleted water. It was checked by deuterium NMR that the lipid chain's quadrupolar splittings had values in good accordance with those reported previously in the literature at equivalent sterol/lipid molar ratio (25Go). Spectra similar to those of Hsueh et al. (30Go) were obtained confirming the presence of Ld phase at 16 mol % ergosterol, Lo phase at 30 mol %, and Ld + Lo phase mixture at 23 mol % at 313 K (data not shown).

Solid-state NMR spectroscopy
NMR experiments were carried out on a Bruker Avance 500 spectrometer (Bruker Optics, Billerica, MA) equipped with a narrow-bore magnet operating at a resonance frequency of 125 MHz for 13C and 500.13 MHz for 1H. All solid-state NMR experiments were carried out at 313 K.

Static deuterium NMR spectra were recorded on a Bruker 7-mm double-resonance probe, with a solenoid coil oriented at 90° with respect to the magnetic field. The acyl-chain quadrupolar splittings of lipids were recorded by using a standard quadrupolar echo sequence. Deuterium {pi}/2 pulses were equal to 4 µs; and the refocusing delay and the repetition time were set to 30 µs and 1 s, respectively.

MAS experiments were recorded with a double-tuned Doty XC5-MAS probe (Doty Scientific, Columbia, SC) equipped with a 5-mm spinning module. The 1H radio-frequency field-strength for heteronuclear two-pulse phase-modulation decoupling was 66 kHz for all experiments. All spectra were acquired at MAS spinning rate of 9 kHz with a repetition delay of 3.5 s to avoid sample heating.

Cross-polarization magic-angle spinning (CP-MAS) and cross-polarization off-magic-angle spinning (CP-OMAS) spectra were acquired using a 1H excitation pulse length of 4 µs and a CP spin-lock field-strength of ~50 kHz. The CP contact time was 2 ms. The relaxation delay was set to 2 s for all experiments.

Two-dimensional non-refocused INADEQUATE (43Go) was obtained by recording 100 t1 increments with 320 scans each. The {tau}-delay was set equal to 0.86 ms. Contact time for the CP period and 1H radio-frequency field were identical to one-dimensional CP spectra.

Dipolar HETCOR spectra were obtained at a 9-kHz spinning rate. CP contact times were set to 2.5 ms or 250 µs. A total of 128 t1 increments with 50-µs dwell-time and 128 scans each, were recorded for each experiment.

Quantum chemistry calculation
The three ß-hydroxyl rotameric ergosterol structures (gauche–, gauche+, and anti–; see Fig. 1) were optimized at the Hartree-Fock level using the STO-3G basis set. The carbon-shielding tensor was then computed by using the gauge-including-atomic-orbital approach within Hartree-Fock theory (6-31G(d,p) basis set). All calculations have been performed using the Gaussian 98 package (44Go). It has been previously shown that this computational strategy is sufficient to reproduce the experimental 1H and 13C isotropic NMR spectra of ergosterol in chloroform solution (45Go). The comparison of theoretical chemical shifts with those experimentally obtained has unambiguously revealed that one must take the theoretical average values among three isoenergetic ß-hydroxyl rotamers to reproduce the 1H and 13C NMR spectra of ergosterol. Averaged carbon chemical shift tensors are thus calculated in a similar way. Diagonalization of the symmetric part of this tensor provides the principal values ({delta}11, {delta}22, and {delta}33) and their associated principal vectors {1}. The principal values are ordered as where the isotropic chemical shift is equal to The static chemical shift anisotropy and the asymmetry parameter are defined as {Delta}{delta}static = {delta}33 – {delta}iso and respectively. The principal vectors define the orientation of the chemical shift tensor principal axis system (denoted PASCS) in the molecular frame (denoted M and defined as axis along the {6} bond, along {7}, and along ). The molecular frame M is independent of the rotameric states and was thus preferred to the C3 molecular frame defined in Marsan et al. (22Go). The ( and ) Euler angles characterize the rotation from PASCS to M (Fig. 2).



View larger version (12K):
[in this window]
[in a new window]
  		FIGURE 1  Newman representation of the three rotational isomers around the C3–O bond. These isomers are discriminated according to the C2–C3–O–H dihedral angle.

 


View larger version (25K):
[in this window]
[in a new window]
  		FIGURE 2  Schematic representation of frame transformation and Euler angles definition. The values are the Euler angles associated to the rotation matrix (RPM) between the chemical shift tensor principal axis system PASCS {1} and the molecular frame M (). Columns of the RPM matrix correspond to chemical shift tensor principal axis {1} defined in the molecular frame. Note that ({alpha}, ß, and 0) are the Euler angles associated to the rotation matrix (RMN) between M and the reorientation frame (N). The values {alpha} and ß correspond to polar coordinates of the ergosterol diffusion axis () in the molecular frame. The values () are the Euler angles corresponding to the rotation from PASCS to N. The values () correspond to the polar coordinates of the ergosterol diffusion axis in PASCS. These angles are extracted from the rotation matrix RPN calculated as the product of RPM and RMN.

 
Static chemical shift tensor parameters are listed in Table 1 for ergosterol carbons 1–19 (see ergosterol structure in Fig. 3).


View this table:
[in this window]
[in a new window]
  		TABLE 1  Theoretical static chemical shift tensor parameters for ergosterol ring carbons

 


View larger version (10K):
[in this window]
[in a new window]
  		FIGURE 3  Structure of ergosterol. Carbons atoms are labeled with a number; hydrogen atoms, with the exception of the hydrogen of the hydroxyl group, are not shown.

 
Simulation of dynamically averaged chemical shift anisotropy
The carbon chemical shift value is given, in the case of fast axial diffusion along the normal bilayer by (46Go)

(1)
where ßlayer is the angle between the bilayer normal and the B0 magnetic field axis, Smol the molecular order parameter (order parameter describing the amplitude of wobbling of the entire rigid moiety), and Sloc is a local order parameter to account for internal motions. It can be set equal to 1 for every ring atom as was shown previously (22Go). The values {Delta}{delta}static and {eta} represent the static chemical shift anisotropy and the asymmetry parameters, respectively. The values ( and ) are the Euler angles associated to the rotation between the chemical shift tensor principal axis system (PASCS) and the diffusion frame (denoted N). The z' axis of the N frame is along the ergosterol axis of diffusion (). The values {alpha}, ß are the polar coordinates of in the molecular frame M. The angles ( and ) are then depending on both ( and ), i.e., rotation from PASCS to M, and ({alpha}, ß), i.e., orientation of N in M (see Fig. 2).

In this article, theoretical carbon chemical shift anisotropy {Delta}{delta}calc in the presence of fast axial diffusion is defined as

(2)
where {delta}90° and {delta}0° are the chemical shifts calculated for ßlayer = 90° and ßlayer = 0°, following Eq. 1.

Using static chemical shift tensor parameters obtained by quantum chemical computations, sets of {Delta}{delta}calc values have been calculated for each ergosterol ring carbon (Sloc = 1) by varying {alpha}, ß and Smol from 0° to 360°, 0° to 90°, and 0 to 1, respectively. The increments were 1° for angles and 0.1 for Smol. The obtained data were filtered according to the root mean-square difference (RMSD) calculated as

(3)
with N{sigma} the number of experimental CSA data (N{sigma} = 18), and and as the theoretical and experimental CSAs for the ith carbon, respectively. The values of were determined as described in the next section.


   RESULTS AND DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
 RESULTS AND DISCUSSION
CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
Production of uniformly 10% and 100% 13C-labeled ergosterol
Methylotrophic yeast proved to be an excellent mean of producing randomly 13C-labeled ergosterol at any extent. Since it can be grown on methanol as a sole carbon source, the extent of labeling can be chosen simply by mixing adequate volumes of labeled and natural abundance methanol. Very high biomass (78 g) could be obtained in a few days, using 300 g of methanol. The purification of ergosterol from total lipids is straightforward, and it was found to be the major sterol, with a small contamination by a closely related sterol possessing one more double bond on its side chain (in position 24). Approximately 1 mg/g of cells (dry weight) could be purified. Two cultures were performed to obtain ergosterol with 10.7% and 99.5% 13C labeling extent.

Complete assignment of ergosterol in DMPC membrane
We have used a combination of one-dimensional and two-dimensional heteronuclear and homonuclear magic-angle spinning (MAS) solid-state NMR techniques—namely, scalar HETCOR, dipolar HETCOR, and INADEQUATE—to perform the complete assignment of ergosterol in membrane without any reference to the known liquid-state assignment (Table 2; structure of ergosterol in Fig. 3). The samples were 8.4:1.6 mol/mol DMPC-d54/ergosterol–11%13C multilamellar vesicles except for the INADEQUATE experiment for which ergosterol, 100% 13C-labeled, was used. Apart from the INADEQUATE experiment, the assignment strategy was identical to that recently described in the case of natural-abundance cholesterol in DMPC membrane (47Go). Two-dimensional pulse sequences used in this work are shown in Fig. 4. Carbon multiplicities have been determined by comparing one-dimensional carbon spectra recorded either via a CP scheme or via an insensitive-nuclei-enhanced-by-polarization-transfer (INEPT) scheme by varying the second delay of the refocusing period (Fig. 5). 1H–13C one-bond correlations were obtained via a liquidlike HETCOR experiment by using a refocused INEPT scheme. Skeletal information has been recovered by combining homonuclear 13C–13C correlations from the INADEQUATE experiment and long-range 1H-13C correlation from the dipolar HETCOR experiment. Fig. 6 shows the INADEQUATE spectrum of 8.4:1.6 perdeuterated acyl-chain DMPC/100% 13C-labeled ergosterol recorded in a total experimental time of 16 h. In the two-dimensional map, two directly bonded carbons share a common frequency in the double-quantum dimension. Starting from the easily identifiable resonance of the hydroxyl carbon C3 at 69.5 ppm, it is straightforward to sequentially assign the complete carbon spectrum. Two sets of correlation peaks are possible for the C3 connectivities (C3–C2, C3–C4), one at 31.5 ppm and another at 39.65 ppm. On the basis of chemical shift considerations, the upfield carbon signal at 39.65 ppm is identified as the methylene C4 carbon. Another way is to remark that each of them has a second clear correlation in the double-quantum frequency dimension with either a methylene for C2 (hence C1 at 38.32 ppm) or a quaternary C-signal (hence C5 at 141.18 ppm, data not shown) for the C4. Interestingly, our assignment differs in one point from what was obtained in the liquid state. In deuterated chloroform solution, carbon C8 was found to resonate at a higher frequency than carbon C5. This inversion was confirmed by analyzing the olefinic region of the dipolar HETCOR experiment (Fig. 7). This discrepancy underlines the importance of assigning compounds in their actual environment. These differences are currently under investigation by analyzing the environment effect (hydrogen bonding and rotameric state of the hydroxyl group) by combining solid-state NMR assignments with quantum chemistry calculations (47Go).


View this table:
[in this window]
[in a new window]
  		TABLE 2  Assignment of carbon and proton chemical shifts for ergosterol inserted in a 8.4:1.6 (mol/mol) DMPC-d54/ergosterol lipid mixture

 


View larger version (17K):
[in this window]
[in a new window]
  		FIGURE 4  Pulse sequences used in this work. (A) Dipolar HETCOR. During the evolution period, magnetization evolves under the 1H chemical shift and the residual proton dipolar couplings. 1H to 13C coherence transfer was then mediated via H–C dipolar couplings. Carbon chemical shift is probed during t2. (B) The CP transfer scheme is replaced by a refocused INEPT scheme and the transfer is mediated via 1J isotropic scalar couplings in the case of scalar HETCOR. (C) Two-dimensional INADEQUATE pulse sequence. Apart from the CP period, the experiment is completely analogous to liquid-state versions. After cross-polarization from 1H, the 13C magnetization evolves during the delay 2{tau} under the isotropic homonuclear J coupling Hamiltonian. The C–C homonuclear dipolar couplings are removed by fast MAS, and the chemical shift evolution is refocused by the 180° pulse. The double-quantum coherence created by the first carbon 90° pulse evolves during t1 at the sum frequency of the two spins, and is converted back into an antiphase transverse coherence by the last 90° pulse, which is detected during t2. (D) Lee-Golburg experiment was used to determine proton-carbon dipolar couplings and permit proton stereospecific assignments of methylene protons.

 


View larger version (30K):
[in this window]
[in a new window]
  		FIGURE 5  Identification of carbon multiplicity. Comparison of one-dimensional carbon MAS spectra obtained with a CP transfer scheme (a) and a refocused INEPT scheme (b) allows us to identify missing resonances as quaternary carbon signals (denoted Q). Comparison of the refocused INEPT spectrum and the spectral edited refocused INEPT (c) shows the expected negative, but also the missing carbon resonances, which can be explained by fast transverse relaxation of the tightly coupled methylene protons of the ergosterol rings during the overall INEPT delay. Hence, negative and missing resonances can be assigned to methylene carbons. Lipids resonances are marked with an asterisk (*). These spectra were recorded on 8.4:1.6 DMPC-d54/ergosterol–11%13C multilamellar vesicles at 313 K.

 


View larger version (26K):
[in this window]
[in a new window]
  		FIGURE 6  Assignment of the ring A of the ergosterol molecule. As in solution-state NMR, directly bonded carbon-13 resonances are identified by the fact that they generate a common double-quantum frequency. The connectivities between antiphase line shapes are indicated by solid lines. The spectrum was obtained in 18 h using 1.5 mg of 100% 13C-labeled ergosterol in DMPC-d54 membrane at 16% molar ratio at 313 K.

 


View larger version (14K):
[in this window]
[in a new window]
  		FIGURE 7  Comparison of 1H–13C chemical shifts correlation extracted from scalar HETCOR (A) and dipolar HETCOR (B) centered onto the olefinic region of ergosterol. From the scalar HETCOR, the proton chemical shift directly bonded to carbon can be easily identified (C6–H6 119.05 ppm, 5.50 ppm, C7–H7 117.05 ppm, 5.27 ppm). The dipolar HETCOR carbon slice extracted at 141.18 ppm shows three cross-peaks: at 2.25, 2.49, and 5.50 ppm—with two correlations corresponding to long-range transfer from H4 protons and H6 proton, respectively. The resonance at 141.18 ppm can be thus identified as the C5 of ergosterol. Signal at 139.2 ppm corresponds then to C8 with a long-range proton transfer from H7 (5.27 ppm) and H9 (1.91 ppm). These spectra were recorded on 8.4:1.6 DMPC-d54/ergosterol–11%13C multilamellar vesicles at 313 K.

 
The ergosterol carbon assignments for the other two lipid sterol mixtures (7.7:2.3 and 7.0:3.0 mol/mol DMPC-d54/ergosterol–11%13C) show differences of <0.2 ppm (data not shown).

Orientation and dynamics of ergosterol in DMPC membrane: chemical shift anisotropy restraints
The determination of sterol's orientation and dynamics in lipid bilayers can be performed using several strategies. One can make use of the residual second rank tensor interactions, which are due to the incomplete motional averaging, by measuring 2H quadrupolar splittings (Ref. 22Go and references cited herein), 1H–13C dipolar couplings (40Go), or chemical shift anisotropies on static liposomes (which give powder spectra) or mechanically oriented bilayer samples (48Go). One can also use MAS techniques to average out these tensor interactions in one dimension and selectively reintroduce them in another dimension (49Go–52Go). The present work was focused on the analysis of 13C CSA residual interactions on liposome samples, under MAS and OMAS conditions.

Critical chemical shift parameters
Molecular reorientations, on a timescale shorter than 10–5 s, will result in a reduction of the chemical shift anisotropy (CSA). The carbon chemical shift tensor can be characterized by six parameters: namely, three principal values ({delta}11, {delta}22, and {delta}33) and three Euler angles specifying the orientation of the principal axis system {1} in the molecular frame. Fig. 8 shows simulated spectra for the cholesterol's carbon C4 resonance in the presence of fast axial diffusion around three axes. Fig. 8 a is for a static sample and Fig. 8, b–d, show the simulated line shape when the diffusion axis is parallel to {2} respectively. One can see that the chemical shift anisotropy is reduced in the presence of fast axial diffusion and depends on the orientation of the diffusion axis. Hence, given the principal values and orientation of at least three static chemical shift tensors within a rigid molecule, one can obtain the diffusion axis orientation and the molecular order parameter, by comparing the experimental motionally averaged anisotropy ({Delta}{delta}exp) to the theoretical one ({Delta}{delta}calc).



View larger version (14K):
[in this window]
[in a new window]
  		FIGURE 8  Simulated 13C spectra of cholesterol–13C4, showing the effect of rapid axially symmetric reorientation on the powder pattern line shape: (a) the static case, with principal values and asymmetry parameter extracted from quantum chemistry calculation ({delta}11 = 58.3 ppm, {delta}22 = 40.8 ppm, {delta}33 = 20.5 ppm, {Delta}{delta}static = –19.4 ppm, and {eta} = 0.9), and rapid axial reorientation about an axis () oriented (b) parallel to {3} ({Delta}{delta}calc = –19.4 ppm), (c) parallel to {4} ({Delta}{delta}calc = 18.4 ppm), and (d) parallel to {5} ({Delta}{delta}calc = 0.9 ppm).

 
Various experimental methods have been designed to estimate the principal values of the chemical shift tensor (see below). By contrast, the experimental determination of the chemical shift tensor's orientations is more difficult since it requires either single crystals (Ref. 53Go and references therein) or chemical synthesis of specifically doubly labeled compounds, which allows us to orient the CSA tensor with respect to a dipolar tensor of known orientation (54Go–59Go). Fortunately, recent developments in quantum chemistry permit the accurate treatment of theoretical chemical shielding (for reviews see Refs. 41Go, 42Go, and 60Go–62Go). In this context, we have recently shown that isotropic carbon chemical shifts of cholesterol and ergosterol can be reproduced accurately, even though these calculations were performed for an isolated molecule, neglecting inner molecular vibrations and using fairly simple basis sets (45Go). However, isotropic chemical shift calculation only requires accurate principal values and does not depend on tensor orientations. To check the quality of eigen-values and eigen-vectors of the theoretical tensors we compared calculated and experimental motionally averaged chemical shift anisotropies of cholesterol's carbons C3 and C4, for which dynamical parameters in a DMPC membrane (7.0:3.0 mol/mol DMPC/cholesterol) have been previously determined (22Go). The carbon theoretical chemical shift tensors have been calculated following the same procedure as for ergosterol (Table 3). As shown in Table 3, there is a reasonably good agreement for each of the two carbon-averaged anisotropies showing that both the principal values and the orientation of the chemical tensors are efficiently provided by quantum chemistry calculations, and that the approximations associated with this approach lead to negligible effects in sterol's orientational analysis. It may be worth noting that the errors on the predicted {Delta}{delta} (1.3 and 1.6 ppm for carbons C3 and C4, respectively) are comparable to the average error in isotropic chemical shift calculations—at ~2 ppm, depending on the computation method and molecule, in Jolibois et al. (45Go). Therefore, orientation effects do not introduce a significantly larger error in chemical shift computation than principal values inaccuracies.


View this table:
[in this window]
[in a new window]
  		TABLE 3  Theoretical ({Delta}{delta}calc) and experimental ({Delta}{delta}exp) motionally averaged CSA of C3 and C4 carbons of cholesterol

 
Motionally averaged chemical shift anisotropy measurements
The chemical shift anisotropy is a powerful probe of molecular structure and dynamics, and its measurement has lead to the design of various one-dimensional and two-dimensional static or MAS NMR methods. As shown above for assignment experiments, under fast rotation at the magic angle, all of the second rank interactions are effectively averaged to zero and high-resolution isotropic spectra can be obtained with high sensitivity and resolution in the case of magnetically dilute spins. In a slower spinning regime, when the spinning frequency {nu}R is comparable to the spread in Larmor frequencies caused by the CSA, the NMR spectrum contains peaks at frequency coordinates viso and (viso+ k vr), where viso is the isotropic shift frequency of site, k an integer (called the side-band order), and vr the spinning frequency. The intensity distribution of the spinning side-bands is characteristic of the CSA tensor, and it may be analyzed to obtain the CSA parameters. Carbon 13 NMR spectra of DMPC-d54/ergosterol–11%13C multilamellar vesicles at three different ergosterol molar ratio (16Go, 23Go, and 30 mol %) were recorded under slow MAS conditions, at 2, 2.5, and 3 kHz. For each spectra, the intensities of the observed side-bands have been used in a Hertzfeld-Berger analysis to extract the components of the motionally averaged chemical shift tensor ({Delta}{delta}exp) associated to five carbon sites of the ergosterol molecule, namely C3 and C5–8, which are characterized by large anisotropies and well-resolved spinning side-bands (Table 4).


View this table:
[in this window]
[in a new window]
  		TABLE 4  The motionally averaged chemical shift tensor ({Delta}{delta}exp) for ergosterol ring carbons of the three lipid mixtures at 313 K

 
In the upfield region of the ergosterol carbon spectrum, spectral overlap, small CSAs, and perturbations from dipolar interactions prevented the extraction of reliable CSA parameters from side-band intensity analyses. To overcome this problem a series of two-dimensional experiments have been proposed, at low (Ref. 63Go and references therein; and Refs. 64Go–66Go) and fast MAS rate (52Go,67Go,68Go), the common principle of these experiments being to use a second dimension to separate overlapping resonances according to their isotropic chemical shifts or side-band order. These experiments are widely used to determine CSAs for numerous nuclei, but they become technically difficult to implement for small CSAs (except for Ref. 68Go, where this point was addressed).

In this article, we have used a one-dimensional approach closely related to the idea of the switched-angle-sample-spinning (49Go,51Go) experiment: in a spectrum arising from fast-rotation OMAS, the CSA is scaled by a reduction factor (which is as long as ßRL, the angle between the rotation axis and the static magnetic field, is kept close to magic angle; see Ref. 69Go). Two of the one-dimensional OMAS spectra at two different angles were then recorded to follow the deviations between {delta}iso of the MAS spectrum and {delta}90° of the OMAS spectra (Fig. 9). Chemical shifts were extracted after calibration on the choline polar-head methyl's resonance, which is known to exhibit a small anisotropy (~2 ppm) (48Go). OMAS scaling factors were calibrated by 31P NMR using the CSA of the DMPC phosphate. Table 4 summarizes the different values of carbon CSA with their uncertainty obtained for ergosterol signals extracted either from the slow MAS or from the OMAS experiments.



View larger version (24K):
[in this window]
[in a new window]
  		FIGURE 9  Ergosterol–11%13C CP-MAS (a), and two CP-OMAS (b and c) spectra in DMPC-d54 multilamellar vesicles at 313 K (8.4:1.6 mol/mol DMPC/ergosterol). The OMAS angles were ß90°,1 = 51.5° and ß90°,2 = 49.5°. The motionally averaged CSAs ({Delta}{delta}exp) for ergosterol ring carbons are calculated using (see Eq. 2 for {Delta}{delta}exp, {delta}iso, and {delta}90° definitions). The values {delta}iso, {delta}90°,1, and {delta}90°,2 are extracted for each carbon from the MAS (a), and the two OMAS (b and c) spectra, respectively. The spectra were calibrated on the choline polar-head methyl resonance, which has a very small CSA (<2 ppm; see Ref. 48Go). For each OMAS experiment, the scaling factor is determined from the 31P NMR spectra acquired on the same sample and same condition, and is equal to The experimental error on {Delta}{delta}exp is estimated at ±3 ppm for each carbon and takes into account the accuracies of {delta}90°,1, {delta}90°,2 measurements in 31P and 13C spectra. Following this procedure we could calculate, for example, {Delta}{delta}exp (C4) = 7.5 ± 3 ppm and {Delta}{delta}exp (C9) = 0 ± 3 ppm.

 
Data analysis
The data analysis was performed in a similar way for three DMPC/ergosterol mixtures (8.4:1.6, 7.7:2.3, and 7.0:3.0 mol/mol). Using static chemical shift tensor parameters obtained by quantum chemical computations, sets of theoretical, motionally averaged, carbon chemical shift anisotropies ({Delta}{delta}calc, Eq. 2) have been calculated by varying {alpha}, ß, and Smol from 0° to 360°, 0° to 90°, and 0 to 1, respectively. For each ({alpha}, ß, Smol) triplet, the 18 carbon's {Delta}{delta}calc values were compared to the experimental ones and their RMSD was evaluated following Eq. 3. Over 32.4 106 possibilities, 386, 199, and 200 triplets correspond to an RMSD smaller than 3.3 ppm (upper limit of the experimental errors) for 8.4:1.6, 7.7:2.3, and 7.0:3.0 (mol/mol) DMPC/ergosterol samples, respectively. For each mixture, these solutions defined one single region in {alpha} and ß, for which the average and the standard deviation values are given in Table 5. The linear regression between experimental anisotropies and theoretical ones obtained with these averaged values shows a good agreement, giving a Pearson coefficient equal to 0.99 and a slope of 1.00 for each of lipid mixture (data not shown).


View this table:
[in this window]
[in a new window]
  		TABLE 5  Diffusion parameters for ergosterol or cholesterol in DMPC bilayers at 313 K

 
The 18 experimental motionally averaged CSA restraints allow us to define orientational and dynamical ergosterol parameters with the same precision as one obtained on cholesterol using quadrupole splittings (the high number of CSA restraints compensating for their lower accuracy, as compared with quadrupole splittings). One important advantage of the CSA strategy is that it does not require chemical synthesis of specifically labeled sterols. One can note that the diffusion axis orientation is different for cholesterol and ergosterol mixtures but is independent of the DMPC/sterol molar ratio. The diffusion axis is indeed approximately parallel to the principal inertia axis (defined by sterol tetra-cycle atoms) in the case of cholesterol ({gamma}I {approx} 2°), whereas it is tilted by ~14° for the ergosterol mixture (Table 5). Various explanations can be proposed to understand this difference. The effective inertia axis of the whole molecule can be more tilted from orientation for ergosterol than for cholesterol, due to the differences between the two sterols' side-chain structures and conformations. Alternatively, specific interactions between adjacent molecules (water or DMPC) and the sterol's hydroxyl group and/or the {pi} system of the second ring could also be different for cholesterol and ergosterol and influence the relative orientation of and (45Go,47Go). Molecular dynamics simulations performed on a 8:1 DMPC/sterol (mol/mol) bilayer system have also suggested some difference in ergosterol and cholesterol orientation (70Go). In these simulations, the average tilt of the sterol molecules, defined as the angle between bilayer normal and the vector connecting carbon atoms C3 and C17, was found slightly higher for ergosterol than for cholesterol.

The two sterols also present significantly different molecular order parameters (Table 5), depending on DMPC/sterol ratio. At low concentration (16 mol %, Ld phase) ergosterol is highly mobile (Smol = 0.76 ± 0.02), and at high concentration (30 mol %, Lo phase) the dynamics of ergosterol and cholesterol are similar (Smol = 0.89 ± 0.02 compared to 0.94 ± 0.01). At intermediate concentration (23 mol %, coexistence of the Ld and Lo phases), ergosterol's molecular order parameter is identical to its value in the Lo phase.

Finally, our experiment gives an estimate for the upper limit of domain size at an intermediate molar ratio. The difference of isotopic chemical shift of ergosterol resonances in the Ld phase (16 mol % ergosterol) and the Lo phase (30 mol % ergosterol) is comprised between 15 and 20 Hz depending on carbon atom. Since for 23 mol % ergosterol sample (coexistence of the Ld and Lo phases) one single sharp resonance is observed for each carbon (fast averaging limit on the chemical-shift scale), the exchange rate of ergosterol between the Ld and Lo domains has to be higher than 102 Hz. Following the equation <{Delta}x2>1/2 = (4D{tau})1/2 (30Go) and using the same diffusion coefficient for phospholipid and sterol, 10–11 m2/s (15Go), we thus determine an upper limit for a domain size of 600 nm (taking {tau} = 10–2 s). This information is complementary to the lower limit of 20 nm obtained from deuterium NMR spectra for which slow exchange was observed on the deuterium-NMR timescale (30Go).


   CONCLUSIONS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
ACKNOWLEDGEMENTS
 REFERENCES
 
The chemical shift anisotropy is routinely used to determine membrane-peptide's orientation (71Go,72Go). It is also a powerful probe of molecular dynamics since molecular reorientations on a timescale shorter than 10–5 s reduce its amplitude. 31P CSA has been used to study phospholipid headgroup's dynamics (73Go) or to reveal specific protein-lipid interactions (74Go,75Go). Fast axial diffusion of transmembrane peptides has been demonstrated using 15N CSA (34Go) and 13C CSA (76Go) of peptide bonds. Nevertheless, to our knowledge the present work is the first really quantitative assessment of a membrane component's orientation and dynamics using theoretical static 13C CSA tensors and experimental motionally averaged values. For this study, we have first produced uniformly 100% and 11% 13C-labeled ergosterol from Pichia pastoris methylotrophically and performed complete proton and carbon assignment of ergosterol in membrane through a combination of one- and two-dimensional 13C–1H heteronuclear (J- and dipolar-HETCOR) and 13C–13C homonuclear (INADEQUATE) MAS experiments. Motionally averaged chemical shift anisotropies for all ring carbons have then been determined by spinning side-band intensity analyses and OMAS experiments. The experimental CSAs have finally been combined with theoretical static CSA tensors to give access to the orientation and molecular order parameters of the diffusion axis. These parameters have been obtained for ergosterol at three different DMPC/ergosterol molar ratios and compared to data obtained previously on cholesterol by deuterium NMR on oriented bilayers (22Go). It has been shown that, in DMPC bilayers, the inertial axis of ergosterol's four-ring system is tilted by ~14° compared to the diffusion axis. Furthermore, ergosterol was shown to be highly mobile in the Ld phase (16 mol % ergosterol) and to present a dynamics similar to cholesterol in the Lo phase (30 mol % sterol).


   ACKNOWLEDGEMENTS
TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
REFERENCES
 
We acknowledge Dr. N. Lindley, Institut National des Sciences Appliquées Toulouse, for his help with P. pastoris cultures.

The nuclear magnetic resonance spectrometers were financed by Centre National de la Recherche Scientifique (Ingénierie des MAcromolécules BIOlogiques program), the Région Midi-Pyrénées, and European Structural funds.

Submitted on February 7, 2005; accepted for publication May 19, 2005.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 MATERIAL AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 ACKNOWLEDGEMENTS
 REFERENCES
 
1. Dahl, C., and J. Dahl. 1988. Biology of Cholesterol. P. L. Yeagle, editor. Academic Publishers, Boca Raton, FL. 147–172.

2. Yeagle, P. L. 1985. Cholesterol and the cell membrane. Biochim. Biophys. Acta. 822:267–287.[Medline]

3. Parks, L. W., S. J. Smith, and J. H. Crowley. 1995. Biochemical and physiological effects of sterol alterations in yeast—a review. Lipids. 30:227–230.[Medline]

4. Parks, L. W., J. H. Crowley, F. W. Leak, S. J. Smith, and M. E. Tomeo. 1999. Use of sterol mutants as probes for sterol functions in the yeast, Saccharomyces cerevisiae. Crit. Rev. Biochem. Mol. Biol. 34:399–404.[CrossRef][Medline]

5. Chapman, D., N. F. Owens, M. C. Phillips, and D. A. Walker. 1969. Mixed monolayers of phospholipids and cholesterol. Biochim. Biophys. Acta. 183:458–465.[Medline]

6. Joos, P., and R. A. Demel. 1969. The interaction energies of cholesterol and lecithin in spread mixed monolayers at the air-water interface. Biochim. Biophys. Acta. 183:447–457.[Medline]

7. Marsan, M. P., E. Bellet-Amalric, I. Muller, G. Zaccai, and A. Milon. 1998. Plant sterols: a neutron diffraction study of sitosterol and stigmasterol in soybean phosphatidylcholine membranes. Biophys. Chem. 75:45–55.[CrossRef][Medline]

8. Smaby, J. M., H. L. Brockman, and R. E. Brown. 1994. Cholesterol's interfacial interactions with sphingomyelins and phosphatidylcholines: hydrocarbon chain structure determines the magnitude of condensation. Biochemistry. 33:9135–9142.[CrossRef][Medline]

9. Stockton, G. W., and I. C. Smith. 1976. A deuterium nuclear magnetic resonance study of the condensing effect of cholesterol on egg phosphatidylcholine bilayer membranes. I. Perdeuterated fatty acid probes. Chem. Phys. Lipids. 17:251–263.[CrossRef][Medline]

10. Demel, R. A., and B. De Kruyff. 1976. The function of sterols in membranes. Biochim. Biophys. Acta. 457:109–132.[Medline]

11. Clejan, S., R. Bittman, P. W. Deroo, Y. A. Isaacson, and A. F. Rosenthal. 1979. Permeability properties of sterol-containing liposomes from analogues of phosphatidylcholine-lacking acyl groups. Biochemistry. 18:2118–2125.[CrossRef][Medline]

12. Schuler, I., A. Milon, Y. Nakatani, G. Ourisson, A. M. Albrecht, P. Benveniste, and M. A. Hartmann. 1991. Differential effects of plant sterols on water permeability and on acyl chain ordering of soybean phosphatidylcholine bilayers. Proc. Natl. Acad. Sci. USA. 88:6926–6930.[Abstract/Free Full Text]

13. Bloom, M., and J. L. Thewalt. 1995. Time and distance scales of membrane domain organization. Mol. Membr. Biol. 12:9–13.[Medline]

14. Kuo, A. L., and C. G. Wade. 1979. Lipid lateral diffusion by pulsed nuclear magnetic resonance. Biochemistry. 18:2300–2308.[CrossRef][Medline]

15. Oradd, G., G. Lindblom, and P. W. Westerman. 2002. Lateral diffusion of cholesterol and dimyristoylphosphatidylcholine in a lipid bilayer measured by pulsed field gradient NMR spectroscopy. Biophys. J. 83:2702–2704.[Abstract/Free Full Text]

16. Filippov, A., G. Oradd, and G. Lindblom. 2003. The effect of cholesterol on the lateral diffusion of phospholipids in oriented bilayers. Biophys. J. 84:3079–3086.[Abstract/Free Full Text]

17. Rubenstein, J. L., B. A. Smith, and H. M. McConnell. 1979. Lateral diffusion in binary mixtures of cholesterol and phosphatidylcholines. Proc. Natl. Acad. Sci. USA. 76:15–18.[Abstract/Free Full Text]

18. Vist, M. R., and J. H. Davis. 1990. Phase equilibria of cholesterol/dipalmitoylphosphatidylcholine mixtures: 2H nuclear magnetic resonance and differential scanning calorimetry. Biochemistry. 29:451–464.[CrossRef][Medline]

19. Huang, T. H., C. W. Lee, S. K. Das Gupta, A. Blume, and R. G. Griffin. 1993. A 13C and 2H nuclear magnetic resonance study of phosphatidylcholine/cholesterol interactions: characterization of liquid-gel phases. Biochemistry. 32:13277–13287.[CrossRef][Medline]

20. McMullen, T. P., R. N. Lewis, and R. N. McElhaney. 1993. Differential scanning calorimetric study of the effect of cholesterol on the thermotropic phase behavior of a homologous series of linear saturated phosphatidylcholines. Biochemistry. 32:516–522.[CrossRef][Medline]

21. McMullen, T. P., and R. N. McElhaney. 1997. Differential scanning calorimetric studies of the interaction of cholesterol with distearoyl and dielaidoyl molecular species of phosphatidylcholine, phosphatidylethanolamine, and phosphatidylserine. Biochemistry. 36:4979–4986.[CrossRef][Medline]

22. Marsan, M. P., I. Muller, C. Ramos, F. Rodriguez, E. J. Dufourc, J. Czaplicki, and A. Milon. 1999. Cholesterol orientation and dynamics in dimyristoylphosphatidylcholine bilayers: a solid state deuterium NMR analysis. Biophys. J. 76:351–359.[Abstract/Free Full Text]

23. Sankaram, M. B., and T. E. Thompson. 1990. Modulation of phospholipid acyl chain order by cholesterol. A solid-state 2H nuclear magnetic resonance study. Biochemistry. 29:10676–10684.[CrossRef][Medline]

24. Semer, R., and E. Gelerinter. 1979. A spin label study of the effects of sterols on egg lecithin bilayers. Chem. Phys. Lipids. 23:201–211.[CrossRef]

25. Urbina, J. A., S. Pekerar, H. B. Le, J. Patterson, B. Montez, and E. Oldfield. 1995. Molecular order and dynamics of phosphatidylcholine bilayer membranes in the presence of cholesterol, ergosterol and lanosterol: a comparative study using 2H-, 13C- and 31P-NMR spectroscopy. Biochim. Biophys. Acta. 1238:163–176.[Medline]

26. McMullen, T. P. W., and R. N. McElhaney. 1996. Physical studies of cholesterol-phospholipid interactions. Curr. Opin. Colloid Interface Sci. 1:83–90.

27. Ohvo-Rekila, H., B. Ramstedt, P. Leppimaki, and J. P. Slotte. 2002. Cholesterol interactions with phospholipids in membranes. Prog. Lipid Res. 41:66–97.[CrossRef][Medline]

28. Ghosh, D., and J. Tinoco. 1972. Monolayer interactions of individual lecithins with natural sterols. Biochim. Biophys. Acta. 266:41–49.[Medline]

29. Schuler, I., G. Duportail, N. Glasser, P. Benveniste, and M. A. Hartmann. 1990. Soybean phosphatidylcholine vesicles containing plant sterols: a fluorescence anisotropy study. Biochim. Biophys. Acta. 1028:82–88.[Medline]

30. Hsueh, Y., K. Gilbert, C. Trandum, M. Zuckermann, and J. Thewalt. 2005. The effect of ergosterol on dipalmitoylphosphatidylcholine bilayers: a deuterium NMR and calorimetric study. Biophys. J. 88:1799–1808.[Abstract/Free Full Text]

31. Kovacs, F. A., J. K. Denny, Z. Song, J. R. Quine, and T. A. Cross. 2000. Helix tilt of the M2 transmembrane peptide from Influenza A virus: an intrinsic property. J. Mol. Biol. 295:117–125.[CrossRef][Medline]

32. Bak, M., R. P. Bywater, M. Hohwy, J. K. Thomsen, K. Adelhorst, H. J. Jakobsen, O. W. Sorensen, and N. C. Nielsen. 2001. Conformation of alamethicin in oriented phospholipid bilayers determined by (15)N solid-state nuclear magnetic resonance. Biophys. J. 81:1684–1698.[Abstract/Free Full Text]

33. Marassi, F. M., C. Ma, J. J. Gesell, and S. J. Opella. 2000. Three-dimensional solid-state NMR spectroscopy is essential for resolution of resonances from in-plane residues in uniformly N-15-labeled helical membrane proteins in oriented lipid bilayers. J. Magn. Reson. 144:156–161.[CrossRef][Medline]

34. Yamaguchi, S., D. Huster, A. Waring, R. I. Lehrer, W. Kearney, B. F. Tack, and M. Hong. 2001. Orientation and dynamics of an antimicrobial peptide in the lipid bilayer by solid-state NMR spectroscopy. Biophys. J. 81:2203–2214.[Abstract/Free Full Text]

35. Park, S. H., A. A. Morse, A. A. Nevzorov, M. F. Mesleh, M. Oblatt-Montal, M. Montal, and S. J. Opella. 2003. Three-dimensional structure of the channel-forming trans-membrane domain of virus protein u (Vpu) from HIV-1. J. Mol. Biol. 333:409–424.[CrossRef][Medline]

36. Glaubitz, C., G. Grobner, and A. Watts. 2000. Structural and orientational information of the membrane-embedded M13 coat protein by C-13-MAS NMR spectroscopy. Biochim. Biophys. Acta. 1463:151–161.[Medline]

37. Creemers, A. F., S. Kiihne, P. H. Bovee-Geurts, W. J. DeGrip, J. Lugtenburg, and H. J. de Groot. 2002. (1)H and (13)C MAS NMR evidence for pronounced ligand-protein interactions involving the ionone ring of the retinylidene chromophore in rhodopsin. Proc. Natl. Acad. Sci. USA. 99:9101–9106.[Abstract/Free Full Text]

38. Haon, S., S. Augé, M. Tropis, A. Milon, and N. Lindley. 1993. Low-cost production of perdeuterated biomass using methylotrophic yeasts. J. Label. Comp. Rad. 33:1053–1063.[CrossRef]

39. Massou, S., S. Auge, M. Tropis, N. D. Lindley, and A. Milon. 1998. NMR analyses of deuterated phospholipids isolated from Pichia angusta. J. Chim. Phys. Phys. Chim. Biol. 95:406–411.

40. Massou, S., V. Puech, F. Talmont, P. Demange, N. D. Lindley, M. Tropis, and A. Milon. 1999. Heterologous expression of a deuterated membrane-integrated receptor and partial deuteration in methylotrophic yeasts. J. Biomol. NMR. 14:231–239.[CrossRef][Medline]

41. de Dios, A. C. 1996. Ab initio calculations of the NMR chemical shift. Prog. Nucl. Magn. Spectrosc. 29:229–278.[CrossRef]

42. de Dios, A. C., and E. Oldfield. 1996. Recent progress in understanding chemical shifts. Solid State Nucl. Magn. Reson. 6:101–125.[CrossRef][Medline]

43. Lesage, A., M. Bardet, and L. Emsley. 1999. Through-bond carbon-carbon connectivities in disordered solids by NMR. J. Am. Chem. Soc. 121:10987–10993.[CrossRef]

44. Frisch, M. J., G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, V. G. Zakrzewski, J. A. Montgomery, R. E. Stratmann, Jr., J. C. Burant, S. Dapprich, J. M. Millam, et al. 1998. Gaussian 98, Rev. A9. Gaussian, Pittsburgh, PA.

45. Jolibois, F., O. Soubias, V. Réat, and A. Milon. 2004. Understanding sterol-membrane interactions. I. Hartree-Fock versus DFT calculations of 13C and 1H NMR isotropic chemical shifts of sterols in solution and hydrogen bonding effects analysis. Chem. Eur. J. 10:5996–6004.[CrossRef]

46. Davis, J. H., and M. Auger. 1999. Static and magic angle spinning NMR of membrane peptides and proteins. Progr. Nucl. Magn. Res. Spectrosc. 35:1–84.[CrossRef]

47. Soubias, O., F. Jolibois, V. Réat, an