Drude Model Based Polarizable Force Field for Lipids

Abstract

Incorporation of electronic polarization of atoms induced by local electric fields in the empirical potential function used to generate molecular dynamics simulations is important for accurate studies of various biological processes such as ion transfer, binding and selectivity. In this work, we present new polarizable force fields for the simulations of zwitterionic lipids membranes based on the classical Drude oscillator model. The lipids considered include: dipalmitoyl phosphatidylcholine (DPPC), dipalmitoyl phosphatidylethanolamine (DPPE), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC).The parameterization strategy relies on the properties of small model compound analogs of functional groups in the molecules. The polar head group parameters are based on the molecular ions dimethyl phosphate and tetramethyl ammonium, methyl acetate for the glycerol linker region and linear alkanes as well as alkenes for the aliphatic tail of the lipid. The small molecule parameters are optimized and validated against a suite of experimental properties. The resulting lipid force field is applied to the simulation of lipid bilayers and monolayers for each lipid. For the lipid bilayer, the area per lipid, deuterium order parameters and electron density profile, and for the monolayers, the surface pressure and dipole potentials are all found to be in good agreement with experimental data when available. Thus, we have successfully validated our force field for a number of lipids.