Modification of the CHARMM force field for DMPC lipid bilayer

The CHARMM force field for DMPC lipids was modified in order to improve agreement with experiment for a number of important properties of hydrated lipid bilayer. The modification consists in introduction of a scaling factor 0.83 for 1–4 electrostatic interactions (between atoms separated by three covalent bonds), which provides correct transgauche ratio in the alkane tails, and recalculation of the headgroup charges on the basis of HF/6‐311(d,p) ab‐initio computations. Both rigid TIP3P and flexible SPC water models were used with the new lipid model, showing similar results. The new model in a 75 ns simulation has shown a correct value of the area per lipid at zero surface tension, as well as good agreement with the experiment for the electron density, structure factor, and order parameters, including those in the headgroup part of lipids. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2008

[1]  Ilpo Vattulainen,et al.  Lessons of slicing membranes: interplay of packing, free area, and lateral diffusion in phospholipid/cholesterol bilayers. , 2004, Biophysical journal.

[2]  Jeffery B. Klauda,et al.  Simulation-based methods for interpreting x-ray data from lipid bilayers. , 2006, Biophysical journal.

[3]  D. J. Price,et al.  A modified TIP3P water potential for simulation with Ewald summation. , 2004, The Journal of chemical physics.

[4]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[5]  Stefan Boresch,et al.  RATIONALIZATION OF THE DIELECTRIC PROPERTIES OF COMMON THREE-SITE WATER MODELS IN TERMS OF THEIR FORCE FIELD PARAMETERS , 1998 .

[6]  J. Callis,et al.  Conformation of the hydrocarbon chains of sodium dodecyl sulfate molecules in micelles: an FTIR study , 1989 .

[7]  P. Hobza,et al.  Lipid hydration: headgroup CH moieties are involved in water binding. , 2004, Biopolymers.

[8]  M. Laguerre,et al.  Detailed structure and dynamics of bicelle phospholipids using selectively deuterated and perdeuterated labels. 2H NMR and molecular mechanics study , 2003 .

[9]  Alexander D. MacKerell,et al.  An all-atom empirical energy function for the simulation of nucleic acids , 1995 .

[10]  E. Lindahl,et al.  Mesoscopic undulations and thickness fluctuations in lipid bilayers from molecular dynamics simulations. , 2000, Biophysical journal.

[11]  Ulf R. Pedersen,et al.  Molecular packing in 1-hexanol-DMPC bilayers studied by molecular dynamics simulation. , 2007, Biophysical chemistry.

[12]  P. P. Ewald.  Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[13]  Siewert J. Marrink,et al.  Methodological issues in lipid bilayer simulations , 2003 .

[14]  H. Casal,et al.  Quantitative determination of hydrocarbon chain conformational order in bilayers of saturated phosphatidylcholines of various chain lengths by Fourier transform infrared spectroscopy. , 1990, Biochemistry.

[15]  L. Sobczyk,et al.  Complexation equilibria in systems composed of phenols and oxygen bases in carbon tetrachloride , 1985 .

[16]  Herman J. C. Berendsen,et al.  Simulation of Water Transport through a Lipid Membrane , 1994 .

[17]  Rahman,et al.  Molecular-dynamics study of atomic motions in water. , 1985, Physical review. B, Condensed matter.

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

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

[20]  P. Kollman,et al.  Application of RESP charges to calculate conformational energies, hydrogen bond energies, and free energies of solvation , 1993 .

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

[22]  Ryan W. Benz,et al.  Experimental validation of molecular dynamics simulations of lipid bilayers: a new approach. , 2005, Biophysical journal.

[23]  J. Nagle,et al.  Structure of fully hydrated fluid phase DMPC and DLPC lipid bilayers using X-ray scattering from oriented multilamellar arrays and from unilamellar vesicles. , 2005, Biophysical journal.

[24]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[25]  W. V. van Gunsteren,et al.  Molecular dynamics simulations of phospholipid bilayers: Influence of artificial periodicity, system size, and simulation time. , 2005, The journal of physical chemistry. B.

[26]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

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

[28]  Carl-Johan Högberg,et al.  A molecular dynamics investigation of the influence of hydration and temperature on structural and dynamical properties of a dimyristoylphosphatidylcholine bilayer. , 2006, The journal of physical chemistry. B.

[29]  Alexei M. Nikitin,et al.  New six‐site acetonitrile model for simulations of liquid acetonitrile and its aqueous mixtures , 2007, J. Comput. Chem..

[30]  R. Griffin,et al.  Dipolar Recoupling in MAS NMR: A Probe for Segmental Order in Lipid Bilayers , 1997 .

[31]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[32]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[33]  M. Hyvönen,et al.  Molecular dynamics simulations of unsaturated lipid bilayers: effects of varying the numbers of double bonds , 2005, European Biophysics Journal.

[34]  Alexander P. Lyubartsev,et al.  M.DynaMix – a scalable portable parallel MD simulation package for arbitrary molecular mixtures , 2000 .

[35]  M. Jensen,et al.  Reparameterization of all-atom dipalmitoylphosphatidylcholine lipid parameters enables simulation of fluid bilayers at zero tension. , 2007, Biophysical journal.

[36]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[37]  G. Peters,et al.  Methodological problems in pressure profile calculations for lipid bilayers. , 2005, The Journal of chemical physics.

[38]  R. C. Weast.  CRC Handbook of Chemistry and Physics , 1973 .

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

[40]  Alexander D. MacKerell,et al.  An ab initio study on the torsional surface of alkanes and its effect on molecular simulations of alkanes and a DPPC bilayer. , 2005, The journal of physical chemistry. B.

[41]  S. Yip,et al.  Molecular dynamics simulation of dielectric properties of water , 1987 .

[42]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[43]  O. Edholm,et al.  The range and shielding of dipole-dipole interactions in phospholipid bilayers. , 2004, Biophysical journal.

[44]  R. Mendelsohn,et al.  Quantitative determination of conformational disorder in the acyl chains of phospholipid bilayers by infrared spectroscopy. , 1989, Biochemistry.