Горный перевал соединяет две долины. В одной белок свёрнут, в другой — развёрнут. Чтобы перейти из одной в другую, нужно пересечь седловину — самую высокую точку на пути. В науке о белках эту точку называют переходным состоянием, а координату, которая показывает расстояние до неё, — коммиттором. Но у реального белка тысячи измерений — углы поворота атомов, — и найти эту седловину почти невозможно.
Авторы работы переписали старый вариационный принцип так, что больше не нужно задавать жёсткие условия на границах. Вместо этого они ввели одну нормировку — так называемую «верность» (fidelity). Теперь можно использовать любую пробную функцию, даже ту, что не обращается в ноль на старте и в единицу на финише. Это похоже на работу с примерным эскизом перевала вместо точной карты — главное, чтобы масштаб был верным.
Метод проверили на трёх белках: AIB9, виллине HP-35 и хигнолине. Для первых двух получили коммитторы в полном пространстве углов, а для хигнолина — скорости сворачивания и разворачивания, используя только данные из зонтичной выборки. Всё считается быстро, по замкнутой формуле. Это открывает путь к изучению редких событий в больших системах — от сворачивания белков до образования кристаллов.
Горный перевал больше не нужно штурмовать вслепую. Достаточно взглянуть на простую карту высот.
как её написали авторы, на английском arXiv:2608.02536v1
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.
Белок виллин HP-35 сворачивается всего за несколько микросекунд — это один из самых быстрых известных процессов сворачивания.