A Polarizable Force Field of Dipalmitoylphosphatidylcholine Based on the Classical Drude Model for Molecular Dynamics Simulations of Lipids

Citations Over TimeTop 10% of 2013 papers

Abstract

A polarizable force field of saturated phosphatidylcholine-containing lipids based on the classical Drude oscillator model is optimized and used in molecular dynamics simulations of bilayer and monolayer membranes. The hierarchical parametrization strategy involves the optimization of parameters for small molecules representative of lipid functional groups, followed by their application in larger model compounds and full lipids. The polar headgroup is based on molecular ions tetramethyl ammonium and dimethyl phosphate, the esterified glycerol backbone is based on methyl acetate, and the aliphatic lipid hydrocarbon tails are based on linear alkanes. Parameters, optimized to best represent a collection of gas and liquid properties for these compounds, are assembled into a complete model of dipalmitoylphosphatidylcholine (DPPC) lipids that is tested against the experimental properties of bilayer and monolayer membranes. The polarizable model yields average structural properties that are in broad accord with experimental data. The area per lipid of the model is 60 Å(2), slightly smaller than the experimental value of 63 Å(2). The order parameters from nuclear magnetic resonance deuterium quadrupolar splitting measures, the electron density profile, and the monolayer dipole potential are in reasonable agreement with experimental data, and with the nonpolarizable CHARMM C36 lipid force field.

Related Papers

• → Molecular dynamics simulations of a DMSO/water mixture using the AMBER force field(2018)20 cited
• → Development and Optimization of a New Force Field for Flexible Aluminosilicates, Enabling Fast Molecular Dynamics Simulations on Parallel Architectures(2012)28 cited
• → Secondary structure dependence of amyloid‐β(1–40) on simulation techniques and force field parameters(2021)24 cited
• → Molecular dynamics simulation of the dipalmitoylphosphatidylcholine (DPPC) lipid bilayer in the fluid phase using the Nosé-Parrinello-Rahman NPT ensemble(1995)37 cited
• → Molecular dynamics simulations of an apoliprotein A–I derived peptide in explicit water(2008)7 cited