Several enhanced sampling techniques rely on the definition of collective variables to effectively explore free energy landscapes. The existing variables that describe the progression along a reactive pathway offer an elegant solution but face a number of limitations. In this paper, we address these challenges by introducing a new path-like collective variable called the "deep-locally non-linear-embedding," which is inspired by principles of the locally linear embedding technique and is trained on a reactive trajectory. The variable mimics the ideal reaction coordinate by automatically generating a non-linear combination of features through a differentiable generalized autoencoder that combines a neural network with a continuous k-nearest neighbor selection. Among the key advantages of this method is its capability to automatically choose the metric for searching neighbors and to learn the path from state A to state B without the need to handpick landmarks a priori. We demonstrate the effectiveness of DeepLNE by showing that the progression along the path variable closely approximates the ideal reaction coordinate in toy models, such as the Müller-Brown potential and alanine dipeptide. Then, we use it in the molecular dynamics simulations of an RNA tetraloop, where we highlight its capability to accelerate transitions and estimate the free energy of folding.
Download full-text PDF |
Source |
---|---|
http://dx.doi.org/10.1063/5.0202156 | DOI Listing |
J Chem Phys
September 2024
School of Pharmaceutical Sciences, University of Geneva, Rue Michel Servet 1, 1206 Genève, Switzerland.
Path-like collective variables (CVs) can be very effective for accurately modeling complex biomolecular processes in molecular dynamics simulations. Recently, we have introduced DeepLNE (deep-locally non-linear-embedding), a machine learning-based path-like CV that provides a progression variable s along the path as a non-linear combination of several descriptors. We have demonstrated the effectiveness of DeepLNE by showing that for simple models such as the Müller-Brown potential and alanine dipeptide, the progression along the path variable closely approximates the ideal reaction coordinate.
View Article and Find Full Text PDFJ Chem Phys
May 2024
School of Pharmaceutical Sciences, University of Geneva, Rue Michel Servet 1, 1206 Genève, Switzerland.
Several enhanced sampling techniques rely on the definition of collective variables to effectively explore free energy landscapes. The existing variables that describe the progression along a reactive pathway offer an elegant solution but face a number of limitations. In this paper, we address these challenges by introducing a new path-like collective variable called the "deep-locally non-linear-embedding," which is inspired by principles of the locally linear embedding technique and is trained on a reactive trajectory.
View Article and Find Full Text PDFJ Phys Chem Lett
February 2024
Research Center for Analytical Sciences, Frontiers Science Center for New Organic Matter, College of Chemistry, Nankai University, Tianjin Key Laboratory of Biosensing and Molecular Recognition, State Key Laboratory of Medicinal Chemical Biology, Tianjin 300071, China.
Enhanced-sampling algorithms relying on collective variables (CVs) are extensively employed to study complex (bio)chemical processes that are not amenable to brute-force molecular simulations. The selection of appropriate CVs characterizing the slow movement modes is of paramount importance for reliable and efficient enhanced-sampling simulations. In this Perspective, we first review the application and limitations of CVs obtained from chemical and geometrical intuition.
View Article and Find Full Text PDFJ Chem Theory Comput
January 2019
Department of Chemistry , University College London, London WC1E 6BT , United Kingdom.
Path Collective Variables (PCVs) are a set of path-like variables that have been successfully used to investigate complex chemical and biological processes and compute their associated free energy surfaces and kinetics. Their current implementation relies on general, but at times inefficient, metrics (such as RMSD or DRMSD) to evaluate the distance between the instantaneous conformational state during the simulation and the reference coordinates defining the path. In this work, we present a new algorithm to construct optimal PCVs metrics as linear combinations of different CVs weighted through a spectral gap optimization procedure.
View Article and Find Full Text PDFEnter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!