Journal
JOURNAL OF CHEMICAL THEORY AND COMPUTATION
Volume 18, Issue 3, Pages 1406-1422Publisher
AMER CHEMICAL SOC
DOI: 10.1021/acs.jctc.1c01049
Keywords
-
Funding
- National Natural Science Foundation of China [22073050]
- US National Institutes of Health [R15-GM139140, P41-GM104601]
- Na-tional Science Foundation [MCB-1517221, CHE-1945465]
- France and Chicago Collaborating in The Sciences (FACCTS) program
- Agence Nationale de la Recherche
- Beckman Institute
Ask authors/readers for more resources
The SMwST algorithm is a method for identifying transition pathways, which is conceptually simple but computationally intensive and requires careful parameter selection. This article presents the theoretical basis, applicability, limitations, and discusses the use of the algorithm in free-energy calculations and kinetics modeling. Through multiple simulations of a prototypical polypeptide, the methodology is examined to optimize computational effort while maintaining accuracy.
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 costeffectiveness 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.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available