mdRNAcat: A Molecular Dynamics RNA Catalog of REST2 Conformational Ensembles
Authors: Anton E. Dorn, Emile de Bruyn, Fabrice von der Lehr, Philipp Knechtges, Stefan Kesselheim, Alexander Schug
Organizations: Jülich Supercomputing Centre, Forschungszentrum Jülich (FZJ) · Faculty of Physics, Georg August University of Göttingen · Institute of Software Technology, Germany Aerospace Center (DLR) · Department of Physics, University of Cologne · Scientific Computing Center, Karlsruhe Institute of Technology (KIT)
Learning molecular models of RNA requires simulation data that capture conformational variability and provide consistent physical training signals. Here, we introduce mdRNAcat, a curated dataset of all-atom molecular dynamics (MD) simulations of 48 RNA systems selected to promote sequence diversity. Each system is simulated using replica-exchange with solute scaling (REST2), with 32 Hamiltonian states and production lengths of 200-300 ns per state. The dataset provides atomic coordinates, instantaneous forces, and forces reevaluated under the unscaled reference Hamiltonian, together with structural annotations and replica-exchange metadata. Coordinates and forces are available at 1 ps intervals for the reference state and 100 ps intervals across all states in a standardized HDF5 format. We characterize conformational sampling, exchange behavior, base pairing, and sugar-pucker distributions to document the properties of the simulated ensembles. mdRNAcat provides a reusable resource for developing and evaluating data-driven RNA models, including coarse-grained force fields.
Figures & tables
Figure 1: ( a ) Distribution of the number of residues per RNA system in the dataset. ( b ) Distribution of the percentage of paired residues in the dataset; most systems display similar behavior, apart from one RNA system in an extended state. Colors indicate different RNA domains.
Total RNA systems
48
Systems with 200ns simulation time
45
Systems with 300ns simulation time
3
Replicas
32
State trajectories
1,536
Total sampled time
316.8
\SIUnitSymbolMicros
Total sampled time at 300K
9.9
\SIUnitSymbolMicros
Table 1: Statistical overview of the dataset.
Field
Size
Type
Unit
Description
<pdbid>/
chain
N
string
Chain ID
element
N
string
Chemical element
resid
N
integer
Residue number
resname
N
string
Residue name
z
N
integer
Atomic number
Table 2: mdRNAcat is made available as HDF5 files with an hierarchical organization as shown in this table. For each system in the data set, there are 2 files both labeled with the respective PDB ID. The first contains only the lowest trajectory with 1ps time steps. The second file contains all replicas and additional derived metrics indicated with a blue background in the table with 100ps time steps. Key: <T> , REST2 effective temperature (K); N , number of atoms; R , number of residues; F , trajectory length in frames ( = frameStridePs ); M , number of replicas; K , energy terms; S , energy samples; X , exchange records; P , log segments. 2 2 2 A comment on the units used in mdRNAcat can be found in the Appendix A.1
Figure 2: Sequence and structural-motif coverage of mdRNAcat against filtered X-ray PDB reference sets. ( a ) Coverage versus locally normalized nearest-representative distance. ( b , c ) UMAP visualizations McInnes et al. (2018a) ; McInnes et al. (2018b) ; stars denote mdRNAcat representatives and colors show smoothed distance ratios on separate scales, with structural values clipped for display. ( d ) Mean structural-descriptor differences in units of the reference-set standard deviation. Coverage is calculated before UMAP projection.
Figure 3: Free energy landscape (FEL) of the lowest temperature ( 300K ) trajectory of the 6qn3 system, the axes are determined with time-lagged independent component analysis (TICA) ( Pérez-Hernández et al. (2013) ), marked with the experimentally resolved structure (exp) and the two lowest minima (1 and 2); experimentally resolved structure is displayed to the left and structures from the two minima to the right; in yellow FEL of a plain single replica 300K simulation of 6qn3 (not part of mdRNAcat) mapped to the same TICA coordinates.
Figure 4: ( a ) Transition probability between two thermodynamic states and energy overlap between two thermodynamic states against the 31 transitions of the states for the 3r4f system, with dotted lines indicating the average over all transitions. ( b ) Density of the thermodynamic states against the energy of the same system. ( c ) Density of the replicas against the energy of the same system.
Avg. transition probability
0.183
±0.065
Avg. energy overlap
0.269
±0.070
Table 3: Replica exchange statistics averaged over all 48 systems with standard deviation between systems.
Figure 5: Free energy landscape (FEL) of the 2qus, 6cu1 and 7elq systems, the axes are determined with time-lagged independent component analysis (TICA) ( Pérez-Hernández et al. (2013) ) at 300K , 354K , 424K and 500K . Note that the axes are capped at the extension of the 424K -plot, so the 500K -plot has values outside of that range not displayed. The conformation space gets explored more with higher temperature.
Figure 6: ( a ) Wheel histogram over the sugar-pucker angle of all residues of all systems; north ( 0\SIUnitSymbolDegree ) indicates C3’-endo, so A-form helical pucker; south ( 180\SIUnitSymbolDegree ) indicates C2’-endo, so loops and bulges. ( b ) Distribution of the residue pair likelihood against the sequence separation between the two paired residue for all systems in the data set; thick is the set mean. ( c ) Same as ( b ) but the x-axis normalized to 1.
Figure 7: Probability of a paired residue in residue and temperature dimensions for the 1xjr, 6lax, 7elq and 9obm systems. On the left of each plot is a bow plot indicating the paired residues in the experimentally resolved structure.
Appendix figures & tables8 assets
Supplementary material from the paper’s appendix.
Appendix
PDB ID
Residues
Simulation length
Status
1duh
44
200ns
Part of mdRNAcat
1ehz
76
Modified residues not parameterized
1hr2
158
Too large
1kh6
53
Modified residues not parameterized
1kxk
70
300ns
Part of mdRNAcat
1nbs
155
Missing residues in the pdb
Appendix
Table 4: RNA systems considered or used in mdRNAcat.
Parameter
Value
integrator
steep
nsteps
100000
emtol
10
emstep
0.01
nstlog
1000
nstenergy
100
Appendix
Table 5: Parameters used for energy minimization.
Parameter
Value
define
-DPOSRES
integrator
md
nsteps
50000
dt
0.002
nstxout
500
nstvout
500
Appendix
Table 6: Parameters used in the NVT equilibration.
Parameter
Value
define
-DPOSRES
integrator
md
nsteps
50000
dt
0.002
nstxout
500
nstvout
500
Appendix
Table 7: Parameters used for NPT equilibration.
Parameter
Value
integrator
md
nsteps
100000000
dt
0.002
nstxout
250
nstvout
250
nstfout
250
Appendix
Table 8: Parameters used for production simulation runs.
Replica
Ti [ K ]
λ
λ
1
300.000000
1.000000
1.000000
2
304.980000
0.983700
0.991817
3
310.050000
0.967600
0.983667
4
315.200000
0.951800
0.975602
5
320.440000
0.936200
0.967574
6
325.760000
0.920900
0.959635
Appendix
Table 9: REST2 effective-temperature ladder and λ scaling factors.
Figure 8: Density of the thermodynamic states against the energy of the first 24 systems in mdRNAcat.
Figure 9: Density of the thermodynamic states against the energy of the second 24 systems in mdRNAcat.
Ribonucleic acid (RNA) functions through conformational changes that are not fully captured by static structures. However, large-scale standardized RNA dynamics data remain limited, and existing approaches typically treat trajectory generation and dynamics understanding as separate objectives. Here, we introduce RNADynBench, a standardized RNA molecular dynamics (MD) benchmark with 2585 quality-controlled 100-ns all-atom trajectories and leakage-controlled splits. Building on RNADynBench, we develop RNADynNet, a unified model for RNA dynamics learning that uses a shared backbone for both trajectory generation and dynamics fingerprint extraction from a single conformer. It combines coordinate denoising, single-frame-to-trajectory alignment, and physical grounding to connect all-atom trajectory generation with dynamics representation learning. Physical grounding improves both generated dynamics and the physical information recoverable from these fingerprints. Across both test sets, including the high-flexibility challenge set, the generated trajectories achieve RMSF correlations of 0.875 and 0.766, while single-conformer predictions show comparable agreement with MD-derived dynamics. RNADynBench and RNADynNet together establish a benchmark and unified modeling framework for generating and understanding RNA dynamics.
Yiming Huang, Lennart Bastian, Hanqun Cao +2
Imperial College London · The Chinese University of Hong Kong · Technical University of Munich
Accurate RNA structure modeling remains difficult because RNA backbones are highly flexible, non-canonical interactions are prevalent, and experimentally determined 3D structures are comparatively scarce. We introduce RiboSphere, a framework that learns discrete geometric representations of RNA by combining vector quantization with flow matching. Our design is motivated by the modular organization of RNA architecture: complex folds are composed from recurring structural motifs. RiboSphere uses a geometric transformer encoder trained using mean-centered coordinates and random rotation augmentation to produce geometry-aware features, which are discretized with finite scalar quantization (FSQ) into a finite vocabulary of latent codes. Conditioned on these discrete codes, a flow-matching decoder reconstructs atomic coordinates, enabling high-fidelity structure generation. We find that the learned code indices are enriched for specific RNA motifs, suggesting that the model captures motif-level compositional structure rather than acting as a purely compressive bottleneck. Across benchmarks, RiboSphere achieves strong performance in structure reconstruction (RMSD 1.25,Å, TM-score 0.84), and its pretrained discrete representations transfer effectively to inverse folding and RNA--ligand binding prediction, with robust generalization in data-scarce regimes. Code is available at https://github.com/Zhangz312/RiboSphere.
Zhou Zhang, Hanqun Cao, Cheng Tan +3
State Key Laboratory for Novel Software Technology at Nanjing University, School of Computer Science, Nanjing University · The Chinese University of Hong Kong · Stanford University
RNA molecules often function through multiple experimentally observed conformational states. A central challenge is to capture ensemble diversity without sacrificing structural fidelity: samples may collapse toward averaged representative structures or drift into unsupported conformational regions. We introduce REFLEX (RNA Ensemble generation calibrated by FLEXibility), a sequence-conditioned generator for RNA backbone-frame ensembles based on a heteroscedastic stochastic bridge on SE(3). During training, experimental B-factors calibrate residue-wise bridge widths, exposing each observed conformer through local perturbation neighborhoods with residue-specific scales while keeping constrained regions sharper. A learned flexibility condition guides the bridge-field network to recover target conformer-specific frames from these calibrated noisy states. At inference, REFLEX uses only the input sequence, requiring no MSA-subsampling heuristics. On a multi-conformer RNA benchmark, REFLEX achieves the best precision-recall trade-off for ensemble coverage across evaluation thresholds, while recovery-stress ablations indicate gains in target-frame recovery under stronger perturbations.