Ojha Anupam Anand, Thakur Saumya, Ahn Surl-Hee, Amaro Rommie E
Department of Chemistry, University of California San Diego, La Jolla, California92093, United States.
Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, Maharashtra400076, India.
J Chem Theory Comput. 2023 Feb 28;19(4):1342-1359. doi: 10.1021/acs.jctc.2c00282. Epub 2023 Jan 31.
Recent advances in computational power and algorithms have enabled molecular dynamics (MD) simulations to reach greater time scales. However, for observing conformational transitions associated with biomolecular processes, MD simulations still have limitations. Several enhanced sampling techniques seek to address this challenge, including the weighted ensemble (WE) method, which samples transitions between metastable states using many weighted trajectories to estimate kinetic rate constants. However, initial sampling of the potential energy surface has a significant impact on the performance of WE, i.e., convergence and efficiency. We therefore introduce deep-learned kinetic modeling approaches that extract statistically relevant information from short MD trajectories to provide a well-sampled initial state distribution for WE simulations. This hybrid approach overcomes any statistical bias to the system, as it runs short unbiased MD trajectories and identifies meaningful metastable states of the system. It is shown to provide a more refined free energy landscape closer to the steady state that could efficiently sample kinetic properties such as rate constants.
计算能力和算法的最新进展使分子动力学(MD)模拟能够达到更长的时间尺度。然而,对于观察与生物分子过程相关的构象转变,MD模拟仍然存在局限性。几种增强采样技术试图应对这一挑战,包括加权系综(WE)方法,该方法使用许多加权轨迹对亚稳态之间的转变进行采样,以估计动力学速率常数。然而,势能面的初始采样对WE的性能有重大影响,即收敛性和效率。因此,我们引入了深度学习动力学建模方法,该方法从短MD轨迹中提取统计相关信息,为WE模拟提供采样良好的初始状态分布。这种混合方法克服了对系统的任何统计偏差,因为它运行短的无偏MD轨迹并识别系统中有意义的亚稳态。结果表明,它能提供更接近稳态的更精细的自由能景观,从而可以有效地采样诸如速率常数等动力学性质。