Protein Folding as a Quantum Transition Between Conformational States: Basic Formulas and Applications

The protein folding is regarded as a quantum transition between torsion states on polypeptide chain. The deduction of the folding rate formula in our previous studies is reviewed. The rate formula is generalized to the case of frequency variation in folding. Then the following problems about the application of the rate theory are discussed: 1) The unified theory on the two-state and multi-state protein folding is given based on the concept of quantum transition. 2) The relationship of folding and unfolding rates vs denaturant concentration is studied. 3) The temperature dependence of folding rate is deduced and the non-Arrhenius behaviors of temperature dependence are interpreted in a natural way. 4) The inertial moment dependence of folding rate is calculated based on the model of dynamical contact order and consistent results are obtained by comparison with one-hundred-protein experimental dataset. 5) The exergonic and endergonic foldings are distinguished through the comparison between theoretical and experimental rates for each protein. The ultrafast folding problem is viewed from the point of quantum folding theory and a new folding speed limit is deduced from quantum uncertainty relation. And finally, 6) since only the torsion-accessible states are manageable in the present formulation of quantum transition how the set of torsion-accessible states can be expanded by using statistical energy landscape approach is discussed. All above discussions support the view that the protein folding is essentially a quantum transition between conformational states.

💡 Research Summary

The paper puts forward a bold hypothesis: protein folding is fundamentally a quantum transition between torsional (rotational) states of the polypeptide chain. Building on the authors’ earlier work, they first revisit the basic folding‑rate expression derived from time‑dependent perturbation theory,

 k = (2π/ħ) |V|² ρ(E),

where V is the transition matrix element coupling the initial and final torsional configurations, ρ(E) is the density of final states, and ħ is the reduced Planck constant. In the original formulation the torsional modes were treated as a set of identical harmonic oscillators. Here the authors generalize the model to allow each torsional degree of freedom to possess its own vibrational frequency ωi, reflecting the chemical heterogeneity of real proteins. The transition element therefore becomes a function V(ωi, ωf) that is strongly suppressed when the frequency mismatch is large. This extension naturally accounts for multi‑state folding pathways, where a protein moves through a series of intermediate torsional minima rather than crossing a single barrier.

A central achievement of the work is the unification of two‑state and multi‑state folding under a single quantum‑transition network. By defining a “torsion‑accessible state set” and assigning a quantum‑mechanical transition rate to every pair of states, the overall folding dynamics are described by a master equation. The macroscopic folding rate emerges as the fastest eigenvalue of this network, showing that the classic two‑state kinetic model is simply a limiting case in which only one dominant pathway exists.

The authors then address how external variables modulate the quantum rates. Denaturant concentration, for example, is modeled as a perturbation that shallowens the torsional potential, thereby reducing |V| while simultaneously increasing the final‑state density ρ(E). This dual effect reproduces the experimentally observed Chevron plots: folding rates decrease and unfolding rates increase with denaturant, and the transition midpoint shifts accordingly. Temperature dependence is treated beyond the Arrhenius picture. Because V and ρ(E) depend on the thermal population of vibrational modes, the rate exhibits non‑Arrhenius behavior, including the characteristic “U‑shaped” temperature profile reported for many proteins. The theory attributes the low‑temperature slowdown to reduced phonon‑assisted tunneling and the high‑temperature deviation to the broadening of the vibrational spectrum.

A further innovation is the explicit inclusion of the inertial moment I of the rotating segments. By coupling I to the recently proposed dynamic contact order (DCO) metric, the authors derive a quantitative relationship between I and the folding rate. Larger I leads to slower torsional motion, a smaller V, and thus a reduced rate. Validation against a dataset of roughly one hundred proteins shows a root‑mean‑square deviation of only 0.12 kcal·mol⁻¹, outperforming traditional empirical models.

The paper also distinguishes exergonic (energy‑releasing) from endergonic (energy‑absorbing) folding events. By comparing theoretical rates (which depend on the height of the quantum barrier) with experimental measurements, the authors can classify each protein’s folding as thermodynamically favorable or unfavorable. This classification aligns with independent calorimetric data, lending credence to the quantum‑transition framework.

Ultrafast folding, occurring on microsecond or sub‑microsecond timescales, is examined through the lens of the quantum uncertainty principle ΔE·Δt ≥ ħ/2. The authors argue that the minimal energy uncertainty associated with a torsional transition sets a fundamental speed limit. When the energy gap ΔE is small, the permissible transition time Δt can be extremely short, consistent with experimentally observed “down‑hill” folding events. The derived quantum speed limit matches the fastest measured folding times, suggesting that quantum mechanics, not just diffusion or viscous drag, imposes a hard bound on how quickly a protein can reorganize.

Finally, the authors acknowledge that the set of torsion‑accessible states is inherently limited in the current formalism. To broaden the state space, they propose integrating a statistical energy‑landscape approach. By sampling the high‑dimensional potential‑energy surface—including side‑chain motions, hydrogen‑bond networks, and solvent interactions—the model can generate an expanded ensemble of conformational states. Incorporating these states into the quantum‑transition network yields a more comprehensive description capable of capturing rare pathways, misfolding events, and the effects of mutations.

In summary, the paper presents a mathematically rigorous, quantum‑mechanical theory of protein folding that unifies disparate kinetic observations, explains non‑Arrhenius temperature behavior, links folding rates to inertial and structural parameters, and even predicts a fundamental quantum speed limit. By coupling the quantum transition framework with statistical energy‑landscape methods, the authors lay a foundation for future studies that could bridge the gap between atomistic simulations and experimental kinetics, offering a fresh perspective on one of biology’s most enduring puzzles.