Involvement of a GNRA tetraloop in long-range RNA tertiary interactions

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.

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.

Molecular Anatomy of the Class I Ligase Ribozyme for Elucidation of the Activity-Generating Unit

The class I ligase ribozyme consists of 121 nucleotides and shows a high catalytic rate comparable to that found in natural proteinaceous polymerases. In this study, we aimed to identify the smaller active unit of the class I ligase ribozyme comprising ~50 nucleotides, comparable to the estimated length of prebiotically synthesized RNA. Based on the three-dimensional structure of the class I ligase ribozyme, mutants were prepared and their ligation activities were analyzed. Sufficient ligation activity was maintained even when shortening to 94 nucleotides. However, because it would be difficult to approach the target of ~50 nucleotides by removing only the partial structure, the class I ligase ribozyme was then split into two molecules. The ligation activity was maintained even when splitting into two molecules of 55 and 39 nucleotides. Using a system with similar split ribozymes, we analyzed the ligation activity of mutants C30, C47, and A71, which have been previously identified as the positions that contribute to catalytic activity, and discussed the structural basis of the activity of these bases. Our findings suggest the rationale for the class I ligase ribozyme's assembling from multiple fragments that would be achievable with prebiotic synthesis.

A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States

An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. We isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.

Cooperativity and Interdependency between RNA Structure and RNA-RNA Interactions

Complex RNA–RNA interactions are increasingly known to play key roles in numerous biological processes from gene expression control to ribonucleoprotein granule formation. By contrast, the nature of these interactions and characteristics of their interfaces, especially those that involve partially or wholly structured RNAs, remain elusive. Herein, we discuss different modalities of RNA–RNA interactions with an emphasis on those that depend on secondary, tertiary, or quaternary structure. We dissect recently structurally elucidated RNA–RNA complexes including RNA triplexes, riboswitches, ribozymes, and reverse transcription complexes. These analyses highlight a reciprocal relationship that intimately links RNA structure formation with RNA–RNA interactions. The interactions not only shape and sculpt RNA structures but also are enabled and modulated by the structures they create. Understanding this two-way relationship between RNA structure and interactions provides mechanistic insights into the expanding repertoire of noncoding RNA functions, and may inform the design of novel therapeutics that target RNA structures or interactions.

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.

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.

Ageritin from Pioppino Mushroom: The Prototype of Ribotoxin-Like Proteins, a Novel Family of Specific Ribonucleases in Edible Mushrooms

Ageritin is a specific ribonuclease, extracted from the edible mushroom Cyclocybe aegerita (synonym Agrocybe aegerita), which cleaves a single phosphodiester bond located within the universally conserved alpha-sarcin loop (SRL) of 23-28S rRNAs. This cleavage leads to the inhibition of protein biosynthesis, followed by cellular death through apoptosis. The structural and enzymatic properties show that Ageritin is the prototype of a novel specific ribonucleases family named 'ribotoxin-like proteins', recently found in fruiting bodies of other edible basidiomycetes mushrooms (e.g., Ostreatin from Pleurotus ostreatus, Edulitins from Boletus edulis, and Gambositin from Calocybe gambosa). Although the putative role of this toxin, present in high amount in fruiting body (>2.5 mg per 100 g) of C. aegerita, is unknown, its antifungal and insecticidal actions strongly support a role in defense mechanisms. Thus, in this review, we focus on structural, biological, antipathogenic, and enzymatic characteristics of this ribotoxin-like protein. We also highlight its biological relevance and potential biotechnological applications in agriculture as a bio-pesticide and in biomedicine as a therapeutic and diagnostic agent.

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.

Further Characterization of the Pseudo-Symmetrical Ribosomal Region

The peptidyl transferase center of the modern ribosome has been found to encompass an area of twofold pseudosymmetry (SymR). This observation strongly suggests that the very core of the ribosome arose from a dimerization event between two modest-sized RNAs. It was previously shown that at least four non-standard interactions exist between the two halves of SymR. Herein, we verify that the structure of the SymR is highly conserved with respect to both ribosome transition state and phylogenetic diversity. These comparisons also reveal two additional sites of interaction between the two halves of SymR and refine our understanding of the previously known interactions. In addition, the possible role that magnesium may have in the coordination, stabilization, association, and evolutionary history of the two halves (A-region and P-region) was examined. Together, the results identify a likely site where structural elements and Mg²⁺ ions may have facilitated the ligation of two aboriginal RNAs into a single unit.

Synthetic Antibody Binding to a Preorganized RNA Domain of Hepatitis C Virus Internal Ribosome Entry Site Inhibits Translation

Structured RNA elements within the internal ribosome entry site (IRES) of hepatitis C virus (HCV) genome hijack host cell machinery for translation initiation through a cap-independent mechanism. Here, using a phage display selection, we obtained two antibody fragments (Fabs), HCV2 and HCV3, against HCV IRES that bind the RNA with dissociation constants of 32 ± 7 nM and 37 ± 8 nM respectively, specifically recognizing the so-called junction IIIabc (JIIIabc). We used these Fabs as crystallization chaperones and determined the high-resolution crystal structures of JIIIabc-HCV2 and -HCV3 complexes at 1.81 Å and 2.75 Å resolution respectively, revealing an antiparallel four-way junction with the IIIa and IIIc subdomains brought together through tertiary interactions. The RNA conformation observed in the structures supports the structural model for this region derived from cryo-EM data for the HCV IRES-40S ribosome complex, suggesting that the tertiary fold of the RNA preorganizes the domain for interactions with the 40S ribosome. Strikingly, both Fabs and the ribosomal protein eS27 not only interact with a common subset of nucleotides within the JIIIabc but also use physiochemically similar sets of protein residues to do so, suggesting that the RNA surface is well-suited for interactions with proteins, perhaps analogous to the "hot spot" concept elaborated for protein-protein interactions. Using a rabbit reticulocyte lysate-based translation assay with a bicistronic reporter construct, we further demonstrated that Fabs HCV2 and HCV3 specifically inhibit the HCV IRES-directed translation, implicating disruption of the JIIIabc-ribosome interaction as a potential therapeutic strategy against HCV.

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.

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.

Ageritin, a Ribotoxin from Poplar Mushroom ( Agrocybe aegerita) with Defensive and Antiproliferative Activities

Ribotoxins make up a group of extracellular rRNA endoribonucleases produced by ascomycetes that display cytotoxicity toward animal cells, having been proposed as insecticidal agents. Recently, the ribotoxin Ageritin has been isolated from the basidiomycetes Agrocybe aegerita (poplar mushroom), suggesting that ribotoxins are widely distributed among fungi. To gain insights into the protective properties of Ageritin against pathogens and its putative biotechnological applications, we have tested several biological activities of Ageritin, comparing them with those of the well-known ribotoxin α-sarcin, and we found that Ageritin displayed, in addition to the already reported activities, (i) antibacterial activity against Micrococcus lysodeikticus, (ii) activity against the tobacco mosaic virus RNA, (iii) endonuclease activity against a supercoiled plasmid, (iv) nuclease activity against genomic DNA, (v) cytotoxicity to COLO 320, HeLa, and Raji cells by promoting apoptosis, and (vi) antifungal activity against the green mold Penicillium digitatum. Therefore, Ageritin and α-sarcin can induce resistance not only to insects but also to viruses, bacteria, and fungi. The multiple biological activities of Ageritin could be exploited to improve resistance to different pathogens by engineering transgenic plants. Furthermore, the induction of cell death by different mechanisms turns these ribotoxins into useful tools for cancer therapy.

Using Rosetta for RNA homology modeling

The three-dimensional structures of RNA molecules provide rich and often critical information for understanding their functions, including how they recognize small molecule and protein partners. Computational modeling of RNA 3D structure is becoming increasingly accurate, particularly with the availability of growing numbers of template structures already solved experimentally and the development of sequence alignment and 3D modeling tools to take advantage of this database. For several recent "RNA puzzle" blind modeling challenges, we have successfully identified useful template structures and achieved accurate structure predictions through homology modeling tools developed in the Rosetta software suite. We describe our semi-automated methodology here and walk through two illustrative examples: an adenine riboswitch aptamer, modeled from a template guanine riboswitch structure, and a SAM I/IV riboswitch aptamer, modeled from a template SAM I riboswitch structure.

Increased Flexibility between Stems of Intramolecular Three-Way Junctions by the Insertion of Bulges

Intramolecular junctions are a ubiquitous structure within DNA and RNA; three-way junctions in particular have high strain around the junction because of the lack of flexibility, preventing the junctions from adopting conformations that would allow for optimal folding. In this work, we used a combination of calorimetric and spectroscopic techniques to study the unfolding of four intramolecular three-way junctions. The control three-way junction, 3H, has the sequence d(GAAATTGCGCT5GCGCGTGCT5GCACAATTTC), which has three arms of different sequences. We studied three other three-way junctions in which one (2HS1H), two (HS12HS1), and three (HS1HS1HS1) cytosine bulges were placed at the junction to allow the arms to adopt a wider range of conformations that may potentially relieve strain. Through calorimetric studies, it was concluded that bulges produce only minor effects on the enthalpic and thermal stability at physiological salt concentrations for 2HS1H and HS1HS1HS1. HS12HS1 displays the strongest effect, with the GTGC stem lacking a defined transition. In addition to unfolding thermodynamics, the differential binding of counterions, water, and protons was determined. It was found that with each bulge, there was a large increase in the binding of counterions; this correlated with a decrease in the immobilization of structural water molecules. The increase in counterion uptake upon folding likely displaces binding of structural water, which is measured by the osmotic stress method, in favor of electrostricted waters. The cytosine bulges do not affect the binding of protons; this finding indicates that the bulges are not forming base-triplet stacks. These results indicate that bulges in junctions do not affect the unfolding profile or the enthalpy of oligonucleotides but do affect the number and amount of molecules immobilized by the junction.

Effect of GCAA stabilizing loops on three- and four-way intramolecular junctions

Tetraloops are a common way of changing the melting behavior of a DNA or RNA structure without changing the sequence of the stem. Because of the ubiquitous nature of tetraloops, our goal is to understand the effect a GCAA tetraloop, which belongs to the GNRA family of tetraloops, has on the unfolding thermodynamics of intramolecular junctions. Specifically, we have described the melting behavior of intramolecular three-way and four-way junctions where a T₅ loop has been replaced with a GCAA tetraloops in different positions. Their thermodynamic profiles, including Δn_Na⁺ and Δn_W, were analyzed based on the position of the tetraloop. We obtained between -16.7 and -27.5 kcal mol^-1 for all junctions studied. The experimental data indicates the influence of the GCAA tetraloop is primarily dictated by the native unfolding of the junction; if the tetraloop is placed on a stem that unfolds as a single domain when the tetraloop is not present, it will unfold as a single domain when the tetraloop is present but with a higher thermal stability. Conversely, if the tetraloop is placed on a stem which unfolds cooperatively with other stems when the tetraloop is not present, the tetraloop will increase the thermal stability of all the stems in the melting domain. The oligonucleotide structure and not the tetraloop itself affects ion uptake; three-way junctions do not gain an increase in ion uptake, but four-way junctions do. This is not the case for water immobilization, where the position of the tetraloop dictates the amount of water immobilized.

Corroboration of Zn(ii)–Mg(ii)-tertiary structure interplays essential for the optimal catalysis of a phosphorothiolate thiolesterase ribozyme

The TW17 ribozyme, a catalytic RNA selected from a pool of artificial RNA, is specific for the Zn²⁺-dependent hydrolysis of a phosphorothiolate thiolester bond. Here, we describe the organic synthesis of both guanosine α-thio-monophosphate and the substrates required for selecting and characterizing the TW17 ribozyme, and for deciphering the catalytic mechanism of the ribozyme. By successively substituting the substrate originally conjugated to the RNA pool with structurally modified substrates, we demonstrated that the TW17 ribozyme specifically catalyzes phosphorothiolate thiolester hydrolysis. Metal titration studies of TW17 ribozyme catalysis in the presence of Zn²⁺ alone, Zn²⁺ and Mg²⁺, and Zn²⁺ and [Co(NH₃)₆]³⁺ supported our findings that Zn²⁺ is absolutely required for ribozyme catalysis, and indicated that optimal ribozyme catalysis involves the presence of outer-sphere and one inner-sphere Mg²⁺. A survey of the TW17 ribozyme activity at various pHs revealed that the activity of the ribozyme critically depends on the alkaline conditions. Moreover, a GNRA tetraloop-containing ribozyme constructed with active catalysis in trans provided catalysis and multiple substrate turnover efficiencies significantly higher than ribozymes lacking a GNRA tetraloop. This research supports the essential roles of Zn²⁺, Mg²⁺, and a GNRA tetraloop in modulating the TW17 ribozyme structure for optimal ribozyme catalysis, leading also to the formulation of a proposed reaction mechanism for TW17 ribozyme catalysis.

Zn(ii) and Mg(ii) and GAGA tetraloop in the ion atmosphere of the TW17 ribozyme is critical to optimal ribozyme catalysis at alkaline pH.

Effect of Loop Length and Sequence on the Stability of DNA Pyrimidine Triplexes with TAT Base Triplets

We report the thermodynamic contributions of loop length and loop sequence to the overall stability of DNA intramolecular pyrimidine triplexes. Two sets of triplexes were designed: in the first set, the C₅ loop closing the triplex stem was replaced with ^5'-CT_nC loops (n = 1-5), whereas in the second set, both the duplex and triplex loops were replaced with a ^5'-GCAA or ^5'-AACG tetraloop. For the triplexes with a ^5'-CT_nC loop, the triplex with five bases in the loop has the highest stability relative to the control. A loop length lower than five compromises the strength of the base-pair stacks without decreasing the thermal stability, leading to a decreased enthalpy, whereas an increase in the loop length leads to a decreased enthalpy and a higher entropic penalty. The incorporation of the GCAA loop yielded more stable triplexes, whereas the incorporation of AACG in the triplex loop yielded a less stable triplex due to an unfavorable enthalpy term. Thus, addition of the GCAA tetraloop can cause an increase in the thermodynamics of the triplex without affecting the sequence or melting behavior and may result in an additional layer of genetic regulation.

Group II Intron Self-Splicing

Group II introns are large, autocatalytic ribozymes that catalyze RNA splicing and retrotransposition. Splicing by group II introns plays a major role in the metabolism of plants, fungi, and yeast and contributes to genetic variation in many bacteria. Group II introns have played a major role in genome evolution, as they are likely progenitors of spliceosomal introns, retroelements, and other machinery that controls genetic variation and stability. The structure and catalytic mechanism of group II introns have recently been elucidated through a combination of genetics, chemical biology, solution biochemistry, and crystallography. These studies reveal a dynamic machine that cycles progressively through multiple conformations as it stimulates the various stages of splicing. A central active site, containing a reactive metal ion cluster, catalyzes both steps of self-splicing. These studies provide insights into RNA structure, folding, and catalysis, as they raise new questions about the behavior of RNA machines.

Study of RNA-A Initiation Translation of The Infectious Pancreatic Necrosis Virus

The infectious pancreatic necrosis virus (IPNV) is a salmonid pathogen that causes significant economic losses to the aquaculture industry. IPNV is a non-enveloped virus containing two uncapped and non-polyadenylated double strand RNA genomic segments, RNA-A and RNA-B. The viral protein Vpg is covalently attached to the 5' end of both segments. There is little knowledge about its viral cycle, particularly about the translation of the RNAs. Through experiments using mono and bicistronic reporters, in this work we show that the 120-nucleotide-long 5'-UTR of RNA-A contains an internal ribosome entry site (IRES) that functions efficiently both in vitro and in salmon cells. IRES activity is strongly dependent on temperature. Also, the IRES structure is confined to the 5'UTR and is not affected by the viral coding sequence. This is the first report of IRES activity in a fish virus and can give us tools to generate antivirals to attack the virus without affecting fish directly.

Structural basis for high-affinity fluorophore binding and activation by RNA Mango

Genetically encoded fluorescent protein tags revolutionized proteome studies, while the lack of intrinsically fluorescent RNAs has hindered transcriptome exploration. Among several RNA-fluorophore complexes that potentially address this problem, RNA Mango has an exceptionally high affinity for its thiazole orange (TO)-derived fluorophore, TO1-Biotin (K_d ~3 nM), and in complex with related ligands, is one of the most red-shifted fluorescent macromolecular tags known. To elucidate how this small aptamer exhibits such properties, which make it well suited for studying low-copy cellular RNAs, we determined its 1.7 Å resolution co-crystal structure. Unexpectedly, the entire ligand, including TO, biotin, and the linker connecting them, abuts one of the near-planar faces of the three-tiered G-quadruplex. The two heterocycles of TO are held in place by two loop adenines and make a 45° angle with respect to each other. Minimizing this angle would increase quantum yield and further improve this tool for in vivo RNA visualization.

A New Method to Predict Ion Effects in RNA Folding

The strong interaction between metal ions in solution and highly charged RNA molecules is critical for RNA structure formation and stabilization. Metal ions binding to RNA can induce RNA structural changes that are important for RNA cellular functions. Therefore, quantitative modeling of the ion effects is essential for RNA structure prediction and RNA-based molecular design. Recently, inspired by the increasing experimental evidence that supports the importance of ion correlation and fluctuation in ion-RNA interactions, we developed a new computational model, Monte Carlo Tightly Bound Ion (MCTBI) model. The validity of the model is shown by the improved accuracy in the predictions for ion binding properties and ion-dependent free energies for RNA structures. In this chapter, using homodimeric tetraloop-receptor docking as an illustrative example, we showcase the MCTBI method for the computational prediction of the ion effects in RNA folding.

Identification of a GUAAY Pentaloop Sequence Involved in a Novel RNA Loop-Helix Interaction

Large RNAs often utilize GNRA tetraloops as structural elements to stabilize the overall tertiary fold. These tetraloop-receptor (TR) interactions have a conserved geometry in which the tetraloop docks into the receptor at an angle of ~15° from the helix containing the receptor. Here, we show that the conserved GUAAY pentaloop found in domain III of group IIB1 introns participates in a novel class of RNA tertiary interaction with a geometry and mode of binding that are significantly different from that found in GNRA TR interactions. This pentaloop is highly conserved within the IIB1 class and interacts with the minor groove of the catalytic domain V. The base planes of the loop and receptor nucleotides are not coplanar and greatly deviate from standard A-minor motifs. The helical axis of the GUAAY stem loop diverges ~70° from the angle of insertion found in a typical GNRA TR interaction. Therefore, the loop architecture and insertion orientation are distinctive, with in vitro splicing data indicating that a GNRA tetraloop is incompatible at this position. The GUAAY pentaloop-receptor motif is also found in the structure of the eukaryotic thiamine pyrophosphate riboswitch in the context of a hexanucleotide loop sequence. We therefore propose, based on phylogenetic, structural, and biochemical data, that the GUAAY pentaloop-receptor interaction represents a novel structural motif that is present in multiple structured RNAs.

Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies

The computer-aided folding of biomolecules, particularly RNAs, is one of the most difficult challenges in computational structural biology. RNA tetraloops are fundamental RNA motifs playing key roles in RNA folding and RNA-RNA and RNA-protein interactions. Although state-of-the-art Molecular Dynamics (MD) force fields correctly describe the native state of these tetraloops as a stable free-energy basin on the microsecond time scale, enhanced sampling techniques reveal that the native state is not the global free energy minimum, suggesting yet unidentified significant imbalances in the force fields. Here, we tested our ability to fold the RNA tetraloops in various force fields and simulation settings. We employed three different enhanced sampling techniques, namely, temperature replica exchange MD (T-REMD), replica exchange with solute tempering (REST2), and well-tempered metadynamics (WT-MetaD). We aimed to separate problems caused by limited sampling from those due to force-field inaccuracies. We found that none of the contemporary force fields is able to correctly describe folding of the 5'-GAGA-3' tetraloop over a range of simulation conditions. We thus aimed to identify which terms of the force field are responsible for this poor description of TL folding. We showed that at least two different imbalances contribute to this behavior, namely, overstabilization of base-phosphate and/or sugar-phosphate interactions and underestimated stability of the hydrogen bonding interaction in base pairing. The first artifact stabilizes the unfolded ensemble, while the second one destabilizes the folded state. The former problem might be partially alleviated by reparametrization of the van der Waals parameters of the phosphate oxygens suggested by Case et al., while in order to overcome the latter effect we suggest local potentials to better capture hydrogen bonding interactions.

Insights into Stability and Folding of GNRA and UNCG Tetraloops Revealed by Microsecond Molecular Dynamics and Well-Tempered Metadynamics

RNA hairpins capped by 5'-GNRA-3' or 5'-UNCG-3' tetraloops (TLs) are prominent RNA structural motifs. Despite their small size, a wealth of experimental data, and recent progress in theoretical simulations of their structural dynamics and folding, our understanding of the folding and unfolding processes of these small RNA elements is still limited. Theoretical description of the folding and unfolding processes requires robust sampling, which can be achieved by either an exhaustive time scale in standard molecular dynamics simulations or sophisticated enhanced sampling methods, using temperature acceleration or biasing potentials. Here, we study structural dynamics of 5'-GNRA-3' and 5'-UNCG-3' TLs by 15-μs-long standard simulations and a series of well-tempered metadynamics, attempting to accelerate sampling by bias in a few chosen collective variables (CVs). Both methods provide useful insights. The unfolding and refolding mechanisms of the GNRA TL observed by well-tempered metadynamics agree with the (reverse) folding mechanism suggested by recent replica exchange molecular dynamics simulations. The orientation of the glycosidic bond of the GL4 nucleobase is critical for the UUCG TL folding pathway, and our data strongly support the hypothesis that GL4-anti forms a kinetic trap along the folding pathway. Along with giving useful insight, our study also demonstrates that using only a few CVs apparently does not capture the full folding landscape of the RNA TLs. Despite using several sophisticated selections of the CVs, formation of the loop appears to remain a hidden variable, preventing a full convergence of the metadynamics. Finally, our data suggest that the unfolded state might be overstabilized by the force fields used.

Hydroxylated histidine of human ribosomal protein uL2 is involved in maintaining the local structure of 28S rRNA in the ribosomal peptidyl transferase center

Protein uL2 is essential for the catalytic activity of the ribosome and has a conserved shape in ribosomes from all domains of life. However, the sequence of its unstructured C-terminal loop apex that contacts the conserved 23S/28S rRNA helix (H) 93 near the ribosomal peptidyl transferase center differs in bacteria, archaea and eukaryotes. Eukaryote-specific residue His216 located in this loop in mammalian uL2 is hydroxylated in ribosomes. We used a set of chemical probes to explore the structure of an RNA that mimicked a segment of 28S rRNA domain V containing part of the uL2 binding site including H93, complexed with either natural (hydroxylated) or recombinant (unmodified) human uL2. It was found that both protein forms engage H93 during binding, but only natural uL2 (uL2n) protects it from hydroxyl radicals. The association of uL2n with RNA leads to changes in its structure at U4532 adjacent to the universally conserved U4531 (U2585, Escherichia coli numbering) involved in peptidyl transferase center formation, and at the universally conserved C4447 (2501) located in the ribosome near A4397 (2451) and C3909 (2063) belonging to the peptidyl transferase center. As a result, both nucleotides become strongly exposed to hydroxyl radicals. Our data argue that the hydroxyl group at His216 in the C-terminal loop apex of mammalian uL2 contributes to stabilization of a protein conformation that is favorable for binding to H93 of 28S rRNA and that this binding induces structural rearrangement in the regions close to the peptidyl transferase center in the mature ribosome.

From static to dynamic: the need for structural ensembles and a predictive model of RNA folding and function

To understand RNA, it is necessary to move beyond a descriptive categorization towards quantitative predictions of its molecular conformations and functional behavior. An incisive approach to understanding the function and folding of biological RNA systems involves characterizing small, simple components that are largely responsible for the behavior of complex systems including helix-junction-helix elements and tertiary motifs. State-of-the-art methods have permitted unprecedented insight into the conformational ensembles of these elements revealing, for example, that conformations of helix-junction-helix elements are confined to a small region of the ensemble, that this region is highly dependent on the junction's topology, and that the correct alignment of tertiary motifs may be a rare conformation on the overall folding landscape. Further characterization of RNA components and continued development of experimental and computational methods with the goal of quantitatively predicting RNA folding and functional behavior will be critical to understanding biological RNA systems.

Tetraloop-like geometries could form the basis of the catalytic activity of the most ancient ribooligonucleotides

The origin of the catalytic activity of ancient oligonucleotides is a largely unexplored field of contemporary science. In the current work we use molecular dynamics simulations to investigate the plausibility of tetraloop-like overhang geometries to initiate transphosphorylation reactions that lead to ligation and terminal cleavage in simple, Watson-Crick (WC) complementary oligoC/oligoG sequences observed experimentally. We show a series of examples of known tetraloop architectures, which can be adopted by the unpaired overhangs of short oligonucleotide sequences for a sufficiently long time to enable chemical reactions that lead to simple ribozyme-like catalytic activity. Thus, our computations demonstrate that the role of non-WC interactions at the emergence of the most ancient catalytic oligonucleotides could be more significant than ever believed.

Annotating RNA motifs in sequences and alignments

RNA performs a diverse array of important functions across all cellular life. These functions include important roles in translation, building translational machinery and maturing messenger RNA. More recent discoveries include the miRNAs and bacterial sRNAs that regulate gene expression, the thermosensors, riboswitches and other cis-regulatory elements that help prokaryotes sense their environment and eukaryotic piRNAs that suppress transposition. However, there can be a long period between the initial discovery of a RNA and determining its function. We present a bioinformatic approach to characterize RNA motifs, which are critical components of many RNA structure-function relationships. These motifs can, in some instances, provide researchers with functional hypotheses for uncharacterized RNAs. Moreover, we introduce a new profile-based database of RNA motifs--RMfam--and illustrate some applications for investigating the evolution and functional characterization of RNA. All the data and scripts associated with this work are available from: https://github.com/ppgardne/RMfam.

Molecular crowding overcomes the destabilizing effects of mutations in a bacterial ribozyme

The native structure of the Azoarcus group I ribozyme is stabilized by the cooperative formation of tertiary interactions between double helical domains. Thus, even single mutations that break this network of tertiary interactions reduce ribozyme activity in physiological Mg(2+) concentrations. Here, we report that molecular crowding comparable to that in the cell compensates for destabilizing mutations in the Azoarcus ribozyme. Small angle X-ray scattering, native polyacrylamide gel electrophoresis and activity assays were used to compare folding free energies in dilute and crowded solutions containing 18% PEG1000. Crowder molecules allowed the wild-type and mutant ribozymes to fold at similarly low Mg(2+) concentrations and stabilized the active structure of the mutant ribozymes under physiological conditions. This compensation helps explains why ribozyme mutations are often less deleterious in the cell than in the test tube. Nevertheless, crowding did not rescue the high fraction of folded but less active structures formed by double and triple mutants. We conclude that crowding broadens the fitness landscape by stabilizing compact RNA structures without improving the specificity of self-assembly.

Computational identification of RNA functional determinants by three-dimensional quantitative structure-activity relationships

Anti-infection drugs target vital functions of infectious agents, including their ribosome and other essential non-coding RNAs. One of the reasons infectious agents become resistant to drugs is due to mutations that eliminate drug-binding affinity while maintaining vital elements. Identifying these elements is based on the determination of viable and lethal mutants and associated structures. However, determining the structure of enough mutants at high resolution is not always possible. Here, we introduce a new computational method, MC-3DQSAR, to determine the vital elements of target RNA structure from mutagenesis and available high-resolution data. We applied the method to further characterize the structural determinants of the bacterial 23S ribosomal RNA sarcin–ricin loop (SRL), as well as those of the lead-activated and hammerhead ribozymes. The method was accurate in confirming experimentally determined essential structural elements and predicting the viability of new SRL variants, which were either observed in bacteria or validated in bacterial growth assays. Our results indicate that MC-3DQSAR could be used systematically to evaluate the drug-target potentials of any RNA sites using current high-resolution structural data.

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.

RNA tectonics (tectoRNA) for RNA nanostructure design and its application in synthetic biology

RNA molecules are versatile biomaterials that act not only as DNA-like genetic materials but also have diverse functions in regulation of cellular biosystems. RNA is capable of regulating gene expression by sequence-specific hybridization. This feature allows the design of RNA-based artificial gene regulators (riboregulators). RNA can also build complex two-dimensional (2D) and 3D nanostructures, which afford protein-like functions and make RNA an attractive material for nanobiotechnology. RNA tectonics is a methodology in RNA nanobiotechnology for the design and construction of RNA nanostructures/nanoobjects through controlled self-assembly of modular RNA units (tectoRNAs). RNA nanostructures designed according to the concept of RNA tectonics are also attractive as tools in synthetic biology, but in vivo RNA tectonics is still in the early stages. This review presents a summary of the achievements of RNA tectonics and its related researches in vitro, and also introduces recent developments that facilitated the use of RNA nanostructures in bacterial cells.

Structural complexity of Dengue virus untranslated regions: cis-acting RNA motifs and pseudoknot interactions modulating functionality of the viral genome

The Dengue virus (DENV) genome contains multiple cis-acting elements required for translation and replication. Previous studies indicated that a 719-nt subgenomic minigenome (DENV-MINI) is an efficient template for translation and (−) strand RNA synthesis in vitro. We performed a detailed structural analysis of DENV-MINI RNA, combining chemical acylation techniques, Pb²⁺ ion-induced hydrolysis and site-directed mutagenesis. Our results highlight protein-independent 5′–3′ terminal interactions involving hybridization between recognized cis-acting motifs. Probing analyses identified tandem dumbbell structures (DBs) within the 3′ terminus spaced by single-stranded regions, internal loops and hairpins with embedded GNRA-like motifs. Analysis of conserved motifs and top loops (TLs) of these dumbbells, and their proposed interactions with downstream pseudoknot (PK) regions, predicted an H-type pseudoknot involving TL1 of the 5′ DB and the complementary region, PK2. As disrupting the TL1/PK2 interaction, via ‘flipping’ mutations of PK2, previously attenuated DENV replication, this pseudoknot may participate in regulation of RNA synthesis. Computer modeling implied that this motif might function as autonomous structural/regulatory element. In addition, our studies targeting elements of the 3′ DB and its complementary region PK1 indicated that communication between 5′–3′ terminal regions strongly depends on structure and sequence composition of the 5′ cyclization region.

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.

Functional replacement of two highly conserved tetraloops in the bacterial ribosome

Ribosomes are RNA-protein complexes responsible for protein synthesis. A dominant structural motif in the rRNAs is an RNA helix capped with a four-nucleotide loop, called a tetraloop. The sequence of the tetraloop is invariant at some positions in the rRNAs but is highly variable at other positions. The biological reason for the conservation of the tetraloop sequence at specific positions in the rRNAs is not clear. In the 16S rRNA, the GAAA tetraloop in helix 8 and the UACG tetraloop in helix 14 are highly conserved and located near the binding site for EF-Tu and EF-G. To investigate whether the structural stability of the tetraloop or the precise sequence of the tetraloop is important for function, we separately changed the GAAA tetraloop in helix 8 to a UACG tetraloop and the UACG tetraloop in helix 14 to a GAAA tetraloop. The effects of the tetraloop replacements on protein synthesis were analyzed in vivo and in vitro. Replacement of the tetraloops in helices 8 and 14 did not significantly affect the growth rate of the Escherichia coli (Δ7rrn) strain. However, the mutant ribosomes showed a slightly reduced rate of protein synthesis in vitro. In addition, we observed a 2-fold increase in the error rate of translation with the mutant ribosomes, which is consistent with an earlier report. Our results suggest that the tetraloops in helices 8 and 14 are highly conserved mainly for their structural stability and the precise sequences of these tetraloops are not critical for protein synthesis.

UNAC tetraloops: to what extent do they mimic GNRA tetraloops?

The structures of four small RNAs each containing a different version of the UNAC loop were determined in solution using NMR spectroscopy and restrained molecular dynamics. The UMAC tetraloops (where M is A or C) exhibited a typical GNRA fold including at least one hydrogen bond between the first U and fourth C. In contrast, UGAC and UUAC tetraloops have a different orientation of the first and fourth residues, such that they do not closely mimic the GNRA fold. Although the UMAC tetraloops are excellent structural mimics of the GNRA tetraloop backbone, sequence comparisons typically do not reveal co-variation between the two loop types. The limited covariation is attributed to differences in the location of potential hydrogen bond donors and acceptors as a result of the replacement of the terminal A of GNRA with C in the UMAC version. Thus, UMAC loops do not readily form the common GNRA tetraloop-receptor interaction. The loop at positions 863-866 in E. coli 16S ribosomal RNA appears to be a major exception. However, in this case the GNRA loop does not in fact engage in the usual base to backbone tertiary interactions. In summary, UMAC loops are not just an alternative sequence version of the GNRA loop family, but instead they expand the types of interactions, or lack thereof, that are possible. From a synthetic biology perspective their inclusion in an artificial RNA may allow the establishment of a stable loop structure while minimizing unwanted long range interactions or permitting alternative long-range interactions. © 2012 Wiley Periodicals, Inc. Biopolymers 97: 617-628, 2012.

Cooperative tertiary interaction network guides RNA folding

Noncoding RNAs form unique 3D structures, which perform many regulatory functions. To understand how RNAs fold uniquely despite a small number of tertiary interaction motifs, we mutated the major tertiary interactions in a group I ribozyme by single-base substitutions. The resulting perturbations to the folding energy landscape were measured using SAXS, ribozyme activity, hydroxyl radical footprinting, and native PAGE. Double- and triple-mutant cycles show that most tertiary interactions have a small effect on the stability of the native state. Instead, the formation of core and peripheral structural motifs is cooperatively linked in near-native folding intermediates, and this cooperativity depends on the native helix orientation. The emergence of a cooperative interaction network at an early stage of folding suppresses nonnative structures and guides the search for the native state. We suggest that cooperativity in noncoding RNAs arose from natural selection of architectures conducive to forming a unique, stable fold.

Determining RNA three-dimensional structures using low-resolution data

Knowing the 3-D structure of an RNA is fundamental to understand its biological function. Nowadays X-ray crystallography and NMR spectroscopy are systematically applied to newly discovered RNAs. However, the application of these high-resolution techniques is not always possible, and thus scientists must turn to lower resolution alternatives. Here, we introduce a pipeline to systematically generate atomic resolution 3-D structures that are consistent with low-resolution data sets. We compare and evaluate the discriminative power of a number of low-resolution experimental techniques to reproduce the structure of the Escherichia coli tRNA(VAL) and P4-P6 domain of the Tetrahymena thermophila group I intron. We test single and combinations of the most accessible low-resolution techniques, i.e. hydroxyl radical footprinting (OH), methidiumpropyl-EDTA (MPE), multiplexed hydroxyl radical cleavage (MOHCA), and small-angle X-ray scattering (SAXS). We show that OH-derived constraints are accurate to discriminate structures at the atomic level, whereas EDTA-based constraints apply to global shape determination. We provide a guide for choosing which experimental techniques or combination of thereof is best in which context. The pipeline represents an important step towards high-throughput low-resolution RNA structure determination.

In vitro selection and characterization of a novel Zn(II)-dependent phosphorothiolate thiolesterase ribozyme

Here we present the in vitro selection of a novel ribozyme specific for Zn2+-dependent catalysis on hydrolysis of a phosphorothiolate thiolester bond. The ribozyme, called the TW17 ribozyme, was evolved and selected from an artificial RNA pool covalently linked to a biotin-containing substrate through the phosphorothiolate thiolester bond. The secondary structure for the evolved ribozyme consisted of three major helices and three loops. Biochemical and chemical studies of ribozyme-catalyzed reaction products provided evidence that the ribozyme specifically catalyzes hydrolysis of the phosphorothiolate thiolester linkage. A successful ribozyme construct with active catalysis in trans further supported the determined ribozyme structure and indicated the potential of the ribozyme for multiple-substrate turnover. The ribozyme also requires Zn2+ and Mg2+ for maximal catalysis. The TW17 ribozyme, in the presence of Zn2+ and Mg2+, conferred a rate enhancement of at least 5 orders of magnitude when compared to the estimated rate of the uncatalyzed reaction. The ribozyme completely lost catalytic activity in the absence of Zn2+, like Zn2+-dependent protein hydrolases. The discovery and characterization of the TW17 ribozyme suggest additional roles for Zn2+ in ribozyme catalysts.

Attenuation of loop-receptor interactions with pseudoknot formation

RNA tetraloops can recognize receptors to mediate long-range interactions in stable natural RNAs. In vitro selected GNRA tetraloop/receptor interactions are usually more ‘G/C-rich’ than their ‘A/U-rich’ natural counterparts. They are not as widespread in nature despite comparable biophysical and chemical properties. Moreover, while AA, AC and GU dinucleotide platforms occur in natural GAAA/11 nt receptors, the AA platform is somewhat preferred to the others. The apparent preference for ‘A/U-rich’ GNRA/receptor interactions in nature might stem from an evolutionary adaptation to avoid folding traps at the level of the larger molecular context. To provide evidences in favor of this hypothesis, several riboswitches based on natural and artificial GNRA receptors were investigated in vitro for their ability to prevent inter-molecular GNRA/receptor interactions by trapping the receptor sequence into an alternative intra-molecular pseudoknot. Extent of attenuation determined by native gel-shift assays and co-transcriptional assembly is correlated to the G/C content of the GNRA receptor. Our results shed light on the structural evolution of natural long-range interactions and provide design principles for RNA-based attenuator devices to be used in synthetic biology and RNA nanobiotechnology.

Structural determinants and mechanism of HIV-1 genome packaging

Like all retroviruses, the human immunodeficiency virus selectively packages two copies of its unspliced RNA genome, both of which are utilized for strand-transfer-mediated recombination during reverse transcription-a process that enables rapid evolution under environmental and chemotherapeutic pressures. The viral RNA appears to be selected for packaging as a dimer, and there is evidence that dimerization and packaging are mechanistically coupled. Both processes are mediated by interactions between the nucleocapsid domains of a small number of assembling viral Gag polyproteins and RNA elements within the 5'-untranslated region of the genome. A number of secondary structures have been predicted for regions of the genome that are responsible for packaging, and high-resolution structures have been determined for a few small RNA fragments and protein-RNA complexes. However, major questions regarding the RNA structures (and potentially the structural changes) that are responsible for dimeric genome selection remain unanswered. Here, we review efforts that have been made to identify the molecular determinants and mechanism of human immunodeficiency virus type 1 genome packaging.

Feature-based 3D motif filtering for ribosomal RNA

MOTIVATION

RNA 3D motifs are recurrent substructures in an RNA subunit and are building blocks of the RNA architecture. They play an important role in binding proteins and consolidating RNA tertiary structures. RNA 3D motif searching consists of two steps: candidate generation and candidate filtering. We proposed a novel method, known as Feature-based RNA Motif Filtering (FRMF), for identifying motifs based on a set of moment invariants and the Earth Mover's Distance in the second step.

RESULTS

A positive set of RNA motifs belonging to six characteristic types, with eight subtypes occurring in HM 50S, is compiled by us. The proposed method is validated on this representative set. FRMF successfully finds most of the positive fragments. Besides the proposed new method and the compiled positive set, we also recognize some new motifs, in particular a π-turn and some non-standard A-minor motifs are found. These newly discovered motifs provide more information about RNA structure conformation.

AVAILABILITY

Matlab code can be downloaded from www.cs.cityu.edu.hk/~yingshen/FRMF.html

CONTACT

cshswong@cityu.edu.hk

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

RNA dynamics: perspectives from spin labels

Dynamics are important and indispensible physical attributes that play essential roles in RNA function. RNA dynamics are complex, spanning vast timescales, and encompassing a large number of physical modes. The technique of site-directed spin labeling (SDSL), which derives information on local structural and dynamic features of a macromolecule by monitoring a chemically stable nitroxide radical using electron paramagnetic resonance spectroscopy, has been applied to monitor intrinsic dynamics at defined structural states as well as to probe conformational transition dynamics of RNAs. The current state of SDSL studies of RNA dynamics is summarized here. Further development and application of SDSL promise to open up many more opportunities for probing RNA dynamics and connecting dynamics to structure and function.

Riboswitches: structures and mechanisms

A critical feature of the hypothesized RNA world would have been the ability to control chemical processes in response to environmental cues. Riboswitches present themselves as viable candidates for a sophisticated mechanism of regulatory control in RNA-based life. These regulatory elements in the modern world are most commonly found in the 5'-untranslated regions of bacterial mRNAs, directly interacting with metabolites as a means of regulating expression of the coding region via a secondary structural switch. In this review, we focus on recent insights into how these RNAs fold into complex architectures capable of both recognizing a specific small molecule compound and exerting regulatory control over downstream sequences, with an emphasis on transcriptional regulation.