4.7 Article

A Companion Guide to the String Method with Swarms of Trajectories: Characterization, Performance, and Pitfalls

Journal

JOURNAL OF CHEMICAL THEORY AND COMPUTATION
Volume 18, Issue 3, Pages 1406-1422

Publisher

AMER CHEMICAL SOC
DOI: 10.1021/acs.jctc.1c01049

Keywords

-

Funding

  1. National Natural Science Foundation of China [22073050]
  2. US National Institutes of Health [R15-GM139140, P41-GM104601]
  3. Na-tional Science Foundation [MCB-1517221, CHE-1945465]
  4. France and Chicago Collaborating in The Sciences (FACCTS) program
  5. Agence Nationale de la Recherche
  6. 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

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available