Open AccessRemoving Thermostat Distortions of Protein Dynamics in Constant-Temperature Molecular Dynamics Simulations
Author(s) -
Alan Hicks,
Matthew MacAinsh,
HuanXiang Zhou
Publication year - 2021
Publication title -
journal of chemical theory and computation
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.001
H-Index - 185
eISSN - 1549-9626
pISSN - 1549-9618
DOI - 10.1021/acs.jctc.1c00448
Subject(s) - thermostat , molecular dynamics , langevin dynamics , time constant , protein dynamics , constant (computer programming) , relaxation (psychology) , statistical physics , fick's laws of diffusion , rotational diffusion , langevin equation , diffusion , chemistry , physics , thermodynamics , computational chemistry , computer science , molecule , quantum mechanics , psychology , social psychology , electrical engineering , programming language , engineering
Molecular dynamics simulations are widely used to determine equilibrium and dynamic properties of proteins. Nearly all simulations, currently, are carried out at constant temperature, with a Langevin thermostat among the most widely used. Thermostats distort protein dynamics, but whether or how such distortions can be corrected has long been an open question. Here, we show that constant-temperature simulations with a Langevin thermostat dilate protein dynamics and present a correction scheme to remove the dynamic distortions. Specifically, ns-scale time constants for overall rotation are dilated significantly but sub-ns time constants for internal motions are dilated modestly, while all motional amplitudes are unaffected. The correction scheme involves contraction of the time constants, with the contraction factor a linear function of the time constant to be corrected. The corrected dynamics of eight proteins are validated by NMR data for rotational diffusion and for backbone amide and side-chain methyl relaxation. The present work demonstrates that even for complex systems like proteins with dynamics spanning multiple timescales, one can predict how thermostats distort protein dynamics and remove such distortions. The correction scheme will have wide applications, facilitating force-field parameterization and propelling simulations to be on par with NMR and other experimental techniques in determining dynamic properties of proteins.