Comprehensive features of natural and in vitro selected GNRA tetraloop-binding receptors

Enhanced hammerhead ribozyme turnover rates: Reevaluating therapeutic space for small catalytic RNAs

We discovered enhanced hammerhead ribozyme (EhhRz) kinetic performance during therapeutic optimization of a lead agent against the human rhodopsin mRNA target. Turnover activity in trans was greater than 300 nMmin-1 under 10-fold substrate excess and physiological Mg2+ (1 mM). We validated a moderate-throughput fluorescence quantitative hhRz kinetic assay, which is linear with substrate and product moles. The EhhRz targets a CUC↓ cleavage site in a substrate with no predicted secondary/tertiary structure and demonstrates classical Michaelis-Menten turnover behavior with efficiency (V max/K m ) up to 3.2 × 106 min-1M-1. EhhRzs show cooperative Mg2+ titration (K Mg1/2  = 0.7 ± 0.02 mM; Hill = 1.7 ± 0.1). Structure-function assays showed that the upstream EhhRz antisense flank (substrate bound) interacts with stem-loop II. Tetraloop sequence variation reveals a marked effect on turnover rate. Downstream substrate U7 is not essential for enhanced activity. Single-turnover reaction rates show substantial improvements in hhRz rates at physiological Mg2+. Embedded within a target-reporter fusion mRNA, EhhRzs cleave regional target elements under intracellular conditions promoting strong knockdown at target mRNA and protein levels. EhhRzs have potential as druggable nucleic acid therapeutics against arbitrary targets, or in the design of improved aptazymes.

Exploring the energetic and conformational properties of the sequence space connecting naturally occurring RNA tetraloop receptor motifs

Folded RNAs contain tertiary contact motifs whose structures and energetics are conserved across different RNAs. The transferable properties of RNA motifs simplify the RNA folding problem, but measuring energetic and conformational properties of many motifs remains a challenge. Here, we use a high-throughput thermodynamic approach to investigate how sequence changes alter the binding properties of naturally occurring motifs, the GAAA tetraloop • tetraloop receptor (TLR) interactions. We measured the binding energies and conformational preferences of TLR sequences that span mutational pathways from the canonical 11ntR to two other natural TLRs, the IC3R and Vc2R. While the IC3R and Vc2R share highly similar energetic and conformational properties, the landscapes that map the sequence changes for their conversion from the 11ntR to changes in these properties differ dramatically. Differences in the energetic landscapes stem from the mutations needed to convert the 11ntR to the IC3R and Vc2R rather than a difference in the intrinsic energetic architectures of these TLRs. The conformational landscapes feature several nonnative TLR variants with conformational preferences that differ from both the initial and final TLRs; these species represent potential branching points along the multidimensional sequence space to sequences with greater fitness in other RNA contexts with alternative conformational preferences. Our high-throughput, quantitative approach reveals the complex nature of sequence-fitness landscapes and leads to models for their molecular origins. Systematic and quantitative molecular approaches provide critical insights into understanding the evolution of natural RNAs as they traverse complex landscapes in response to selective pressures.

A comprehensive survey of long-range tertiary interactions and motifs in non-coding RNA structures

Understanding the 3D structure of RNA is key to understanding RNA function. RNA 3D structure is modular and can be seen as a composition of building blocks of various sizes called tertiary motifs. Currently, long-range motifs formed between distant loops and helical regions are largely less studied than the local motifs determined by the RNA secondary structure. We surveyed long-range tertiary interactions and motifs in a non-redundant set of non-coding RNA 3D structures. A new dataset of annotated LOng-RAnge RNA 3D modules (LORA) was built using an approach that does not rely on the automatic annotations of non-canonical interactions. An original algorithm, ARTEM, was developed for annotation-, sequence- and topology-independent superposition of two arbitrary RNA 3D modules. The proposed methods allowed us to identify and describe the most common long-range RNA tertiary motifs. Along with the prevalent canonical A-minor interactions, a large number of previously undescribed staple interactions were observed. The most frequent long-range motifs were found to belong to three main motif families: planar staples, tilted staples, and helical packing motifs.

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.

RNA 3D Modeling with FARFAR2, Online

Understanding the three-dimensional structure of an RNA molecule is often essential to understanding its function. Sampling algorithms and energy functions for RNA structure prediction are improving, due to the increasing diversity of structural data available for training statistical potentials and testing structural data, along with a steady supply of blind challenges through the RNA-Puzzles initiative. The recent FARFAR2 algorithm enables near-native structure predictions on fairly complex RNA structures, including automated selection of final candidate models and estimation of model accuracy. Here, we describe the use of a publicly available webserver for RNA modeling for realistic scenarios using FARFAR2, available at https://rosie.rosettacommons.org/farfar2 . We walk through two cases in some detail: a simple model pseudoknot from the frameshifting element of beet western yellows virus modeled using the "basic interface" to the webserver and a replication of RNA-Puzzle 20, a metagenomic twister sister ribozyme, using the "advanced interface." We also describe example runs of FARFAR2 modeling including two kinds of experimental data: a c-di-GMP riboswitch modeled with low-resolution restraints from MOHCA-seq experiments and a tandem GA motif modeled with 1H NMR chemical shifts.

Importance of Anion-π Interactions in RNA GAAA and GGAG Tetraloops: A Combined MD and QM Study

In this study, we demonstrate that anion−π interactions (an attractive noncovalent force between electron deficient π-systems and anions) are involved in the stabilization of GAAA and GGAG RNA tetraloops. Using the single recognition particle (SRP)–RNA complexes as a case of study, we combined molecular dynamics (MD) and quantum mechanics (QM) calculations to shed light on the structural influence of phosphate–G anion−π interactions and hydrogen bonds (HBs) involving K⁺/Mg²⁺ water clusters. In addition, the RNA assemblies herein were further characterized by means of the “atoms in molecules” (AIM) and noncovalent interactions plot (NCIplot) methodologies. We believe the results derived from this study might be important in the fields of chemical biology (RNA folding and engineering) and supramolecular chemistry (anion−π interactions) as well as to further expand the current knowledge regarding RNA structural motifs.

Features and Functions of the A-Minor Motif, the Most Common Motif in RNA Structure

A-minor motifs are RNA tertiary structure motifs that generally involve a canonical base pair and an adenine base forming hydrogen bonds with the minor groove of the base pair. Such motifs are among the most common tertiary interactions in known RNA structures, comparable in number with the non-canonical base pairs. They are often found in functionally important regions of non-coding RNAs and, in particular, play a central role in protein synthesis. Here, we review local variations of the A-minor geometry and discuss difficulties associated with their annotation, as well as various structural contexts and common A-minor co-motifs, and diverse functions of A-minors in various processes in a living cell.

High-throughput dissection of the thermodynamic and conformational properties of a ubiquitous class of RNA tertiary contact motifs

Despite RNA's diverse secondary and tertiary structures and its complex conformational changes, nature utilizes a limited set of structural "motifs"-helices, junctions, and tertiary contact modules-to build diverse functional RNAs. Thus, in-depth descriptions of a relatively small universe of RNA motifs may lead to predictive models of RNA tertiary conformational landscapes. Motifs may have different properties depending on sequence and secondary structure, giving rise to subclasses that expand the universe of RNA building blocks. Yet we know very little about motif subclasses, given the challenges in mapping conformational properties in high throughput. Previously, we used "RNA on a massively parallel array" (RNA-MaP), a quantitative, high-throughput technique, to study thousands of helices and two-way junctions. Here, we adapt RNA-MaP to study the thermodynamic and conformational properties of tetraloop/tetraloop receptor (TL/TLR) tertiary contact motifs, analyzing 1,493 TLR sequences from different classes. Clustering analyses revealed variability in TL specificity, stability, and conformational behavior. Nevertheless, natural GAAA/11ntR TL/TLRs, while varying in tertiary stability by ∼2.5 kcal/mol, exhibited conserved TL specificity and conformational properties. Thus, RNAs may tune stability without altering the overall structure of these TL/TLRs. Furthermore, their stability correlated with natural frequency, suggesting thermodynamics as the dominant selection pressure. In contrast, other TL/TLRs displayed heterogenous conformational behavior and appear to not be under strong thermodynamic selection. Our results build toward a generalizable model of RNA-folding thermodynamics based on the properties of isolated motifs, and our characterized TL/TLR library can be used to engineer RNAs with predictable thermodynamic and conformational behavior.

Engineering Crystal Packing in RNA Structures I: Past and Future Strategies for Engineering RNA Packing in Crystals

X-ray crystallography remains a powerful method to gain atomistic insights into the catalytic and regulatory functions of RNA molecules. However, the technique requires the preparation of diffraction-quality crystals. This is often a resource- and time-consuming venture because RNA crystallization is hindered by the conformational heterogeneity of RNA, as well as the limited opportunities for stereospecific intermolecular interactions between RNA molecules. The limited success at crystallization explains in part the smaller number of RNA-only structures in the Protein Data Bank. Several approaches have been developed to aid the formation of well-ordered RNA crystals. The majority of these are construct-engineering techniques that aim to introduce crystal contacts to favor the formation of well-diffracting crystals. A typical example is the insertion of tetraloop-tetraloop receptor pairs into non-essential RNA segments to promote intermolecular association. Other methods of promoting crystallization involve chaperones and crystallization-friendly molecules that increase RNA stability and improve crystal packing. In this review, we discuss the various techniques that have been successfully used to facilitate crystal packing of RNA molecules, recent advances in construct engineering, and directions for future research in this vital aspect of RNA crystallography.

Annotation of the local context of the RNA secondary structure improves the classification and prediction of A-minors

Non-coding RNAs play a crucial role in various cellular processes in living organisms, and RNA functions heavily depend on molecule structures composed of stems, loops, and various tertiary motifs. Among those, the most frequent are A-minor interactions, which are often involved in the formation of more complex motifs such as kink-turns and pseudoknots. We present a novel classification of A-minors in terms of RNA secondary structure where each nucleotide of an A-minor is attributed to the stem or loop, and each pair of nucleotides is attributed to their relative position within the secondary structure. By analyzing classes of A-minors in known RNA structures, we found that the largest classes are mostly homogeneous and preferably localize with known A-minor co-motifs, e.g. tetraloop-tetraloop receptor and coaxial stacking. Detailed analysis of local A-minors within internal loops revealed a novel recurrent RNA tertiary motif, the across-bulged motif. Interestingly, the motif resembles the previously known GAAA/11nt motif but with the local adenines performing the role of the GAAA-tetraloop. By using machine learning, we show that particular classes of local A-minors can be predicted from sequence and secondary structure. The proposed classification is the first step toward automatic annotation of not only A-minors and their co-motifs but various types of RNA tertiary motifs as well.

In vitro selected GUAA tetraloop-binding receptors with structural plasticity and evolvability towards natural RNA structural modules

GNRA tetraloop-binding receptor interactions are key components in the macromolecular assembly of a variety of functional RNAs. In nature, there is an apparent bias for GAAA/11nt receptor and GYRA/helix interactions, with the former interaction being thermodynamically more stable than the latter. While past in vitro selections allowed isolation of novel GGAA and GUGA receptors, we report herein an in vitro selection that revealed several novel classes of specific GUAA receptors with binding affinities comparable to those from natural GAAA/11nt interactions. These GUAA receptors have structural homology with double-locked bulge RNA modules naturally occurring in ribosomal RNAs. They display mutational robustness that enables exploration of the sequence/phenotypic space associated to GNRA/receptor interactions through epistasis. Their thermodynamic self-assembly fitness landscape is characterized by a rugged neutral network with possible evolutionary trajectories toward natural GNRA/receptor interactions. High throughput sequencing analysis revealed synergetic mutations located away from the tertiary interactions that positively contribute to assembly fitness. Our study suggests that the repertoire of GNRA/receptor interactions is much larger than initially thought from the analysis of natural stable RNA molecules and also provides clues for their evolution towards natural GNRA/receptors.

Structural distinctions between NAD+ riboswitch domains 1 and 2 determine differential folding and ligand binding

Riboswitches are important gene regulatory elements frequently encountered in bacterial mRNAs. The recently discovered nadA riboswitch contains two similar, tandemly arrayed aptamer domains, with the first domain possessing high affinity for nicotinamide adenine dinucleotide (NAD⁺). The second domain which comprises the ribosomal binding site in a putative regulatory helix, however, has withdrawn from detection of ligand-induced structural modulation thus far, and therefore, the identity of the cognate ligand and the regulation mechanism have remained unclear. Here, we report crystal structures of both riboswitch domains, each bound to NAD⁺. Furthermore, we demonstrate that ligand binding to domain 2 requires significantly higher concentrations of NAD⁺ (or ADP retaining analogs) compared to domain 1. Using a fluorescence spectroscopic approach, we further shed light on the structural features which are responsible for the different ligand affinities, and describe the Mg²⁺-dependent, distinct folding and pre-organization of their binding pockets. Finally, we speculate about possible scenarios for nadA RNA gene regulation as a putative two-concentration sensor module for a time-controlled signal that is primed and stalled by the gene regulation machinery at low ligand concentrations (domain 1), and finally triggers repression of translation as soon as high ligand concentrations are reached in the cell (domain 2).

Responsive self-assembly of tectoRNAs with loop-receptor interactions from the tetrahydrofolate (THF) riboswitch

Naturally occurring RNAs are known to exhibit a high degree of modularity, whereby specific structural modules (or motifs) can be mixed and matched to create new molecular architectures. The modular nature of RNA also affords researchers the ability to characterize individual structural elements in controlled synthetic contexts in order to gain new and critical insights into their particular structural features and overall performance. Here, we characterized the binding affinity of a unique loop–receptor interaction found in the tetrahydrofolate (THF) riboswitch using rationally designed self-assembling tectoRNAs. Our work suggests that the THF loop–receptor interaction has been fine-tuned for its particular role as a riboswitch component. We also demonstrate that the thermodynamic stability of this interaction can be modulated by the presence of folinic acid, which induces a local structural change at the level of the loop–receptor. This corroborates the existence of a THF binding site within this tertiary module and paves the way for its potential use as a THF responsive module for RNA nanotechnology and synthetic biology.

Conformational adaptation of UNCG loops upon crowding

If the A-form helix is the major structural motif found in RNA, the loops that cap them constitute the second most important family of motifs. Among those, two are overrepresented, GNRA and UNCG tetraloops. Recent surveys of RNA structures deposited in the PDB show that GNRA and UNCG tetraloops can adopt tertiary folds that are very different from their canonical conformations, characterized by the presence of a U-turn of a Z-turn, respectively. Crystallographic data from both a lariat-capping (LC) ribozyme and a group II intron ribozyme reveal that a given UUCG tetraloop can adopt a distinct fold depending on its structural environment. Specifically, when the crystal packing applies relaxed constraints on the loop, the canonical Z-turn conformation is observed. In contrast, a highly packed environment induces "squashing" of the tetraloop by distorting its sugar-phosphate backbone in a specific way that expels the first and fourth nucleobases out of the loop, and falls in van der Waals distance of the last base pair of the helix, taking the place of the pair formed between the first and fourth residues in Z-turn loops. The biological relevance of our observations is supported by the presence of similarly deformed loops in the highly packed environment of the ribosome and in a complex between a dsRNA and a RNase III. The finding that Z-turn loops change conformation under higher molecular packing suggests that, in addition to their demonstrated role in stabilizing RNA folding, they may contribute to the three-dimensional structure of RNA by mediating tertiary interactions with distal residues.

Linking RNA Sequence, Structure, and Function on Massively Parallel High-Throughput Sequencers

High-throughput sequencing methods have revolutionized our ability to catalog the diversity of RNAs and RNA-protein interactions that can exist in our cells. However, the relationship between RNA sequence, structure, and function is enormously complex, demonstrating the need for methods that can provide quantitative thermodynamic and kinetic measurements of macromolecular interaction with RNA, at a scale commensurate with the sequence diversity of RNA. Here, we discuss a class of methods that extend the core functionality of DNA sequencers to enable high-throughput measurements of RNA folding and RNA-protein interactions. Topics discussed include a description of the method and multiple applications to RNA-binding proteins, riboswitch design and engineering, and RNA tertiary structure energetics.

Rational Design of an Orthogonal Pair of Bimolecular RNase P Ribozymes through Heterologous Assembly of Their Modular Domains

The modular structural domains of multidomain RNA enzymes can often be dissected into separate domain RNAs and their noncovalent assembly can often reconstitute active enzymes. These properties are important to understand their basic characteristics and are useful for their application to RNA-based nanostructures. Bimolecular forms of bacterial RNase P ribozymes consisting of S-domain and C-domain RNAs are attractive as platforms for catalytic RNA nanostructures, but their S-domain/C-domain assembly was not optimized for this purpose. Through analysis and engineering of bimolecular forms of the two bacterial RNase P ribozymes, we constructed a chimeric ribozyme with improved catalytic ability and S-domain/C-domain assembly and developed a pair of bimolecular RNase P ribozymes the assembly of which was considerably orthogonal to each other.

Deducing putative ancestral forms of GNRA/receptor interactions from the ribosome

Stable RNAs rely on a vast repertoire of long-range interactions to assist in the folding of complex cellular machineries such as the ribosome. The universally conserved L39/H89 interaction is a long-range GNRA-like/receptor interaction localized in proximity to the peptidyl transferase center of the large subunit of the ribosome. Because of its central location, L39/H89 likely originated at an early evolutionary stage of the ribosome and played a significant role in its early function. However, L39/H89 self-assembly is impaired outside the ribosomal context. Herein, we demonstrate that structural modularity principles can be used to re-engineer L39/H89 to self-assemble in vitro. The new versions of L39/H89 improve affinity and loop selectivity by several orders of magnitude and retain the structural and functional features of their natural counterparts. These versions of L39/H89 are proposed to be ancestral forms of L39/H89 that were capable of assembling and folding independently from proteins and post-transcriptional modifications. This work demonstrates that novel RNA modules can be rationally designed by taking advantage of the modular syntax of RNA. It offers the prospect of creating new biochemical models of the ancestral ribosome and increases the tool kit for RNA nanotechnology and synthetic biology.

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.

Identification of receptors for UNCG and GNRA Z-turns and their occurrence in rRNA

In contrast to GNRA tetraloop receptors that are common in RNA, receptors for the more thermostable UNCG loops have remained elusive for almost three decades. An analysis of all RNA structures with resolution ≤3.0 Å from the PDB allowed us to identify three previously unnoticed receptors for UNCG and GNRA tetraloops that adopt a common UNCG fold, named ‘Z-turn’ in agreement with our previously published nomenclature. These receptors recognize the solvent accessible second Z-turn nucleotide in different but specific ways. Two receptors participating in a complex network of tertiary interactions are associated with the rRNA UUCG and GAAA Z-turns capping helices H62 and H35a in rRNA large subunits. Structural comparison of fully assembled ribosomes and comparative sequence analysis of >6500 rRNA sequences helped us recognize that these motifs are almost universally conserved in rRNA, where they may contribute to organize the large subunit around the subdomain-IV four-way junction. The third UCCG receptor was identified in a rRNA/protein construct crystallized at acidic pH. These three non-redundant Z-turn receptors are relevant for our understanding of the assembly of rRNA and other long-non-coding RNAs, as well as for the design of novel folding motifs for synthetic biology.

How the Conformations of an Internal Junction Contribute to Fold an RNA Domain

Like proteins, some RNAs fold to compact structures. We can model functional RNAs as a series of short, rigid, base-paired elements, connected by non-base-paired nucleotides that serve as junctions. These connecting regions bend and twist, facilitating the formation of tertiary contacts that stabilize compact states. Here, we explore the roles of salt and junction sequence in determining the structures of a ubiquitous connector: an asymmetric internal loop. We focus on the J5/5a junction from the widely studied P4-P6 domain of the Tetrahymena ribozyme. Following the addition of magnesium ions to fold P4-P6, this junction bends dramatically, bringing the two halves of the RNA domain together for tertiary contact engagement. Using single-molecule fluorescence resonance energy transfer (smFRET), we examine the role of sequence and salt on model RNA constructs that contain these junction regions. We explore the wild-type J5/5a junction as well as two sequence variants. These junctions display distinct, salt-dependent conformations. Small-angle X-ray scattering (SAXS) measurements verify that these effects persist in the full-length P4-P6 domain. These measurements underscore the importance of junction sequence and interactions with ions in facilitating RNA folding.

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.

High-Throughput Investigation of Diverse Junction Elements in RNA Tertiary Folding

RNAs fold into defined tertiary structures to function in critical biological processes. While quantitative models can predict RNA secondary structure stability, we are still unable to predict the thermodynamic stability of RNA tertiary structure. Here, we probe conformational preferences of diverse RNA two-way junctions to develop a predictive model for the formation of RNA tertiary structure. We quantitatively measured tertiary assembly energetics of >1,000 of RNA junctions inserted in multiple structural scaffolds to generate a "thermodynamic fingerprint" for each junction. Thermodynamic fingerprints enabled comparison of junction conformational preferences, revealing principles for how sequence influences 3-dimensional conformations. Utilizing fingerprints of junctions with known crystal structures, we generated ensembles for related junctions that predicted their thermodynamic effects on assembly formation. This work reveals sequence-structure-energetic relationships in RNA, demonstrates the capacity for diverse compensation strategies within tertiary structures, and provides a path to quantitative modeling of RNA folding energetics based on "ensemble modularity."

Blind prediction of noncanonical RNA structure at atomic accuracy

Prediction of RNA structure from nucleotide sequence remains an unsolved grand challenge of biochemistry and requires distinct concepts from protein structure prediction. Despite extensive algorithmic development in recent years, modeling of noncanonical base pairs of new RNA structural motifs has not been achieved in blind challenges. We report a stepwise Monte Carlo (SWM) method with a unique add-and-delete move set that enables predictions of noncanonical base pairs of complex RNA structures. A benchmark of 82 diverse motifs establishes the method's general ability to recover noncanonical pairs ab initio, including multistrand motifs that have been refractory to prior approaches. In a blind challenge, SWM models predicted nucleotide-resolution chemical mapping and compensatory mutagenesis experiments for three in vitro selected tetraloop/receptors with previously unsolved structures (C7.2, C7.10, and R1). As a final test, SWM blindly and correctly predicted all noncanonical pairs of a Zika virus double pseudoknot during a recent community-wide RNA-Puzzle. Stepwise structure formation, as encoded in the SWM method, enables modeling of noncanonical RNA structure in a variety of previously intractable problems.

Single-Molecule Fluorescence Reveals Commonalities and Distinctions among Natural and in Vitro-Selected RNA Tertiary Motifs in a Multistep Folding Pathway

Decades of study of the RNA folding problem have revealed that diverse and complex structured RNAs are built from a common set of recurring structural motifs, leading to the perspective that a generalizable model of RNA folding may be developed from understanding of the folding properties of individual structural motifs. We used single-molecule fluorescence to dissect the kinetic and thermodynamic properties of a set of variants of a common tertiary structural motif, the tetraloop/tetraloop-receptor (TL/TLR). Our results revealed a multistep TL/TLR folding pathway in which preorganization of the ubiquitous AA-platform submotif precedes the formation of the docking transition state and tertiary A-minor hydrogen bond interactions form after the docking transition state. Differences in ion dependences between TL/TLR variants indicated the occurrence of sequence-dependent conformational rearrangements prior to and after the formation of the docking transition state. Nevertheless, varying the junction connecting the TL/TLR produced a common kinetic and ionic effect for all variants, suggesting that the global conformational search and compaction electrostatics are energetically independent from the formation of the tertiary motif contacts. We also found that in vitro-selected variants, despite their similar stability at high Mg²⁺ concentrations, are considerably less stable than natural variants under near-physiological ionic conditions, and the occurrence of the TL/TLR sequence variants in Nature correlates with their thermodynamic stability in isolation. Overall, our findings are consistent with modular but complex energetic properties of RNA structural motifs and will aid in the eventual quantitative description of RNA folding from its secondary and tertiary structural elements.

Artificial RNA Motifs Expand the Programmable Assembly between RNA Modules of a Bimolecular Ribozyme Leading to Application to RNA Nanostructure Design

A bimolecular ribozyme consisting of a core ribozyme (ΔP5 RNA) and an activator module (P5abc RNA) has been used as a platform to design assembled RNA nanostructures. The tight and specific assembly between the P5abc and ΔP5 modules depends on two sets of intermodule interactions. The interface between P5abc and ΔP5 must be controlled when designing RNA nanostructures. To expand the repertoire of molecular recognition in the P5abc/ΔP5 interface, we modified the interface by replacing the parent tertiary interactions in the interface with artificial interactions. The engineered P5abc/ΔP5 interfaces were characterized biochemically to identify those suitable for nanostructure design. The new interfaces were used to construct 2D-square and 1D-array RNA nanostructures.

Structural Principles of Fluorescent RNA Aptamers

Several aptamer RNAs have been selected in vitro that bind to otherwise weakly fluorescent small molecules and enhance their fluorescence several thousand-fold. By genetically tagging cellular RNAs of interest with these aptamers and soaking cells in their cell-permeable cognate small-molecule fluorophores, it is possible to use them to study RNA localization and trafficking. These aptamers have also been fused to metabolite-binding RNAs to generate fluorescent biosensors. The 3D structures of three unrelated fluorogenic RNAs have been determined, and reveal a shared reliance on base quadruples (tetrads) to constrain the photo-excited chromophore. The structural diversity of fluorogenic RNAs and the chemical diversity of potential fluorophores to be activated are likely to yield a variety of future fluorogenic RNA tags that are optimized for different applications in RNA imaging and in the design of fluorescent RNA biosensors.

Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway

The past decade has seen a wealth of 3D structural information about complex structured RNAs and identification of functional intermediates. Nevertheless, developing a complete and predictive understanding of the folding and function of these RNAs in biology will require connection of individual rate and equilibrium constants to structural changes that occur in individual folding steps and further relating these steps to the properties and behavior of isolated, simplified systems. To accomplish these goals we used the considerable structural knowledge of the folded, unfolded, and intermediate states of P4-P6 RNA. We enumerated structural states and possible folding transitions and determined rate and equilibrium constants for the transitions between these states using single-molecule FRET with a series of mutant P4-P6 variants. Comparisons with simplified constructs containing an isolated tertiary contact suggest that a given tertiary interaction has a stereotyped rate for breaking that may help identify structural transitions within complex RNAs and simplify the prediction of folding kinetics and thermodynamics for structured RNAs from their parts. The preferred folding pathway involves initial formation of the proximal tertiary contact. However, this preference was only ∼10 fold and could be reversed by a single point mutation, indicating that a model akin to a protein-folding contact order model will not suffice to describe RNA folding. Instead, our results suggest a strong analogy with a modified RNA diffusion-collision model in which tertiary elements within preformed secondary structures collide, with the success of these collisions dependent on whether the tertiary elements are in their rare binding-competent conformations.

Mutational analysis of structural elements in a class-I cyclic di-GMP riboswitch to elucidate its regulatory mechanism

The Vc2 riboswitch possesses an aptamer domain belonging to the class-I c-di-GMP riboswitch family. This domain has been analysed and the molecular mechanism by which it recognizes the c-di-GMP ligand has been elucidated. On the other hand, the regulatory mechanism of the full-length Vc2 riboswitch to control its downstream open reading frame (ORF) remains largely unknown. In this study, we performed in vivo reporter assays and in vitro biochemical analyses of the full-length riboswitch and its aptamer domain. We evaluated the results of in vivo and in vitro analyses to elucidate the regulatory mechanism of the Vc2 riboswitch. The present results suggest that recognition of c-di-GMP ligand by the Vc2 riboswitch aptamer domain downregulates expression of its downstream ORF primarily at the translational level.

Characterization of an RNA receptor motif that recognizes a GCGA tetraloop

Tertiary interactions between a new RNA motif and RNA tetraloops were analyzed to determine whether this new motif shows preference for a GCGA tetraloop. In the structural context of a ligase ribozyme, this motif discriminated GCGA loop from 3 other tetraloops. The affinity between the GCGA loop and its receptor is strong enough to carry out the ribozyme activity.

Optimization of RNA-based c-di-GMP fluorescent sensors through tuning their structural modules

Cyclic diguanylate (c-di-GMP) is a second messenger of bacteria and its detection is an important issue in basic and applied microbiology. As c-di-GMP riboswitch ligand-binding domains (aptamer domains) capture c-di-GMP with high affinity and selectivity, they are promising platforms for the development of RNA-based c-di-GMP sensors. We analyzed two previously reported c-di-GMP sensor RNAs derived from the Vc2 riboswitch. We also designed and tested their variants, some of which showed improved properties as RNA-based c-di-GMP sensors.

RNA-Redesign: a web server for fixed-backbone 3D design of RNA

RNA is rising in importance as a design medium for interrogating fundamental biology and for developing therapeutic and bioengineering applications. While there are several online servers for design of RNA secondary structure, there are no tools available for the rational design of 3D RNA structure. Here we present RNA-Redesign (http://rnaredesign.stanford.edu), an online 3D design tool for RNA. This resource utilizes fixed-backbone design to optimize the sequence identity and nucleobase conformations of an RNA to match a desired backbone, analogous to fundamental tools that underlie rational protein engineering. The resulting sequences suggest thermostabilizing mutations that can be experimentally verified. Further, sequence preferences that differ between natural and computationally designed sequences can suggest whether natural sequences possess functional constraints besides folding stability, such as cofactor binding or conformational switching. Finally, for biochemical studies, the designed sequences can suggest experimental tests of 3D models, including concomitant mutation of base triples. In addition to the designs generated, detailed graphical analysis is presented through an integrated and user-friendly environment.

Triggering RNAi with multifunctional RNA nanoparticles and their delivery

Proteins are considered to be the key players in structure, function, and metabolic regulation of our bodies. The mechanisms used in conventional therapies often rely on inhibition of proteins with small molecules, but another promising method to treat disease is by targeting the corresponding mRNAs. In 1998, Craig Mellow and Andrew Fire discovered dsRNA-mediated gene silencing via RNA interference or RNAi. This discovery introduced almost unlimited possibilities for new gene silencing methods, thus opening new doors to clinical medicine. RNAi is a biological process that inhibits gene expression by targeting the mRNA. RNAi-based therapeutics have several potential advantages (i) a priori ability to target any gene, (ii) relatively simple design process, (iii) site-specificity, (iv) potency, and (v) a potentially safe and selective knockdown of the targeted cells. However, the problem lies within the formulation and delivery of RNAi therapeutics including rapid excretion, instability in the bloodstream, poor cellular uptake, and inefficient intracellular release. In an attempt to solve these issues, different types of RNAi therapeutic delivery strategies including multifunctional RNA nanoparticles are being developed. In this mini-review, we will briefly describe some of the current approaches.

Engineered RNA Nanodesigns for Applications in RNA Nanotechnology

Nucleic acids have emerged as an extremely promising platform for nanotechnological applications because of their unique biochemical properties and functions. RNA, in particular, is characterized by relatively high thermal stability, diverse structural flexibility, and its capacity to perform a variety of functions in nature. These properties make RNA a valuable platform for bio-nanotechnology, specifically RNA Nanotechnology, that can create de novo nanostructures with unique functionalities through the design, integration, and re-engineering of powerful mechanisms based on a variety of existing RNA structures and their fundamental biochemical properties. This review highlights the principles that underlie the rational design of RNA nanostructures, describes the main strategies used to construct self-assembling nanoparticles, and discusses the challenges and possibilities facing the application of RNA Nanotechnology in the future.

Role of RNA structure motifs in IRES-dependent translation initiation of the coxsackievirus B3: new insights for developing live-attenuated strains for vaccines and gene therapy

Internal ribosome entry site (IRES) elements are highly structured RNA sequences that function to recruit ribosomes for the initiation of translation. In contrast to the canonical cap-binding, the mechanism of IRES-mediated translation initiation is still poorly understood. Translation initiation of the coxsackievirus B3 (CVB3), a causative agent of viral myocarditis, has been shown to be mediated by a highly ordered structure of the 5' untranslated region (5'UTR), which harbors an IRES. Taking into account that efficient initiation of mRNA translation depends on temporally and spatially orchestrated sequence of RNA-protein and RNA-RNA interactions, and that, at present, little is known about these interactions, we aimed to describe recent advances in our understanding of molecular structures and biochemical functions of the translation initiation process. Thus, this review will explore the IRES elements as important RNA structures and the significance of these structures in providing an alternative mechanism of translation initiation of the CVB3 RNA. Since translation initiation is the first intracellular step during the CVB3 infection cycle, the IRES region provides an ideal target for antiviral therapies. Interestingly, the 5' and 3'UTRs represent promising candidates for the study of CVB3 cardiovirulence and provide new insights for developing live-attenuated vaccines.

RNA self-assembly and RNA nanotechnology

CONSPECTUS: Nanotechnology's central goal involves the direct control of matter at the molecular nanometer scale to build nanofactories, nanomachines, and other devices for potential applications including electronics, alternative fuels, and medicine. In this regard, the nascent use of nucleic acids as a material to coordinate the precise arrangements of specific molecules marked an important milestone in the relatively recent history of nanotechnology. While DNA served as the pioneer building material in nucleic acid nanotechnology, RNA continues to emerge as viable alternative material with its own distinct advantages for nanoconstruction. Several complementary assembly strategies have been used to build a diverse set of RNA nanostructures having unique structural attributes and the ability to self-assemble in a highly programmable and controlled manner. Of the different strategies, the architectonics approach uniquely endeavors to understand integrated structural RNA architectures through the arrangement of their characteristic structural building blocks. Viewed through this lens, it becomes apparent that nature routinely uses thermodynamically stable, recurrent modular motifs from natural RNA molecules to generate unique and more complex programmable structures. With the design principles found in natural structures, a number of synthetic RNAs have been constructed. The synthetic nanostructures constructed to date have provided, in addition to affording essential insights into RNA design, important platforms to characterize and validate the structural self-folding and assembly properties of RNA modules or building blocks. Furthermore, RNA nanoparticles have shown great promise for applications in nanomedicine and RNA-based therapeutics. Nevertheless, the synthetic RNA architectures achieved thus far consist largely of static, rigid particles that are still far from matching the structural and functional complexity of natural responsive structural elements such as the ribosome, large ribozymes, and riboswitches. Thus, the next step in synthetic RNA design will involve new ways to implement these same types of dynamic and responsive architectures into nanostructures functioning as real nanomachines in and outside the cell. RNA nanotechnology will likely garner broader utility and influence with a greater focus on the interplay between thermodynamic and kinetic influences on RNA self-assembly and using natural RNAs as guiding principles.

The art of editing RNA structural alignments

Manual editing of RNA structural alignments may be considered more art than science, since it still requires an expert biologist to take multiple levels of information into account and be slightly creative when constructing high-quality alignments. Even though the task is rather tedious, it is rewarded by great insight into the evolution of structure and function of your favorite RNA molecule. In this chapter I will review the methods and considerations that go into constructing RNA structural alignments at the secondary and tertiary structure level; introduce software, databases, and algorithms that have proven useful in semiautomating the work process; and suggest future directions towards full automatization.

RNA modularity for synthetic biology

RNA molecules are highly modular components that can be used in a variety of contexts for building new metabolic, regulatory and genetic circuits in cells. The majority of synthetic RNA systems to date predominately rely on two-dimensional modularity. However, a better understanding and integration of three-dimensional RNA modularity at structural and functional levels is critical to the development of more complex, functional bio-systems and molecular machines for synthetic biology applications.

Comparative sequence and structure analysis reveals the conservation and diversity of nucleotide positions and their associated tertiary interactions in the riboswitches

The tertiary motifs in complex RNA molecules play vital roles to either stabilize the formation of RNA 3D structure or to provide important biological functionality to the molecule. In order to better understand the roles of these tertiary motifs in riboswitches, we examined 11 representative riboswitch PDB structures for potential agreement of both motif occurrences and conservations. A total of 61 unique tertiary interactions were found in the reference structures. In addition to the expected common A-minor motifs and base-triples mainly involved in linking distant regions the riboswitch structures three highly conserved variants of A-minor interactions called G-minors were found in the SAM-I and FMN riboswitches where they appear to be involved in the recognition of the respective ligand’s functional groups. From our structural survey as well as corresponding structure and sequence alignments, the agreement between motif occurrences and conservations are very prominent across the representative riboswitches. Our analysis provide evidence that some of these tertiary interactions are essential components to form the structure where their sequence positions are conserved despite a high degree of diversity in other parts of the respective riboswitches sequences. This is indicative of a vital role for these tertiary interactions in determining the specific biological function of riboswitch.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

An in vitro-selected RNA receptor for the GAAC loop: modular receptor for non-GNRA-type tetraloop

Although artificial RNA motifs that can functionally replace the GNRA/receptor interaction, a class of RNA-RNA interacting motifs, were isolated from RNA libraries and used to generate designer RNA structures, receptors for non-GNRA tetraloops have not been found in nature or selected from RNA libraries. In this study, we report successful isolation of a receptor motif interacting with GAAC, a non-GNRA tetraloop, from randomized sequences embedded in a catalytic RNA. Biochemical characterization of the GAAC/receptor interacting motif within three structural contexts showed its binding affinity, selectivity and structural autonomy. The motif has binding affinity comparable with that of a GNRA/receptor, selectivity orthogonal to GNRA/receptors and structural autonomy even in a large RNA context. These features would be advantageous for usage of the motif as a building block for designer RNAs. The isolated motif can also be used as a query sequence to search for unidentified naturally occurring GANC receptor motifs.

The GA-minor submotif as a case study of RNA modularity, prediction, and design

Complex natural RNAs such as the ribosome, group I and group II introns, and RNase P exemplify the fact that three-dimensional (3D) RNA structures are highly modular and hierarchical in nature. Tertiary RNA folding typically takes advantage of a rather limited set of recurrent structural motifs that are responsible for controlling bends or stacks between adjacent helices. Herein, the GA minor and related structural motifs are presented as a case study to highlight several structural and folding principles, to gain further insight into the structural evolution of naturally occurring RNAs, as well as to assist the rational design of artificial RNAs.

Fixation and accumulation of thermotolerant catalytic competence of a pair of ligase ribozymes through complex formation and cross ligation

In the early stages of the hypothetical RNA world, some primitive RNA catalysts (ribozymes) may have emerged through self-assembly of short RNA oligomers. Although they may be unstable against temperature fluctuations and other environmental changes, ligase ribozymes (ribozymes with RNA strand-joining activity) may resolve structural instability of self-assembling RNAs by converting them to the corresponding unimolecular formats. To investigate this possibility, we constructed a model system using a cross-ligation system composed of a pair of self-assembling ligase ribozymes. Their abilities to act as catalysts, substrates, and a cross-ligation system were analyzed with or without thermal pretreatment before the reactions. A pair of self-assembling ligase ribozymes, each of which can form multiple conformations, demonstrated that thermotolerance was acquired and accumulated through complex-formation that stabilized the active forms of the bimolecular ribozymes and also cross-ligation that produced the unimolecular ribozymes.

Sequence-dependent structural dynamics of primate adenosine-to-inosine editing substrates

Humans have the highest level of adenosine-to-inosine (A-to-I) editing amongst primates, yet the reasons for this difference remain unclear. Sequence analysis of the Alu Sg elements (A-to-I RNA substrates) corresponding to the Nup50 gene in human, chimp, and rhesus reveals subtle sequence variations surrounding the edit sites. We have developed three constructs that represent human (HuAp5), chimp (ChAp5), and rhesus (RhAp5) Nup50 Alu Sg A-to-I editing substrates. Here, 2-aminopurine (2-Ap) was substituted for edited adenosine (A5) so as to monitor the fluorescence intensity with respect to temperature. UV and steady-state fluorescence (SSF) T(M) plots indicate that local and global unfolding are coincident, with the human construct displaying a T(M) of approximately 70°C, compared to 60°C for chimp and 54°C for rhesus. However, time-resolved fluorescence (TRF) resolves three different fluorescence lifetimes that we assign to folded, intermediate(s), and unfolded states. The TRF data fit well to a two-intermediate model, whereby both intermediates (M, J) are in equilibrium with each other, and the folded/unfolded states. Our model suggests that, at 37°C, human state J and the folded state will be the most heavily populated in comparison to the other primate constructs. In order for adenosine deaminase acting on RNA (ADAR) to efficiently dock, a stable duplex must be present that corresponds to the human construct, globally. Next, the enzyme must "flip out" the base of interest to facilitate the A-to-I conversion; a nucleotide in an intermediate-like position would enhance this conformational change. Our experiments demonstrate that subtle variations in RNA sequence might contribute to the high A-to-I editing levels found in humans.

Structural variation and uniformity among tetraloop-receptor interactions and other loop-helix interactions in RNA crystal structures

Tetraloop-receptor interactions are prevalent structural units in RNAs, and include the GAAA/11-nt and GNRA-minor groove interactions. In this study, we have compiled a set of 78 nonredundant loop-helix interactions from X-ray crystal structures, and examined them for the extent of their sequence and structural variation. Of the 78 interactions in the set, only four were classical GAAA/11-nt motifs, while over half (48) were GNRA-minor groove interactions. The GNRA-minor groove interactions were not a homogeneous set, but were divided into five subclasses. The most predominant subclass is characterized by two triple base pair interactions in the minor groove, flanked by two ribose zipper contacts. This geometry may be considered the "standard" GNRA-minor groove interaction, while the other four subclasses are alternative ways to form interfaces between a minor groove and tetraloop. The remaining 26 structures in the set of 78 have loops interacting with mostly idiosyncratic receptors. Among the entire set, a number of sequence-structure correlations can be identified, which may be used as initial hypotheses in predicting three-dimensional structures from primary sequences. Conversely, other sequence patterns are not predictive; for example, GAAA loop sequences and GG/CC receptors bind to each other with three distinct geometries. Finally, we observe an example of structural evolution in group II introns, in which loop-receptor motifs are substituted for each other while maintaining the larger three-dimensional geometry. Overall, the study gives a more complete view of RNA loop-helix interactions that exist in nature.

The right angle (RA) motif: a prevalent ribosomal RNA structural pattern found in group I introns

The right angle (RA) motif, previously identified in the ribosome and used as a structural module for nano-construction, is a recurrent structural motif of 13 nucleotides that establishes a 90° bend between two adjacent helices. Comparative sequence analysis was used to explore the sequence space of the RA motif within ribosomal RNAs in order to define its canonical sequence space signature. We investigated the sequence constraints associated with the RA signature using several artificial self-assembly systems. Thermodynamic and topological investigations of sequence variants associated with the RA motif in both minimal and expanded structural contexts reveal that the presence of a helix at the 3' end of the RA motif increases the thermodynamic stability and rigidity of the resulting three-helix junction domain. A search for the RA in naturally occurring RNAs as well as its experimental characterization led to the identification of the RA in groups IC1 and ID intron ribozymes, where it is suggested to play an integral role in stabilizing peripheral structural domains. The present study exemplifies the need of empirical analysis of RNA structural motifs for facilitating the rational design and structure prediction of RNAs.