Combined Monte Carlo and molecular dynamics simulation of fully hydrated dioleyl and palmitoyl-oleyl phosphatidylcholine lipid bilayers.

We have applied a new equilibration procedure for the atomic level simulation of a hydrated lipid bilayer to hydrated bilayers of dioleyl-phosphatidylcholine (DOPC) and palmitoyl-oleyl phosphatidylcholine (POPC). The procedure consists of alternating molecular dynamics trajectory calculations in a constant surface tension and temperature ensemble with configurational bias Monte Carlo moves to different regions of the configuration space of the bilayer in a constant volume and temperature ensemble. The procedure is applied to bilayers of 128 molecules of POPC with 4628 water molecules, and 128 molecules of DOPC with 4825 water molecules. Progress toward equilibration is almost three times as fast in central processing unit (CPU) time compared with a purely molecular dynamics (MD) simulation. Equilibration is complete, as judged by the lack of energy drift in 200-ps runs of continuous MD. After the equilibrium state was reached, as determined by agreement between the simulation volume per lipid molecule with experiment, continuous MD was run in an ensemble in which the lateral area was restrained to fluctuate about a mean value and a pressure of 1 atm applied normal to the bilayer surface. Three separate continuous MD runs, 200 ps in duration each, separated by 10,000 CBMC steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of POPC and DOPC. Analysis of the hydration environment in the headgroups supports a mechanism by which unsaturation contributes to reduced transition temperatures. In this view, the relatively horizontal orientation of the unsaturated bond increases the area per lipid, resulting in increased water penetration between the headgroups. As a result the headgroup-headgroup interactions are attenuated and shielded, and this contributes to the lowered transition temperature.

[1]  Jean-Paul Ryckaert,et al.  Molecular dynamics of liquid n-butane near its boiling point , 1975 .

[2]  P. Huang,et al.  Molecular dynamics simulations of phospholipid bilayers. , 1994, Journal of biomolecular structure & dynamics.

[3]  Ching-Hsien Huang,et al.  Molecular mechanics simulation studies of dienoic hydrocarbons: From alkenes to 1‐Palmitoyl‐2‐linoleoyl‐phosphatidylcholines , 1996 .

[4]  Shusen Li,et al.  Molecular mechanics simulation studies of dienoic hydrocarbons: From alkenes to 1-Palmitoyl-2-linoleoyl-phosphatidylcholines , 1996, J. Comput. Chem..

[5]  D P Tieleman,et al.  A computer perspective of membranes: molecular dynamics studies of lipid bilayer systems. , 1997, Biochimica et biophysica acta.

[6]  Terry R. Stouch,et al.  Lipid Membrane Structure and Dynamics Studied by All-Atom Molecular Dynamics Simulations of Hydrated Phospholipid Bilayers , 1993 .

[7]  S H White,et al.  Determination of the hydrocarbon core structure of fluid dioleoylphosphocholine (DOPC) bilayers by x-ray diffraction using specific bromination of the double-bonds: effect of hydration. , 1998, Biophysical journal.

[8]  R M Venable,et al.  Molecular dynamics simulations of a lipid bilayer and of hexadecane: an investigation of membrane fluidity. , 1993, Science.

[9]  K. Schulten,et al.  Molecular dynamics simulation of a bilayer of 200 lipids in the gel and in the liquid crystal phase , 1993 .

[10]  Alexander D. MacKerell,et al.  Molecular dynamics simulation of unsaturated lipid bilayers at low hydration: parameterization and comparison with diffraction studies. , 1997, Biophysical journal.

[11]  J. Nagle,et al.  Structure of fully hydrated bilayer dispersions. , 1988, Biochimica et biophysica acta.

[12]  O. Berger,et al.  Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. , 1997, Biophysical journal.

[13]  S H White,et al.  Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction data. III. Complete structure. , 1992, Biophysical journal.

[14]  Michael M. Clark,et al.  Application of combined Monte Carlo and molecular dynamics method to simulation of dipalmitoyl phosphatidylcholine lipid bilayer , 1999, J. Comput. Chem..

[15]  Daan Frenkel,et al.  Configurational bias Monte Carlo: a new sampling scheme for flexible chains , 1992 .

[16]  Benoît Roux,et al.  Biological membranes : a molecular perspective from computation and experiment , 1996 .

[17]  S. Feller,et al.  Determination of component volumes of lipid bilayers from simulations. , 1997, Biophysical journal.

[18]  J. Nagle,et al.  Structure and interactions of fully hydrated dioleoylphosphatidylcholine bilayers. , 1998, Biophysical journal.

[19]  E. Jakobsson,et al.  Application of combined Monte Carlo and molecular dynamics method to simulation of dipalmitoyl phosphatidylcholine lipid bilayer , 1999, J. Comput. Chem..

[20]  S. White,et al.  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. , 1992, Biophysical journal.

[21]  Eric Jakobsson,et al.  Optimization of Hydrocarbon Chain Interaction Parameters: Application to the Simulation of Fluid Phase Lipid Bilayers , 1999 .

[22]  E Jakobsson,et al.  Computer simulation studies of biological membranes: progress, promise and pitfalls. , 1997, Trends in biochemical sciences.

[23]  J. Seelig,et al.  Molecular order in cis and trans unsaturated phospholipid bilayers. , 1978, Biochemistry.

[24]  Richard W. Pastor,et al.  Molecular dynamics and Monte Carlo simulations of lipid bilayers , 1994 .

[25]  E. Jakobsson,et al.  Incorporation of surface tension into molecular dynamics simulation of an interface: a fluid phase lipid bilayer membrane. , 1995, Biophysical journal.

[26]  H. Berendsen,et al.  Molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine bilayer with different macroscopic boundary conditions and parameters , 1996 .

[27]  R. Pastor,et al.  Time Scales of Lipid Dynamics and Molecular Dynamics , 1996 .

[28]  K V Damodaran,et al.  A comparison of DMPC- and DLPE-based lipid bilayers. , 1994, Biophysical journal.

[29]  H. C. Andersen,et al.  A theory of the effects of head-group structure and chain unsaturation on the chain melting transition of phospholipid dispersions. , 1980, Biochemistry.

[30]  R. Pastor,et al.  Computer simulation of liquid/liquid interfaces. II. Surface tension-area dependence of a bilayer and monolayer , 1995 .

[31]  M. Klein,et al.  Constant pressure and temperature molecular dynamics simulation of a fully hydrated liquid crystal phase dipalmitoylphosphatidylcholine bilayer. , 1995, Biophysical journal.

[32]  R. Suter,et al.  X-ray structure determination of fully hydrated L alpha phase dipalmitoylphosphatidylcholine bilayers. , 1996, Biophysical journal.

[33]  M. Klein,et al.  Constant pressure and temperature molecular-dynamics simulation of the hydrated diphytanolphosphatidylcholine lipid bilayer , 1998 .

[34]  R. Armen,et al.  Phospholipid component volumes: determination and application to bilayer structure calculations. , 1998, Biophysical journal.

[35]  K. Merz,et al.  Head group-water interactions in lipid bilayers: a comparison between DMPC- and DLPE-based lipid bilayers , 1993 .

[36]  Eric Jakobsson,et al.  Simulations of lipid membranes with atomic resolution , 1998 .