TAPS: A traveling-salesman based automated path searching method for functional conformational changes of biological macromolecules.

Abstract

Locating the minimum free energy paths (MFEPs) between two conformational states is among the most important tasks of biomolecular simulations. For example, knowledge of the MFEP is critical for focusing the effort of unbiased simulations that are used for the construction of Markov state models to the biologically relevant regions of the system. Typically, existing path searching methods perform local sampling around the path nodes in a pre-selected collective variable (CV) space to allow a gradual downhill evolution of the path toward the MFEP. Despite the wide application of such a strategy, the gradual path evolution and the non-trivial a priori choice of CVs are also limiting its overall efficiency and automation. Here we demonstrate that non-local perpendicular sampling can be pursued to accelerate the search, provided that all nodes are reordered thereafter via a traveling-salesman scheme. Moreover, path-CVs can be computed on-the-fly and used as a coordinate system, minimizing the necessary prior knowledge about the system. Our traveling-salesman based automated path searching method achieves a 5-8 times speedup over the string method with swarms-of-trajectories for two peptide systems in vacuum and solution, making it a promising method for obtaining initial pathways when investigating functional conformational changes between a pair of structures.