Improved Angle Potentials for Coarse-Grained Molecular Dynamics Simulations

Monica Bulacu*, Nicolae Goga, Wei Zhao, Giulia Rossi, Luca Monticelli, Xavier Periole, D. Peter Tieleman, Siewert J. Marrink

*Corresponding author for this work

Research output: Contribution to journalArticleAcademicpeer-review

41 Citations (Scopus)
915 Downloads (Pure)

Abstract

Potentials routinely used in atomistic molecular dynamics simulations are not always suitable for modeling systems at coarse-grained resolution. For example, in the calculation of traditional torsion angle potentials, numerical instability is often encountered in the case of very flexible molecules. To improve the stability and accuracy of coarse-grained molecular dynamics simulations, we propose two approaches. The first makes use of improved forms for the angle potentials: the restricted bending (ReB) potential prevents torsion angles from visiting unstable or unphysical configurations and the combined bending-torsion (CBT) potential smoothly flattens the interactions when such configurations are sampled. In the second approach, dummy-assisted dihedral (DAD), the torsion potential is applied differently: instead of acting directly on the beads, it acts on virtual beads, bound to the real ones. For simple geometrical reasons, the unstable region is excluded from the accessible conformational space. The benefits of the new approaches are demonstrated in simulations of polyethylene glycol (PEG), polystyrene (PS), and polypeptide molecules described by the MARTINI coarse-grained force field. The new potentials are implemented in an in-house version of the Gromacs package, publicly available.

Original languageEnglish
Pages (from-to)3282-3292
Number of pages11
JournalJournal of Chemical Theory and Computation
Volume9
Issue number8
DOIs
Publication statusPublished - Aug-2013

Keywords

  • FORCE-FIELD
  • PULL-OUT
  • HYDROCARBONS
  • BUTANE
  • MODEL

Cite this