Efficient free energy calculations by combining two complementary tempering sampling methods

• Published in: The Journal of chemical physics 2017-01, Vol.146 (2), p.024103-024103
• Main Authors: Xie, Liangxu, Shen, Lin, Chen, Zhe-Ning, Yang, Mingjun
• Format: Article
• Language: English
• Subjects: Barriers; Computer simulation; Dipeptides - chemistry; Energy Transfer; Free energy; Isomerism; Mathematical Concepts; Molecular dynamics; Molecular Dynamics Simulation; Protein Conformation; Proteins - chemistry; Sampling methods; Temperature; Tempering; Thermodynamic properties; Thermodynamics
• ISSN: 0021-9606; 1089-7690
• DOI: 10.1063/1.4973607

Although energy barriers can be efficiently crossed in the reaction coordinate (RC) guided sampling, this type of method suffers from identification of the correct RCs or requirements of high dimensionality of the defined RCs for a given system. If only the approximate RCs with significant barriers are used in the simulations, hidden energy barriers with small to medium height would exist in other degrees of freedom (DOFs) relevant to the target process and consequently cause the problem of insufficient sampling. To address the sampling in this so-called hidden barrier situation, here we propose an effective approach to combine temperature accelerated molecular dynamics (TAMD), an efficient RC-guided sampling method, with the integrated tempering sampling (ITS), a generalized ensemble sampling method. In this combined ITS-TAMD method, the sampling along the major RCs with high energy barriers is guided by TAMD and the sampling of the rest of the DOFs with lower but not negligible barriers is enhanced by ITS. The performance of ITS-TAMD to three systems in the processes with hidden barriers has been examined. In comparison to the standalone TAMD or ITS approach, the present hybrid method shows three main improvements. (1) Sampling efficiency can be improved at least five times even if in the presence of hidden energy barriers. (2) The canonical distribution can be more accurately recovered, from which the thermodynamic properties along other collective variables can be computed correctly. (3) The robustness of the selection of major RCs suggests that the dimensionality of necessary RCs can be reduced. Our work shows more potential applications of the ITS-TAMD method as the efficient and powerful tool for the investigation of a broad range of interesting cases.