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.

[1]  H. G. Petersen,et al.  Error estimates on averages of correlated data , 1989 .

[2]  William H. Press,et al.  Numerical recipes , 1990 .

[3]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[4]  S. White,et al.  Fluid bilayer structure determination by the combined use of x-ray and neutron diffraction. I. Fluid bilayer models and the limits of resolution. , 1991, Biophysical journal.

[5]  B. Roux.  Commentary: surface tension of biomembranes. , 1996, Biophysical journal.

[6]  S. White,et al.  Fluid bilayer structure determination by the combined use of x-ray and neutron diffraction. II. "Composition-space" refinement method. , 1991, Biophysical journal.

[7]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

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

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

[10]  Alexander D. MacKerell,et al.  An Improved Empirical Potential Energy Function for Molecular Simulations of Phospholipids , 2000 .

[11]  Richard H. Sherman,et al.  Chaotic communications in the presence of noise , 1993, Optics & Photonics.

[12]  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.

[13]  Alexander D. MacKerell,et al.  An Empirical Potential Energy Function for Phospholipids: Criteria for Parameter Optimization and Applications , 1996 .

[14]  D. T. Cromer,et al.  X-ray scattering factors computed from numerical Hartree–Fock wave functions , 1968 .

[15]  S. White,et al.  Transbilayer distribution of bromine in fluid bilayers containing a specifically brominated analogue of dioleoylphosphatidylcholine. , 1991, Biochemistry.

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

[17]  S. Lowen.  The Biophysical Journal , 1960, Nature.

[18]  B. Roux,et al.  Molecular dynamics simulation of the gramicidin channel in a phospholipid bilayer. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[19]  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.

[20]  S H White,et al.  Structure of a fluid dioleoylphosphatidylcholine bilayer determined by joint refinement of x-ray and neutron diffraction data. I. Scaling of neutron data and the distributions of double bonds and water. , 1991, Biophysical journal.

[21]  Alexander D. MacKerell.  Molecular Dynamics Simulation Analysis of a Sodium Dodecyl Sulfate Micelle in Aqueous Solution: Decreased Fluidity of the Micelle Hydrocarbon Interior , 1995 .

[22]  R. Pastor,et al.  On simulating lipid bilayers with an applied surface tension: periodic boundary conditions and undulations. , 1996, Biophysical journal.

[23]  K. Lonsdale.  X-Ray Diffraction , 1971, Nature.

[24]  F. Jähnig,et al.  What is the surface tension of a lipid bilayer membrane? , 1996, Biophysical journal.

[25]  Klaus Schulten,et al.  Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-range Interactions , 1991 .

[26]  E Jakobsson,et al.  Molecular simulation of dioleoylphosphatidylcholine lipid bilayers at differing levels of hydration. , 2001, Biophysical journal.

[27]  E Jakobsson,et al.  Combined Monte Carlo and molecular dynamics simulation of fully hydrated dioleyl and palmitoyl-oleyl phosphatidylcholine lipid bilayers. , 1999, Biophysical journal.

[28]  J. E. Glynn,et al.  Numerical Recipes: The Art of Scientific Computing , 1989 .

[29]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[30]  T. Woolf,et al.  Interpretation of small angle X-ray measurements guided by molecular dynamics simulations of lipid bilayers. , 2003, Chemistry and physics of lipids.

[31]  Douglas J. Tobias,et al.  Constant Pressure and Temperature Molecular Dynamics Simulations of Crystals of the Lecithin Fragments: Glycerylphosphorylcholine and Dilauroylglycerol , 1995 .

[32]  J. Nagle,et al.  Structure of lipid bilayers. , 2000, Biochimica et biophysica acta.

[33]  J. Nagle,et al.  Specific volumes of lipids in fully hydrated bilayer dispersions. , 1988, Biochimica et biophysica acta.

[34]  C. Brooks.  Computer simulation of liquids , 1989 .

[35]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[36]  S. Simon,et al.  Permeability and Stability of Lipid Bilayers , 1995 .

[37]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[38]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[39]  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.

[40]  S H White,et al.  The nature of the hydrophobic binding of small peptides at the bilayer interface: implications for the insertion of transbilayer helices. , 1989, Biochemistry.

[41]  Richard W. Pastor,et al.  Computer Simulation of a DPPC Phospholipid Bilayer: Structural Changes as a Function of Molecular Surface Area , 1997 .

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

[43]  E. Eanes.  LOW-ANGLE X-RAY DIFFRACTION. , 1964, Journal. Royal Microscopical Society.

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