Understanding Reaction Mechanisms from Start to Finish

Abstract

Understanding mechanisms of rare but important events in complex molecular systems, such as protein folding or ligand (un)binding, requires accurately mapping transition paths from an initial to a final state. The committor is the ideal reaction coordinate for this purpose, but calculating it for high-dimensional, nonlinear systems has long been considered intractable. Here, we introduce an iterative path sampling strategy for computing the committor function for systems with high free energy barriers. We start with an initial guess to define isocommittor interfaces for transition interface sampling. The resulting path ensemble is then reweighted and used to train a neural network, yielding a more accurate committor model. This process is repeated until convergence, effectively solving the long-standing circular problem in enhanced sampling where a good reaction coordinate is needed to generate efficient sampling, and vice-versa. The final, converged committor model can be interrogated to extract mechanistic insights. We demonstrate the power of our method on a benchmark 2D potential and a more complex host-guest (un)binding process in explicit solvent.

References in corpus (9)

• Statistically optimal analysis of samples from multiple equilibrium states
• A Novel Path Sampling Method for the Calculation of Rate Constants
• Rational Construction of Stochastic Numerical Methods for Molecular Sampling
• Reaction rate calculation by parallel path swapping
• Autonomous artificial intelligence discovers mechanisms of molecular self-organization in virtual experiments
• Computing the Committor with the Committor: an Anatomy of the Transition State Ensemble
• Learning reaction coordinates via cross-entropy minimization: Application to alanine dipeptide
• Molecular free energy profiles from force spectroscopy experiments by inversion of observed committors
• Exchanging replicas with unequal cost, infinitely and permanently