Slow-motion ESR study of order and dynamics in oriented lipid multibilayers: effects of unsaturation and hydration

Electron spin resonance (ESR) experiments were carried out on 3-doxyl-5 alpha-cholestane spin-label (CSL) molecules embedded in macroscopically oriented multibilayers of dimyristoylphosphatidylcholine (DMPC), palmitoyloleoylphosphatidylcholine (POPC), dioleoylphosphatidylcholine (DOPC) and dilinoleoylphosphatidylcholine (DLPC). For these lipids we studied the effects of temperature, hydration and unsaturation on the orientational order parameters and rotational motions of the probe molecules in the liquid crystalline phase. The experimental ESR spectra were simulated by a numerical solution of the stochastic Liouville equation (SLE) for the density matrix of a spin-label molecule. This allows extraction of detailed information about both molecular order and rotational dynamics. The data show that, in our temperature range, the lipid systems are in the slow-motion regime, thereby precluding a motional narrowing interpretation. This is illustrated by a simple model calculation which shows that a fast-motion interpretation seriously overestimates the order parameters. We have compared our results with data obtained independently from angle-resolved fluorescence depolarization (AFD) experiments on oriented bilayers in which 1-[4-(trimethylammonio)phenyl]-6-phenyl-1,3,5-hexatriene (TMA-DPH) molecules were used as fluorescent probes (Deinum et al., (1988) Biochemistry 27, 852-860). It is found that the orientational order and the rotational dynamics obtained with both techniques agree well. This shows that the probe molecules do not perturb the local bilayer structure to any large extent and that they indeed reflect the intrinsic behaviour of the lipid molecules. Upon increase in temperature or hydration, we observe faster reorientational motion and lower molecular ordering. In contrast, we do not find any systematic effect of unsaturation on molecular reorientational motion. Our results indicate that changes in membrane molecular order and reorientational dynamics have to be considered separately and are not necessarily correlated as implied by the common concept of membrane fluidity.     

Effects of steroid molecules on the dynamical structure of dioleoylphosphatidylcholine and digalactosyldiacylglycerol bilayers

The ESR spectra of cholestane spin labels (CSL) in dioleoylphosphatidylcholine (DOPC) bilayers containing 20 wt% of cholesterol, 7-dehydrocholesterol, beta-sitosterol, stigmasterol and lanosterol exhibit a marked similarity, thus indicating that these steroids induced the same effects on the lipid bilayer over the temperature range 21-55 degrees C. The incorporation of these steroids into the DOPC bilayers enhances the orientational order of the CSL molecules at every temperature studied, but only induces a pronounced slow-down in their rotational motions at temperatures above 35 degrees C. Similar results were obtained in DOPC/ergosterol multilamellar liposomes, but the changes are now less pronounced than in the other five DOPC/steroid systems. In contrast, the addition of stigmasterol to digalactosyldiacylglycerol (DGDG) bilayers appears to increase the order parameter mean value of P2, without affecting the diffusion coefficients. Furthermore, the incorporation of 7-dehydrocholesterol to DGDG bilayers causes a large enhancement in the orientational order, but has only a small effect on D perpendicular of the CSL molecules. Importantly, this latter effect appears to be independent of temperature. The marked changes in the rates of the rotational motion brought about by the addition of steroids, contrasts with the lack of a significant effect of unsaturation on the bilayer dynamics reported by us previously (Korstanje et al. (1989), Biochim. Biophys. Acta 980, 225-233, and 982, 196-204).     

Molecular dynamics of the cyclic lipodepsipeptides' action on model membranes: effects of syringopeptin22A, syringomycin E, and syringotoxin studied by EPR technique

Interaction of pore-forming toxins, syringopeptin22A (SP22A), syringomycin E (SRE) and syringotoxin (ST), with model membranes were investigated. Liposomes were prepared from saturated phospholipids (DPPC or DMPC) or from binary mixtures of DPPC with varying amount of DOPC or cholesterol. The effects of the three toxins on the molecular order and dynamics of the lipids were studied using electron paramagnetic resonance (EPR) techniques. SP22A was the most-, SRE less-, and ST the least effective to increase the ordering and to decrease the rotational correlation time of the lipid molecules. The effects were more pronounced: (a) on small unilamellar vesicles (SUVs) than on multilamellar vesicles (MUVs); (b) on pure DPPC than on DPPC-cholesterol or DPPC-DOPC mixtures. Fluidity changes, determined from EPR spectra at different concentrations of the toxin, suggested the shell structure of the lipid molecules in pore formation. EPR spectra observed at different depth of the hydrocarbon chain of the lipid molecules implied an active role of the lipid molecules in the architecture of the pores created in the presence of the three toxins. Temperature dependence of the fluidity of the SUVs treated with toxins has shown an abrupt and irreversible change in the molecular dynamics of the lipid molecules at a temperature close to the pretransition, depending on the toxin species and the lipid composition. Coalescence and aggregation of the SUVs were proposed as the origin of this irreversible change.     

Structural and dynamic effects of oxidatively modified phospholipids in unsaturated lipid membranes.

Phospholipid hydroperoxides and phospholipid alcohols are two of the major forms of oxidatively modified phospholipids produced during oxidant stress and lipid peroxidation. The process of lipid peroxidation is known to affect the physiological function of membranes. We, therefore, investigated the effects of lipid peroxidation products on the molecular interactions in membranes. Our study was specifically focused on the effects of lipid peroxidation products on static membrane structure (molecular orientational order) and on the reorientational dynamics of the probe molecules in lipid bilayers. The study was done by performing angle-resolved fluorescence depolarization measurements (AFD) on the fluorescent probe diphenylhexatriene (DPH) and by performing angle-resolved electron spin resonance (A-ESR) measurements on cholestane (CSL) nitroxide spin probes embedded in macroscopically oriented planar bilayers consisting of 2-10% 1-palmitoyl-2-(9/13-hydroperoxylinoleoyl)phosphatidylcholine (PLPC-OOH) or 1-palmitoyl-2-(9/13-hydroxylinoleoyl)phosphatidylcholine (PLPC-OH) in 1-palmitoyl-2-linoleoylphosphatidylcholine (PLPC) or dilinoleoylphosphatidylcholine (DLPC). Both probe molecules have rigid cylindrical geometries and report on the overall molecular order and dynamics. However, being more polar, the nitroxide spin probe CSL is preferentially located near the surface of the membrane, while the less polar fluorescent probe DPH reports preferentially near the central hydrophobic region of the lipid bilayers. The results show that the presence of relatively small amounts of oxidatively modified phospholipids within the PLPC or DLPC membranes causes pronounced structural effects as the molecular orientational order of the probe molecules is strongly decreased. In contrast, the effect on membrane reorientational dynamics is minimal.     

Atomistic simulation studies of cholesteryl oleates: model for the core of lipoprotein particles

We have conducted molecular dynamics simulations to gain insight into the atomic-scale properties of an isotropic system of cholesteryl oleate (CO) molecules. Cholesteryl esters are major constituents of low density lipoprotein particles, the key players in the formation of atherosclerosis, as well as the storage form of cholesterol. Here the aim is to clarify structural and dynamical properties of CO molecules under conditions, which are suggestive of those in the core of low density lipoprotein particles. The simulations in the fluid phase indicate that the system of CO molecules is characterized by an absence of translational order, as expected, while the orientational order between distinct CO molecules is significant at short distances, persisting over a molecular size. As for intramolecular properties, the bonds along the oleate chain are observed to be weakly ordered with respect to the sterol structure, unlike the bonds along the short hydrocarbon chain of cholesterol where the ordering is significant. The orientational distribution of the oleate chain as a whole with respect to the sterol moiety is of broad nature, having a major amount of extended and a less considerable proportion of bended structures. Distinct transient peaks at specific angles also appear. The diffusion of CO molecules is found to be a slow process and characterized by a diffusion coefficient of the order of 2x10(-9) cm2/s. This is considerably slower than diffusion, e.g., in ordered domains of lipid membranes rich in sphingomyelin and cholesterol. Analysis of the rotational diffusion rates and trans-to-gauche transition rates yield results consistent with experiments.     

Dynamics and ordering in mixed model membranes of dimyristoylphosphatidylcholine and dimyristoylphosphatidylserine: a 250-GHz electron spin resonance study using cholestane.

We report here on a 250-GHz electron spin resonance (ESR) study of macroscopically aligned model membranes composed of mixtures of dimyristoylphosphatidylcholine (DMPC) and dimyristoylphosphatidylserine (DMPS), utilizing the nixtroxide-labeled cholesterol analog cholestane (CSL). Two clearly resolved spectral components, distinct in both their ordering and dynamics, are resolved. The major component in membranes composed mostly of DMPC shows typical characteristics, with the long axis of CSL parallel to the bilayer normal with slow (10(6) </= R </= 10(7) s-1) rotational diffusion rates, as expected for cholesterol. The second component grows in as the mole fraction of DMPS increases. A detailed analysis shows that CSL senses a local, strongly biaxial environment. Our results imply that the inefficient packing between cholesterol and DMPS occurs probably because of the strong interactions between the PS headgroups, which provide the local biaxiality. Such a packing of the headgroups has been predicted by molecular dynamics simulations but had not been observed experimentally. The analysis of these spectral components was greatly aided by the excellent orientational resolution provided by the 250-GHz spectra. This enabled the key qualitative features of this interpretation to be "read" off the spectra before the detailed analysis.     

Molecular dynamics investigation of the structure of a fully hydrated gel-phase dipalmitoylphosphatidylcholine bilayer.

We report the results of a constant pressure and temperature molecular dynamics simulation of a gel-phase dipalmitoylphosphatidylcholine bilayer with nw = 11.8 water molecules/lipid at 19 degrees C. The results of the simulation were compared in detail with a variety of x-ray and neutron diffraction data. The average positions of specific carbon atoms along the bilayer normal and the interlamellar spacing and electron density profile were in very good agreement with neutron and x-ray diffraction results. The area per lipid and the details of the in-plane hydrocarbon chain structure were in excellent agreement with wide-angle x-ray diffraction results. The only significant deviation is that the chains met in a pleated arrangement at the bilayer center, although they should be parallel. Novel discoveries made in the present work include the observation of a bimodal headgroup orientational distribution. Furthermore, we found that there are a significant number of gauche conformations near the ends of the hydrocarbon chains and, in addition to verifying a previous suggestion that there is partial rotational ordering in the hydrocarbon chains, that the two chains in a given molecule are inequivalent with respect to rotations. Finally, we have investigated the lipid/water interface and found that the water penetrates beneath the headgroups, but not as far as the carbonyl groups, that the phosphates are strongly hydrated almost exclusively at the nonesterified oxygen atoms, and that the hydration of the ammonium groups is more diffuse, with some water molecules concentrated in the grooves between the methyl groups.     

Interactions of pyrethroids with phosphatidylcholine bilayers: comparisons in liposomal systems exhibiting large or small radii of curvature

Pyrethroid interactions with dipalmitoyl phosphatidylcholine (DPPC) vesicles have been characterized in bilayers having large and small radii of curvature. The abilities of pyrethroids to alter the gel-fluid phase transition profiles were determined by steady state fluorescence anisotropy and phase-modulation lifetime techniques using the fluorescent probes cis- and trans-parinaric acid. Using the geometric isomers of parinaric acid as membrane probes, pyrethroids were found to lower the phase transition temperature (Tc) of DPPC large multilamellar vesicles with the same order of comparative effectiveness as previously reported using the fluorescent probe 1,6-diphenyl-1,3,5-hexatriene (DPH). Permethrin had a greater depressive effect upon the Tc of DPPC in the small unilamellar vesicle (SUV) system than in the large multilamellar system. Conversely, allethrin was less effective in reducing the Tc of DPPC SUVs. The enhanced effect of permethrin in decreasing the Tc of DPPC SUVs was greatest in regions of more rigid lipid packing, as determined by trans-parinaric acid fluorescence parameters. The results indicate that changes in lipid packing configuration caused by differing bilayer radii of curvature may alter the interactive characteristics of pyrethroids with lipid membranes.     

Hydration properties and structure of phosphatidylcholine membranes in the presence of n-nonyl bromide

Interaction of chemical fusogen n-nonyl bromide with a model membrane formed from phosphatidylcholine was studied using 2D-NMR spectra of heavy water and 31P-NMR proton decoupled spectra of the lipid phosphate group in multilamellar lipid dispersions. n-Nonyl bromide was found to influence the hydration layer of the model membrane. No participation of phosphatidylcholine molecules in non-bilayer configurations of the membrane was observed.     

The dynamics of protein hydration water: a quantitative comparison of molecular dynamics simulations and neutron-scattering experiments

We present results from an extensive molecular dynamics simulation study of water hydrating the protein Ribonuclease A, at a series of temperatures in cluster, crystal, and powder environments. The dynamics of protein hydration water appear to be very similar in crystal and powder environments at moderate to high hydration levels. Thus, we contend that experiments performed on powder samples are appropriate for discussing hydration water dynamics in native protein environments. Our analysis reveals that simulations performed on cluster models consisting of proteins surrounded by a finite water shell with free boundaries are not appropriate for the study of the solvent dynamics. Detailed comparison to available x-ray diffraction and inelastic neutron-scattering data shows that current generation force fields are capable of accurately reproducing the structural and dynamical observables. On the time scale of tens of picoseconds, at room temperature and high hydration, significant water translational diffusion and rotational motion occur. At low hydration, the water molecules are translationally confined but display appreciable rotational motion. Below the protein dynamical transition temperature, both translational and rotational motions of the water molecules are essentially arrested. Taken together, these results suggest that water translational motion is necessary for the structural relaxation that permits anharmonic and diffusive motions in proteins. Furthermore, it appears that the exchange of protein-water hydrogen bonds by water rotational/librational motion is not sufficient to permit protein structural relaxation. Rather, the complete exchange of protein-bound water molecules by translational displacement seems to be required.     

ESR studies of spin-labeled membranes aligned by isopotential spin-dry ultracentrifugation: lipid-protein interactions.

Electron spin resonance (ESR) studies have been performed on spin-labeled model membranes aligned using the isopotential spin-dry ultracentrifugation (ISDU) method of Clark and Rothschild. This method relies on sedimentation of the membrane fragments onto a gravitational isopotential surface with simultaneous evaporation of the solvent in a vacuum ultracentrifuge to promote alignment. The degree of alignment obtainable using ISDU, as monitored by ESR measurements of molecular ordering for both lipid (16-PC) and cholestane spin labels (CSL), in dipalmitoylphosphatidylcholine (DPPC) model membranes compares favorably with that obtainable by pressure-annealing. The much gentler conditions under which membranes may be aligned by ISDU greatly extends the range of macroscopically aligned membrane samples that may be investigated by ESR. We report the first ESR study of an integral membrane protein, bacteriorhodopsin (BR) in well-aligned multilayers. We have also examined ISDU-aligned DPPC multilayers incorporating a short peptide gramicidin A' (GA), with higher water content than previously studied. 0.24 mol% BR/DPPC membranes with CSL probe show two distinct components, primarily in the gel phase, which can be attributed to bulk and boundary regions of the bilayer. The boundary regions show sharply decreased molecular ordering and spectral effects comparable to those observed from 2 mol% GA/DPPC membranes. The boundary regions for both BR and GA also exhibit increased fluidity as monitored by the rotational diffusion rates. The high water content of the GA/DPPC membranes reduces the disordering effect as evidenced by the reduced populations of the disordered components. The ESR spectra obtained slightly below the main phase transition of DPPC from both the peptide- and protein-containing membranes reveals a new component with increased ordering of the lipids associated with the peptide or protein. This increase coincides with a broad endothermic peak in the DSC, suggesting a disaggregation of both the peptide and the protein before the main phase transition of the lipid. Detailed simulations of the multicomponent ESR spectra have been performed by the latest nonlinear least-squares methods, which have helped to clarify the spectral interpretations. It is found that the simulations of ESR spectra from CSL in the gel phase for all the lipid membranes studied could be significantly improved by utilizing a model with CSL molecules existing as both hydrogen-bonded to the bilayer interface and non-hydrogen-bonded within the bilayer.     

Liposome complexation efficiency monitored by FRET: effect of charge ratio,helper lipid and plasmid size

Cationic lipid/DNA complexes (lipoplexes) are promising vehicles for DNA vaccines or gene therapy. In these systems, transfection efficiency is highly related to lipoplex charge ratio, since lipoplexes with charge ratios (±) lower than electroneutrality have most DNA uncovered by the liposomes, and thus are unprotected from enzyme degradation. However, a large excess of cationic lipids is undesirable because of eventual cytotoxicity. The aim of this work was to determine the minimum charge ratio from which all DNA molecules are complexed by the liposomes varying the lipid formulation and plasmid size, using a new FRET (fluorescence resonance energy transfer) methodology. The similarity of FRET results, fluorescence intensity data and fluorescence decays of several charge ratios above (±) ≥ 4 or 5 confirmed that once all DNA is covered by the liposomes, additional lipid molecules do not affect the lipoplex multilamellar repeat distance. It was also verified by FRET that the presence of helper lipid reduces the amount of cationic lipid required for DNA protection but does not affect the lipoplex multilamellar repeat distance. This distance varies with the plasmid size when supercoiled plasmid is used, being apparently larger when longer plasmids are used. Our study indicates that, despite the complexity of these systems not being totally described by our model, FRET is an informative technique in lipoplex characterization.     

Electron spin resonance and electron-spin-echo study of oriented multilayers of L alpha-dipalmitoylphosphatidylcholine water systems.

A detailed electron spin resonance (ESR) study of spin-labeled-oriented multilayers of L alpha-dipalmitoylphosphatidylcholine (DPPC) water systems for low water content (2-10% by weight) is reported with the purpose of characterizing the dynamical and structural properties of model membrane systems. Emphasis is placed on the value of combining such experiments with detailed simulations based on current slow-motional theories. Information is obtained regarding ordering and anisotropic rotational diffusion rates via ESR lineshape analysis over the entire motional range, from the fast motional region through the moderately slow and slow to the rigid limit. This includes the low-temperature gel phase, the liquid crystalline L alpha (1) phase and what appears to be a third high-temperature phase above the L alpha phase. Cholestane (CSL) and spin-labeled DPPC (5-PC, 8-PC, and 16-PC) have been used to probe different depths of the bilayer. While CSL and 5-PC both reflect the high ordering of the bilayer close to the lipid-water interface, CSL appears to be located close enough to the water for the nitroxide to be involved in hydrogen bonding with water molecules. 16-PC reflects the relatively low ordering near the tail of the hydrocarbon chain in the bilayer. Quantitative estimates of ordering and motion are obtained for these cases. The results from CSL indicate that close to the lipid-water interface the DPPC molecule is oriented approximately perpendicular to the bilayer in these low water-content systems. However, all three labeled lipid probes indicate that the hydrocarbon chain of DPPC may be bent away from the bilayer normal by as much as 30 degrees and this evidence is stronger at low temperatures. When cholesterol is added to the DPPC-water system at a concentration greater than or equal to 2.5 mol %, the ordering is greatly increased although the rotational diffusion rate remains almost unaffected in the gel phase. Electron spin echoes (ESE) are observed for the first time from oriented lipid-water multilayers. Results obtained from cw ESR lineshape analysis are correlated with data from ESE experiments, which give a more direct measurement of relaxation times. These results indicate that for detection of very slow motions (close to the rigid limit) ESE experiments are more sensitive to dynamics than continuous wave ESR for which inhomogeneous broadening becomes a major problem.     

NMR investigations of interactions between anesthetics and lipid bilayers

Interactions between anesthetics (lidocaine and short chain alcohols) and lipid membranes formed by dimyristoylphosphatidylcholine (DMPC) were studied using NMR spectroscopy. The orientational order of lidocaine was investigated using deuterium NMR on a selectively labelled compound whereas segmental ordering in the lipids was probed by two-dimensional ¹H-¹³C separated local field experiments under magic-angle spinning conditions. In addition, trajectories generated in molecular dynamics (MD) computer simulations were used for interpretation of the experimental results. Separate simulations were carried out with charged and uncharged lidocaine molecules. Reasonable agreement between experimental dipolar interactions and the calculated counterparts was observed. Our results clearly show that charged lidocaine affects significantly the lipid headgroup. In particular the ordering of the lipids is increased accompanied by drastic changes in the orientation of the P-N vector in the choline group.     

Molecular order and dynamics in planar lipid bilayers: effects of unsaturation and sterols

Angle-resolved fluorescence depolarization experiments were carried out on 1,6-diphenyl-1,3,5-hexatriene (DPH) and 1-[4-(trimethylammonio)phenyl]-6-phenyl-1,3,5-hexatriene (TMA-DPH) molecules embedded in macroscopically oriented multilayers of saturated [dimyristoylphosphatidylcholine (DMPC)] and unsaturated [palmitoyloleoylphosphatidylcholine (POPC), dioleoylphosphatidylcholine (DOPC), dilineoylphosphatidylcholine (DLPC), plant digalactosyldiglyceride (DGDG)] lipids with and without cholesterol. In all the lipid systems studied the order parameter (P2) of TMA-DPH molecules was found to be higher than that for DPH. Considerations of the order parameter (P4), however, indicate that DPH molecules have a heterogeneous distribution in bilayers of unsaturated lipids, with a significant fraction of the molecules lying with their long axes parallel to the bilayer planes. Both the DPH and TMA-DPH molecules exhibit a decrease in the molecular order as well as a decrease in their rates of motion on increasing the unsaturation of the hydrocarbon chains. The addition of cholesterol tends to reverse this effect, with an increase in both the order and dynamics. Bilayers of DOPC, however, exhibit a somewhat different result. It is suggested that the discrepancies between these observations and findings with lipid vesicle systems simply reflect the effects of curvature on the behavior of the probe molecules. The results indicate that the concept of membrane fluidity must be used with great caution.     

Molecular dynamics simulations of unsaturated lipid bilayers: effects of varying the numbers of double bonds

The importance of unsaturated, and especially polyunsaturated phosphatidylcholine molecules for the functional properties of biological membranes is widely accepted. Here, the effects of unsaturation on the nanosecond-scale structural and dynamic properties of the phosphatidylcholine bilayer were elucidated by performance of multinanosecond molecular dynamics simulations of all-atom bilayer models. Bilayers of dipalmitoylphosphatidylcholine and its mono-, di-, and tetraunsaturated counterparts were simulated, containing, respectively, oleoyl, linoleoyl, or arachidonoyl chains in the sn-2 position. Analysis of the simulations focused on comparison of the structural properties, especially the ordering of the chains in the membranes. Although the results suggest some problems in the CHARMM force field of the lipids when applied in a constant pressure ensemble, the features appearing in the ordering of the unsaturated chains are consistent with the behaviour known from ²H NMR experiments. The rigidity of the double bonds is compensated by the flexibility of skew state single bonds juxtaposed with double bonds. The presence of double bonds in the sn-2 chains considerably reduces the order parameters of the CH bonds. Moreover, the double bond region of tetraunsaturated chains is shown to span all the way from the bilayer centre to the head group region.     

Orientational and rotational velocity correlation functions for water-hydrating B- and Z-DNA double helices

The molecular dynamics simulations reported earlier for the structure and dynamics of water molecules hydrating B- and Z-DNA double helices are analyzed for the orientational correlation functions and the proton rotational velocity autocorrelation functions. The spectra of the rotational velocity autocorrelation functions obtained from the simulation results are compared with the neutron inelastic scattering experiments on hydrated Na-DNA samples. The results predict a small frequency component associated with water molecules bound to the double helices that disappears for waters away from the double helix.     

Molecular dynamics simulations and 2H NMR study of the GalCer/DPPG lipid bilayer

Molecular dynamics simulations were performed on a two-component lipid bilayer system in the liquid crystalline phase at constant pressure and constant temperature. The lipid bilayers were composed of a mixture of neutral galactosylceramide (GalCer) and charged dipalmitoylphosphatidylglycerol (DPPG) lipid molecules. Two lipid bilayer systems were prepared with GalCer:DPPG ratio 9:1 (10%-DPPG system) and 3:1 (25%-DPPG system). The 10%-DPPG system represents a collapsed state lipid bilayer, with a narrow water space between the bilayers, and the 25%-DPPG system represents an expanded state with a fluid space of approximately 10 nm. The number of lipid molecules used in each simulation was 1024, and the length of the production run simulation was 10 ns. The simulations were validated by comparing the results with experimental data for several important aspects of the bilayer structure and dynamics. Deuterium order parameters obtained from (2)H NMR experiments for DPPG chains are in a very good agreement with those obtained from molecular dynamics simulations. The surface area per GalCer lipid molecule was estimated to be 0.608 +/- 0.011 nm(2). From the simulated electron density profiles, the bilayer thickness defined as the distance between the phosphorus peaks across the bilayer was calculated to be 4.21 nm. Both simulation systems revealed a tendency for cooperative bilayer undulations, as expected in the liquid crystalline phase. The interaction of water with the GalCer and DPPG oxygen atoms results in a strong water ordering in a spherical hydration shell and the formation of hydrogen bonds (H-bonds). Each GalCer lipid molecule makes 8.6 +/- 0.1 H-bonds with the surrounding water, whereas each DPPG lipid molecule makes 8.3 +/- 0.1 H-bonds. The number of water molecules per GalCer or DPPG in the hydration shell was estimated to be 10-11 from an analysis of the radial distribution functions. The formation of the intermolecular hydrogen bonds was observed between hydroxyl groups from the opposing GalCer sugar headgroups, giving an energy of adhesion in the range between -1.0 and -3.4 erg/cm(2). We suggest that this value is the contribution of the hydrogen-bond component to the net adhesion energy between GalCer bilayers in the liquid crystalline phase.