Massively Parallel Implementation of Steered Molecular Dynamics in Tinker-HP: Comparisons of Polarizable and Non-Polarizable Simulations of Realistic Systems

Steered molecular dynamic (SMD) is a powerful technique able to accelerate rare event sampling in Molecular Dynamics (MD) simulations by applying an external force to a set of chosen atoms. Despite generating nonequilibrium simulations, SMD remains capable of reconstructing equilibrium properties such as the Potential of Mean Force (PMF). Of course, one would like to use all types of force fields (FF) ranging from classical ones to more advanced polarizable models using point induced dipoles and distributed multipoles such as AMOEBA. To enable such studies, the SMD methodology has been implemented in the framework of the massively parallel Tinker-HP software allowing for both long polarizable and non-polarizable MD simulations of large proteins. To validate this new implementation, we first compared the Tinker-HP SMD results to the literature. Tests have been performed on three different benchmark systems: the M–A deca-alanine (112 atoms), the ubiquitin (9737 atoms), and the CD2CD58 complex (97594 atoms)....