bioRxivBiophysicsPreprint
Enhanced-Sampling Molecular Dynamics Recovers Rare Functional RNA Conformations Across Diverse Structural Contexts
Abstract
Accurate determination of RNA conformational ensembles is essential for understanding RNA function and advancing RNA-targeted drug discovery, yet lowly-populated alternative states remain difficult to resolve with atomistic detail. A central constraint is that experimental refinement can only select conformations already present in the starting library, making library generation the limiting step. Using the HIV-1 trans-activation response element (TAR) as a model system, we benchmarked conventional MD (cMD) against the enhanced sampling methods Gaussian-accelerated MD (GaMD), replica-exchange Gaussian-accelerated MD (Rex-GaMD), replica-exchange with solute tempering (REST2), and temperature replica-exchange MD (T-REMD), as well as the structure-prediction based methods FARFAR2 and AlphaFold 3. Each library was refined against experimental residual dipolar couplings (RDC) and validated independently using ensemble-averaged QM/MM chemical shifts. We showed that T-REMD produced the most accurate ensemble by both measures, and its advantage tracked with broader, more continuous coverage of the interhelical conformational landscape. Broad temperature-range T-REMD also sampled conformations resembling excited state 1 (ES1) and the U23-A27-U38 base-triple, without requiring these states to be specified during library generation. More accurate ensembles further improved coverage of experimentally observed ligand-bound TAR conformations and enhanced ensemble-based virtual screening, linking structural accuracy to functional utility. The same workflow applied to the preQ1 class I riboswitch and the UUCG tetraloop improved agreement with experimental data in both cases. Together, these results establish replica-exchange enhanced sampling, particularly T-REMD, as an effective strategy for constructing experimentally validated RNA ensembles and accessing conformations corresponding to rare functional substates.
- assay: structural biology
- research focus: protein structure prediction
- approach: deep learning
- approach: statistical model
- approach: simulation
- Citations
- 0
Reproducibility ledger
Nobody has reported trying to reproduce this paper yet. Did you replicate an experiment or run the code? Your report saves the next lab a week.
Sign in to file a reproducibility report.
Related papers
- Integrative Ensemble Modeling reveals RNA conformations targetable by small molecules
Stefano Bosio, Vincent Schnapka, … Matteo Masetti·bioRxiv·Jul 29, 2026
- Toward Accurate RNA Folding Thermodynamics: Evaluation of Enhanced Sampling Methods for Force Field Benchmarking
Petra Kührová, Vojtěch Mlýnský, … Pavel Banáš·Journal of Chemical Theory and ComputationPeer-reviewed·May 14, 2026
- An RNA dynamic ensemble at atomic resolution
H. Shi, A. Rangadurai, … H. M. Al-Hashimi·bioRxiv·May 13, 2020
- Towards Convergence in Folding Simulations of RNA Tetraloops: Comparison of Enhanced Sampling Techniques and Effects of Force Field Corrections
V. Mlynsky, M. Janecek, … J. Sponer·bioRxiv·Nov 30, 2021
- DynaRNA: Dynamic RNA Conformation Ensemble Generation with Diffusion Model
Zhengxin Li, Junjie Zhu, … Haifeng Chen·bioRxiv·May 22, 2025
- Topology, Energy Landscapes and Folding Kinetics of RNA Hairpins
Anjali R. Verma, Pratyush Tiwary, Yunrui Qiu·bioRxiv·Sep 22, 2026