Abstract:The study of phenomena such as protein folding and conformational changes in molecules is a central theme in chemical physics. Molecular dynamics (MD) simulation is the primary tool for the study of transition processes in biomolecules, but it is hampered by a huge timescale gap between the processes of interest and atomic vibrations that dictate the time step size. Therefore, it is imperative to combine MD simulations with other techniques in order to quantify the transition processes taking place on large timescales. In this work, the diffusion map with Mahalanobis kernel, a meshless approach for approximating the Backward Kolmogorov Operator (BKO) in collective variables, is upgraded to incorporate standard enhanced sampling techniques, such as metadynamics. The resulting algorithm, which we call the target measure Mahalanobis diffusion map ( tm-mmap ), is suitable for a moderate number of collective variables in which one can approximate the diffusion tensor and free energy. Imposing appropriate boundary conditions allows use of the approximated BKO to solve for the committor function and utilization of transition path theory to find the reactive current delineating the transition channels and the transition rate. The proposed algorithm, tm-mmap , is tested on the two-dimensional Moro–Cardin two-well system with position-dependent diffusion coefficient and on alanine dipeptide in two collective variables where the committor, the reactive current, and the transition rate are compared to those computed by the finite element method (FEM). Finally, tm-mmap is applied to alanine dipeptide in four collective variables where the use of finite elements is infeasible.

Committor Guided Estimates of Molecular Transition Rates

Committor guided estimates of molecular transition rates

Computing the Committor with the Committor: an Anatomy of the Transition State Ensemble

Computing the committor with the committor to study the transition state ensemble

Computing Committor Functions for the Study of Rare Events Using Deep Learning

From high-dimensional committors to reactive insights

Iterative variational learning of committor-consistent transition pathways using artificial neural networks

On committor functions in milestoning

Everything everywhere all at once, a probability-based enhanced sampling approach to rare events

Deep Learning Method for Computing Committor Functions with Adaptive Sampling

Data-driven Computation of Molecular Reaction Coordinates

Splitting probabilities as optimal controllers of rare reactive events

Biomolecular transitions: efficient computation of pathways, free energies, and rates

Computing committors via Mahalanobis diffusion maps with enhanced sampling data

Transition pathways in complex systems: Reaction coordinates, isocommittor surfaces, and transition tubes

Coupling rare event algorithms with data-based learned committor functions using the analogue Markov chain

Coherent control from quantum committment probabilities