An Efficient GaMD Multi-Level Enhanced Sampling Strategy: Application to Polarizable Force Fields Simulations of Large Biological Systems



We introduce a novel multi-level enhanced sampling strategy grounded on Gaussian accelerated Molecular Dynamics (GaMD). First, we propose a GaMD multi-GPUs-accelerated implementation within the Tinker-HP molecular dynamics package. We introduce the new "dual-water" mode and its use with the flexible AMOEBA polarizable force field.By adding harmonic boosts to the water stretching and bonding terms, it accelerates the solvent-solute interactions while enabling speedups thanks to the use of fast multiple--timestep integrators. To further reduce time-to-solution, we couple GaMD to Umbrella Sampling (US). The GaMD—US/dual--water approach is tested on the 1D Potential of Mean Force (PMF) of the solvated CD2--CD58 system (168000 atoms) allowing the AMOEBA PMF to converge within 1 kcal/mol of the experimental value. Finally, Adaptive Sampling (AS) is added enabling AS-GaMD capabilities but also the introduction of the new Adaptive Sampling--US--GaMD (ASUS-GaMD) scheme. The highly parallel ASUS--GaMD setup decreases time to convergence by respectively 10 and 20 times compared to GaMD-US and US. Overall, beside the acceleration of PMF computations, Tinker-HP now allows for the simultaneous use of Adaptive Sampling and GaMD-"dual water" enhanced sampling approaches increasing the applicability of polarizable force fields to large scale simulations of biological systems.

Version notes

Updated abstract


Supplementary material

Supplementary Materials
See main text for a detailed description of SI.