Determination of component volumes of lipid bilayers from simulations.

An efficient method for extracting volumetric data from simulations is developed. The method is illustrated using a recent atomic-level molecular dynamics simulation of L alpha phase 1,2-dipalmitoyl-sn-glycero-3-phosphocholine bilayer. Results from this simulation are obtained for the volumes of water (VW), lipid (V1), chain methylenes (V2), chain terminal methyls (V3), and lipid headgroups (VH), including separate volumes for carboxyl (Vcoo), glyceryl (Vgl), phosphoryl (VPO4), and choline (Vchol) groups. The method assumes only that each group has the same average volume regardless of its location in the bilayer, and this assumption is then tested with the current simulation. The volumes obtained agree well with the values VW and VL that have been obtained directly from experiment, as well as with the volumes VH, V2, and V3 that require certain assumptions in addition to the experimental data. This method should help to support and refine some assumptions that are necessary when interpreting experimental data.     

Mesoscopic simulation of cell membrane damage, morphology change and rupture by nonionic surfactants.

A new simulation method, dissipative particle dynamics, is applied to model biological membranes. In this method, several atoms are united into a single simulation particle. The solubility and compressibility of the various liquid components are reproduced by the simulation model. When applied to a bilayer of phosphatidylethanolamine, the membrane structure obtained matches quantitatively with full atomistic simulations and with experiments reported in the literature. The method is applied to investigate the cause of cell death when bacteria are exposed to nonionic surfactants. Mixed bilayers of lipid and nonionic surfactant were studied, and the diffusion of water through the bilayer was monitored. Small transient holes are seen to appear at 40% mole-fraction C(9)E(8), which become permanent holes between 60 and 70% surfactant. When C(12)E(6) is applied, permanent holes only arise at 90% mole-fraction surfactant. Some simulations have been carried out to determine the rupture properties of mixed bilayers of phosphatidylethanolamine and C(12)E(6). These simulations indicate that the area of a pure lipid bilayer can be increased by a factor 2. The inclusion of surfactant considerably reduces both the extensibility and the maximum stress that the bilayer can withstand. This may explain why dividing cells are more at risk than static cells.     

Diffraction-based density restraints for membrane and membrane-peptide molecular dynamics simulations

We have recently shown that current molecular dynamics (MD) atomic force fields are not yet able to produce lipid bilayer structures that agree with experimentally-determined structures within experimental errors. Because of the many advantages offered by experimentally validated simulations, we have developed a novel restraint method for membrane MD simulations that uses experimental diffraction data. The restraints, introduced into the MD force field, act upon specified groups of atoms to restrain their mean positions and widths to values determined experimentally. The method was first tested using a simple liquid argon system, and then applied to a neat dioleoylphosphatidylcholine (DOPC) bilayer at 66% relative humidity and to the same bilayer containing the peptide melittin. Application of experiment-based restraints to the transbilayer double-bond and water distributions of neat DOPC bilayers led to distributions that agreed with the experimental values. Based upon the experimental structure, the restraints improved the simulated structure in some regions while introducing larger differences in others, as might be expected from imperfect force fields. For the DOPC-melittin system, the experimental transbilayer distribution of melittin was used as a restraint. The addition of the peptide caused perturbations of the simulated bilayer structure, but which were larger than observed experimentally. The melittin distribution of the simulation could be fit accurately to a Gaussian with parameters close to the observed ones, indicating that the restraints can be used to produce an ensemble of membrane-bound peptide conformations that are consistent with experiments. Such ensembles pave the way for understanding peptide-bilayer interactions at the atomic level.     

Eutectic interactions between saturated and unsaturated chain cholesteryl esters: comparison of calculated and observed phase diagrams

Binary phase behavior of saturated chain with unsaturated chain cholesteryl esters is evaluated by analysis of the phase diagrams in terms of ideal solution theory. Cholesteryl palmitate, which crystallizes in the bilayer structure, forms a eutectic with either cholesteryl oleate or cholesteryl linoleate and, as indicated by low angle X-ray data, the components are nearly totally fractionated in the solid state. The fit of the two experimental liquidus curves by a calculation of freezing point depression for an ideal solution indicates that the molecular interactions are nonspecific in the binary liquid state. Cholesteryl caprylate and cholesteryl oleate, both of which crystallize as the monolayer II form, also form a eutectic. X-ray data again indicate nearly total fractionation. The liquidus curve is reasonably well matched by calculation of ideal freezing point depression. However, dissimilar molecular volumes can cause the melt-cholesteric transition line to deviate from an ideal concentration dependence. Possible fractionation mechanisms for cholesteryl esters in arterial lesions are thereby indicated. For example, when the molecules have greatly different volumes, clustering can occur in the liquid crystalline state. Even when the molecular volumes are similar, the saturated component can solidify in regions where it is relatively abundant, because of the incompatibility of two crystal structures with greatly different layer structures.     

1-Alkanols and membranes: A story of attraction

Although 1-alkanols have long been known to act as penetration enhancers and anesthetics, the mode of operation is not yet understood. In this study, long-time molecular dynamics simulations have been performed to investigate the effect of 1-alkanols of various carbon chain lengths onto the structure and dynamics of dimyristoylphosphatidylcholine bilayers. The simulations were complemented by microcalorimetry, continuous bleaching and film balance experiments. In the simulations, all investigated 1-alkanols assembled inside the lipid bilayer within tens of nanoseconds. Their hydroxyl groups bound preferentially to the lipid carbonyl group and the hydrocarbon chains stretched into the hydrophobic core of the bilayer. Both molecular dynamics simulations and experiments showed that all 1-alkanols drastically affected the bilayer properties. Insertion of long-chain 1-alkanols decreased the area per lipid while increasing the thickness of the bilayer and the order of the lipids. The bilayer elasticity was reduced and the diffusive motion of the lipids within the bilayer plane was suppressed. On the other hand, integration of ethanol into the bilayer enlarged the area per lipid. The bilayer became softer and lipid diffusion was enhanced.     

1-Alkanols and membranes: a story of attraction

Although 1-alkanols have long been known to act as penetration enhancers and anesthetics, the mode of operation is not yet understood. In this study, long-time molecular dynamics simulations have been performed to investigate the effect of 1-alkanols of various carbon chain lengths onto the structure and dynamics of dimyristoylphosphatidylcholine bilayers. The simulations were complemented by microcalorimetry, continuous bleaching and film balance experiments. In the simulations, all investigated 1-alkanols assembled inside the lipid bilayer within tens of nanoseconds. Their hydroxyl groups bound preferentially to the lipid carbonyl group and the hydrocarbon chains stretched into the hydrophobic core of the bilayer. Both molecular dynamics simulations and experiments showed that all 1-alkanols drastically affected the bilayer properties. Insertion of long-chain 1-alkanols decreased the area per lipid while increasing the thickness of the bilayer and the order of the lipids. The bilayer elasticity was reduced and the diffusive motion of the lipids within the bilayer plane was suppressed. On the other hand, integration of ethanol into the bilayer enlarged the area per lipid. The bilayer became softer and lipid diffusion was enhanced.     

Structure determination of membrane proteins in five easy pieces

Rotational Alignment (RA) solid-state NMR provides the basis for a general method for determining the structures of membrane proteins in phospholipid bilayers under physiological conditions. Membrane proteins are high priority targets for structure determination, and are challenging for existing experimental methods. Because membrane proteins reside in liquid crystalline phospholipid bilayer membranes it is important to study them in this type of environment. The RA solid-state NMR approach we have developed can be summarized in five steps, and incorporates methods of molecular biology, biochemistry, sample preparation, the implementation of NMR experiments, and structure calculations. It relies on solid-state NMR spectroscopy to obtain high-resolution spectra and residue-specific structural restraints for membrane proteins that undergo rotational diffusion around the membrane normal, but whose mobility is otherwise restricted by interactions with the membrane phospholipids. High resolution spectra of membrane proteins alone and in complex with other proteins and ligands set the stage for structure determination and functional studies of these proteins in their native, functional environment.     

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.     

The water permeability of the monoolein/triolein bilayer membrane

The water permeability of the lipid bilayer can be used as a probe of membrane structure. A simple model of the bilayer, the liquid hydrocarbon model, views the membrane as a thin slice of bulk hydrocarbon liquid. A previous study (Petersen, D. (1980) Biochim. Biophys. Acta 600, 666–677) showed that this model does not accurately predict the water permeability of the monoolein/n-hexadecane bilayer: the measured activation energy for water permeation is 50% above the predicted value. From this it was inferred that the hydrocarbon chains in the lipid bilayer are more ordered than in the bulk hydrocarbon liquid. The present study tests the liquid hydrocarbon model for the monoolein/triolein bilayer, which has been shown to contain very little triolein in the plane of the membrane (Waldbillig, R.C. and Szabo, G. (1979) Biochim. Biophys. Acta 557, 295–305). Measurements of the water permeability coefficient of the bilayer are compared with predictions of the liquid hydrocarbon model based on measurements of the water permeability coefficient of bulk 8-heptadecene. The predicted and measured values agree quite closely over the temperature range studied (15–35°C): the predicted activation energy is 11.1±0.2 kcal/mol, whereas the measured activation energy for the bilayer is 9.8±0.7 kcal/mol. This close agreement is in contrast with the monoolein/n-hexadecane results and suggests that, insofar as water permeation is concerned, the liquid hydrocarbon model quite closely represents the monoolein/triolein bilayer.     

Molecular dynamics simulation of lipid packing and mobility in bilayer membranes

A model of lipid bilayer membrane in water has been developed. Parameters have been selected that allow molecular dynamics simulation of lipid bilayers in the all-atom approximation. The calculated indices of packing and mobility of lipid molecules for the liquid crystalline state of the bilayer agree well with the experimental data. Based on the model of the liquid crystalline state of the membrane, a system in the gel-like state has been constructed. The gel-state model reproduces well the packing of lipids in real bilayers, whereas the mobility of molecules proves to be overestimated.     

Monolayers at the oil/water interface as a proper model for bilayer membranes

In order to clarify the structural relationship between lipid monolayer and bilayer membranes, physical states of these membranes are discussed from their energetic points of view. It is concluded that the monolayer formed at the oil/water interface is a proper model system to represent the physical state of half of a bilayer in its liquid crystalline state. The theoretical prediction is confirmed by the monolayer surface tension measurements and the bilayer conductance experiments with water soluble (extrinsic) proteins. It is also deduced that the surface pressure of the bilayer in the liquid crystalline state is quite high, about 45 dyn/cm, and the interaction of cytochrome c with the bilayer is mainly electrostatic at the bilayer membrane periphery.     

Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction data. II. Distribution and packing of terminal methyl groups.

We continue in this paper the presentation of theoretical and experimental methods for the joint refinement of neutron and x-ray lamellar diffraction data for the analysis of fluid (L alpha phase) bilayer structure (Wiener, M. C., and S. H. White. 1991 a, b, c. Biophys. J. 59:162-173 and 174-185; Biochemistry. 30:6997-7008; Wiener, M. C., G. I. King, and S. H. White. Biophys. J. 60: 568-576). We show how to obtain the distribution and packing of the terminal methyls in the interior of a fluid dioleoylphosphatidylcholine bilayer (66% RH) by combining x-ray and neutron scattering-length transbilayer profiles with no a priori assumptions about the functional form of the distribution. We find that the methyls can be represented by a Gaussian function with 1/e-halfwidth of 2.95 +/- 0.28 A situated at the bilayer center. There is substantial mixing of the methyls and methylenes in the bilayer center. The Gaussian representation of the methyl distribution is narrower and has a different shape than predicted by several simulations of fluid bilayers (Gruen, D. W. R., and E. H. B. de Lacey. 1984. Surfactants in Solution, Vol. 1. Plenum Publishing Corp., New York. 279-306; de Loof, H., et al. 1991. Biochemistry. 30:2099-2133) but this may be due to the smaller area/lipid of our experiments and the presence of the double-bonds. Determination of the absolute specific volume of DOPC and an analysis of bulk alkane volumetric data over a range of hydrostatic pressures lead to estimates of methylene and methyl volumes at the bilayer center of 27 +/- 1 A3 and 57.2 +/- 3.6 A3, respectively. This result provides direct confirmation of the common assumption that the molecular packing of methyl and methylene groups in bilayers is the same as in bulk liquid alkanes.     

Structure, topology, and tilt of cell-signaling peptides containing nuclear localization sequences in membrane bilayers determined by solid-state NMR and molecular dynamics simulation studies

Cell-signaling peptides have been extensively used to transport functional molecules across the plasma membrane into living cells. These peptides consist of a hydrophobic sequence and a cationic nuclear localization sequence (NLS). It has been assumed that the hydrophobic region penetrates the hydrophobic lipid bilayer and delivers the NLS inside the cell. To better understand the transport mechanism of these peptides, in this study, we investigated the structure, orientation, tilt of the peptide relative to the bilayer normal, and the membrane interaction of two cell-signaling peptides, SA and SKP. Results from CD and solid-state NMR experiments combined with molecular dynamics simulations suggest that the hydrophobic region is helical and has a transmembrane orientation with the helical axis tilted away from the bilayer normal. The influence of the hydrophobic mismatch, between the hydrophobic length of the peptide and the hydrophobic thickness of the bilayer, on the tilt angle of the peptides was investigated using thicker POPC and thinner DMPC bilayers. NMR experiments showed that the hydrophobic domain of each peptide has a tilt angle of 15 +/- 3 degrees in POPC, whereas in DMPC, 25 +/- 3 degree and 30 +/- 3 degree tilts were observed for SA and SKP peptides, respectively. These results are in good agreement with molecular dynamics simulations, which predict a tilt angle of 13.3 degrees (SA in POPC), 16.4 degrees (SKP in POPC), 22.3 degrees (SA in DMPC), and 31.7 degrees (SKP in DMPC). These results and simulations on the hydrophobic fragment of SA or SKP suggest that the tilt of helices increases with a decrease in bilayer thickness without changing the phase, order, and structure of the lipid bilayers.     

A thermodynamic study of the effects of cholesterol on the interaction between liposomes and ethanol

The association of ethanol with unilamellar dimyristoyl phosphatidylcholine (DMPC) liposomes of varying cholesterol content has been investigated by isothermal titration calorimetry over a wide temperature range (8-45 degrees C). The calorimetric data show that the interaction of ethanol with the lipid membranes is endothermic and strongly dependent on the phase behavior of the mixed lipid bilayer, specifically whether the lipid bilayer is in the solid ordered (so), liquid disordered (ld), or liquid ordered (lo) phase. In the low concentration regime (<10 mol%), cholesterol enhances the affinity of ethanol for the lipid bilayer compared to pure DMPC bilayers, whereas higher levels of cholesterol (>10 mol%) reduce affinity of ethanol for the lipid bilayer. Moreover, the experimental data reveal that the affinity of ethanol for the DMPC bilayers containing small amounts of cholesterol is enhanced in the region around the main phase transition. The results suggest the existence of a close relationship between the physical structure of the lipid bilayer and the association of ethanol with the bilayer. In particular, the existence of dynamically coexisting domains of gel and fluid lipids in the transition temperature region may play an important role for association of ethanol with the lipid bilayers. Finally, the relation between cholesterol content and the affinity of ethanol for the lipid bilayer provides some support for the in vivo observation that cholesterol acts as a natural antagonist against alcohol intoxication.     

Experimental validation of molecular dynamics simulations of lipid bilayers: a new approach

A novel protocol has been developed for comparing the structural properties of lipid bilayers determined by simulation with those determined by diffraction experiments, which makes it possible to test critically the ability of molecular dynamics simulations to reproduce experimental data. This model-independent method consists of analyzing data from molecular dynamics bilayer simulations in the same way as experimental data by determining the structure factors of the system and, via Fourier reconstruction, the overall transbilayer scattering-density profiles. Multi-nanosecond molecular dynamics simulations of a dioleoylphosphatidylcholine bilayer at 66% RH (5.4 waters/lipid) were performed in the constant pressure and temperature ensemble using the united-atom GROMACS and the all-atom CHARMM22/27 force fields with the GROMACS and NAMD software packages, respectively. The quality of the simulated bilayer structures was evaluated by comparing simulation with experimental results for bilayer thickness, area/lipid, individual molecular-component distributions, continuous and discrete structure factors, and overall scattering-density profiles. Neither the GROMACS nor the CHARMM22/27 simulations reproduced experimental data within experimental error. The widths of the simulated terminal methyl distributions showed a particularly strong disagreement with the experimentally observed distributions. A comparison of the older CHARMM22 with the newer CHARMM27 force fields shows that significant progress is being made in the development of atomic force fields for describing lipid bilayer systems empirically.     

Structural properties of a highly polyunsaturated lipid bilayer from molecular dynamics simulations.

The structure of a fully hydrated mixed (saturated/polyunsaturated) chain lipid bilayer in the biologically relevant liquid crystalline phase has been examined by performing a molecular dynamics study. The model membrane, a 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphocholine (SDPC, 18:0/22:6 PC) lipid bilayer, was investigated at constant (room) temperature and (ambient) pressure, and the results obtained in the nanosecond time scale reproduced quite well the available experimental data. Polyunsaturated fatty acids are found in high concentrations in neuronal and retinal tissues and are essential for the development of human brain function. The docosahexaenoic fatty acid, in particular, is fundamental for the proper function of the visual receptor rhodopsin. The lipid bilayer order has been investigated through the orientational order parameters. The water-lipid interface has been explored thoroughly in terms of its dimensions and the organization of the different components. Several types of interactions occurring in the system have been analyzed, specifically, the water-hydrocarbon chain, lipid-lipid and lipid-water interactions. The distribution of dihedral angles along the chains and the molecular conformations of the polyunsaturated chain of the lipids have also been studied. Special attention has been focused on the microscopic (molecular) origin of the effects of polyunsaturations on the different physical properties of membranes.     

Temperature dependence of the repulsive pressure between phosphatidylcholine bilayers.

Bilayer structure and interbilayer repulsive pressure were measured from 5 to 50 degrees C by the osmotic stress/x-ray diffraction method for both gel and liquid crystalline phase lipid bilayers. For gel phase dibehenoylphosphatidylcholine (DBPC) the bilayer thickness and pressure-distance relations were nearly temperature-independent, and at full hydration the equilibrium fluid spacing increased approximately 1 A, from 10 A at 5 degrees C to 11 A at 50 degrees C. In contrast, for liquid crystalline phase egg phosphatidylcholine (EPC), the bilayer thickness, equilibrium fluid spacing, and pressure-distance relation were all markedly temperature-dependent. As the temperature was increased from 5 to 50 degrees C the EPC bilayer thickness decreased approximately 4 A, and the equilibrium fluid spacing increased from 14 to 21 A. Over this temperature range there was little change in the pressure-distance relation for fluid spacings less than approximately 10 A, but a substantial increase in the total pressure for fluid spacings greater than 10 A. These data show that for both gel and liquid crystalline bilayers there is a short-range repulsive pressure that is nearly temperature-independent, whereas for liquid crystalline bilayers there is also a longer-range pressure that increases with temperature. From analysis of the energetics of dehydration we argue that the temperature-independent short-range pressure is consistent with a hydration pressure due to polarization or electrostriction of water molecules by the phosphorylcholine moiety. For the liquid crystalline phase, the 7 A increase in equilibrium fluid spacing with increasing temperature can be predicted by an increase in the undulation pressure as a consequence of a temperature-dependent decrease in bilayer bending modulus.     

Interactions of liquid crystal-forming molecules with phospholipid bilayers studied by molecular dynamics simulations

Recent experiments have shown that liquid crystals can be used to image mammalian cell membranes and to amplify structural reorganization in phospholipid-laden liquid crystal-aqueous interfaces. In this work, molecular dynamics simulations were employed to explore the interactions between commonly used liquid crystal-forming molecules and phospholipid bilayers. In particular, umbrella sampling was used to obtain the potential of mean force of 4-cyano-4'-pentylbiphenyl (5CB) and 4'-(3,4-difluor-phenyl)-4-pentyl-bicylohexyl (5CF) molecules partitioning into a dipalmitoylphosphatidylcholine bilayer. In addition, results of simulations are presented for systems consisting of a fully hydrated bilayer with 5CB or 5CF molecules at the lowest (4.5 mol %) and highest (20 mol %) concentrations used in recent laboratory experiments. It is found that mesogens preferentially partition from the aqueous phase into the membrane; the potential of mean force exhibits highly favorable free energy differences for partitioning (-18 k(B)T for 5CB and -26 k(B)T for 5CF). The location and orientation of mesogens associated with the most stable free energies in umbrella sampling simulations of dilute systems were found to be consistent with those observed in liquid-crystal-rich bilayers. It is found that the presence of mesogens in the bilayer enhances the order of lipid acyl tails, and changes the spatial and orientational arrangement of lipid headgroup atoms. These effects are more pronounced at higher liquid-crystal concentrations. In comparing the behavior of 5CB and 5CF, a stronger spatial correlation (i.e., possibly leading to aggregation) is observed between 5CB molecules within a bilayer than between 5CF molecules. Also, the range of molecular orientations and positions along the bilayer normal is larger for 5CB molecules. At the same time, 5CF molecules were found to bind more strongly to lipid headgroups, thereby slowing the lateral motion of lipid molecules.     

蟾毒灵与脂膜相互作用的研究

应用＾１３Ｃ－ＣＰ／ＭＡＳ和ＤＳＣ方法研究蟾毒灵与磷脂膜相互作用的执致相变特性及动力学特性。ＤＳＣ曲线表明蟾毒灵使磷脂膜相变温度降低,吸热峰变宽。＾１３Ｃ－ＣＰ／ＭＡＳ谱表明磷脂膜的ＮＭＲ信号峰化学位移随温度稍有变化,提示磷脂膜在液晶态脂肪烃链有不同程序的反－旁式异构化。含蟾毒灵的ＥＰＣ脂双层ＮＭＲ谱,随温度升高有蟾毒灵信号峰出现,ＥＰＣ脂双层分子内各部分的信号峰强度和峰形变化明显,说明脂双层分子     

Molecular structure of membrane tethers

Membrane tethers are nanotubes formed by a lipid bilayer. They play important functional roles in cell biology and provide an experimental window on lipid properties. Tethers have been studied extensively in experiments and described by theoretical models, but their molecular structure remains unknown due to their small diameters and dynamic nature. We used molecular dynamics simulations to obtain molecular-level insight into tether formation. Tethers were pulled from single-component lipid bilayers by application of an external force to a lipid patch along the bilayer normal or by lateral compression of a confined bilayer. Tether development under external force proceeded by viscoelastic protrusion followed by viscous lipid flow. Weak forces below a threshold value produced only a protrusion. Larger forces led to a crossover to tether elongation, which was linear at a constant force. Under lateral compression, tethers formed from undulations of unrestrained bilayer area. We characterized in detail the tether structure and its formation process, and obtained the material properties of the membrane. To our knowledge, these results provide the first molecular view of membrane tethers.