Structure of POPC lipid bilayers in OPLS3e force field

12 April 2022, Version 1
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

It is crucial for molecular dynamics simulations of biomembranes that the force field parameters give a realistic model of the membrane behavior. In this study, we examined the OPLS3e force field for the carbon–hydrogen order parameters SCH of POPC (1-palmitoyl-2-oleylphosphatidylcholine) lipid bilayers at varying hydration conditions and ion concentrations. The results show that OPLS3e behaves similarly to the CHARMM36 force field and relatively accurately follows the experimentally measured SCH for the lipid headgroup, the glycerol backbone, and the acyl tails. Thus, OPLS3e is a good choice for simulations of most membrane systems. The exception are systems with an abundancy of ions, as similarly to most other force fields OPLS3e strongly overestimates the membrane-binding of cations, especially Ca2+. This leads to undesirable positive charge of the membrane surface and drastically lowers the concentration of Ca2+ in the surrounding solvent, which might cause issues in systems sensitive to correct charge distribution profiles across the membrane.

Keywords

molecular dynamics simulations
force field
NMR
biomembrane structure
lipid bilayer

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.