Melting of Branched RNA Molecules

Stability of the branching structure of an RNA molecule is an important condition for its function. In this letter we show that the melting thermodynamics of RNA molecules is very sensitive to their branching geometry for the case of a molecule whose groundstate has the branching geometry of a Cayley Tree and whose pairing interactions are described by the Go model. Whereas RNA molecules with a linear geometry melt via a conventional continuous phase transition with classical exponents, molecules with a Cayley Tree geometry are found to have a free energy that seems smooth, at least within our precision. Yet, we show analytically that this free energy in fact has a mathematical singularity at the stability limit of the ordered structure. The correlation length appears to diverge on the high-temperature side of this singularity.

💡 Research Summary

The paper investigates how the branching geometry of an RNA molecule influences its thermal denaturation (melting) behavior, focusing on a system whose native (ground‑state) secondary structure is a Cayley tree and whose base‑pair interactions are modeled by the Go model. The Go model simplifies the energetics by assigning favorable interaction energy only to base pairs that appear in the native structure, thereby eliminating non‑native contacts and allowing a clean theoretical analysis of geometry‑driven effects.

First, the authors recapitulate the well‑known melting transition of linear (one‑dimensional) RNA. Using the same Go Hamiltonian, they show that as temperature increases, native base pairs break progressively, leading to a continuous second‑order phase transition. The free‑energy curvature (heat capacity) exhibits a sharp peak at a critical temperature (T_c), and the correlation length (\xi) diverges as (\xi\propto|T-T_c|^{-1}), consistent with mean‑field critical exponents for a one‑dimensional system with short‑range interactions.

The core of the study then turns to a Cayley tree—a regular, infinite‑branching lattice where each node has a fixed coordination number (k). This topology dramatically increases the number of branches with depth, creating a “high‑dimensional” effective space. By exploiting the recursive self‑similarity of the tree, the authors derive an exact recursion relation for the partition function of a subtree of depth (n). The recursion distinguishes configurations where the root bond is intact from those where it is broken, leading to a nonlinear map (P_{n+1}=f(P_n)). Iterating this map numerically yields the free energy (F_N(T)) and heat capacity (C_N(T)) for trees of finite depth (N).

Numerical results reveal two striking features. (1) The free‑energy curve appears smooth over the entire temperature range; no obvious cusp or discontinuity is visible, suggesting the absence of a conventional thermodynamic singularity. (2) The heat capacity does not develop a pronounced peak even for large (N); instead it shows a broad, modest variation. At first glance this would imply that a Cayley‑tree RNA melts without a phase transition, in stark contrast to the linear case.

To resolve this apparent paradox, the authors perform an analytical study of the recursion’s fixed points in the complex temperature plane. They demonstrate that a branch‑point singularity exists at a temperature (T_c^\ast) that marks the stability limit of the ordered (fully paired) state. While the singularity is hidden on the real axis—making the free energy look analytic—it manifests as a non‑analyticity in higher derivatives when the function is analytically continued. Importantly, the correlation length diverges on the high‑temperature side of this point, following (\xi\propto|T-T_c^\ast|^{-\nu}) with an exponent (\nu) that depends on the branching factor (k). Larger (k) yields smaller (\nu), indicating a sharper divergence. This divergence reflects the fact that, even though individual bonds break, the extensive branching sustains long‑range correlations throughout the tree.

The paper therefore concludes that Cayley‑tree RNA exhibits a “hidden” critical behavior: the thermodynamic observables (free energy, heat capacity) appear smooth, yet the underlying statistical mechanics contains a genuine singularity and a diverging correlation length. This contrasts with the conventional continuous transition of linear RNA and highlights the profound impact of topology on RNA stability.

From an experimental perspective, the authors suggest that high‑resolution single‑molecule techniques such as temperature‑dependent FRET, optical tweezers, or small‑angle X‑ray scattering could detect subtle asymmetries in the heat capacity or the rapid growth of correlation length near (T_c^\ast). Synthetic RNA nanostructures engineered to adopt a Cayley‑tree secondary structure would provide a direct test of the theoretical predictions.

In broader terms, the work implies that biological RNAs with complex branching (e.g., ribozymes, viral genomes) may possess intrinsic thermal robustness or sensitivity that is not captured by simple linear models. The interplay between branching geometry and energetic frustration could be a design principle exploited by nature to fine‑tune functional stability under physiological temperature fluctuations. Future extensions could explore random or heterogeneous branching, incorporate sequence‑dependent stacking energies, or apply the framework to other polymeric systems (e.g., dendrimers) where similar hidden criticalities may arise.