Visualizing RNA Structures by SAXS-Driven MD Simulations

Sequence-dependent conformational preferences of disordered single-stranded RNA

Disordered single-stranded RNA (ssRNA) molecules, like their well-folded counterparts, have crucial functions that depend on their structures. However, since native ssRNAs constitute a highly heterogeneous conformer population, their structural characterization poses challenges. One important question regards the role of sequence in influencing ssRNA structure. Here, we adopt an integrated approach that combines solution-based measurements, including small-angle X-ray scattering (SAXS) and Förster resonance energy transfer (FRET), with experimentally guided all-atom molecular dynamics (MD) simulations, to construct structural ensembles of a 30-nucleotide RNA homopolymer (rU30) and a 30-nucleotide RNA heteropolymer with an A-/C-rich sequence. We compare the size, shape, and flexibility of the two different ssRNAs. While the average properties align with polymer-physics descriptions of flexible polymers, we discern distinct, sequence-dependent conformations at the molecular level that demand a more detailed representation than provided by polymer models. These findings emphasize the role of sequence in shaping the overall properties of ssRNA.

Comparison of Three Computational Tools for the Prediction of RNA Tertiary Structures

Understanding the structures of noncoding RNAs (ncRNAs) is important for the development of RNA-based therapeutics. There are inherent challenges in employing current experimental techniques to determine the tertiary (3D) structures of RNAs with high complexity and flexibility in folding, which makes computational methods indispensable. In this study, we compared the utilities of three advanced computational tools, namely RNAComposer, Rosetta FARFAR2, and the latest AlphaFold 3, to predict the 3D structures of various forms of RNAs, including the small interfering RNA drug, nedosiran, and the novel bioengineered RNA (BioRNA) molecule showing therapeutic potential. Our results showed that, while RNAComposer offered a malachite green aptamer 3D structure closer to its crystal structure, the performances of RNAComposer and Rosetta FARFAR2 largely depend upon the secondary structures inputted, and Rosetta FARFAR2 predictions might not even recapitulate the typical, inverted "L" shape tRNA 3D structure. Overall, AlphaFold 3, integrating molecular dynamics principles into its deep learning framework, directly predicted RNA 3D structures from RNA primary sequence inputs, even accepting several common post-transcriptional modifications, which closely aligned with the experimentally determined structures. However, there were significant discrepancies among three computational tools in predicting the distal loop of human pre-microRNA and larger BioRNA (tRNA fused pre-miRNA) molecules whose 3D structures have not been characterized experimentally. While computational predictions show considerable promise, their notable strengths and limitations emphasize the needs for experimental validation of predictions besides characterization of more RNA 3D structures.

Insight into structural biophysics from solution X-ray scattering

The current challenges of structural biophysics include determining the structure of large self-assembled complexes, resolving the structure of ensembles of complex structures and their mass fraction, and unraveling the dynamic pathways and mechanisms leading to the formation of complex structures from their subunits. Modern synchrotron solution X-ray scattering data enable simultaneous high-spatial and high-temporal structural data required to address the current challenges of structural biophysics. These data are complementary to crystallography, NMR, and cryo-TEM data. However, the analysis of solution scattering data is challenging; hence many different analysis tools, listed in the SAS Portal (http://smallangle.org/), were developed. In this review, we start by briefly summarizing classical X-ray scattering analyses providing insight into fundamental structural and interaction parameters. We then describe recent developments, integrating simulations, theory, and advanced X-ray scattering modeling, providing unique insights into the structure, energetics, and dynamics of self-assembled complexes. The structural information is essential for understanding the underlying physical chemistry principles leading to self-assembled supramolecular architectures and computational structural refinement.

A generalizable scaffold-based approach for structure determination of RNAs by cryo-EM

Single-particle cryo-electron microscopy (cryo-EM) can reveal the structures of large and often dynamic molecules, but smaller biomolecules (≤50 kDa) remain challenging targets due to their intrinsic low signal to noise ratio. Methods to help resolve small proteins have been applied but development of similar approaches to aid in structural determination of small, structured RNA elements have lagged. Here, we present a scaffold-based approach that we used to recover maps of sub-25 kDa RNA domains to 4.5-5.0 Å. While lacking the detail of true high-resolution maps, these maps are suitable for model building and preliminary structure determination. We demonstrate this method helped faithfully recover the structure of several RNA elements of known structure, and that it promises to be generalized to other RNAs without disturbing their native fold. This approach may streamline the sample preparation process and reduce the optimization required for data collection. This first-generation scaffold approach provides a robust system to aid in RNA structure determination by cryo-EM and lays the groundwork for further scaffold optimization to achieve higher resolution.

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

Advances, Applications, and Perspectives in Small-Angle X-ray Scattering of RNA

RNAs exhibit a plethora of functions far beyond transmitting genetic information. Often, RNA functions are entailed in their structure, be it as a regulatory switch, protein binding site, or providing catalytic activity. Structural information is a prerequisite for a full understanding of RNA-regulatory mechanisms. Owing to the inherent dynamics, size, and instability of RNA, its structure determination remains challenging. Methods such as NMR spectroscopy, X-ray crystallography, and cryo-electron microscopy can provide high-resolution structures; however, their limitations make structure determination, even for small RNAs, cumbersome, if at all possible. Although at a low resolution, small-angle X-ray scattering (SAXS) has proven valuable in advancing structure determination of RNAs as a complementary method, which is also applicable to large-sized RNAs. Here, we review the technological and methodological advancements of RNA SAXS. We provide examples of the powerful inclusion of SAXS in structural biology and discuss possible future applications to large RNAs.

Magnesium ions mediate ligand binding and conformational transition of the SAM/SAH riboswitch

The SAM/SAH riboswitch binds S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH) with similar affinities. Mg2+ is generally known to stabilize RNA structures by neutralizing phosphates, but how it contributes to ligand binding and conformational transition is understudied. Here, extensive molecular dynamics simulations (totaling 120 μs) predicted over 10 inner-shell Mg2+ ions in the SAM/SAH riboswitch. Six of them line the two sides of a groove to widen it and thereby pre-organize the riboswitch for ligand entry. They also form outer-shell coordination with the ligands and stabilize an RNA-ligand hydrogen bond, which effectively diminishes the selectivity between SAM and SAH. One Mg2+ ion unique to the apo form maintains the Shine-Dalgarno sequence in an autonomous mode and thereby facilitates its release for ribosome binding. Mg2+ thus plays vital roles in SAM/SAH riboswitch function.

A Generalizable Scaffold-Based Approach for Structure Determination of RNAs by Cryo-EM

Single-particle cryo-electron microscopy (cryo-EM) can reveal the structures of large and often dynamic molecules, but smaller biomolecules remain challenging targets due to their intrinsic low signal to noise ratio. Methods to resolve small proteins have been applied but development of similar approaches for small structured RNA elements have lagged. Here, we present a scaffold-based approach that we used to recover maps of sub-25 kDa RNA domains to 4.5 - 5.0 Å. While lacking the detail of true high-resolution maps, these are suitable for model building and preliminary structure determination. We demonstrate this method faithfully recovers the structure of several RNA elements of known structure, and it promises to be generalized to other RNAs without disturbing their native fold. This approach may streamline the sample preparation process and reduce the optimization required for data collection. This first-generation scaffold approach provides a system for RNA structure determination by cryo-EM and lays the groundwork for further scaffold optimization to achieve higher resolution.

Recognizing the power of machine learning and other computational methods to accelerate progress in small molecule targeting of RNA

RNA structures regulate a wide range of processes in biology and disease, yet small molecule chemical probes or drugs that can modulate these functions are rare. Machine learning and other computational methods are well poised to fill gaps in knowledge and overcome the inherent challenges in RNA targeting, such as the dynamic nature of RNA and the difficulty of obtaining RNA high-resolution structures. Successful tools to date include principal component analysis, linear discriminate analysis, k-nearest neighbor, artificial neural networks, multiple linear regression, and many others. Employment of these tools has revealed critical factors for selective recognition in RNA:small molecule complexes, predictable differences in RNA- and protein-binding ligands, and quantitative structure activity relationships that allow the rational design of small molecules for a given RNA target. Herein we present our perspective on the value of using machine learning and other computation methods to advance RNA:small molecule targeting, including select examples and their validation as well as necessary and promising future directions that will be key to accelerate discoveries in this important field.

Magnesium ions mediate ligand binding and conformational transition of the SAM/SAH riboswitch

The SAM/SAH riboswitch binds S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH) with similar affinities. Mg 2+ is generally known to stabilize RNA structures by neutralizing phosphates, but how it contributes to ligand binding and conformational transition is understudied. Here, extensive molecular dynamics simulations (totaling 120 μs) identified over 10 inner-shell Mg 2+ ions in the SAM/SAH riboswitch. Six of them line the two sides of a groove to widen it and thereby pre-organize the riboswitch for ligand entry. They also form outer-shell coordination with the ligands and stabilize an RNA-ligand hydrogen bond, which effectively diminish the selectivity between SAM and SAH. One Mg 2+ ion unique to the apo form maintains the Shine-Dalgarno sequence in an autonomous mode and thereby facilitates its release for ribosome binding. Mg 2+ thus plays vital roles in SAM/SAH riboswitch function.

Structure and ensemble refinement against SAXS data: Combining MD simulations with Bayesian inference or with the maximum entropy principle

Small-angle X-ray scattering (SAXS) is a powerful method for tracking conformational transitions of proteins or soft-matter complexes in solution. However, the interpretation of the experimental data is challenged by the low spatial resolution and the low information content of the data, which lead to a high risk of overinterpreting the data. Here, we illustrate how SAXS data can be integrated into all-atom molecular dynamics (MD) simulation to derive atomic structures or heterogeneous ensembles that are compatible with the data. Besides providing atomistic insight, the MD simulation adds physicochemical information, as encoded in the MD force fields, which greatly reduces the risk of overinterpretation. We present an introduction into the theory of SAXS-driven MD simulations as implemented in GROMACS-SWAXS, a modified version of the GROMACS simulation software. We discuss SAXS-driven parallel-replica ensemble refinement with commitment to the maximum entropy principle as well as a Bayesian formulation of SAXS-driven structure refinement. Practical considerations for running and interpreting the simulations are presented. The methods are freely available via GitLab at https://gitlab.com/cbjh/gromacs-swaxs.

An overview of structural approaches to study therapeutic RNAs

RNAs provide considerable opportunities as therapeutic agent to expand the plethora of classical therapeutic targets, from extracellular and surface proteins to intracellular nucleic acids and its regulators, in a wide range of diseases. RNA versatility can be exploited to recognize cell types, perform cell therapy, and develop new vaccine classes. Therapeutic RNAs (aptamers, antisense nucleotides, siRNA, miRNA, mRNA and CRISPR-Cas9) can modulate or induce protein expression, inhibit molecular interactions, achieve genome editing as well as exon-skipping. A common RNA thread, which makes it very promising for therapeutic applications, is its structure, flexibility, and binding specificity. Moreover, RNA displays peculiar structural plasticity compared to proteins as well as to DNA. Here we summarize the recent advances and applications of therapeutic RNAs, and the experimental and computational methods to analyze their structure, by biophysical techniques (liquid-state NMR, scattering, reactivity, and computational simulations), with a focus on dynamic and flexibility aspects and to binding analysis. This will provide insights on the currently available RNA therapeutic applications and on the best techniques to evaluate its dynamics and reactivity.

Insights into the structural stability of major groove RNA triplexes by WAXS-guided MD simulations

RNA triple helices are commonly observed tertiary motifs that are associated with critical biological functions, including signal transduction. Because the recognition of their biological importance is relatively recent, their full range of structural properties has not yet been elucidated. The integration of solution wide-angle X-ray scattering (WAXS) with molecular dynamics (MD) simulations, described here, provides a new way to capture the structures of major-groove RNA triplexes that evade crystallographic characterization. This method yields excellent agreement between measured and computed WAXS profiles and allows for an atomically detailed visualization of these motifs. Using correlation maps, the relationship between well-defined features in the scattering profiles and real space characteristics of RNA molecules is defined, including the subtle conformational variations in the double-stranded RNA upon the incorporation of a third strand by base triples. This readily applicable approach has the potential to provide insight into interactions that stabilize RNA tertiary structure that enables function.