Massively parallel implementation of Steered Molecular Dynamics in Tinker-HP: comparisons of polarizable and non-polarizable simulations of realistic systems

International audienceSteered Molecular Dynamic (SMD) is a powerful technique able to accelerate rare events sampling in Molecular Dynamics (MD) simulations by applying an external force to a set of chosen atoms. Despite generating non-equilibrium 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 dis- tributed multipoles such as AMOEBA. To enable such studies, the SMD methodology has been implemented in the framework of the massively parallel Tinker–HP softwareallowing 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 benchmarks: the M–A deca–alanine (112 atoms), the Ubiquitin (9737 atoms) and the CD2CD58 com- plex (97594 atoms). Non–polarizable (AMBER99, AMBER99SB, CHARMM22CMAP and OPLSAAL) and polarizable (AMOEBABIO09, AMOEBAPRO13 and AMOEBABIO18) force fields have been used. For each one of them, PMFs have been re-constructed and compared in terms of free energy barrier results and hydrogen bonding behavior fluctuations over time. Using a SMD velocity of 0.01 Å/ps applied to a set of 20 trajectories, we show that polarizable and non–polarizable force fields do notalways agree. As it could be anticipated, strong differences are noticed between polarizable and non-polarizable models when considered in vacuum, whereas results are morecomparable when a water environment is added. However, for the largest system, i.e. the CD2CD58 complex, strong differences related to the modeling of a salt bridge are noticed exhibiting some potential issues with classical FFs. Overall, such simulations highlight the importance of the inclusion of polarization effects as PMFs free energy barriers computed with AMOEBA always decrease compared to non–polarizable force fields