Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA’s diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

RNA is involved in a broad range of biological processes that extend far beyond translation. Many of RNA's recently discovered functions rely on folding to a specific conformation or transitioning between conformations. The RNA structure contains rigid, short basepaired regions connected by more flexible linkers. Studies of model constructs such as small helix-junction-helix (HJH) motifs are useful in understanding how these elements work together to determine RNA conformation. Here, we reveal the full ensemble of solution structures assumed by a model RNA HJH. We apply small-angle x-ray scattering and an ensemble optimization method to selectively refine models generated by all-atom molecular dynamics simulations. The expectation of a broad distribution of helix orientations, at and above physiological ionic strength, is not met. Instead, this analysis shows that the HJH structures are dominated by two distinct conformations at moderate to high ionic strength. Atomic structures, selected from the molecular dynamics simulations, reveal strong base-base interactions in the junction that critically constrain the conformational space available to the HJH molecule and lead to a surprising re-extension at high salt. These results are corroborated by comparison with previous single-molecule fluorescence resonance energy transfer experiments on the same constructs.

No abstract

The biological functions of RNA range from gene regulation through catalysis and depend critically on its structure and flexibility. Conformational variations of flexible, non-base-paired components, including RNA hinges, bulges, or single-stranded tails, are well documented. Recent work has also identified variations in the structure of ubiquitous, base-paired duplexes found in almost all functional RNAs. Duplexes anchor the structures of folded RNAs, and their surface features are recognized by partner molecules. To date, no consistent picture has been obtained that describes the range of conformations assumed by RNA duplexes. Here, we apply wide angle, solution X-ray scattering (WAXS) to quantify these variations, by sampling length scales characteristic of helical geometries under different solution conditions. To identify the radius, helical rise, twist, and length of dsRNA helices, we exploit molecular dynamics generated structures, explicit solvent models, and ensemble optimization methods. Our results quantify the substantial and salt-dependent deviations of double-stranded (ds) RNA duplexes from the assumed canonical A-form conformation. Recent experiments underscore the need to properly describe the structures of RNA duplexes when interpreting the salt dependence of RNA conformations.

Macromolecular structures can be determined from solution X-ray scattering. Small-angle X-ray scattering (SAXS) provides global structural information on length scales of 10s to 100s of Ångstroms, and many algorithms are available to convert SAXS data into low-resolution structural envelopes. Extension of measurements to wider scattering angles (WAXS or wide-angle X-ray scattering) can sharpen the resolution to below 10 Å, filling in structural details that can be critical for biological function. These WAXS profiles are especially challenging to interpret because of the significant contribution of solvent in addition to solute on these smaller length scales. Based on training with molecular dynamics generated models, the application of extreme gradient boosting (XGBoost) is discussed, which is a supervised machine learning (ML) approach to interpret features in solution scattering profiles. These ML methods are applied to predict key structural parameters of double-stranded ribonucleic acid (dsRNA) duplexes. Duplex conformations vary with salt and sequence and directly impact the foldability of functional RNA molecules. The strong structural periodicities in these duplexes yield scattering profiles with rich sets of features at intermediate-to-wide scattering angles. In the ML models, these profiles are treated as 1D images or features. These ML models identify specific scattering angles, or regions of scattering angles, which correspond with and successfully predict distinct structural parameters. Thus, this work demonstrates that ML strategies can integrate theoretical molecular models with experimental solution scattering data, providing a new framework for extracting highly relevant structural information from solution experiments on biological macromolecules.