Protein Folding as a Quantum Transition Between Conformational States

The importance of torsion vibration in the transmission of life information is indicated. The localization of quantum torsion state is proved. Following these analyses a formalism on the quantum theory of conformation-electron system is proposed. The conformational-electronic transition is calculated by non-adiabatic operator method. The protein folding is viewed from conformational quantum transition and the folding rate is calculated. The time-scale of microsecond to millisecond for the fundamental folding event (nucleation, collapse, etc) is deduced. The dependence of transition rate W on N inertial moments is given. It indicates how W increases with the number N of torsion angles and decreases with the inertial moment I of atomic group in cooperative transition. The temperature dependence is also deduced which is different from chemical reaction in high-temperature region. It is demonstrated that the conformational dynamics gives deep insights into the folding mechanism and provides a useful tool for analyzing and explaining experimental facts on the rate of protein folding.

💡 Research Summary

The paper proposes a quantum‑mechanical framework for protein folding, treating the process as a transition between discrete torsional conformational states rather than a classical diffusion over a rugged energy landscape. The authors begin by identifying torsional vibrations of peptide bonds and side‑chain rotations as the primary carriers of biological information. Each torsional angle φ_i (i = 1…N) is quantized, leading to a wavefunction ψ(φ_1,…,φ_N) governed by the Hamiltonian Ĥ_torsion = ∑_i (L_i²/2I_i) + V(φ_1,…,φ_N), where L_i is the angular momentum operator, I_i the moment of inertia of the rotating atomic group, and V the torsional potential. By analyzing the potential barriers, they prove that the torsional wavefunction can become localized in specific minima, forming well‑defined quantum “normal modes”.

Next, the authors construct a combined conformational‑electronic system. The total Hamiltonian is Ĥ_total = Ĥ_torsion + Ĥ_electron + Ĥ_coupling, where Ĥ_coupling describes the interaction between the electronic subsystem and the torsional coordinates. Instead of invoking the usual Born‑Oppenheimer (adiabatic) approximation, they retain non‑adiabatic effects by introducing a transition operator V̂_non‑adiabatic = ∑_i (∂Ĥ/∂φ_i)(∂/∂L_i). This operator directly couples changes in torsional angles to changes in electronic angular momentum, allowing the calculation of transition matrix elements between initial (ψ_i) and final (ψ_f) conformational‑electronic states.

Using Fermi’s golden rule, the transition rate is derived as
W = (2π/ħ) |⟨ψ_f|V̂_non‑adiabatic|ψ_i⟩|² ρ(E),
where ρ(E) is the density of final states. The matrix element scales with the number of torsional degrees of freedom N and inversely with the square root of the effective moment of inertia I_eff of the rotating groups. Consequently, the authors obtain a compact scaling law:
W ∝ N · (k_BT/ħ) · exp(−ΔE/k_BT) / √I_eff.
This expression predicts that adding more rotatable bonds accelerates folding, while increasing the inertia of a cooperative segment slows it down.

Temperature dependence emerges from two competing contributions. At low temperatures, quantum tunneling through torsional barriers dominates, leading to a steep decline of W with decreasing T. At higher temperatures, thermal activation over the barriers adds a milder, Arrhenius‑like component, but the overall T‑dependence remains distinct from that of ordinary chemical reactions because the non‑adiabatic coupling modifies the effective activation free energy.

To estimate the intrinsic time scale, the authors use typical torsional frequencies ν ≈ 10¹²–10¹³ Hz and moments of inertia I ≈ 10⁻⁴⁶–10⁻⁴⁵ kg·m², yielding a fundamental transition time τ ≈ 1/ν·(I/ħ) on the order of microseconds to milliseconds. When realistic protein sizes (N ≈ 100–300) are inserted, the predicted τ matches experimentally observed nucleation and collapse phases of folding.

The paper further discusses how mutations that alter side‑chain mass or rigidity change I_eff, thereby modulating W. External factors such as pressure, solvent polarity, or temperature also reshape the torsional potential V(φ) and the inertia landscape, offering a quantitative route to predict folding rate changes under diverse conditions.

In conclusion, the authors argue that the quantum transition model complements classical energy‑landscape theories by providing a mechanistic link between torsional dynamics, electronic coupling, and folding kinetics. It explains non‑linear temperature effects, cooperative transitions, and the observed micro‑ to millisecond folding times without invoking ad hoc diffusion coefficients. The framework opens avenues for integrating quantum‑chemical calculations with molecular dynamics simulations and for designing experiments (e.g., ultrafast spectroscopy) that directly probe the non‑adiabatic transition matrix elements, ultimately advancing our ability to predict and control protein folding rates.