Abstract
AbstractOwing 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 available polarizable lipid models, CHARMM-Drude and AMOEBA-based parameters, against high-fidelity experimental data and compare them to top-performing non-polarizable models. While some improvement in 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, best non-polarizable models tend to outperformed 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.
Publisher
Cold Spring Harbor Laboratory
Reference108 articles.
1. The ugly, bad, and good stories of large-scale biomolecular simulations
2. Polarisable force fields: what do they add in biomolecular simulations?;Current Opinion in Structural Biology,2020
3. Nochebuena, J. ; Piquemal, J.-P. ; Liu, S. ; Cisneros, G. A . Cooperativity and Frustration Effects (or Lack Thereof) in Polarizable and Non-polarizable Force Fields. Journal of Chemical Theory and Computation 0, 0, null.
4. Ion selectivity in channels and transporters
5. Induced polarization in molecular dynamics simulations of the 5-HT3 receptor channel;Journal of the American Chemical Society,2020