Structural dynamics of a single-stranded RNA-helix junction using NMR

Nanosecond chain dynamics of single-stranded nucleic acids

The conformational dynamics of single-stranded nucleic acids are fundamental for nucleic acid folding and function. However, their elementary chain dynamics have been difficult to resolve experimentally. Here we employ a combination of single-molecule Förster resonance energy transfer, nanosecond fluorescence correlation spectroscopy, and nanophotonic enhancement to determine the conformational ensembles and rapid chain dynamics of short single-stranded nucleic acids in solution. To interpret the experimental results in terms of end-to-end distance dynamics, we utilize the hierarchical chain growth approach, simple polymer models, and refinement with Bayesian inference to generate structural ensembles that closely align with the experimental data. The resulting chain reconfiguration times are exceedingly rapid, in the 10-ns range. Solvent viscosity-dependent measurements indicate that these dynamics of single-stranded nucleic acids exhibit negligible internal friction and are thus dominated by solvent friction. Our results provide a detailed view of the conformational distributions and rapid dynamics of single-stranded nucleic acids.

Sarcin/Ricin Domain RNA Retains Its Structure Better Than A-RNA in Adaptively Biased Molecular Dynamics Simulations

Less than one in thirty of the RNA sequences transcribed in humans are translated into protein. The noncoding RNA (ncRNA) functions in catalysis, structure, regulation, and more. However, for the most part, these functions are poorly characterized. RNA is modular and described by motifs that include helical A-RNA with canonical Watson-Crick base-pairing as well as structures with only noncanonical base pairs. Understanding the structure and dynamics of motifs will aid in deciphering functions of specific ncRNAs. We present computational studies on a standard sarcin/ricin domain (SRD), citrus bark cracking viroid SRD, as well as A-RNA. We have applied enhanced molecular dynamics techniques that construct an inverse free-energy surface (iFES) determined by collective variables that monitor base-pairing and backbone conformation. Each SRD RNA is flanked on each side by A-RNA, allowing comparison of the behavior of these motifs in the same molecule. The RNA iFESs have single peaks, indicating that the combined motifs should denature as a single cohesive unit, rather than by regional melting. Local root-mean-square deviation (RMSD) analysis and communication propensity (CProp, variance in distances between residue pairs) reveal distinct motif properties. Our analysis indicates that the standard SRD is more stable than the viroid SRD, which is more stable than A-RNA. Base pairs at SRD to A-RNA transitions have limited flexibility. Application of CProp reveals extraordinary stiffness of the SRD, allowing residues on opposite sides of the motif to sense each other's motions.

Influence of position of hydroxyl group of flavonoids on their binding with single stranded polyriboadenylic acid: A spectroscopic evaluation

Single stranded polyriboadenylic acid [poly (rA)] has been accepted widely as a suitable drug target owing to its vital role in the development of cancer since it controls gene expression during cell growth and differentiation. The biological properties of poly (rA) depend on its structural morphology. Pharmacologically active flavonoids can act as suitable binders to poly (rA) and significantly change its biophysical properties. Different factors favour flavonoid-poly (rA) binding. In our present work we have explored the role played by the position of hydroxyl groups in the flavonoids namely 3, 5, 6 and 7 hydroxyflavones in their course of interaction with poly (rA). A range of spectroscopic experiments reveal that 3HF binds best to poly (rA) among the four chosen flavonoids. This is probably due to the presence of a hydroxyl group in '3' position that enables it to exhibit ESIPT phenomenon which is missing for the other used flavonoids.

Anomalous Reverse Transcription through Chemical Modifications in Polyadenosine Stretches

Thermostable reverse transcriptases are workhorse enzymes underlying nearly all modern techniques for RNA structure mapping and for the transcriptome-wide discovery of RNA chemical modifications. Despite their wide use, these enzymes' behaviors at chemical modified nucleotides remain poorly understood. Wellington-Oguri et al. recently reported an apparent loss of chemical modification within putatively unstructured polyadenosine stretches modified by dimethyl sulfate or 2' hydroxyl acylation, as probed by reverse transcription. Here, reanalysis of these and other publicly available data, capillary electrophoresis experiments on chemically modified RNAs, and nuclear magnetic resonance spectroscopy on (A)₁₂ and variants show that this effect is unlikely to arise from an unusual structure of polyadenosine. Instead, tests of different reverse transcriptases on chemically modified RNAs and molecules synthesized with single 1-methyladenosines implicate a previously uncharacterized reverse transcriptase behavior: near-quantitative bypass through chemical modifications within polyadenosine stretches. All tested natural and engineered reverse transcriptases (MMLV; SuperScript II, III, and IV; TGIRT-III; and MarathonRT) exhibit this anomalous bypass behavior. Accurate DMS-guided structure modeling of the polyadenylated HIV-1 3' untranslated region requires taking into account this anomaly. Our results suggest that poly(rA-dT) hybrid duplexes can trigger an unexpectedly effective reverse transcriptase bypass and that chemical modifications in mRNA poly(A) tails may be generally undercounted.

Visualizing Disordered Single-Stranded RNA: Connecting Sequence, Structure, and Electrostatics

Disordered homopolymeric regions of single-stranded RNA, such as U or A tracts, are found within functional RNAs where they play distinct roles in defining molecular structure and facilitating recognition by partners. Despite this prominence, details of conformational and biophysical properties of these regions have not yet been resolved. We apply a number of experimental techniques to investigate the conformations of these biologically important motifs and provide quantitative measurements of their ion atmospheres. Single strands of RNA display pronounced sequence-dependent conformations that relate to the unique ion atmospheres each attracts. Chains of rU bases are relatively unstructured under all conditions, while chains of rA bases display distinct ordering through stacking or clustering motifs, depending on the composition of the surrounding solution. These dramatic structural differences are consistent with the measured disparity in ion composition and atmospheres around each homopolymer, revealing a complex interplay of base, ion, and single-strand ordering. The unique structural and ionic signatures of homopolymer ssRNAs explains their role(s) in folding structured RNAs and may explain their distinct recognition by protein partners.

RNA Dynamics by NMR Spectroscopy

An ever-increasing number of functional RNAs require a mechanistic understanding. RNA function relies on changes in its structure, so-called dynamics. To reveal dynamic processes and higher energy structures, new NMR methods have been developed to elucidate these dynamics in RNA with atomic resolution. In this Review, we provide an introduction to dynamics novices and an overview of methods that access most dynamic timescales, from picoseconds to hours. Examples are provided as well as insight into theory, data acquisition and analysis for these different methods. Using this broad spectrum of methodology, unprecedented detail and invisible structures have been obtained and are reviewed here. RNA, though often more complicated and therefore neglected, also provides a great system to study structural changes, as these RNA structural changes are more easily defined-Lego like-than in proteins, hence the numerous revelations of RNA excited states.

Sequence-dependent RNA helix conformational preferences predictably impact tertiary structure formation

Structured RNAs and RNA complexes underlie biological processes ranging from control of gene expression to protein translation. Approximately 50% of nucleotides within known structured RNAs are folded into Watson-Crick (WC) base pairs, and sequence changes that preserve these pairs are typically assumed to preserve higher-order RNA structure and binding of macromolecule partners. Here, we report that indirect effects of the helix sequence on RNA tertiary stability are, in fact, significant but are nevertheless predictable from a simple computational model called RNAMake-∆∆G. When tested through the RNA on a massively parallel array (RNA-MaP) experimental platform, blind predictions for >1500 variants of the tectoRNA heterodimer model system achieve high accuracy (rmsd 0.34 and 0.77 kcal/mol for sequence and length changes, respectively). Detailed comparison of predictions to experiments support a microscopic picture of how helix sequence changes subtly modulate conformational fluctuations at each base-pair step, which accumulate to impact RNA tertiary structure stability. Our study reveals a previously overlooked phenomenon in RNA structure formation and provides a framework of computation and experiment for understanding helix conformational preferences and their impact across biological RNA and RNA-protein assemblies.

Unique crystallographic signatures of Boc-Gly-Phe-Phe-OMe and Boc-Gly-Phg-Phe-OMe and their self-association

The crystal structures and morphology of Boc-Gly-Phe-Phe-OMe (1) and Boc-Gly-Phg-Phe-OMe (2) are reported. While 1 forms a rare open turn conformation, 2 forms a β-sheet conformation.

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

Structure of RNA Stem Loop B from the Picornavirus Replication Platform

The presumptive RNA cloverleaf at the start of the 5'-untranslated region of the picornavirus genome is an essential element in replication. Stem loop B (SLB) of the cloverleaf is a recognition site for the host polyC-binding protein, which initiates a switch from translation to replication. Here we present the solution structure of human rhinovirus isotype 14 SLB using nuclear magnetic resonance spectroscopy. SLB adopts a predominantly A-form helical structure. The stem contains five Watson-Crick base pairs and one wobble base pair and is capped by an eight-nucleotide loop. The wobble base pair introduces perturbations into the helical parameters but does not appear to introduce flexibility. However, the helix major groove appears to be accessible. Flexibility is seen throughout the loop and in the terminal nucleotides. The pyrimidine-rich region of the loop, the apparent recognition site for the polyC-binding protein, is the most disordered region of the structure.

Direct NOE simulation from long MD trajectories

A software package, MD2NOE, is presented which calculates Nuclear Overhauser Effect (NOE) build-up curves directly from molecular dynamics (MD) trajectories. It differs from traditional approaches in that it calculates correlation functions directly from the trajectory instead of extracting inverse sixth power distance terms as an intermediate step in calculating NOEs. This is particularly important for molecules that sample conformational states on a timescale similar to molecular reorientation. The package is tested on sucrose and results are shown to differ in small but significant ways from those calculated using an inverse sixth power assumption. Results are also compared to experiment and found to be in reasonable agreement despite an expected underestimation of water viscosity by the water model selected.

Efficient large-scale preparation and purification of short single-stranded RNA oligonucleotides

Sequence-specific RNA recognition by RNA-binding proteins plays a crucial role in the post-translational regulation of gene expression. Biophysical and biochemical studies help to unravel the principles of sequence-specific RNA recognition, but the methods used require large amounts of single-stranded RNA (ssRNA). Here we present a fast and robust method for large-scale preparation and purification of short ssRNA oligonucleotides for biochemical, biophysical, and structural studies. We designed an efficiently folding, self-cleaving hammerhead (HH) ribozyme to prepare ssRNA oligonucleotides. Hammerhead ribozyme RNAs self-cleave with over 95% efficiency during in vitro transcription as a function of magnesium concentration to produce high yields of the desired ssRNA products. The resulting ssRNAs can be purified from crude transcription reactions by denaturing anion-exchange chromatography and then desalted by weak anion-exchange chromatography using volatile ammonium bicarbonate buffer solutions. The ssRNA oligonucleotides produced this way are homogenous, as judged by mass spectrometry (MS), and are suitable for biochemical and biophysical studies. Moreover, for high-resolution NMR structure determination of RNA-protein complexes, our protocol enables efficient preparation of ssRNA oligonucleotides with various isotope-labeling schemes which are not commercially available.

Ligand-detected relaxation dispersion NMR spectroscopy: dynamics of preQ1-RNA binding

An NMR-based approach to characterizing the binding kinetics of ligand molecules to biomolecules, like RNA or proteins, by ligand-detected Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments is described. A (15)N-modified preQ1 ligand is used to acquire relaxation dispersion experiments in the presence of low amounts of the Fsu class I preQ1 aptamer RNA, and increasing ligand concentrations to probe the RNA small molecule interaction. Our experimental data strongly support the conformational selection mechanism postulated. The approach gives direct access to two parameters of a ligand-receptor interaction: the off rate and the population of the small molecule-receptor complex. A detailed description of the kinetics underlying the ligand binding process is of crucial importance to fully understanding a riboswitch's function and to evaluate potential new antibiotics candidates targeting the noncoding RNA species. Ligand-detected NMR relaxation dispersion experiments represent a valuable diagnostic tool for the characterization of binding mechanisms.