Elastic energy of proteins and the stages of protein folding

We propose a universal elastic energy for proteins, which depends only on the radius of gyration $R_{g}$ and the residue number $N$. It is constructed using physical arguments based on the hydrophobic effect and hydrogen bonding. Adjustable parameters are fitted to data from the computer simulation of the folding of a set of proteins using the CSAW (conditioned self-avoiding walk) model. The elastic energy gives rise to scaling relations of the form $R_{g}\sim N^{\nu}$ in different regions. It shows three folding stages characterized by the progression with exponents $\nu = 3/5, 3/7, 2/5$, which we identify as the unfolded stage, pre-globule, and molten globule, respectively. The pre-globule goes over to the molten globule via a break in behavior akin to a first-order phase transition, which is initiated by a sudden acceleration of hydrogen bonding.

💡 Research Summary

The authors present a minimalist, universal elastic‑energy framework for protein folding that depends only on two macroscopic variables: the radius of gyration (R g) and the number of residues (N). Their construction is grounded in two dominant physical forces that drive protein compaction: the hydrophobic effect, which tends to reduce the solvent‑exposed surface area, and hydrogen bonding, which stabilizes internal secondary‑structure elements. By treating the hydrophobic contribution as a surface‑area term proportional to N 2⁄3 R g⁻¹ and the hydrogen‑bond contribution as a volume‑like term proportional to N R g⁻², they arrive at an elastic‑energy expression

E(R g,N) = A N^{2/3} R g⁻¹ + B N R g⁻²,

where A and B are phenomenological constants. Minimizing this energy with respect to R g yields scaling relations of the form R g ∝ N^{ν}. The authors fit A and B to data generated by the Conditioned Self‑Avoiding Walk (CSAW) model, a coarse‑grained Monte‑Carlo simulation that incorporates excluded‑volume, hydrophobic, and hydrogen‑bonding interactions. The fit reveals three distinct exponent values: ν ≈ 3⁄5, 3⁄7, and 2⁄5, each corresponding to a recognizable stage in the folding pathway.

The ν = 3⁄5 regime matches the classical Flory exponent for a polymer in a good solvent and is identified as the unfolded, random‑coil stage. Here, hydrophobic forces dominate but hydrogen bonds are scarce, so the chain remains relatively expanded. The intermediate ν = 3⁄7 regime is interpreted as a “pre‑globule” phase. In this region the chain has begun to collapse under hydrophobic pressure, yet the hydrogen‑bond network is still incomplete, leading to a compact but not fully ordered structure. Finally, ν = 2⁄5 characterizes the molten‑globule state, where a rapid increase in hydrogen‑bond formation produces a more tightly packed, yet still dynamic, ensemble. The transition from pre‑globule to molten‑globule appears as a sharp kink in the R g‑N plot, reminiscent of a first‑order phase transition: the first derivative of the free‑energy landscape changes discontinuously as hydrogen bonding accelerates, driving the system into a new energetic minimum.

The paper’s strength lies in its ability to capture the essential thermodynamic driving forces of protein folding with a simple analytic form, bypassing the need for atomistic simulations while still reproducing experimentally observed scaling behavior. It also provides a clear physical interpretation of the three folding stages, linking macroscopic exponents to microscopic interactions. However, several limitations are acknowledged. The parameters A and B are calibrated exclusively against CSAW simulations, which, although sophisticated, are still coarse‑grained and may not capture all nuances of real proteins (e.g., electrostatic interactions, salt effects, temperature dependence). Moreover, the model assumes uniform hydrophobic and hydrogen‑bond contributions across all residues, ignoring sequence‑specific propensities that can significantly modulate folding pathways. Finally, the proposed first‑order‑like transition is inferred from static scaling data; dynamic studies would be required to confirm the presence of a true thermodynamic discontinuity.

Future work suggested by the authors includes (i) validation against experimental measurements of R g (e.g., small‑angle X‑ray scattering) across a broad range of protein sizes, (ii) extension of the energy functional to incorporate additional terms for electrostatic screening, solvent viscosity, and temperature, and (iii) kinetic simulations that track the time evolution of hydrogen‑bond formation to elucidate the microscopic mechanism underlying the abrupt pre‑globule‑to‑molten‑globule conversion. By addressing these points, the universal elastic‑energy model could become a powerful tool for predicting folding behavior, guiding protein design, and interpreting experimental data across the diverse landscape of globular proteins.