Capturing the essence of folding and functions of biomolecules using Coarse-Grained Models

The distances over which biological molecules and their complexes can function range from a few nanometres, in the case of folded structures, to millimetres, for example during chromosome organization. Describing phenomena that cover such diverse length, and also time scales, requires models that capture the underlying physics for the particular length scale of interest. Theoretical ideas, in particular, concepts from polymer physics, have guided the development of coarse-grained models to study folding of DNA, RNA, and proteins. More recently, such models and their variants have been applied to the functions of biological nanomachines. Simulations using coarse-grained models are now poised to address a wide range of problems in biology.

💡 Research Summary

The review by Thirumalai and colleagues provides a comprehensive overview of how coarse‑grained (CG) modeling, grounded in polymer physics, can be used to capture the essential physics of biomolecular folding and function across an extraordinary range of length scales—from a few nanometers in folded proteins, RNA and DNA to millimeter‑scale chromosome organization. The authors begin by emphasizing that while atomistic simulations offer detailed insight, they are limited by the enormous time and size scales inherent to biological processes. Instead, CG models, which retain only the degrees of freedom relevant at a given scale, enable the study of phenomena that would otherwise be computationally intractable.

The paper is organized around the concept that the appropriate level of coarse‑graining is dictated by the length scale of interest. At sub‑nanometer resolution (≈0.5 nm), explicit chemical interactions such as hydrogen bonding and base stacking must be represented to describe nucleic‑acid base‑pairing and protein side‑chain contacts. At the persistence‑length scale of double‑stranded DNA (≈50 nm), the polymer can be treated as a semi‑flexible filament, and the Worm‑Like Chain (WLC) model accurately reproduces force‑extension curves, loop‑formation kinetics, and supercoiling behavior. For longer DNA (microns to millimeters), the chain behaves as a flexible self‑avoiding polymer, allowing the use of Flory‑type scaling arguments and lattice models to explore chromosome folding, segregation, and territorial organization.

Specific case studies illustrate the power of this multiscale approach. Loop formation (cyclization) of dsDNA is shown to depend non‑linearly on the ratio L/lp: short chains (L/lp ≈ 1) experience high bending energy, while very long chains suffer entropic penalties, leading to a minimum cyclization time around L/lp ≈ 2–3. This insight helps explain why bacterial promoters often contain DNA segments of ~100 bp that are optimally sized for regulatory looping.

The authors also discuss DNA stretching experiments on λ‑phage DNA, demonstrating that the WLC model maps onto a quantum rotor problem and yields persistence length values (≈53 nm) that match experimental data across a wide force range. This validates the use of simple elastic models for interpreting single‑molecule force spectroscopy.

A particularly striking example is the entropic mechanism of bacterial chromosome segregation. By simulating self‑avoiding polymers confined within cylindrical volumes that mimic the nucleoid, the authors show that newly replicated “daughter” chains spontaneously migrate to the periphery, driven solely by an increase in configurational entropy. This suggests that proteins may play a modulatory rather than a primary driving role in bacterial chromosome partitioning.

In the realm of chromatin organization, the review contrasts equilibrium globules (which predict a contact probability decay I(s) ∝ s⁻¹·⁵) with fractal globules (I(s) ∝ s⁻¹). Monte Carlo simulations of 4000‑bead polymers reproduce the experimentally observed s⁻¹ scaling in Hi‑C data, supporting the fractal globule model as a realistic description of human genome folding that preserves local genomic proximity while avoiding topological entanglements.

RNA folding is examined through the lens of polyelectrolyte theory. Multivalent cations (Z > 1) are far more efficient at neutralizing the negatively charged backbone than monovalent ions, shifting the folding midpoint Cₘ by orders of magnitude and causing the radius of gyration to scale as Rg ∝ 1/Z². Coarse‑grained representations that collapse 5–6 nucleotides into a single bead successfully reproduce SAXS‑derived intermediate structures of the Tetrahymena ribozyme, revealing a rapid early compaction followed by slower, fluid‑like rearrangements.

Finally, the review addresses the complex kinetics of hairpin formation in short RNA or ssDNA. Single‑molecule force spectroscopy combined with temperature quenches uncovers a bistable free‑energy landscape with two basins (folded and unfolded) that become equally populated at a critical (Tm, fm). The folding pathway differs markedly between force‑quench and temperature‑quench protocols, with an initial slow nucleation (loop formation) followed by rapid “zipping” of the stem. This multistep behavior cannot be captured by simple two‑state models but emerges naturally from CG simulations that incorporate both thermal fluctuations and mechanical forces.

Overall, the paper convincingly argues that coarse‑grained models, when guided by polymer physics and calibrated against experimental observables, provide a versatile and powerful framework for dissecting the mechanisms of biomolecular folding and function. The authors advocate for continued development of CG approaches that integrate crowding, post‑translational modifications, and protein–RNA/DNA interactions, thereby extending the reach of computational biophysics into ever more complex cellular contexts.