Abstract

The recently presented LOPLS-AA all-atom force field for long hydrocarbon chains, based on the OPLS-AA force field, was extended to alcohols, esters, and glyceryl monooleate (GMO) lipids as a model lipid. Dihedral angles were fitted against high level ab initio calculations, and ester charges were increased to improve their hydration properties. Additionally, the ester Lennard-Jones parameters were readjusted to reproduce experimental liquid bulk properties, densities, and heats of vaporization. This extension enabled the setup of LOPLS-AA parameters for GMO molecules. The properties of the lipid force field were tested for the liquid-crystalline phase of a GMO bilayer. The obtained area per lipid for GMO is in good agreement with experiment. Additionally, the lipid dynamics on the subpicosecond to the nanosecond time scale is in excellent agreement with results from time-of-flight (TOF) quasielastic neutron scattering (QENS) experiments on a multilamellar monoolein system, enabling here for the first time the critical evaluation of the short-time dynamics obtained from a molecular dynamics simulation of a membrane system.