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.     

Mean-field calculations of chain packing and conformational statistics in lipid bilayers: comparison with experiments and molecular dynamics studies.

A molecular, mean-field theory of chain packing statistics in aggregates of amphiphilic molecules is applied to calculate the conformational properties of the lipid chains comprising the hydrophobic cores of dipalmitoyl-phosphatidylcholine (DPPC), dioleoyl-phosphatidylcholine (DOPC), and palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers in their fluid state. The central quantity in this theory, the probability distribution of chain conformations, is evaluated by minimizing the free energy of the bilayer assuming only that the segment density within the hydrophobic region is uniform (liquidlike). Using this distribution we calculate chain conformational properties such as bond orientational order parameters and spatial distributions of the various chain segments. The lipid chains, both the saturated palmitoyl (-(CH2)14-CH3) and the unsaturated oleoyl (-(CH2)7-CH = CH-(CH2)7-CH3) chains are modeled using rotational isomeric state schemes. All possible chain conformations are enumerated and their statistical weights are determined by the self-consistency equations expressing the condition of uniform density. The hydrophobic core of the DPPC bilayer is treated as composed of single (palmitoyl) chain amphiphiles, i.e., the interactions between chains originating from the same lipid headgroup are assumed to be the same as those between chains belonging to different molecules. Similarly, the DOPC system is treated as a bilayer of oleoyl chains. The POPC bilayer is modeled as an equimolar mixture of palmitoyl and oleoyl chains. Bond orientational order parameter profiles, and segment spatial distributions are calculated for the three systems above, for several values of the bilayer thickness (or, equivalently, average area/headgroup) chosen, where possible, so as to allow for comparisons with available experimental data and/or molecular dynamics simulations. In most cases the agreement between the mean-field calculations, which are relatively easy to perform, and the experimental and simulation data is very good, supporting their use as an efficient tool for analyzing a variety of systems subject to varying conditions (e.g., bilayers of different compositions or thicknesses at different temperatures).     

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.     

Location of the TEMPO moiety of TEMPO-PC in phosphatidylcholine bilayers is membrane phase dependent

The (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) moiety tethered to the headgroup of phosphatidylcholine (PC) lipid is employed in spin labeling electron paramagnetic resonance spectroscopy to probe the water dynamics near lipid bilayer interfaces. Due to its amphiphilic character, however, TEMPO spin label could partition between aqueous and lipid phases, and may even be stabilized in the lipid phase. Accurate assessment of the TEMPO-PC configuration in bilayer membranes is essential for correctly interpreting the data from measurements. Here, we carry out all-atom molecular dynamics (MD) simulations of TEMPO-PC probe in single-component lipid bilayers at varying temperatures, using two standard MD force fields. We find that, for a dipalmitoylphosphatidylcholine (DPPC) membrane whose gel-to-fluid lipid phase transition occurs at 314 K, while the TEMPO spin label is stabilized above the bilayer interface in the gel phase, there is a preferential location of TEMPO below the membrane interface in the fluid phase. For bilayers made of unsaturated lipids, 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), which adopt the fluid phase at ambient temperature, TEMPO is unequivocally stabilized inside the bilayers. Our finding of membrane phase-dependent positioning of the TEMPO moiety highlights the importance of assessing the packing order and fluidity of lipids under a given measurement condition.     

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).     

Effect of hydrostatic pressure on water penetration and rotational dynamics in phospholipid-cholesterol bilayers.

The effect of high hydrostatic pressure on the lipid bilayer hydration, the mean order parameter, and rotational dynamics of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) cholesterol vesicles has been studied by time-resolved fluorescence spectroscopy up to 1500 bar. Whereas the degree of hydration in the lipid headgroup and interfacial region was assessed from fluorescence lifetime data using the probe 1-(4-trimethylammonium-phenyl)-6-phenyl-1,3,5-hexatriene (TMA-DPH), the corresponding information in the upper acyl chain region was estimated from its effect on the fluorescence lifetime of and 3-(diphenylhexatrienyl)propyl-trimethylammonium (TMAP-DPH). The lifetime data indicate a greater level of interfacial hydration for DPPC bilayers than for POPC bilayers, but there is no marked difference in interchain hydration of the two bilayer systems. The addition of cholesterol at levels from 30 to 50 mol% to DPPC has a greater effect on the increase of hydrophobicity in the interfacial region of the bilayer than the application of hydrostatic pressure of several hundred to 1000 bar. Although the same trend is observed in the corresponding system, POPC/30 mol% cholesterol, the observed effects are markedly less pronounced. Whereas the rotational correlation times of the fluorophores decrease in passing the pressure-induced liquid-crystalline to gel phase transition of DPPC, the wobbling diffusion coefficient remains essentially unchanged. The wobbling diffusion constant of the two fluorophores changes markedly upon incorporation of 30 mol% cholesterol, and increases at higher pressures, also in the case of POPC/30 mol% cholesterol. The observed effects are discussed in terms of changes in the rotational characteristics of the fluorophores and the phase-state of the lipid mixture. The results demonstrate the ability of cholesterol to adjust the structural and dynamic properties of membranes composed of different phospholipid components, and to efficiently regulate the motional freedom and hydrophobicity of membranes, so that they can withstand even drastic changes in environmental conditions, such as high external hydrostatic pressure.     

Structure and dynamics of sphingomyelin bilayer: insight gained through systematic comparison to phosphatidylcholine

Sphingomyelin, one of the main lipid components of biological membranes, is actively involved in various cellular processes such as protein trafficking and signal transduction. In particular, specific lateral domains enriched in sphingomyelin and cholesterol have been proposed to play an important functional role in biomembranes, although their precise characteristics have remained unclear. A thorough understanding of the functional role of membranes requires detailed knowledge of their individual lipid components. Here, we employ molecular dynamics simulations to conduct a systematic comparison of a palmitoylsphingomyelin (PSM, 16:0-SM) bilayer with a membrane that comprises dipalmitoylphosphatidylcholine (DPPC) above the main phase transition temperature. We clarify atomic-scale properties that are specific to sphingomyelin due to its sphingosine moiety, and further discuss their implications for SM-rich membranes. We find that PSM bilayers, and in particular the dynamics of PSM systems, are distinctly different from those of a DPPC bilayer. When compared with DPPC, the strong hydrogen bonding properties characteristic to PSM are observed to lead to considerable structural changes in the polar headgroup and interface regions. The strong ordering of PSM acyl chains and specific ordering effects in the vicinity of a PSM-water interface reflect this issue clearly. The sphingosine moiety and related hydrogen bonding further play a crucial role in the dynamics of PSM bilayers, as most dynamic properties, such as lateral and rotational diffusion, are strongly suppressed. This is most evident in the rotational motion characterized by spin-lattice relaxation times and the decay of hydrogen bond autocorrelation functions that are expected to be important in complexation of SM with other lipids in many-component bilayers. A thorough understanding of SM bilayers would greatly benefit from nuclear magnetic resonance experiments for acyl chain ordering and dynamics, allowing full comparison of these simulations to experiments.     

An electron spin resonance study of interactions between gramicidin A'' and phosphatidylcholine bilayers.

The model of microscopic order and macroscopic disorder was used to stimulate electron spin resonance spectra of spin-labeled lipids, 5-PC, 10-PC, and 16-PC in multilamellar vesicles of dipalmitoylphosphatidylcholine (DPPC) containing gramicidin A' (GA) at temperatures above the gel-to-liquid crystal transition of DPPC. The simulations show that at a lower concentration of GA (i.e., molar ratios of DPPC/GA greater than 3), GA has only a slight effect on the acyl chain dynamics. The rotational diffusion rate around the axis parallel to the long hydrocarbon chain remains unchanged or increases slightly, while the rate around the perpendicular axes decreases slightly. These spectra from DPPC/GA mixtures could only be fit successfully with two or more components consistent with the well-known concept of "boundary lipids," that is, the peptide induces structural inhomogeneity in lipid bilayers. However, the spectra were significantly better fit with additional components that exhibit increased local ordering, implying decreased amplitude of rotational motion, rather than immobilized components with sharply a reduced rotational rate. The largest relative effects occur at the end of the acyl chains, where the average local order parameter St of 16-PC increases from 0.06 for pure lipid to 0.66 for 1:1 DPPC/GA. The inhomogeneity in ordering in DPPC bilayers due to GA decreases with increasing temperature. The hyperfine tensor component Azz increases for 10-PC and 16-PC when GA is incorporated into DPPC bilayers, indicating that water has deeply penetrated into the DPPC bilayers. Simulations of published electron spin resonance spectra of 14-PC in dimyristoylphosphatidylcholine/cytochrome oxidase complexes were also better fit by additional components that were more ordered, rather than immobilized. The average local order parameter in this case is found to increase from 0.11 for pure dimyristoylphosphatidylcholine to 0.61 for a lipid/protein ratio of 50. These spectra and their simulations are similar to the results obtained with 16-PC in the DPPC/GA mixtures. The relevance to studies of lipid-protein interactions for other proteins is briefly discussed.     

Phase behavior and nanoscale structure of phospholipid membranes incorporated with acylated C14-peptides

The thermotropic phase behavior and lateral structure of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers containing an acylated peptide has been characterized by differential scanning calorimetry (DSC) on vesicles and atomic force microscopy (AFM) on mica-supported bilayers. The acylated peptide, which is a synthetic decapeptide N-terminally linked to a C14 acyl chain (C14-peptide), is incorporated into DPPC bilayers in amounts ranging from 0-20 mol %. The calorimetric scans of the two-component system demonstrate a distinct influence of the C14-peptide on the lipid bilayer thermodynamics. This is manifested as a concentration-dependent downshift of both the main phase transition and the pretransition. In addition, the main phase transition peak is significantly broadened, indicating phase coexistence. In the AFM imaging scans we found that the C14-peptide, when added to supported gel phase DPPC bilayers, inserts preferentially into preexisting defect regions and has a noticeable influence on the organization of the surrounding lipids. The presence of the C14-peptide gives rise to a laterally heterogeneous bilayer structure with coexisting lipid domains characterized by a 10 A height difference. The AFM images also show that the appearance of the ripple phase of the DPPC lipid bilayers is unaffected by the C14-peptide. The experimental results are supported by molecular dynamics simulations, which show that the C14-peptide has a disordering effect on the lipid acyl chains and causes a lateral expansion of the lipid bilayer. These effects are most pronounced for gel-like bilayer structures and support the observed downshift in the phase-transition temperature. Moreover, the molecular dynamics data indicate a tendency of a tryptophan residue in the peptide sequence to position itself in the bilayer headgroup region.     

Order and fluidity of lipid membranes as determined by fluorescence anisotropy decay

The fluorescence anisotropy decay of four different probes in bilayers of dimyristoylphosphatidylcholine was measured. The probes are diphenylhexatriene, diphenyloctatetraene, trimethylaminodiphenylhexatriene, and trans-parinaric acid. The data for each probe were analyzed in terms of two orientational order parameters, the ordinary order parameter and a higher one, and two rotational diffusion coefficients. The order parameters are largely independent of probe size, but depend on the position of the probes along the membrane normal, thus reflecting the profile of lipid order. If a probe is located in the plateau region of lipid order, its order parameters are interpreted as representing the rigid-body order of lipids. According to this interpretation, the total lipid order in the plateau region originates about equally from rigid-body order and conformational order. The two order parameters obtained for each probe are used to derive approximate angular distributions of the probe molecules. The diffusion coefficient for rotation about the long molecular axis is found to be infinitely large, indicating unhindered rotation about this axis. The diffusion coefficient for rotation about the short molecular axes is evaluated for a viscosity which results as 0.2 poise. This viscosity for rotational diffusion is an order of magnitude smaller than the viscosity for lateral diffusion indicating that at least two viscosities are required to characterize the fluidity of a lipid membrane.Abbreviations FAD fluorescence anisotropy decay - DMR deuterium magnetic resonance - ESR electron spin resonance - DMPC dimyristoylphosphatidylcholine - DPPC dipalmitoylphosphatidylcholine - DPH 1,6-diphenyl-1,3,5-hexatriene - DPO 1,6-diphenyl-1,3,5,7-octatetraene - TMA-DPH 1-[4-(trimethylamino)phenyl]-6-phenyl-1,3,5-hexatriene - tPnA trans-parinaric acid - NPN N-phenyl-1-naphthylamine - BBO 2,5-bis(4-biphenylyl)oxazole     

Reparameterization of all-atom dipalmitoylphosphatidylcholine lipid parameters enables simulation of fluid bilayers at zero tension

Molecular dynamics simulations of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers using the CHARMM27 force field in the tensionless isothermal-isobaric (NPT) ensemble give highly ordered, gel-like bilayers with an area per lipid of approximately 48 A(2). To obtain fluid (L(alpha)) phase properties of DPPC bilayers represented by the CHARMM energy function in this ensemble, we reparameterized the atomic partial charges in the lipid headgroup and upper parts of the acyl chains. The new charges were determined from the electron structure using both the Mulliken method and the restricted electrostatic potential fitting method. We tested the derived charges in molecular dynamics simulations of a fully hydrated DPPC bilayer. Only the simulation with the new restricted electrostatic potential charges shows significant improvements compared with simulations using the original CHARMM27 force field resulting in an area per lipid of 60.4 +/- 0.1 A(2). Compared to the 48 A(2), the new value of 60.4 A(2) is in fair agreement with the experimental value of 64 A(2). In addition, the simulated order parameter profile and electron density profile are in satisfactory agreement with experimental data. Thus, the biologically more interesting fluid phase of DPPC bilayers can now be simulated in all-atom simulations in the NPT ensemble by employing our modified CHARMM27 force field.     

Residue specific partitioning of KL4 into phospholipid bilayers

KL₄, which has demonstrated success in the treatment of respiratory distress, is a synthetic helical, amphipathic peptide mimetic of lung surfactant protein B. The unusual periodicity of charged residues within KL₄ and its relatively high hydrophobicity distinguish it from canonical amphipathic helical peptides. Here we utilized site specific spin labeling of both lipids and the peptide coupled with EPR spectroscopy to discern the effects of KL₄ on lipid dynamics, the residue specific dynamics of hydrophobic regions within KL₄, and the partitioning depths of specific KL₄ residues into the DPPC/POPG and POPC/POPG lipid bilayers under physiologically relevant conditions. KL₄ induces alterations in acyl chain dynamics in a lipid-dependent manner, with the peptide partitioning more deeply into DPPC-rich bilayers. Combined with an earlier NMR study of changes in lipid dynamics on addition of KL₄ (V.C. Antharam et al., 2009), we are able to distinguish how KL₄ affects both collective bilayer motions and intramolecular acyl chain dynamics in a lipid-dependent manner. EPR power saturation results for spin labeled lipids demonstrate that KL₄ also alters the accessibility profiles of paramagnetic colliders in a lipid-dependent manner. Measurements of dynamics and depth parameters for individual spin-labeled residues within KL₄ are consistent with a model where the peptide partitions deeply into the lipid bilayers but lies parallel to the bilayer interface in both lipid environments; the depth of partitioning is dependent on the degree of lipid acyl chain saturation within the bilayer.     

A computer simulation of free-volume distributions and related structural properties in a model lipid bilayer.

A novel combined approach of molecular dynamics (MD) and Monte Carlo simulations is developed to calculate various free-volume distributions as a function of position in a lipid bilayer membrane at 323 K. The model bilayer consists of 2 x 100 chain molecules with each chain molecule having 15 carbon segments and one head group and subject to forces restricting bond stretching, bending, and torsional motions. At a surface density of 30 A2/chain molecule, the probability density of finding effective free volume available to spherical permeants displays a distribution with two exponential components. Both pre-exponential factors, p1 and p2, remain roughly constant in the highly ordered chain region with average values of 0.012 and 0.00039 A-3, respectively, and increase to 0.049 and 0.0067 A-3 at the mid-plane. The first characteristic cavity size V1 is only weakly dependent on position in the bilayer interior with an average value of 3.4 A3, while the second characteristic cavity size V2 varies more dramatically from a plateau value of 12.9 A3 in the highly ordered chain region to 9.0 A3 in the center of the bilayer. The mean cavity shape is described in terms of a probability distribution for the angle at which the test permeant is in contact with one of and does not overlap with anyone of the chain segments in the bilayer. The results show that (a) free volume is elongated in the highly ordered chain region with its long axis normal to the bilayer interface approaching spherical symmetry in the center of the bilayer and (b) small free volume is more elongated than large free volume. The order and conformational structures relevant to the free-volume distributions are also examined. It is found that both overall and internal motions have comparable contributions to local disorder and couple strongly with each other, and the occurrence of kink defects has higher probability than predicted from an independent-transition model.     

Reorientational dynamics in lipid vesicles and liposomes studied with ESR: effects of hydration, curvature and unsaturation

Electron spin resonance experiments were carried out on 3-doxyl-5 alpha-cholestane spin-label (CSL) molecules embedded in multilamellar liposomes and small unilamellar vesicles (SUVs) of palmitoyloleoylphosphatidylcholine (POPC), dioleoylphosphatidylcholine (DOPC) and dilinoleoylphosphatidylcholine (DLPC). The experimental spectra were analyzed by a numerical solution of the stochastic Liouville equation. Effects of temperature, presence of unsaturated bonds and high bilayer curvature on the dynamic behaviour of the lipid molecules were studied. Our results, combined with results from planar multibilayers with a varying hydration rate (Korstanje et al. (1989) Biochim. Biophys. Acta 980, 225-233), give a consistent picture of the orientational order and rotational dynamics of CSL molecules embedded in lipid matrices with various geometrical configurations. Increase of hydration or temperature reduces molecular ordering and increases molecular dynamics. In highly curved vesicle configurations, SUVs, molecular order is found to be lower than in multilamellar liposomes. In contrast, rotational motion is not affected by increase of curvature. In all lipid configurations studied, increase of the number of unsaturated bonds in the fatty acid chains reduces molecular ordering. We find, however, no effect of unsaturation on the rotational mobility of the CSL probe molecules. These results clearly show that changes in molecular orientational order and reorientational dynamics have to be considered separately, and that they are not necessarily correlated as implied by the common concept of membrane fluidity. Comparing our results with data from a motional narrowing analysis shows that the latter approach seriously overestimates the rate of molecular reorientation.     

Structure and dynamics of the dilauroylphosphatidylethanolamine lipid bilayer.

A 200-ps molecular dynamics (MD) simulation trajectory of a model dilauroylphosphatidylethanolamine (DLPE) bilayer in water at 315 K has been generated. Segmental order parameters, electron density profiles, and water pair distribution functions have been calculated. Comparison to experiment is made where possible. The dynamics of the system has been studied by analyzing the velocity autocorrelation functions (VAF) of both water and lipid atoms. Furthermore, the diffusive properties of water have been analyzed by computing the mean square displacement (MSD) and orientational correlation function (OCF) of water in two regions around the bilayer. The calculated order parameters show a behavior similar to the liquid crystalline phase of other bilayers, but the region around C1-C3 does not show the expected behavior. The electron density profile shows features that are characteristic of the liquid crystalline phase. The radial distribution functions suggest ordering of water near the charged head groups, which results in about 15 water molecules solvating each lipid molecule. We find from the VAF, MSD, and OCF calculation that the water molecules near the head groups of the lipid bilayer move more slowly than those further away. The VAF of the hydrocarbon chains have features of low-frequency motions that are probably cooperative nature in addition to the high-frequency motions associated with bond angle and torsional motions.     

Membrane topology of a 14-mer model amphipathic peptide: a solid-state NMR spectroscopy study

We have investigated the interaction between a synthetic amphipathic 14-mer peptide and model membranes by solid-state NMR. The 14-mer peptide is composed of leucines and phenylalanines modified by the addition of crown ethers and forms a helical amphipathic structure in solution and bound to lipid membranes. To shed light on its membrane topology, 31P, 2H, 15N solid-state NMR experiments have been performed on the 14-mer peptide in interaction with mechanically oriented bilayers of dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), and dipalmitoylphosphatidylcholine (DPPC). The 31P, 2H, and 15N NMR results indicate that the 14-mer peptide remains at the surface of the DLPC, DMPC, and DPPC bilayers stacked between glass plates and perturbs the lipid orientation relative to the magnetic field direction. Its membrane topology is similar in DLPC and DMPC bilayers, whereas the peptide seems to be more deeply inserted in DPPC bilayers, as revealed by the greater orientational and motional disorder of the DPPC lipid headgroup and acyl chains. 15N{31P} rotational echo double resonance experiments have also been used to measure the intermolecular dipole-dipole interaction between the 14-mer peptide and the phospholipid headgroup of DMPC multilamellar vesicles, and the results indicate that the 14-mer peptide is in contact with the polar region of the DMPC lipids. On the basis of these studies, the mechanism of membrane perturbation of the 14-mer peptide is associated to the induction of a positive curvature strain induced by the peptide lying on the bilayer surface and seems to be independent of the bilayer hydrophobic thickness.     

The static and dynamic behaviour of fluorescent probe molecules in lipid bilayers

Angle-resolved fluorescence depolarization experiments were carried out on 1,6-diphenyl-1,3,5-hexatriene (DPH) and 1-[4-(trimethylammonium)phenyl]-6-phenyl-1,3,5-hexatriene (TMA-DPH) molecules embedded in multibilayers of dimyristoylphosphatidylcholine (DMPC) and palmitoyloleoylphosphatidylcholine (POPC) above their respective phase transitions. The finding that the order parameter 〈P₂〉 of the absorption moment is significantly higher than that for the emission moment for each probe is shown to arise from a tilt of the emission moment relative to the molecular symmetry axis. It is further shown that while the order parameter 〈P₂〉 is the same for both probes in DMPC bilayers, it is higher for TMA-DPH than for DPH molecules in POPC bilayers. Considerations of the order parameters 〈P₄〉, however, show that this difference can be ascribed solely to the higher fraction of DPH molecules lying with their axes parallel to the bilayer surface. Furthermore it is found that TMA-DPH molecules undergo slower reorientational motions than DPH molecules in the same bilayer system. Nevertheless the motion of both probe molecules is faster in DMPC than in POPC bilayers. The results indicate that TMA-DPH is a more useful probe than DPH in the systems investigated.     

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.     

Molecular dynamics of individual alpha-helices of bacteriorhodopsin in dimyristol phosphatidylocholine. I. Structure and dynamics.

Understanding the role of the lipid bilayer in membrane protein structure and dynamics is needed for tertiary structure determination methods. However, the molecular details are not well understood. Molecular dynamics computer calculations can provide insight into these molecular details of protein:lipid interactions. This paper reports on 10 simulations of individual alpha-helices in explicit lipid bilayers. The 10 helices were selected from the bacteriorhodopsin structure as representative alpha-helical membrane folding components. The bilayer is constructed of dimyristoyl phosphatidylcholine molecules. The only major difference between simulations is the primary sequence of the alpha-helix. The results show dramatic differences in motional behavior between alpha-helices. For example, helix A has much smaller root-mean-squared deviations than does helix D. This can be understood in terms of the presence of aromatic residues at the interface for helix A that are not present in helix D. Additional motions are possible for the helices that contain proline side chains relative to other amino acids. The results thus provide insight into the types of motion and the average structures possible for helices within the bilayer setting and demonstrate the strength of molecular simulations in providing molecular details that are not directly visualized in experiments.