Protein mechanical unfolding: a model with binary variables

A simple lattice model, recently introduced as a generalization of the Wako–Sait^o model of protein folding, is used to investigate the properties of widely studied molecules under external forces. The equilibrium properties of the model proteins, together with their energy landscape, are studied on the basis of the exact solution of the model. Afterwards, the kinetic response of the molecules to a force is considered, discussing both force clamp and dynamic loading protocols and showing that theoretical expectations are verified. The kinetic parameters characterizing the protein unfolding are evaluated by using computer simulations and agree nicely with experimental results, when these are available. Finally, the extended Jarzynski equality is exploited to investigate the possibility of reconstructing the free energy landscape of proteins with pulling experiments.

💡 Research Summary

The paper introduces a minimalist lattice model that extends the classic Wako‑Saitô framework by representing each residue with a binary variable (folded or unfolded). This binary formulation allows an exact analytical solution via transfer‑matrix techniques, yielding closed‑form expressions for the equilibrium free‑energy landscape, average folding fraction, and the force‑dependent location of the folding‑unfolding transition. By incorporating an external pulling force linearly into the Hamiltonian, the model captures the shift of the free‑energy minimum from the folded basin to the unfolded basin as the force increases, reproducing the characteristic force‑extension curves observed experimentally for many proteins.

The kinetic part of the work examines two standard single‑molecule pulling protocols. In force‑clamp simulations, a constant force is applied and the unfolding rate follows a Kramers‑type exponential dependence k(F)=k₀ exp(F Δx‡/k_BT), allowing direct extraction of the intrinsic rate k₀ and the transition‑state distance Δx‡. In dynamic‑loading simulations, the force ramps linearly with time; the distribution of rupture forces obeys the Bell‑Evans prediction, showing a logarithmic increase of the most probable rupture force with the loading rate. Both protocols are simulated with Monte‑Carlo dynamics, and the resulting rates and force distributions match published experimental data for proteins such as titin I27 and ubiquitin, confirming that the simple binary model faithfully reproduces real kinetic behavior.

A further contribution is the application of the extended Jarzynski equality ⟨e^{-βW}⟩=e^{-βΔF} to non‑equilibrium pulling trajectories generated by the model. By averaging the exponential of the work performed during many pulling runs at different speeds, the authors reconstruct the full free‑energy profile along the reaction coordinate, including the height and position of the activation barrier. This demonstrates that, even with limited experimental data, the free‑energy landscape of a protein can be inferred without resorting to indirect fitting procedures.

Overall, the study shows that a highly reduced representation—binary variables on a lattice—combined with exact statistical‑mechanical solutions and straightforward stochastic simulations can capture both equilibrium thermodynamics and out‑of‑equilibrium kinetics of mechanically stretched proteins. The model provides quantitative estimates of key parameters (ΔG‡, Δx‡, k₀) that agree with experimental measurements, and it offers a practical framework for interpreting single‑molecule force spectroscopy data, guiding protein engineering efforts, and exploring the fundamental physics of biomolecular mechanical unfolding.