Martini straight: Boosting performance using a shorter cutoff and GPUs

Abstract In molecular dynamics simulations, sufficient sampling is of key importance and a continuous challenge in the field. The coarse grain Martini force field has been widely used to enhance sampling. In its original implementation, this force field applied a shifted Lennard-Jones potential for the non-bonded van der Waals interactions, to avoid problems related to a relatively short cutoff. Here we investigate the use of a straight cutoff Lennard-Jones potential with potential modifiers. Together with a Verlet neighbor search algorithm, the modified potential allows the use of GPUs to accelerate the computations in Gromacs. We find that this alternative potential has little influence on most of the properties studied, including partitioning free energies, bulk liquid properties and bilayer properties. At the same time, energy conservation is kept within reasonable bounds. We conclude that the newly proposed straight cutoff approach is a viable alternative to the standard shifted potentials used in Martini, offering significant speedup even in the absence of GPUs.

[1]  Wilfred F van Gunsteren,et al.  Biomolecular modeling: Goals, problems, perspectives. , 2006, Angewandte Chemie.

[2]  Vijay S. Pande,et al.  Screen Savers of the World Unite! , 2000, Science.

[3]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

[4]  Helgi I. Ingólfsson,et al.  Computational Lipidomics with insane: A Versatile Tool for Generating Custom Membranes for Molecular Simulations. , 2015, Journal of chemical theory and computation.

[5]  L. Verlet Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[6]  A. Mark,et al.  Coarse grained model for semiquantitative lipid simulations , 2004 .

[7]  R. C. Weast CRC Handbook of Chemistry and Physics , 1973 .

[8]  D. Tieleman,et al.  Perspective on the Martini model. , 2013, Chemical Society reviews.

[9]  Giovanni Bussi,et al.  Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration , 2013, Entropy.

[10]  Wilfred F van Gunsteren,et al.  On using a too large integration time step in molecular dynamics simulations of coarse-grained molecular models. , 2009, Physical chemistry chemical physics : PCCP.

[11]  Helgi I Ingólfsson,et al.  The power of coarse graining in biomolecular simulations , 2013, Wiley interdisciplinary reviews. Computational molecular science.

[12]  Durba Sengupta,et al.  Polarizable Water Model for the Coarse-Grained MARTINI Force Field , 2010, PLoS Comput. Biol..

[13]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[14]  Siewert J Marrink,et al.  Comment on "On using a too large integration time step in molecular dynamics simulations of coarse-grained molecular models" by M. Winger, D. Trzesniak, R. Baron and W. F. van Gunsteren, Phys. Chem. Chem. Phys., 2009, 11, 1934. , 2010, Physical chemistry chemical physics : PCCP.

[15]  D. Tieleman,et al.  The MARTINI force field: coarse grained model for biomolecular simulations. , 2007, The journal of physical chemistry. B.

[16]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[17]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

[18]  D. Rana,et al.  Dynamic Interfacial Tension Behavior of Water/Oil Systems Containing In situ-Formed Surfactants , 2002 .

[19]  John L. Klepeis,et al.  Millisecond-scale molecular dynamics simulations on Anton , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[20]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .