Search papers, labs, and topics across Lattice.
This study reformulates the variational principle for estimating committors and reaction rates by relaxing traditional boundary conditions, allowing for a broader class of trial functions. By employing a normalization of a single boundary observable, the fidelity, the authors successfully derive committors in high-dimensional spaces using only one-dimensional profiles and minimal input data. The approach yields efficient closed-form solutions for committors and rates in complex protein folding scenarios, demonstrating its applicability to systems like AIB9 and villin HP-35, as well as chignolin's folding dynamics.
Relaxing boundary conditions allows for accurate committor and reaction rate estimations in high-dimensional spaces using minimal input data.
The committor is the optimal reaction coordinate for a rare transition: it pinpoints the transition state and fixes the rate, and it governs events from protein folding to crystal nucleation. It minimises a Dirichlet energy, whose value at the minimum is the reactive flux, over functions that vanish on the reactant state and equal one on the product state. In high dimensions such a trial space is very hard to build. Here we rewrite the variational principle so that the boundary conditions are replaced by a normalisation of one boundary observable, the fidelity. Any trial function is then admissible, including functions that cannot satisfy the boundary values at all. On this basis we estimate the high-dimensional committor and rates from one-dimensional profiles along projected coordinates, taking as input only pre-existing equilibrium or reweighted configurations, a diffusion constant estimate, and the two state definitions. The optimum has a closed form that is cheap to evaluate. We obtain committors for AIB9 and villin HP-35 in full torsion space, and folding and unfolding rates for chignolin from umbrella sampling alone.
