Equilibrium properties and force-driven unfolding pathways of RNA molecules

The mechanical unfolding of a simple RNA hairpin and of a 236–bases portion of the Tetrahymena thermophila ribozyme is studied by means of an Ising–like model. Phase diagrams and free energy landscapes are computed exactly and suggest a simple two–state behaviour for the hairpin and the presence of intermediate states for the ribozyme. Nonequilibrium simulations give the possible unfolding pathways for the ribozyme, and the dominant pathway corresponds to the experimentally observed one.

💡 Research Summary

The paper presents a theoretical study of mechanically induced unfolding of RNA using an Ising‑like statistical model. Two systems are examined: a simple hairpin and a 236‑nucleotide fragment of the Tetrahymena thermophila ribozyme. In the model each base‑pair is represented by a binary spin variable (+1 for folded, –1 for unfolded). Nearest‑neighbor interactions encode base‑pair stacking, while an external pulling force contributes a linear term proportional to the end‑to‑end extension. Because the Hamiltonian is a sum of one‑dimensional Ising terms, the partition function can be evaluated exactly, allowing the authors to compute free‑energy landscapes, entropy, and specific heat as functions of temperature and force without resorting to Monte‑Carlo sampling for equilibrium properties.

For the hairpin, the exact solution yields a clear two‑state picture. The free‑energy profile as a function of extension shows two minima separated by a modest barrier. At low force the folded minimum is global; beyond a critical force (~15 pN) the unfolded minimum becomes lower, reproducing the classic two‑state unfolding transition observed in single‑molecule force spectroscopy. The phase diagram in the temperature–force plane displays a single line of first‑order transitions, confirming the simplicity of the system.

The ribozyme fragment, by contrast, exhibits a richer landscape. Exact calculations reveal several shallow minima corresponding to partially unfolded intermediates. In the force range 12–18 pN an intermediate basin becomes thermodynamically significant, indicating that specific helices (notably the P2‑P3‑P4 region) unfold before the rest of the structure. The phase diagram therefore contains multiple transition lines, reflecting a multistate behavior.

To explore kinetic pathways, the authors perform nonequilibrium Monte‑Carlo simulations in which the pulling force is ramped linearly in time. These simulations generate unfolding trajectories that can be classified into distinct pathways. The dominant pathway—observed in the majority of runs—starts with the disruption of the P2‑P3‑P4 helices, followed sequentially by the unfolding of the P5‑P6‑P7 region. A secondary pathway, in which P5‑P6‑P7 unfolds first, appears with much lower probability. The distribution of pathways matches experimental observations from optical‑trap pulling experiments on the same ribozyme fragment, providing strong validation of the model.

The study highlights several advantages of the Ising‑like approach: (1) analytical tractability for equilibrium thermodynamics, (2) minimal parameter set (stacking energy, force coupling, temperature), and (3) computational efficiency that enables systematic scanning of large RNA segments. Limitations are also acknowledged: the model treats base‑pair interactions isotropically, neglects explicit ion screening and solvent effects, and discretizes conformational space, which may miss subtle structural fluctuations captured by all‑atom molecular dynamics.

In summary, the work demonstrates that a simple Ising‑type representation can accurately reproduce both the equilibrium phase behavior and the dominant nonequilibrium unfolding routes of RNA molecules ranging from a simple hairpin to a complex ribozyme fragment. The findings provide a useful theoretical framework for interpreting single‑molecule force spectroscopy data and for guiding the design of RNA‑based nanodevices where controlled mechanical response is essential.