Constant pressure and temperature molecular dynamics simulation of a fully hydrated liquid crystal phase dipalmitoylphosphatidylcholine bilayer.

We report a constant pressure and temperature molecular dynamics simulation of a fully hydrated liquid crystal (L alpha) phase bilayer of dipalmitoylphosphatidylcholine at 50 degrees C and 28 water molecules/lipid. We have shown that the bilayer is stable throughout the 1550-ps simulation and have demonstrated convergence of the system dimensions. Several important aspects of the bilayer structure have been investigated and compared favorably with experimental results. For example, the average positions of specific carbon atoms along the bilayer normal agree well with neutron diffraction data, and the electron density profile is in accord with x-ray diffraction results. The hydrocarbon chain deuterium order parameters agree reasonably well with NMR results for the middles of the chains, but the simulation predicts too much order at the chain ends. In spite of the deviations in the order parameters, the hydrocarbon chain packing density appears to be essentially correct, inasmuch as the area/lipid and bilayer thickness are in agreement with the most refined experimental estimates. The deuterium order parameters for the glycerol and choline groups, as well as the phosphorus chemical shift anisotropy, are in qualitative agreement with those extracted from NMR measurements.     

Constant-pressure molecular dynamics investigation of cholesterol effects in a dipalmitoylphosphatidylcholine bilayer.

We report a 1.4-ns constant-pressure molecular dynamics simulation of cholesterol at 12.5 mol% in a dipalmitoylphosphatidylcholine (DPPC) bilayer at 50 degrees C and compare the results to our previous simulation of a pure DPPC bilayer. The interlamellar spacing was increased by 2.5 A in the cholesterol-containing bilayer, consistent with x-ray diffraction results, whereas the bilayer thickness was increased by only 1 A. The bilayer/water interface was more abrupt because the lipid headgroups lie flatter to fill spaces left by the cholesterol molecules. This leads to less compensation by the lipid headgroups of the oriented water contribution to the membrane dipole potential and could explain the experimentally observed increase in the magnitude of the dipole potential by cholesterol. Our calculations suggested that 12.5 mol% cholesterol does not significantly affect the conformations and packing of the hydrocarbon chains and produces only a slight reduction in the empty free volume. However, cholesterol has a significant influence on the subnanosecond time scale lipid dynamics: the diffusion constant for the center-of-mass "rattling" motion was reduced by a factor of 3, and the reorientational motion of the methylene groups was slowed along the entire length of the hydrocarbon chains.     

Distance measurements using paramagnetic ion-induced relaxation in the saturation transfer electron spin resonance of spin-labeled biomolecules: Application to phospholipid bilayers and interdigitated gel phases

The saturation transfer electron spin resonance (STESR) spectra of spin-labeled phosphatidylcholines in gel phase lipid bilayers are shown to be sensitive to dipolar spin-spin interactions with paramagnetic ions in the aqueous phase. The reciprocal integrated intensity of the STESR spectrum is linearly dependent on aqueous Ni²⁺ ion concentration, hence, confirming the expectation that the STESR intensity is directly proportional to the spin-lattice relaxation time of the spin label. The gradient of the relaxation rate with respect to Ni²⁺ ion concentration decreases strongly with the position of the nitroxide group down the sn-2 chain of the spin-labeled lipid and is consistent with a 1/R³ dependence on the distance, R, from the bilayer surface. The values derived for the dimensions of the bilayer and lipid molecules in the case of dipalmitoyl phosphatidylcholine (DPPC) are in good agreement with those available from x-ray diffraction studies. Allowance for the multibilayer nature of the DPPC dispersions gives an estimate of the water layer thickness that is also consistent with results from x-ray diffraction. The profile of the paramagnetic ion-induced relaxation is drastically changed with DPPC dispersions in glycerol for which the lipid chains are known to be interdigitated in the gel phase. The terminal methyl groups of the lipid chains are located approximately in register with the C-3 atoms of the sn-2 chain of the oppositely oriented lipid molecules in the interdigitated phase. The thickness of the lipid layer and the effective thickness of the lipid polar group are reduced by ~40% in the interdigitated phase as compared with the bilayer phase. The calibrations of the distance dependence established by use of spin labels at defined chain positions should be applicable to STESR measurements on other biological systems.     

Relations for lipid bilayers. Connection of electron density profiles to other structural quantities.

Three relations are derived that connect low angle diffraction/scattering results obtained from lipid bilayers to other structural quantities of interest. The first relates the area along the surface of the bilayer, the measured specific volume, and the zeroth order structure factor, F(0). The second relates the size of the trough in the center of the electron density profile, the volume of the terminal methyl groups, and the volume of the methylene groups in the fatty acid chains. The third relates the size of the headgroup electron density peak, the volume of the headgroup, and the volumes of water and hydrocarbon in the headgroup region. These relations, which are easily modified for neutron diffraction, are useful for obtaining structural quantities from electron density profiles obtained by fitting model profiles to measured low angle x-ray intensities.     

On the position of the hydro-phobic/philic boundary in lipid bilayers

The sensitivity of calculated structural dimensions of hydrated lipids to the position of the hydrophobic/hydrophilic boundary is reviewed. The position of this boundary is critical in determining the extent of hydration and location of water in the bilayer. A calculation of the dimensions of the hydrophilic and hydrophobic parts of the phosphatidylcholine and ethanolamine bilayer from literature values of the x-ray long spacing shows that the choice of boundary in phospholipids is not arbitrary and is best placed at the average position of the first CH₂ group in the hydrocarbon chains. Calculated dimensions of the hydrocarbon core and the water accessible regions agree with neutron and x-ray diffraction measurements. Hydration differences between phosphatidylcholines and phosphatidylethanolamines are readily explained from derived estimates of the layers of water which cover these headgroups.     

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.     

New structural model for mixed-chain phosphatidylcholine bilayers

Multilamellar suspensions of a mixed-chain saturated phosphatidylcholine with 18 carbon atoms in the sn-1 chain and 10 carbon atoms in the sn-2 chain have been analyzed by X-ray diffraction techniques. The structural parameters for this lipid in the gel state are quite different than usual phosphatidylcholine bilayer phases. A symmetric and sharp wide-angle reflection at 4.11 A indicates that the hydrocarbon chains in hydrated C(18):C(10)PC bilayers are more tightly packed than in usual gel-state phosphatidylcholine bilayers and that there is no hydrocarbon chain tilt. The lipid thickness is about 12 A smaller than would be expected in a normal bilayer phase, and the area per molecule is 3 times the area per hydrocarbon chain. In addition, the bilayer thickness increases upon melting to the liquid-crystalline state, whereas normal bilayer phases decrease in thickness upon melting. On the basis of these data, we propose a new lipid packing model for gel-state C(18):C(10)PC bilayers in which the long C(18) chain spans the entire width of the hydrocarbon region of the bilayer and the short C(10) chain aligns or abuts with the C(10) chain from the apposing molecule. This model is novel in that there are three hydrocarbon chains per head group at the lipid-water interface. Calculations show that this phase is energetically favorable for mixed-chain lipids provided the long acyl chain is nearly twice the length of the shorter chain. In the liquid-crystalline state C(18):C(10)PC forms a normal fluid bilayer, with two chains per head group.(ABSTRACT TRUNCATED AT 250 WORDS)     

Acyl-chain methyl distributions of liquid-ordered and -disordered membranes

A central feature of the lipid raft concept is the formation of cholesterol-rich lipid domains. The introduction of relatively rigid cholesterol molecules into fluid liquid-disordered (L_d) phospholipid bilayers can produce liquid-ordered (L_o) mixtures in which the rigidity of cholesterol causes partial ordering of the flexible hydrocarbon acyl chains of the phospholipids. Several lines of evidence support this concept, but direct structural information about L_o membranes is lacking. Here we present the structure of L_o membranes formed from cholesterol and dioleoylphosphatidylcholine (DOPC). Specific deuteration of the DOPC acyl-chain methyl groups and neutron diffraction measurements reveal an extraordinary disorder of the acyl chains of neat L_d DOPC bilayers. The disorder is so great that >20% of the methyl groups are in intimate contact with water in the bilayer interface. The ordering of the DOPC acyl chains by cholesterol leads to retraction of the methyl groups away from the interface. Molecular dynamics simulations based on experimental systems reveal asymmetric transbilayer distributions of the methyl groups associated with each bilayer leaflet.     

Hydration of POPC bilayers studied by 1H-PFG-MAS-NOESY and neutron diffraction

The stability of lipid bilayers is ultimately linked to the hydrophobic effect and the properties of water of hydration. Magic angle spinning (MAS) nuclear Overhauser enhancement spectroscopy (NOESY) with application of pulsed magnetic field gradients (PFG) was used to study the interaction of water with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers in the fluid phase. NOESY cross-relaxation between water and polar groups of lipids, but also with methylene resonances of hydrophobic hydrocarbon chains, has been observed previously. This observation led to speculations that substantial amounts of water may reside in the hydrophobic core of bilayers. Here, the results of a quantitative analysis of cross-relaxation in a lipid 1-palmitoyl-2-oleoyl-sn-glycero-3 phosphocholine (POPC)/water mixture are reported. Coherences were selected via application of pulsed magnetic field gradients. This technique shortens acquisition times of NOESY spectra to 20 min and reduces t ₁-spectral noise, enabling detection of weak crosspeaks, like those between water and lipids, with higher precision than with non-gradient NOESY methods. The analysis showed that water molecules interact almost exclusively with sites of the lipid–water interface, including choline, phosphate, glycerol, and carbonyl groups. The lifetime of lipid–water associations is rather short, on the order of 100 ps, at least one order of magnitude shorter than the lifetime of lipid–lipid associations. The distribution of water molecules over the lipid bilayer was measured at identical water content by neutron diffraction. Water molecules penetrate deep into the interfacial region of bilayers but water concentration in the hydrophobic core is below the detection limit of one water molecule per lipid, in excellent agreement with the cross-relaxation data. Dedicated to Prof. K. Arnold on the occasion of his 65th birthday.     

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.     

Molecular Dynamics Simulations of Phospholipid Bilayers

Abstract

Molecular dynamics (MD) simulations at 37°C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a_t/a_h), correctly predict the aggregation preferences of the three lipids observed experimentally at 37°C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.     

Temperature and compositional dependence of the structure of hydrated dimyristoyl lecithin.

Differential scanning calorimetry and x-ray diffraction techniques have been used to investigate the structure and phase behavior of hydrated dimyristoyl lecithin (DML) in the hydration range 7.5 to 60 weight % water and the temperature range -10 to +60 degrees C. Four different calorimetric transitions have been observed: T1, a low enthalpy transition (deltaH approximately equal to 1 kcal/mol of DML) at 0 degrees C between lamellar phases (L leads to Lbeta); T2, the low enthalpy "pretransition" at water contents greater than 20 weight % corresponding to the transition Lbeta leads to Pbeta; T3, the hydrocarbon chain order-disorder transition (deltaH = 6 to 7 kcal/mol of DML) representing the transition of the more ordered low temperature phases (Lbeta, Pbeta, or crystal C, depending on the water content) to the lamellar Lalpha phase; T4, a transition occurring at 25--27 degrees C at low water contents representing the transition from the lamellar Lbeta phase to a hydrated crystalline phase C. The structures of the Lbeta, Pbeta, C, and Lalpha phases have been examined as a function of temperature and water content. The Lbeta structure has a lamellar bilayer organization with the hydrocarbon chains fully extended and tilted with respect to the normal to the bilayer plane, but packed in a distorted quasihexagonal lattice. The Pbeta structure consists of lipid bilayer lamellae distorted by a periodic "ripple" in the plane of the lamellae; the hydrocarbon chains are tilted but appear to be packed in a regular hexagonal lattice. The diffraction pattern from the crystalline phase C indexes according to an orthorhombic cell with a = 53.8 A, b = 9.33 A, c = 8.82 A. In the lamellae bilayer Lalpha strucure, the hydrocarbon chains adopt a liquid-like conformation. Analysis of the hydration characteristics and bilayer parameters (lipid thickness, surface area/molecule) of synthetic lecithins permits an evaluation of the generalized hydration and structural behavior of this class of lipids.     

Behavior of cholesterol and its effect on head group and chain conformations in lipid bilayers: a molecular dynamics study.

Cholesterol molecules were put into a computer-modeled hydrated bilayer of dimyristoyl phosphatidyl choline molecules, and molecular dynamics simulations were run to characterize the effect of this important molecule on membrane structure and dynamics. The effect was judged by observing differences in order parameters, tilt angles, and the fraction of gauche bonds along the hydrocarbon chains between lipids adjacent to cholesterol molecules and comparing them with those further away. It was observed that cholesterol causes an increase in the fraction of trans dihedrals and motional ordering of chains close to the rigid steroid ring system with a decrease in the kink population. The hydrogen-bonding interactions between cholesterol and lipid molecules were determined from radial distribution calculations and showed the cholesterol hydroxyl groups either solvated by water, or forming hydrogen bond contacts with the oxygens of lipid carbonyl and phosphate groups. The dynamics and conformation of the cholesterol molecules were investigated and it was seen that they had a smaller tilt with respect to the bilayer normal than the lipid chains and furthermore that the hydrocarbon tail of the cholesterol was conformationally flexible.     

The structure of polyunsaturated lipid bilayers important for rhodopsin function: a neutron diffraction study

The structure of oriented 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine bilayers with perdeuterated stearoyl- or docosahexaenoyl hydrocarbon chains was investigated by neutron diffraction. Experiments were conducted at two different relative humidities, 66 and 86%. At both humidities we observed that the polyunsaturated docosahexaenoyl chain has a preference to reside near the lipid water interface. That leaves voids in the bilayer center that are occupied by saturated stearoyl chain segments. This uneven distribution of saturated- and polyunsaturated chain densities is likely to result in membrane elastic stress that modulates function of integral receptor proteins like rhodopsin.     

Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature.

Molecular dynamics simulations of 500 ps were performed on a system consisting of a bilayer of 64 molecules of the lipid dipalmitoylphosphatidylcholine and 23 water molecules per lipid at an isotropic pressure of 1 atm and 50 degrees C. Special attention was devoted to reproduce the correct density of the lipid, because this quantity is known experimentally with a precision better than 1%. For this purpose, the Lennard-Jones parameters of the hydrocarbon chains were adjusted by simulating a system consisting of 128 pentadecane molecules and varying the Lennard-Jones parameters until the experimental density and heat of vaporization were obtained. With these parameters the lipid density resulted in perfect agreement with the experimental density. The orientational order parameter of the hydrocarbon chains agreed perfectly well with the experimental values, which, because of its correlation with the area per lipid, makes it possible to give a proper estimate of the area per lipid of 0.61 +/- 0.01 nm2.     

Effects of phenylamide herbicides on the physical properties of phosphatidylcholine membranes

A number of phenylamide herbicides are observed to uncouple electron transport in isolated chloroplasts and mitochondria and alter the H+ permeability of artificial liposomes. Several of these phenylamides were incorporated into phosphatidylcholine multilamellar and small unilamellar vesicles to measure their effects on the physical properties of membranes. X-ray diffraction analysis of the multilamellar vesicles revealed that the herbicides partitioned into the hydrocarbon chain region of the bilayer, but caused only minimal perturbations on hydrocarbon chain packing. 31P-NMR spectroscopy of these multilamellar vesicles showed both a broadening and lowering of the phase transition temperature of the bilayer lipids upon addition of the herbicides. 13C-NMR spectroscopy of small, unilamellar phosphatidylcholine vesicles was performed to measure the effects of the phenylamides on the chemical shifts and the spin-lattice relaxation times of the individual phosphatidylcholine carbon atoms. None of the added compounds had any measurable effect on the 13C-NMR chemical shifts of the phosphatidylcholine. However, the herbicides significantly modified spin-lattice relaxation times of certain of the lipid carbon atoms. These results generally indicate that the herbicides orient in the lipid bilayers such that the hydrocarbon chains of the phenylamides associate with the hydrocarbon chains of the lipid, whereas the phenyl moiety resides in the polar region of the bilayer.     

Membrane lateral compressibility determined by NMR and x-ray diffraction: effect of acyl chain polyunsaturation.

The elastic area compressibility modulus, Ka, of lamellar liquid crystalline bilayers was determined by a new experimental approach using 2H-NMR order parameters of lipid hydrocarbon chains together with lamellar repeat spacings measured by x-ray diffraction. The combination of NMR and x-ray techniques yields accurate determination of lateral area per lipid molecule. Samples of saturated, monounsaturated, and polyunsaturated phospholipids were equilibrated with polyethylene glycol (PEG) 20,000 solutions in water at concentrations from 0 to 55 wt % PEG at 30 degrees C. This procedure is equivalent to applying 0 to 8 dyn/cm lateral pressure to the bilayers. The resulting reductions in area per lipid were measured with a resolution of +/-0.2 A2 and the fractional area decrease was proportional to applied lateral pressure. For 1,2-dimyristoyl(d54)-sn-glycero-3-phosphocholine, 1-stearoyl(d35)-2-oleoyl-sn-glycero-3-phosphocholine (SOPC-d35), and 1-stearoyl(d35)-2-docosahexaenoyl-sn-glycero-3-phosphocholine (SDPC-d35) cross-sectional areas per molecule in excess water of 59.5, 61.4, and 69.2 A2 and bilayer elastic area compressibility moduli of 141, 221, and 121 dyn/cm were determined, respectively. Combining NMR and x-ray results enables the determination of compressibility differences between saturated and unsaturated hydrocarbon chains. In mixed-chain SOPC-d35 both chains have similar compressibility moduli; however, in mixed-chain polyunsaturated SDPC-d35, the saturated stearic acid chain appears to be far less compressible than the polyunsaturated docosahexaenoic acid chain.     

The structure of lipid bilayers and the effects of general anaesthetics: An X-ray and neutron diffraction study

We have looked for the effects of three clinically used inhalational anaesthetics (nitrous oxide, halothane and cyclopropane) on the structure of lecithin/ cholesterol bilayers. The anaesthetics were delivered to the membranes in the gaseous phase, so that effects at clinical concentrations could be determined.High-resolution X-ray diffraction patterns were recorded out to 4 Å and analyzed using swelling experiments. Parallel neutron diffraction experiments were performed and analyzed using H₂O-²H₂O exchange. Methods were developed which enabled us to obtain confidence limits for the X-ray and neutron structure factors.The resultant X-ray and neutron scattering density profiles clearly define the positions of the principal molecular groups in the unperturbed bilayer. In particular, the high-resolution electron density profiles reveal features directly attributable to the cholesterol molecule. A comparison with the neutron scattering density profiles shows that cholesterol is anchored with its hydroxyl group at the water/hydrocarbon interface, aligned with the fatty acid ester groups of the lecithin molecule. We suggest that this positioning of the cholesterol molecule allows it to act as a thickness buffer for plasma membranes.In the presence of very high concentrations of general anaesthetics, the bilayers show increased disorder while maintaining constant membrane thickness. At surgical concentrations, however, there are no significant changes in bilayer structure at 95% confidence levels. We briefly review the literature previously used to support lipid bilayer hypotheses of general anaesthesia. We conclude that the lipid bilayer per se is not the primary site of action of general anaesthetics.     

Characterization of a quasicrystalline phase in codispersions of phosphatidylethanolamine and glucocerebroside

Synchrotron x-ray diffraction, differential scanning calorimetry, and electron spin resonance spectroscopy have been employed to characterize a quasicrystalline phase formed in aqueous dispersions of binary mixtures of glucocerebroside and palmitoyloleoylphosphatidylethanolamine. Small- and wide-angle x-ray scattering intensity patterns were recorded during temperature scans between 20 degrees and 90 degrees C from mixtures of composition 2, 5, 10, 20, 30, and 40 mol glucocerebroside per 100 mol phospholipid. The quasicrystalline phase was characterized by a broad lamellar d-spacing of 6.06 nm at 40 degrees C and a broad wide-angle x-ray scattering band centered at approximately 0.438 nm, close to the gel phase centered at approximately 0.425 nm and distinct from a broad peak centered at 0.457 nm observed for a liquid-crystal phase at 80 degrees C. The quasicrystalline phase coexisted with gel and fluid phase of the pure phospholipid. An analysis of the small-angle x-ray scattering intensity profiles indicated a stoichiometry of one glucosphingolipid per two phospholipid molecules in the complex. Structural transitions monitored in cooling scans by synchrotron x-ray diffraction indicated that a cubic phase transforms initially into a lamellar gel. Thermal studies showed that the gel phase subsequently relaxes into the quasicrystalline phase in an exothermic transition. Electron spin resonance spectroscopy using spin labels located at positions 7, 12, and 16 carbons of phospholipid hydrocarbon chains indicated that order and motional constraints at the 7 and 12 positions were indistinguishable between gel and quasicrystalline phases but there was a significant decrease in order and increase in rate of motion at the 16 position on transformation to the quasicrystalline phase. The results are interpreted as an arrangement of polar groups of the complex in a crystalline array and a quasicrystalline packing of the hydrocarbon chains predicated by packing problems in the bilayer core requiring disordering of the highly asymmetric chains. The possible involvement of quasicrystalline phases in formation of membrane rafts is considered.     

Mechanism of alamethicin insertion into lipid bilayers.

Alamethicin adsorbs on the membrane surface at low peptide concentrations. However, above a critical peptide-to-lipid ratio (P/L), a fraction of the peptide molecules insert in the membrane. This critical ratio is lipid dependent. For diphytanoyl phosphatidylcholine it is about 1/40. At even higher concentrations P/L > or = 1/15, all of the alamethicin inserts into the membrane and forms well-defined pores as detected by neutron in-plane scattering. A previous x-ray diffraction measurement showed that alamethicin adsorbed on the surface has the effect of thinning the bilayer in proportion to the peptide concentration. A theoretical study showed that the energy cost of membrane thinning can indeed lead to peptide insertion. This paper extends the previous studies to the high-concentration region P/L > 1/40. X-ray diffraction shows that the bilayer thickness increases with the peptide concentration for P/L > 1/23 as the insertion approaches 100%. The thickness change with the percentage of insertion is consistent with the assumption that the hydrocarbon region of the bilayer matches the hydrophobic region of the inserted peptide. The elastic energy of a lipid bilayer including both adsorption and insertion of peptide is discussed. The Gibbs free energy is calculated as a function of P/L and the percentage of insertion phi in a simplified one-dimensional model. The model exhibits an insertion phase transition in qualitative agreement with the data. We conclude that the membrane deformation energy is the major driving force for the alamethicin insertion transition.