Un puerto de montaña conecta dos valles. En uno, la proteína está plegada; en el otro, desplegada. Para pasar de uno a otro, hay que cruzar la silla, el punto más alto del camino. En la ciencia de las proteínas, este punto se llama estado de transición, y la coordenada que indica la distancia hasta él se denomina commitador. Pero una proteína real tiene miles de dimensiones (los ángulos de rotación de los átomos), y encontrar esa silla es casi imposible.
Los autores del trabajo reformularon el antiguo principio variacional de modo que ya no es necesario imponer condiciones rígidas en los límites. En su lugar, introdujeron una única normalización, la llamada "fidelidad". Ahora se puede usar cualquier función de prueba, incluso una que no sea cero al inicio ni uno al final. Es como trabajar con un boceto aproximado del puerto en lugar de un mapa exacto: lo importante es que la escala sea correcta.
El método se probó en tres proteínas: AIB9, vilina HP-35 y higrolina. Para las dos primeras se obtuvieron commitadores en el espacio completo de ángulos, y para la higrolina, las velocidades de plegamiento y despliegue, utilizando solo datos de muestreo de paraguas. Todo se calcula rápidamente, mediante una fórmula cerrada. Esto abre el camino para estudiar eventos raros en sistemas grandes, desde el plegamiento de proteínas hasta la formación de cristales.
Ya no es necesario asaltar el puerto de montaña a ciegas. Basta con mirar un simple mapa de altitudes.
tal como lo escribieron los autores, en inglés 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.
La vilina HP-35 se pliega en solo unos microsegundos: es uno de los procesos de plegamiento más rápidos conocidos.