Free energy of helical transmembrane peptide dimerization in OPLS-AA/Berger force field simulations: inaccuracy and implications for partner-specific Lennard-Jones parameters between peptides and lipids

Abstract

Interactions between transmembrane (TM) peptides are important in biophysical chemistry, but there are few studies assessing atomistic simulation parameters concerning the energetics of interactions of TM helical peptides. Our potential of mean force analysis using OPLS-AA protein/Berger lipid force fields (FFs) shows that the dimerisation energy of (AALALAA)3 helical peptides in the dioleoylphosphatidylcholine bilayer is −4.4 kJ/mol, which was much smaller than the reported experimental value (−12.7 kJ/mol), thus calling for improvement of parameters of the combined FFs. As each of the FFs has been independently developed, we then tested the effects of downscaling the Lennard-Jones (LJ) energy terms between the OPLS-AA atoms and Berger lipid atoms, preserving the parameters within each FF. A 0.9-fold rescaling of the LJ energies was found to render the dimerisation energy close to the experimental value. Solvation of backbone atoms as well as side chain atoms in lipids is crucial for the TM helix interaction. In similar analyses, GROMOS 53A6 FF exhibited as weak dimerisation propensity (~−5.2 kJ/mol) as OPLS-AA/Berger, but CHARMM36 showed relatively accurate propensity (~−9.9 kJ/mol). Challenges and strategies in rendering the TM interaction energy realistic within the framework of current FFs are discussed.