Hidden Structural Modules in a Cooperative RNA Folding Transition

Optimized periphery-core interface increases fitness of the Bacillus subtilis glmS ribozyme

Like other functional RNAs, ribozymes encode a conserved catalytic center supported by peripheral domains that vary among ribozyme sub-families. To understand how core-periphery interactions contribute to ribozyme fitness, we compared the cleavage kinetics of all single base substitutions at 152 sites across the Bacillus subtilis glmS ribozyme by high-throughput sequencing (k-seq). The in vitro activity map mirrored phylogenetic sequence conservation in glmS ribozymes, indicating that biological fitness reports all biochemically important positions. The k-seq results and folding assays showed that most deleterious mutations lower activity by impairing ribozyme self-assembly. All-atom molecular dynamics simulations of the complete ribozyme revealed how individual mutations in the core or the IL4 peripheral loop introduce a non-native tertiary interface that rewires the catalytic center, eliminating activity. We conclude that the need to avoid non-native helix packing powerfully constrains the evolution of tertiary structure motifs in RNA.

Structural and dynamic mechanisms for coupled folding and tRNA recognition of a translational T-box riboswitch

T-box riboswitches are unique riboregulators where gene regulation is mediated through interactions between two highly structured RNAs. Despite extensive structural insights, how RNA-RNA interactions drive the folding and structural transitions of T-box to achieve functional conformations remains unclear. Here, by combining SAXS, single-molecule FRET and computational modeling, we elaborate the folding energy landscape of a translational T-box aptamer consisting of stems I, II and IIA/B, which Mg2+-induced global folding and tRNA binding are cooperatively coupled. smFRET measurements reveal that high Mg2+ stabilizes IIA/B and its stacking on II, which drives the pre-docking of I and II into a competent conformation, subsequent tRNA binding promotes docking of I and II to form a high-affinity tRNA binding groove, of which the essentiality of IIA/B and S-turn in II is substantiated with mutational analysis. We highlight a delicate balance among Mg2+, the intra- and intermolecular RNA-RNA interactions in modulating RNA folding and function.

Impact of Ion-Mixing Entropy on Orientational Preferences of DNA Helices: FRET Measurements and Computer Simulations

Biological processes require DNA and RNA helices to pack together in specific interhelical orientations. While electrostatic repulsion between backbone charges is expected to be maximized when helices are in parallel alignment, such orientations are commonplace in nature. To better understand how the repulsion is overcome, we used experimental and computational approaches to investigate how the orientational preferences of DNA helices depend on the concentration and valence of mobile cations. We used Förster resonance energy transfer (FRET) to probe the relative orientations of two 24-bp helices held together via a freely rotating PEG linker. At low cation concentrations, the helices preferred more "cross"-like orientations over those closer to parallel, and this preference was reduced with increasing salt concentrations. The results were in good quantitative agreement with Poisson-Boltzmann (PB) calculations for monovalent salt (Na+). However, PB underestimated the ability of mixtures of monovalent and divalent ions (Mg2+) to reduce the conformational preference. As a complementary approach, we performed all-atom molecular dynamics (MD) simulations and found better agreement with the experimental results. While MD and PB predict similar electrostatic forces, MD predicts a greater accumulation of Mg2+ in the ion atmosphere surrounding the DNA. Mg2+ occupancy is predicted to be greater in conformations close to the parallel orientation than in conformations close to the crossed orientation, enabling a greater release of Na+ ions and providing an entropic gain (one bound ion for two released). MD predicts an entropy gain larger than that of PB because of the increased Mg2+ occupancy. The entropy changes have a negligible effect at low Mg2+ concentrations because the free energies are dominated by electrostatic repulsion. However, as the Mg2+ concentration increases, charge screening is more effective and the mixing entropy produces readily detectable changes in packing preferences. Our results underline the importance of mixing entropy of counterions in nucleic acid interactions and provide a new understanding on the impact of a mixed ion atmosphere on the packing of DNA helices.

Dissecting the energetic architecture within an RNA tertiary structural motif via high-throughput thermodynamic measurements

Structured RNAs and RNA/protein complexes perform critical cellular functions. They often contain structurally conserved tertiary contact "motifs," whose occurrence simplifies the RNA folding landscape. Prior studies have focused on the conformational and energetic modularity of intact motifs. Here, we turn to the dissection of one common motif, the 11nt receptor (11ntR), using quantitative analysis of RNA on a massively parallel array to measure the binding of all single and double 11ntR mutants to GAAA and GUAA tetraloops, thereby probing the energetic architecture of the motif. While the 11ntR behaves as a motif, its cooperativity is not absolute. Instead, we uncovered a gradient from high cooperativity amongst base-paired and neighboring residues to additivity between distant residues. As expected, substitutions at residues in direct contact with the GAAA tetraloop resulted in the largest decreases to binding, and energetic penalties of mutations were substantially smaller for binding to the alternate GUAA tetraloop, which lacks tertiary contacts present with the canonical GAAA tetraloop. However, we found that the energetic consequences of base partner substitutions are not, in general, simply described by base pair type or isostericity. We also found exceptions to the previously established stability-abundance relationship for 11ntR sequence variants. These findings of "exceptions to the rule" highlight the power of systematic high-throughput approaches to uncover novel variants for future study in addition to providing an energetic map of a functional RNA.

Measurement of the specific and non-specific binding energies of Mg2+ to RNA

Determining the non-specific and specific electrostatic contributions of magnesium binding to RNA is a challenging problem. We introduce a single-molecule method based on measuring the folding energy of a native RNA in magnesium and at its equivalent sodium concentration. The latter is defined so that the folding energy in sodium equals the non-specific electrostatic contribution in magnesium. The sodium equivalent can be estimated according to the empirical 100/1 rule (1 M NaCl is equivalent to 10 mM MgCl2), which is a good approximation for most RNAs. The method is applied to an RNA three-way junction (3WJ) that contains specific Mg2+ binding sites and misfolds into a double hairpin structure without binding sites. We mechanically pull the RNA with optical tweezers and use fluctuation theorems to determine the folding energies of the native and misfolded structures in magnesium (10 mM MgCl2) and at the equivalent sodium condition (1 M NaCl). While the free energies of the misfolded structure are equal in magnesium and sodium, they are not for the native structure, the difference being due to the specific binding energy of magnesium to the 3WJ, which equals ΔG≃ 10 kcal/mol. Besides stabilizing the 3WJ, Mg2+ also kinetically rescues it from the misfolded structure over timescales of tens of seconds in a force-dependent manner. The method should generally be applicable to determine the specific binding energies of divalent cations to other tertiary RNAs.

NMR-derived secondary structure of the full-length Ox40 mRNA 3'UTR and its multivalent binding to the immunoregulatory RBP Roquin

Control of posttranscriptional mRNA decay is a crucial determinant of cell homeostasis and differentiation. mRNA lifetime is governed by cis-regulatory elements in their 3' untranslated regions (UTR). Despite ongoing progress in the identification of cis elements we have little knowledge about the functional and structural integration of multiple elements in 3'UTR regulatory hubs and their recognition by mRNA-binding proteins (RBPs). Structural analyses are complicated by inconsistent mapping and prediction of RNA fold, by dynamics, and size. We here, for the first time, provide the secondary structure of a complete mRNA 3'UTR. We use NMR spectroscopy in a divide-and-conquer strategy complemented with SAXS, In-line probing and SHAPE-seq applied to the 3'UTR of Ox40 mRNA, which encodes a T-cell co-receptor repressed by the protein Roquin. We provide contributions of RNA elements to Roquin-binding. The protein uses its extended bi-modal ROQ domain to sequentially engage in a 2:1 stoichiometry with a 3'UTR core motif. We observe differential binding of Roquin to decay elements depending on their structural embedment. Our data underpins the importance of studying RNA regulation in a full sequence and structural context. This study serves as a paradigm for an approach in analysing structured RNA-regulatory hubs and their binding by RBPs.

Pseudoknot length modulates the folding, conformational dynamics, and robustness of Xrn1 resistance of flaviviral xrRNAs

To understand how RNA dynamics is regulated and connected to its function, we investigate the folding, conformational dynamics and robustness of Xrn1 resistance of a set of flaviviral xrRNAs using SAXS, smFRET and in vitro enzymatic assays. Flaviviral xrRNAs form discrete ring-like 3D structures, in which the length of a conserved long-range pseudoknot (PK2) ranges from 2 bp to 7 bp. We find that xrRNAs' folding, conformational dynamics and Xrn1 resistance are strongly correlated and highly Mg²⁺-dependent, furthermore, the Mg²⁺-dependence is modulated by PK2 length variations. xrRNAs with long PK2 require less Mg²⁺ to stabilize their folding, exhibit reduced conformational dynamics and strong Xrn1 resistance even at low Mg²⁺, and tolerate mutations at key tertiary motifs at high Mg²⁺, which generally are destructive to xrRNAs with short PK2. These results demonstrate an unusual regulatory mechanism of RNA dynamics providing insights into the functions and future biomedical applications of xrRNAs.

Assembly factors chaperone ribosomal RNA folding by isolating helical junctions that are prone to misfolding

While RNAs are known to misfold, the underlying molecular causes have been mainly studied in fragments of biologically relevant larger RNAs. As these small RNAs are dominated by secondary structures, misfolding of these secondary structures remains the most-explored cause for global RNA misfolding. Conversely, how RNA chaperones function in a biological context to promote native folding beyond duplex annealing remains unknown. Here, in a combination of dimethylsulfate mutational profiling with sequencing (DMS-MaPseq), structural analyses, biochemical experiments, and yeast genetics, we show that three-helix junctions are prone to misfolding during assembly of the small ribosomal subunit in vivo. We identify ubiquitous roles for ribosome assembly factors in chaperoning their folding by preventing the formation of premature tertiary interactions, which otherwise kinetically trap misfolded junctions, thereby blocking further progress in the assembly cascade. While these protein chaperones act indirectly by binding the interaction partners of junctions, our analyses also suggest direct roles for small nucleolar RNAs (snoRNAs) in binding and chaperoning helical junctions during transcription. While these assembly factors do not utilize energy to ameliorate misfolding, our data demonstrate how their dissociation renders reversible folding steps irreversible, thereby driving native folding and assembly and setting up a timer that dictates the propensity of misfolded intermediates to escape quality control. Finally, the data demonstrate that RNA chaperones act locally on individual tertiary interactions, in contrast to protein chaperones, which globally unfold misfolded proteins.

RNA architecture influences plant biology

The majority of the genome is transcribed to RNA in living organisms. RNA transcripts can form astonishing arrays of secondary and tertiary structures via Watson-Crick, Hoogsteen, or wobble base pairing. In vivo, RNA folding is not a simple thermodynamic event of minimizing free energy. Instead, the process is constrained by transcription, RNA-binding proteins, steric factors, and the microenvironment. RNA secondary structure (RSS) plays myriad roles in numerous biological processes, such as RNA processing, stability, transportation, and translation in prokaryotes and eukaryotes. Emerging evidence has also implicated RSS in RNA trafficking, liquid-liquid phase separation, and plant responses to environmental variations such as temperature and salinity. At molecular level, RSS is correlated with splicing, polyadenylation, protein synthesis, and miRNA biogenesis and functions. In this review, we summarize newly reported methods for probing RSS in vivo and functions and mechanisms of RSS in plant physiology.

Chemical shift prediction of RNA imino groups: application toward characterizing RNA excited states

NH groups in proteins or nucleic acids are the most challenging target for chemical shift prediction. Here we show that the RNA base pair triplet motif dictates imino chemical shifts in its central base pair. A lookup table is established that links each type of base pair triplet to experimental chemical shifts of the central base pair, and can be used to predict imino chemical shifts of RNAs to remarkable accuracy. Strikingly, the semiempirical method can well interpret the variations of chemical shifts for different base pair triplets, and is even applicable to non-canonical motifs. This finding opens an avenue for predicting chemical shifts of more complicated RNA motifs. Furthermore, we combine the imino chemical shift prediction with NMR relaxation dispersion experiments targeting both ¹⁵N and ¹H^N of the imino group, and verify a previously characterized excited state of P5abc subdomain including an earlier speculated non-native G•G mismatch.

Prediction of chemical shifts is critical for extracting structural and dynamic information from biomolecular NMR data. Here the authors report an RNA imino group chemical shift predictor, showing that the imino chemical shifts of a residue are dictated by the surrounding base pair triplet.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

The roles of structural dynamics in the cellular functions of RNAs

RNAs fold into 3D structures that range from simple helical elements to complex tertiary structures and quaternary ribonucleoprotein assemblies. The functions of many regulatory RNAs depend on how their 3D structure changes in response to a diverse array of cellular conditions. In this Review, we examine how the structural characterization of RNA as dynamic ensembles of conformations, which form with different probabilities and at different timescales, is improving our understanding of RNA function in cells. We discuss the mechanisms of gene regulation by microRNAs, riboswitches, ribozymes, post-transcriptional RNA modifications and RNA-binding proteins, and how the cellular environment and processes such as liquid-liquid phase separation may affect RNA folding and activity. The emerging RNA-ensemble-function paradigm is changing our perspective and understanding of RNA regulation, from in vitro to in vivo and from descriptive to predictive.

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.

Dynamic ensemble of HIV-1 RRE stem IIB reveals non-native conformations that disrupt the Rev-binding site

The HIV-1 Rev response element (RRE) RNA element mediates the nuclear export of intron containing viral RNAs by forming an oligomeric complex with the viral protein Rev. Stem IIB and nearby stem II three-way junction nucleate oligomerization through cooperative binding of two Rev molecules. Conformational flexibility at this RRE region has been shown to be important for Rev binding. However, the nature of the flexibility has remained elusive. Here, using NMR relaxation dispersion, including a new strategy for directly observing transient conformational states in large RNAs, we find that stem IIB alone or when part of the larger RREII three-way junction robustly exists in dynamic equilibrium with non-native excited state (ES) conformations that have a combined population of ∼20%. The ESs disrupt the Rev-binding site by changing local secondary structure, and their stabilization via point substitution mutations decreases the binding affinity to the Rev arginine-rich motif (ARM) by 15- to 80-fold. The ensemble clarifies the conformational flexibility observed in stem IIB, reveals long-range conformational coupling between stem IIB and the three-way junction that may play roles in cooperative Rev binding, and also identifies non-native RRE conformational states as new targets for the development of anti-HIV therapeutics.

Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

Challenges and approaches to predicting RNA with multiple functional structures

The revolution in sequencing technology demands new tools to interpret the genetic code. As in vivo transcriptome-wide chemical probing techniques advance, new challenges emerge in the RNA folding problem. The emphasis on one sequence folding into a single minimum free energy structure is fading as a new focus develops on generating RNA structural ensembles and identifying functional structural features in ensembles. This review describes an efficient combinatorially complete method and three free energy minimization approaches to predicting RNA structures with more than one functional fold, as well as two methods for analysis of a thermodynamics-based Boltzmann ensemble of structures. The review then highlights two examples of viral RNA 3'-UTR regions that fold into more than one conformation and have been characterized by single molecule fluorescence energy resonance transfer or NMR spectroscopy. These examples highlight the different approaches and challenges in predicting structure and function from sequence for RNA with multiple biological roles and folds. More well-defined examples and new metrics for measuring differences in RNA structures will guide future improvements in prediction of RNA structure and function from sequence.

The Story of RNA Folding, as Told in Epochs

The past decades have witnessed tremendous developments in our understanding of RNA biology. At the core of these advances have been studies aimed at discerning RNA structure and at understanding the forces that influence the RNA folding process. It is easy to take the present state of understanding for granted, but there is much to be learned by considering the path to our current understanding, which has been tortuous, with the birth and death of models, the adaptation of experimental tools originally developed for characterization of protein structure and catalysis, and the development of novel tools for probing RNA. In this review we tour the stages of RNA folding studies, considering them as "epochs" that can be generalized across scientific disciplines. These epochs span from the discovery of catalytic RNA, through biophysical insights into the putative primordial RNA World, to characterization of structured RNAs, the building and testing of models, and, finally, to the development of models with the potential to yield generalizable predictive and quantitative models for RNA conformational, thermodynamic, and kinetic behavior. We hope that this accounting will aid others as they navigate the many fascinating questions about RNA and its roles in biology, in the past, present, and future.