Research Center for Analytical Sciences, College of Chemistry, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Nankai University, Tianjin 300071, China.
Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States.
J Chem Theory Comput. 2022 Mar 8;18(3):1406-1422. doi: 10.1021/acs.jctc.1c01049. Epub 2022 Feb 9.
The string method with swarms of trajectories (SMwST) is an algorithm that identifies a physically meaningful transition pathway─a one-dimensional curve, embedded within a high-dimensional space of selected collective variables. The SMwST algorithm leans on a series of short, unbiased molecular dynamics simulations spawned at different locations of the discretized path, from whence an average dynamic drift is determined to evolve the string toward an optimal pathway. However conceptually simple in both its theoretical formulation and practical implementation, the SMwST algorithm is computationally intensive and requires a careful choice of parameters for optimal cost-effectiveness in applications to challenging problems in chemistry and biology. In this contribution, the SMwST algorithm is presented in a self-contained manner, discussing with a critical eye its theoretical underpinnings, applicability, inherent limitations, and use in the context of path-following free-energy calculations and their possible extension to kinetics modeling. Through multiple simulations of a prototypical polypeptide, combining the search of the transition pathway and the computation of the potential of mean force along it, several practical aspects of the methodology are examined with the objective of optimizing the computational effort, yet without sacrificing accuracy. In light of the results reported here, we propose some general guidelines aimed at improving the efficiency and reliability of the computed pathways and free-energy profiles underlying the conformational transitions at hand.
群体轨迹字符串方法 (SMwST) 是一种算法,用于识别物理上有意义的转变途径——一条嵌入在所选集体变量高维空间中的一维曲线。SMwST 算法依赖于一系列短的、无偏的分子动力学模拟,这些模拟在离散路径的不同位置产生,从中确定平均动态漂移,以将字符串向最佳路径演变。尽管在理论公式和实际实现方面都非常简单,但 SMwST 算法计算量很大,并且需要仔细选择参数,以在化学和生物学中的挑战性问题的应用中实现最佳的成本效益。在本贡献中,以自包含的方式呈现了 SMwST 算法,用批判的眼光讨论了其理论基础、适用性、内在局限性以及在路径跟随自由能计算中的应用,以及它们可能扩展到动力学建模。通过对典型多肽的多次模拟,结合转变途径的搜索和沿其计算平均力势,检查了该方法的几个实际方面,目的是优化计算工作量,同时又不牺牲准确性。根据这里报告的结果,我们提出了一些通用准则,旨在提高所计算途径的效率和可靠性,并改善当前构象转变的自由能分布。