Salt Dependence of A-Form RNA Duplexes: Structures and Implications

Molecular Origin of Distinct Hydration Dynamics in Double Helical DNA and RNA Sequences

Water molecules are essential to determine the structure of nucleic acids and mediate their interactions with other biomolecules. Here, we characterize the hydration dynamics of analogous DNA and RNA double helices with unprecedented resolution and elucidate the molecular origin of their differences: first, the localization of the slowest hydration water molecules─in the minor groove in DNA, next to phosphates in RNA─and second, the markedly distinct hydration dynamics of the two phosphate oxygen atoms OR and OS in RNA. Using our Extended Jump Model for water reorientation, we assess the relative importance of previously proposed factors, including the local topography, water bridges, and the presence of ions. We show that the slow hydration dynamics at RNA OR sites is not due to bridging water molecules but is caused by both the larger excluded volume and the stronger initial H-bond next to OR, due to the different phosphate orientations in A-form double helical RNA.

RNA structures and dynamics with Å resolution revealed by x-ray free-electron lasers

RNA macromolecules, like proteins, fold to assume shapes that are intimately connected to their broadly recognized biological functions; however, because of their high charge and dynamic nature, RNA structures are far more challenging to determine. We introduce an approach that exploits the high brilliance of x-ray free-electron laser sources to reveal the formation and ready identification of angstrom-scale features in structured and unstructured RNAs. Previously unrecognized structural signatures of RNA secondary and tertiary structures are identified through wide-angle solution scattering experiments. With millisecond time resolution, we observe an RNA fold from a dynamically varying single strand through a base-paired intermediate to assume a triple-helix conformation. While the backbone orchestrates the folding, the final structure is locked in by base stacking. This method may help to rapidly characterize and identify structural elements in nucleic acids in both equilibrium and time-resolved experiments.

RNA structures and dynamics with Å resolution revealed by x-ray free electron lasers

RNA macromolecules, like proteins, fold to assume shapes that are intimately connected to their broadly recognized biological functions; however, because of their high charge and dynamic nature, RNA structures are far more challenging to determine. We introduce an approach that exploits the high brilliance of x-ray free electron laser sources to reveal the formation and ready identification of Å scale features in structured and unstructured RNAs. New structural signatures of RNA secondary and tertiary structures are identified through wide angle solution scattering experiments. With millisecond time resolution, we observe an RNA fold from a dynamically varying single strand through a base paired intermediate to assume a triple helix conformation. While the backbone orchestrates the folding, the final structure is locked in by base stacking. In addition to understanding how RNA triplexes form and thereby function as dynamic signaling elements, this new method can vastly increase the rate of structure determination for these biologically essential, but mostly uncharacterized macromolecules.

RNA Captures More Cations than DNA: Insights from Molecular Dynamics Simulations

The distribution of cations around nucleic acids is essential for a broad variety of processes ranging from DNA condensation and RNA folding to the detection of biomolecules in biosensors. Predicting the exact distribution of ions remains challenging since the distribution and, hence, a broad variety of nucleic acid properties depend on the salt concentration, the valency of the ions, and the ion type. Despite the importance, a general theory to quantify ion-specific effects for highly charged biomolecules is still lacking. Moreover, recent experiments reveal that despite their similar building blocks, DNA and RNA duplexes can react differently to the same ionic conditions. The aim of our current work is to provide a comprehensive set of molecular dynamics simulations using more than 180 μs of simulation time. For the mono- and divalent cations Li⁺, Na⁺, K⁺, Cs⁺, Ca²⁺, Sr²⁺, and Ba²⁺, the simulations allow us to reveal the ion-specific distributions and binding patterns for DNA and RNA duplexes. The microscopic insights from the simulations display the origin of ion-specificity and shed light on the question of why DNA and RNA show opposing behavior in the same ionic conditions. Finally, the detailed binding patterns from the simulations reveal why RNA can capture more cations than DNA.

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.

Visualizing RNA Structures by SAXS-Driven MD Simulations

The biological role of biomolecules is intimately linked to their structural dynamics. Experimental or computational techniques alone are often insufficient to determine accurate structural ensembles in atomic detail. We use all-atom molecular dynamics (MD) simulations and couple it to small-angle X-ray scattering (SAXS) experiments to resolve the structural dynamics of RNA molecules. To accomplish this task, we utilize a set of re-weighting and biasing techniques tailored for RNA molecules. To showcase our approach, we study two RNA molecules: a riboswitch that shows structural variations upon ligand binding, and a two-way junction RNA that displays structural heterogeneity and sensitivity to salt conditions. Integration of MD simulations and experiments allows the accurate construction of conformational ensembles of RNA molecules. We observe a dynamic change of the SAM-I riboswitch conformations depending on its binding partners. The binding of SAM and Mg2+ cations stabilizes the compact state. The absence of Mg2+ or SAM leads to the loss of tertiary contacts, resulting in a dramatic expansion of the riboswitch conformations. The sensitivity of RNA structures to the ionic strength demonstrates itself in the helix junction helix (HJH). The HJH shows non-monotonic compaction as the ionic strength increases. The physics-based picture derived from the experimentally guided MD simulations allows biophysical characterization of RNA molecules. All in all, SAXS-guided MD simulations offer great prospects for studying RNA structural dynamics.

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

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.

Machine learning deciphers structural features of RNA duplexes measured with solution X-ray scattering

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.

Structural and functional conservation of the programmed -1 ribosomal frameshift signal of SARS coronavirus 2 (SARS-CoV-2)

Approximately 17 years after the severe acute respiratory syndrome coronavirus (SARS-CoV) epidemic, the world is currently facing the COVID-19 pandemic caused by SARS corona virus 2 (SARS-CoV-2). According to the most optimistic projections, it will take more than a year to develop a vaccine, so the best short-term strategy may lie in identifying virus-specific targets for small molecule-based interventions. All coronaviruses utilize a molecular mechanism called programmed -1 ribosomal frameshift (-1 PRF) to control the relative expression of their proteins. Previous analyses of SARS-CoV have revealed that it employs a structurally unique three-stemmed mRNA pseudoknot that stimulates high -1 PRF rates and that it also harbors a -1 PRF attenuation element. Altering -1 PRF activity impairs virus replication, suggesting that this activity may be therapeutically targeted. Here, we comparatively analyzed the SARS-CoV and SARS-CoV-2 frameshift signals. Structural and functional analyses revealed that both elements promote similar -1 PRF rates and that silent coding mutations in the slippery sites and in all three stems of the pseudoknot strongly ablate -1 PRF activity. We noted that the upstream attenuator hairpin activity is also functionally retained in both viruses, despite differences in the primary sequence in this region. Small-angle X-ray scattering analyses indicated that the pseudoknots in SARS-CoV and SARS-CoV-2 have the same conformation. Finally, a small molecule previously shown to bind the SARS-CoV pseudoknot and inhibit -1 PRF was similarly effective against -1 PRF in SARS-CoV-2, suggesting that such frameshift inhibitors may be promising lead compounds to combat the current COVID-19 pandemic.

Structural and functional conservation of the programmed -1 ribosomal frameshift signal of SARS-CoV-2

17 years after the SARS-CoV epidemic, the world is facing the COVID-19 pandemic. COVID-19 is caused by a coronavirus named SARS-CoV-2. Given the most optimistic projections estimating that it will take over a year to develop a vaccine, the best short-term strategy may lie in identifying virus-specific targets for small molecule interventions. All coronaviruses utilize a molecular mechanism called -1 PRF to control the relative expression of their proteins. Prior analyses of SARS-CoV revealed that it employs a structurally unique three-stemmed mRNA pseudoknot to stimulate high rates of -1 PRF, and that it also harbors a -1 PRF attenuation element. Altering -1 PRF activity negatively impacts virus replication, suggesting that this molecular mechanism may be therapeutically targeted. Here we present a comparative analysis of the original SARS-CoV and SARS-CoV-2 frameshift signals. Structural and functional analyses revealed that both elements promote similar rates of -1 PRF and that silent coding mutations in the slippery sites and in all three stems of the pseudoknot strongly ablated -1 PRF activity. The upstream attenuator hairpin activity has also been functionally retained. Small-angle x-ray scattering indicated that the pseudoknots in SARS-CoV and SARS-CoV-2 had the same conformation. Finally, a small molecule previously shown to bind the SARS-CoV pseudoknot and inhibit -1 PRF was similarly effective against -1 PRF in SARS-CoV-2, suggesting that such frameshift inhibitors may provide promising lead compounds to counter the current pandemic.

Salt effect on thermodynamics and kinetics of a single RNA base pair

Due to the polyanionic nature of RNAs, the structural folding of RNAs are sensitive to solution salt conditions, while there is still lack of a deep understanding of the salt effect on the thermodynamics and kinetics of RNAs at a single base-pair level. In this work, the thermodynamic and the kinetic parameters for the base-pair AU closing/opening at different salt concentrations were calculated by 3-µsec all-atom molecular dynamics (MD) simulations at different temperatures. It was found that for the base-pair formation, the enthalpy change [Formula: see text] is nearly independent of salt concentration, while the entropy change [Formula: see text] exhibits a linear dependence on the logarithm of salt concentration, verifying the empirical assumption based on thermodynamic experiments. Our analyses revealed that such salt concentration dependence of the entropy change mainly results from the dependence of ion translational entropy change for the base pair closing/opening on salt concentration. Furthermore, the closing rate increases with the increasing of salt concentration, while the opening rate is nearly independent of salt concentration. Additionally, our analyses revealed that the free energy surface for describing the base-pair opening and closing dynamics becomes more rugged with the decrease of salt concentration.