Evaluating polarizable biomembrane simulations against experiments
Preprint
Abstract
Owing to the increase of available computational capabilities and the potential for providing more accurate description, polarizable molecular dynamics force fields are gaining popularity in modelling biomolecular systems. It is, however, crucial to evaluate how much precision is truly gained with the increased cost and complexity of the simulation. Here, we leverage the NMRlipids open collaboration and Databank to assess the performance of available polarizable lipid models—the CHARMM-Drude and the AMOEBA-based parameters—against high-fidelity experimental data and compare them to the top-performing non-polarizable models. While some improvement in the description of ion binding to membranes is observed in the most recent CHARMM-Drude parameters, and the conformational dynamics of AMOEBA-based parameters are excellent, the best non-polarizable models tend to outperform their polarizable counterparts for each property we explored. The identified shortcomings range from inaccuracies in describing the conformational space of lipids to excessively slow conformational dynamics. Our results provide valuable insights for further refinement of polarizable lipid force fields and for selecting the best simulation parameters for specific applications.Keywords
All Related Versions
- Published version: Journal of Chemical Theory and Computation, 20 (10), 4325.
This publication has 97 references indexed in Scilit:
- Polarizable Atomic Multipole-Based AMOEBA Force Field for ProteinsJournal of Chemical Theory and Computation, 2013
- Ion selectivity in channels and transportersJournal of General Physiology, 2011
- MDAnalysis: A toolkit for the analysis of molecular dynamics simulationsJournal of Computational Chemistry, 2011
- Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid TypesThe Journal of Physical Chemistry B, 2010
- Current Status of the AMOEBA Polarizable Force FieldThe Journal of Physical Chemistry B, 2010
- CHARMM-GUI Membrane Builder for Mixed Bilayers and Its Application to Yeast MembranesBiophysical Journal, 2009
- CHARMM: The biomolecular simulation programJournal of Computational Chemistry, 2009
- Lipid Bilayer Structure Determined by the Simultaneous Analysis of Neutron and X-Ray Scattering DataBiophysical Journal, 2008
- Rotation of Lipids in Membranes: Molecular Dynamics Simulation, 31P Spin-Lattice Relaxation, and Rigid-Body DynamicsBiophysical Journal, 2008
- Automated Builder and Database of Protein/Membrane Complexes for Molecular Dynamics SimulationsPLOS ONE, 2007