RNA Structure: Reading the Ribosome

Positioning of CCA-arms of the A- and the P-tRNAs towards the 28S rRNA in the human ribosome

Nucleotides of 28S rRNA involved in binding of the human 80S ribosome with acceptor ends of the A site and the P site tRNAs were determined using two complementary approaches, namely, cross-linking with application of tRNA(Asp) analogues substituted with 4-thiouridine in position 75 or 76 and hydroxyl radical footprinting with the use of the full sized tRNA and the tRNA deprived of the 3'-terminal trinucleotide CCA. In general, these 28S rRNA nucleotides are located in ribosomal regions homologous to the A, P and E sites of the prokaryotic 50S subunit. However, none of the approaches used discovered interactions of the apex of the large rRNA helix 80 with the acceptor end of the P site tRNA typical with prokaryotic ribosomes. Application of the results obtained to available atomic models of 50S and 60S subunits led us to a conclusion that the A site tRNA is actually present in both A/A and A/P states and the P site tRNA in the P/P and P/E states. Thus, the present study gives a biochemical confirmation of the data on the structure and dynamics of the mammalian ribosomal pretranslocation complex obtained with application of cryo-electron microscopy and single-molecule FRET [Budkevich et al., 2011]. Moreover, in our study, particular sets of 28S rRNA nucleotides involved in oscillations of tRNAs CCA-termini between their alternative locations in the mammalian 80S ribosome are revealed.

Toward a molecular understanding of RNA remodeling by DEAD-box proteins

DEAD-box proteins are superfamily 2 helicases that function in all aspects of RNA metabolism. They employ ATP binding and hydrolysis to generate tight, yet regulated RNA binding, which is used to unwind short RNA helices non-processively and promote structural transitions of RNA and RNA-protein substrates. In the last few years, substantial progress has been made toward a detailed, quantitative understanding of the structural and biochemical properties of DEAD-box proteins. Concurrently, progress has been made toward a physical understanding of the RNA rearrangements and folding steps that are accelerated by DEAD-box proteins in model systems. Here, we review the recent progress on both of these fronts, focusing on the mitochondrial DEAD-box proteins Mss116 and CYT-19 and their mechanisms in promoting the splicing of group I and group II introns.

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.

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.

Interaction of nucleic acids with Coomassie Blue G-250 in the Bradford assay

The Bradford assay has been used reliably for decades to quantify protein in solution. The analyte is incubated in acidic solution of Coomassie Blue G-250 dye, during which reversible ionic and nonionic binding interactions form. Bradford assay color yields were determined for salmon, bovine, shrimp, and kiwi fruit genomic DNA; baker's yeast RNA; bovine serum albumin (BSA); and hen egg lysozyme. Pure DNA and RNA bound the dye, with color yields of 0.0017 mg⁻¹ cm⁻¹ and 0.0018 mg⁻¹ cm⁻¹, respectively. The nucleic acid-Coomassie Blue response was significant, at roughly 9% of that for BSA and 18% of that for lysozyme.

Domain III of the T. thermophilus 23S rRNA folds independently to a near-native state

The three-dimensional structure of the ribosomal large subunit (LSU) reveals a single morphological element, although the 23S rRNA is contained in six secondary structure domains. Based upon maps of inter- and intra-domain interactions and proposed evolutionary pathways of development, we hypothesize that Domain III is a truly independent structural domain of the LSU. Domain III is primarily stabilized by intra-domain interactions, negligibly perturbed by inter-domain interactions, and is not penetrated by ribosomal proteins or other rRNA. We have probed the structure of Domain III rRNA alone and when contained within the intact 23S rRNA using SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension), in the absence and presence of magnesium. The combined results support the hypothesis that Domain III alone folds to a near-native state with secondary structure, intra-domain tertiary interactions, and inter-domain interactions that are independent of whether or not it is embedded in the intact 23S rRNA or within the LSU. The data presented support previous suggestions that Domain III was added relatively late in ribosomal evolution.

Another burst of smoke: atomic resolution structures of RF3 bound to the ribosome

Two recent reports provide atomic resolution information detailing the interaction of the class II release factor, RF3, with the bacterial ribosome. Differences in the composition of the two crystal forms allow us to learn a considerable amount about how translational GTPases engage the ribosome to facilitate and define conformational rearrangements involved in protein synthesis.

Structural architecture of the human long non-coding RNA, steroid receptor RNA activator

While functional roles of several long non-coding RNAs (lncRNAs) have been determined, the molecular mechanisms are not well understood. Here, we report the first experimentally derived secondary structure of a human lncRNA, the steroid receptor RNA activator (SRA), 0.87 kB in size. The SRA RNA is a non-coding RNA that coactivates several human sex hormone receptors and is strongly associated with breast cancer. Coding isoforms of SRA are also expressed to produce proteins, making the SRA gene a unique bifunctional system. Our experimental findings (SHAPE, in-line, DMS and RNase V1 probing) reveal that this lncRNA has a complex structural organization, consisting of four domains, with a variety of secondary structure elements. We examine the coevolution of the SRA gene at the RNA structure and protein structure levels using comparative sequence analysis across vertebrates. Rapid evolutionary stabilization of RNA structure, combined with frame-disrupting mutations in conserved regions, suggests that evolutionary pressure preserves the RNA structural core rather than its translational product. We perform similar experiments on alternatively spliced SRA isoforms to assess their structural features.

The spectroscopic basis of fluorescence triple correlation spectroscopy

We have developed fluorescence triple correlation spectroscopy (F3CS) as an extension of the widely used fluorescence microscopy technique fluorescence correlation spectroscopy. F3CS correlates three signals at once and provides additional capabilities for the study of systems with complex stoichiometry, kinetic processes, and irreversible reactions. A general theory of F3CS was developed to describe the interplay of molecular dynamics and microscope optics, leading to an analytical function to predict experimental triple correlations of molecules that freely diffuse through the tight focus of the microscope. Experimental correlations were calculated from raw fluorescence data using triple correlation integrals that extend multiple-tau correlation theory to delay times in two dimensions. The quality of experimental data was improved by tuning specific spectroscopic parameters and employing multiple independent detectors to minimize optoelectronic artifacts. Experiments with the reversible system of freely diffusing 16S rRNA revealed that triple correlation functions contain symmetries predicted from time-reversal arguments. Irreversible systems are shown to break these symmetries, and correlation strategies were developed to detect time-reversal asymmetries in a comprehensive way with respect to two delay times, each spanning many orders of magnitude in time. The correlation strategies, experimental approaches, and theory developed here enable studies of the composition and dynamics of complex systems using F3CS.

The HIV-1 Rev response element: an RNA scaffold that directs the cooperative assembly of a homo-oligomeric ribonucleoprotein complex

The HIV-1 Rev response element (RRE) is a ~350 nucleotide, highly structured, cis-acting RNA element essential for viral replication. It is located in the env coding region of the viral genome and is extremely well conserved across different HIV-1 isolates. It is present on all partially spliced and unspliced viral mRNA transcripts, and serves as an RNA framework onto which multiple molecules of the viral protein Rev assemble. The Rev-RRE oligomeric complex mediates the export of these messages from the nucleus to the cytoplasm, where they are translated to produce essential viral proteins and/or packaged as genomes for new virions.

The formation of a potential spring in the ribosome

Time-dependent chemical modification and cleavage results have provided intriguing insights into structural changes that occur in the distal loop of helix 11 in 16S ribosomal RNA (rRNA). Located distant from the decoding region, between proteins S17 and S20, the results of this study suggest that this region of rRNA may act as a buffer or a spring between these two proteins during protein biosynthesis. During the assembly process, protein S17 apparently produces the major structural deformations in this region, causing it to be folded in a spring-like structure. Base C264 in this region shows erratic behavior during assembly and also shows time-dependent enhancement when elongation factor G with GTP is added to 70S ribosomes. Evidence is presented to suggest that this region of rRNA may be used to allow relative motion to occur between proteins S17 and S20 during translocation.

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.

Nonparametric clustering for studying RNA conformations

The local conformation of RNA molecules is an important factor in determining their catalytic and binding properties. The analysis of such conformations is particularly difficult due to the large number of degrees of freedom, such as the measured torsion angles per residue and the interatomic distances among interacting residues. In this work, we use a nearest-neighbor search method based on the statistical mechanical Potts model to find clusters in the RNA conformational space. The proposed technique is mostly automatic and may be applied to problems, where there is no prior knowledge on the structure of the data space in contrast to many other clustering techniques. Results are reported for both single residue conformations, where the parameter set of the data space includes four to seven torsional angles, and base pair geometries, where the data space is reduced to two dimensions. Moreover, new results are reported for base stacking geometries. For the first two cases, i.e., single residue conformations and base pair geometries, we get a very good match between the results of the proposed clustering method and the known classifications with only few exceptions. For the case of base stacking geometries, we validate our classification with respect to geometrical constraints and describe the content, and the geometry of the new clusters.

Mechanism of antisense oligonucleotide interaction with natural RNAs

Oligonucleotides find several numbers of applications: as diagnostic probes, RT and PCR primers and antisense agents due to their ability of forming specific interactions with complementary nucleotide sequences within nucleic acids. These interactions are strongly affected by accessibility of the target sequence in the RNA structure. In the present work the mechanism of invasion of RNA structure by oligonucleotide was investigated using a model system: yeast tRNA(Phe) and oligonucleotides complementary to the 3'-part of this molecule. Kinetics of interaction of oligonucleotides with in vitro transcript of yeast tRNAPhe was studied using stopped-flow technique with fluorescence quenching detection, 5'-DABCYL labeled oligonucleotide was hybridized with 3'-fluorescein labeled tRNA(Phe). The results evidence for a four-step invasion process of the oligonucleotide-RNA complex formation. The process is initiated by formation of transition complexes with nucleotides in the T-loop and ACCA sequence. This complex formation is followed by RNA unfolding and formation of an extended heteroduplex with the oligonucleotide via strand displacement process. Computer modeling of oligonucleotide-tRNA(Phe) interaction revealed potential factors that could favor transition complexes formation and confirmed the proposed mechanism, showing the oligonucleotide to be a molecular "wedge". Our data evidence that oligonucleotide invasion into structured RNA is initiated by loop-single strand interactions, similar to the initial step of the antisense RNA-RNA interactions. The obtained results can be used for choosing efficient oligonucleotide probes.

Assembly of bacterial ribosomes

The assembly of ribosomes from a discrete set of components is a key aspect of the highly coordinated process of ribosome biogenesis. In this review, we present a brief history of the early work on ribosome assembly in Escherichia coli, including a description of in vivo and in vitro intermediates. The assembly process is believed to progress through an alternating series of RNA conformational changes and protein-binding events; we explore the effects of ribosomal proteins in driving these events. Ribosome assembly in vivo proceeds much faster than in vitro, and we outline the contributions of several of the assembly cofactors involved, including Era, RbfA, RimJ, RimM, RimP, and RsgA, which associate with the 30S subunit, and CsdA, DbpA, Der, and SrmB, which associate with the 50S subunit.

Tertiary interactions between helices h13 and h44 in 16S RNA contribute to the fidelity of translation

The A-minor interaction, formed between single-stranded adenosines and the minor groove of a receptor helix, is among the most common motifs found in rRNA. Among the A-minors found in 16S rRNA are a set of interactions between adenosines at positions 1433, 1434 and 1468 in helix 44 (h44) and their receptors in the nucleotide 320-340 region of helix 13 (h13). These interactions have been implicated in the maintenance of translational accuracy, because base substitutions at the adjacent C1469 increase miscoding errors. We have tested their functional significance through mutagenesis of h13 and h44. Mutations at the h44 A residues, or the A-minor receptors in h13, increase a variety of translational errors and a subset of the mutants show decreased association between 30S and 50S ribosomal subunits. These results are consistent with the involvement of h13-h44 interactions in the alignment and packing of these helices in the 30S subunit and the importance of this helical alignment for tRNA selection and subunit-subunit interaction.

Understanding the errors of SHAPE-directed RNA structure modeling

Single-nucleotide-resolution chemical mapping for structured RNA is being rapidly advanced by new chemistries, faster readouts, and coupling to computational algorithms. Recent tests have shown that selective 2'-hydroxyl acylation by primer extension (SHAPE) can give near-zero error rates (0-2%) in modeling the helices of RNA secondary structure. Here, we benchmark the method using six molecules for which crystallographic data are available: tRNA(phe) and 5S rRNA from Escherichia coli, the P4-P6 domain of the Tetrahymena group I ribozyme, and ligand-bound domains from riboswitches for adenine, cyclic di-GMP, and glycine. SHAPE-directed modeling of these highly structured RNAs gave an overall false negative rate (FNR) of 17% and a false discovery rate (FDR) of 21%, with at least one helix prediction error in five of the six cases. Extensive variations of data processing, normalization, and modeling parameters did not significantly mitigate modeling errors. Only one varation, filtering out data collected with deoxyinosine triphosphate during primer extension, gave a modest improvement (FNR = 12%, and FDR = 14%). The residual structure modeling errors are explained by the insufficient information content of these RNAs' SHAPE data, as evaluated by a nonparametric bootstrapping analysis. Beyond these benchmark cases, bootstrapping suggests a low level of confidence (<50%) in the majority of helices in a previously proposed SHAPE-directed model for the HIV-1 RNA genome. Thus, SHAPE-directed RNA modeling is not always unambiguous, and helix-by-helix confidence estimates, as described herein, may be critical for interpreting results from this powerful methodology.

ATP-dependent roles of the DEAD-box protein Mss116p in group II intron splicing in vitro and in vivo

The yeast DEAD-box protein Mss116p functions as a general RNA chaperone in splicing mitochondrial group I and group II introns. For most of its functions, Mss116p is thought to use ATP-dependent RNA unwinding to facilitate RNA structural transitions, but it has been suggested to assist in the folding of one group II intron (aI5γ) primarily by stabilizing a folding intermediate. Here we compare three aI5γ constructs: one with long exons, one with short exons, and a ribozyme construct lacking exons. The long exons result in slower splicing, suggesting that they misfold and/or stabilize nonnative intronic structures. Nevertheless, Mss116p acceleration of all three constructs depends on ATP and is inhibited by mutations that compromise RNA unwinding, suggesting similar mechanisms. Results of splicing assays and a new two-stage assay that separates ribozyme folding and catalysis indicate that maximal folding of all three constructs by Mss116p requires ATP-dependent RNA unwinding. ATP-independent activation is appreciable for only a subpopulation of the minimal ribozyme construct and not for constructs containing exons. As expected for a general RNA chaperone, Mss116p can also disrupt the native ribozyme, which can refold after Mss116p removal. Finally, using yeast strains with mitochondrial DNA containing only the single intron aI5γ, we show that Mss116p mutants promote splicing in vivo to degrees that correlate with their residual ATP-dependent RNA-unwinding activities. Together, our results indicate that, although DEAD-box proteins play multiple roles in RNA folding, the physiological function of Mss116p in aI5γ splicing includes a requirement for ATP-dependent local unfolding, allowing the conversion of nonfunctional RNA structure into functional RNA structure.

Prediction of RNA-binding proteins by voting systems

It is important to identify which proteins can interact with RNA for the purpose of protein annotation, since interactions between RNA and proteins influence the structure of the ribosome and play important roles in gene expression. This paper tries to identify proteins that can interact with RNA using voting systems. Firstly through Weka, 34 learning algorithms are chosen for investigation. Then simple majority voting system (SMVS) is used for the prediction of RNA-binding proteins, achieving average ACC (overall prediction accuracy) value of 79.72% and MCC (Matthew's correlation coefficient) value of 59.77% for the independent testing dataset. Then mRMR (minimum redundancy maximum relevance) strategy is used, which is transferred into algorithm selection. In addition, the MCC value of each classifier is assigned to be the weight of the classifier's vote. As a result, best average MCC values are attained when 22 algorithms are selected and integrated through weighted votes, which are 64.70% for the independent testing dataset, and ACC value is 82.04% at this moment.

Classification of pseudo pairs between nucleotide bases and amino acids by analysis of nucleotide-protein complexes

Nucleotide bases are recognized by amino acid residues in a variety of DNA/RNA binding and nucleotide binding proteins. In this study, a total of 446 crystal structures of nucleotide–protein complexes are analyzed manually and pseudo pairs together with single and bifurcated hydrogen bonds observed between bases and amino acids are classified and annotated. Only 5 of the 20 usual amino acid residues, Asn, Gln, Asp, Glu and Arg, are able to orient in a coplanar fashion in order to form pseudo pairs with nucleotide bases through two hydrogen bonds. The peptide backbone can also form pseudo pairs with nucleotide bases and presents a strong bias for binding to the adenine base. The Watson–Crick side of the nucleotide bases is the major interaction edge participating in such pseudo pairs. Pseudo pairs between the Watson–Crick edge of guanine and Asp are frequently observed. The Hoogsteen edge of the purine bases is a good discriminatory element in recognition of nucleotide bases by protein side chains through the pseudo pairing: the Hoogsteen edge of adenine is recognized by various amino acids while the Hoogsteen edge of guanine is only recognized by Arg. The sugar edge is rarely recognized by either the side-chain or peptide backbone of amino acid residues.

A quantitative measure of chirality inside nucleic acid databank

We show the capability of a chirality index (Pietropaolo et al., Proteins 2008;70:667-677) to investigate nucleic acid structures because of its high sensitivity to helical conformations. By analyzing selected structures of DNA and RNA, we have found that sequences rich in cytosine and guanine have a tendency to left-handed chirality, in contrast to regions rich in adenine or thymine which show strong negative, right-handed, chirality values. We also analyze RNA structures, where specific loops and hairpin motifs are characterized by a well-defined chirality value. We find that in nucleosome the chirality is exalted, whereas in ribosome it is reduced. Our results illustrate the sensitivity of this descriptor for nucleic acid conformations.

Removal of covalent heterogeneity reveals simple folding behavior for P4-P6 RNA

RNA folding landscapes have been described alternately as simple and as complex. The limited diversity of RNA residues and the ability of RNA to form stable secondary structures prior to adoption of a tertiary structure would appear to simplify folding relative to proteins. Nevertheless, there is considerable evidence for long-lived misfolded RNA states, and these observations have suggested rugged energy landscapes. Recently, single molecule fluorescence resonance energy transfer (smFRET) studies have exposed heterogeneity in many RNAs, consistent with deeply furrowed rugged landscapes. We turned to an RNA of intermediate complexity, the P4-P6 domain from the Tetrahymena group I intron, to address basic questions in RNA folding. P4-P6 exhibited long-lived heterogeneity in smFRET experiments, but the inability to observe exchange in the behavior of individual molecules led us to probe whether there was a non-conformational origin to this heterogeneity. We determined that routine protocols in RNA preparation and purification, including UV shadowing and heat annealing, cause covalent modifications that alter folding behavior. By taking measures to avoid these treatments and by purifying away damaged P4-P6 molecules, we obtained a population of P4-P6 that gave near-uniform behavior in single molecule studies. Thus, the folding landscape of P4-P6 lacks multiple deep furrows that would trap different P4-P6 molecules in different conformations and contrasts with the molecular heterogeneity that has been seen in many smFRET studies of structured RNAs. The simplicity of P4-P6 allowed us to reliably determine the thermodynamic and kinetic effects of metal ions on folding and to now begin to build more detailed models for RNA folding behavior.

Structural diversity of eukaryotic 18S rRNA and its impact on alignment and phylogenetic reconstruction

Ribosomal RNAs are important because they catalyze the synthesis of peptides and proteins. Comparative studies of the secondary structure of 18S rRNA have revealed the basic locations of its many length-conserved and length-variable regions. In recent years, many more sequences of 18S rDNA with unusual lengths have been documented in GenBank. These data make it possible to recognize the diversity of the secondary and tertiary structures of 18S rRNAs and to identify the length-conserved parts of 18S rDNAs. The longest 18S rDNA sequences of almost every known eukaryotic phylum were included in this study. We illustrated the bioinformatics-based structure to show that, the regions that are more length-variable, regions that are less length-variable, the splicing sites for introns, and the sites of A-minor interactions are mostly distributed in different parts of the 18S rRNA. Additionally, this study revealed that some length-variable regions or insertion positions could be quite close to the functional part of the 18S rRNA of Foraminifera organisms. The tertiary structure as well as the secondary structure of 18S rRNA can be more diverse than what was previously supposed. Besides revealing how this interesting gene evolves, it can help to remove ambiguity from the alignment of eukaryotic 18S rDNAs and to improve the performance of 18S rDNA in phylogenetic reconstruction. Six nucleotides shared by Archaea and Eukaryota but rarely by Bacteria are also reported here for the first time, which might further support the supposed origin of eukaryote from archaeans.

A mutate-and-map strategy accurately infers the base pairs of a 35-nucleotide model RNA

We present a rapid experimental strategy for inferring base pairs in structured RNAs via an information-rich extension of classic chemical mapping approaches. The mutate-and-map method, previously applied to a DNA/RNA helix, systematically searches for single mutations that enhance the chemical accessibility of base-pairing partners distant in sequence. To test this strategy for structured RNAs, we have carried out mutate-and-map measurements for a 35-nt hairpin, called the MedLoop RNA, embedded within an 80-nt sequence. We demonstrate the synthesis of all 105 single mutants of the MedLoop RNA sequence and present high-throughput DMS, CMCT, and SHAPE modification measurements for this library at single-nucleotide resolution. The resulting two-dimensional data reveal visually clear, punctate features corresponding to RNA base pair interactions as well as more complex features; these signals can be qualitatively rationalized by comparison to secondary structure predictions. Finally, we present an automated, sequence-blind analysis that permits the confident identification of nine of the 10 MedLoop RNA base pairs at single-nucleotide resolution, while discriminating against all 1460 false-positive base pairs. These results establish the accuracy and information content of the mutate-and-map strategy and support its feasibility for rapidly characterizing the base-pairing patterns of larger and more complex RNA systems.

A three-dimensional RNA motif in Potato spindle tuber viroid mediates trafficking from palisade mesophyll to spongy mesophyll in Nicotiana benthamiana

Cell-to-cell trafficking of RNA is an emerging biological principle that integrates systemic gene regulation, viral infection, antiviral response, and cell-to-cell communication. A key mechanistic question is how an RNA is specifically selected for trafficking from one type of cell into another type. Here, we report the identification of an RNA motif in Potato spindle tuber viroid (PSTVd) required for trafficking from palisade mesophyll to spongy mesophyll in Nicotiana benthamiana leaves. This motif, called loop 6, has the sequence 5'-CGA-3'...5'-GAC-3' flanked on both sides by cis Watson-Crick G/C and G/U wobble base pairs. We present a three-dimensional (3D) structural model of loop 6 that specifies all non-Watson-Crick base pair interactions, derived by isostericity-based sequence comparisons with 3D RNA motifs from the RNA x-ray crystal structure database. The model is supported by available chemical modification patterns, natural sequence conservation/variations in PSTVd isolates and related species, and functional characterization of all possible mutants for each of the loop 6 base pairs. Our findings and approaches have broad implications for studying the 3D RNA structural motifs mediating trafficking of diverse RNA species across specific cellular boundaries and for studying the structure-function relationships of RNA motifs in other biological processes.

Quantum chemical studies of nucleic acids: can we construct a bridge to the RNA structural biology and bioinformatics communities?

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.

Absolute binding free energy calculations of sparsomycin analogs to the bacterial ribosome

The interactions of the 50S subunit of bacterial ribosome with antibiotic sparsomycin (SPS) and five analogs (AN) are investigated through the calculation of the standard (absolute) binding free energy and the characterization of conformational dynamics. The standard binding free energies of the complexes are computed using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. Restraining potentials are applied and then released during the simulation to efficiently sample the changes in translational, orientational, and conformational freedom of the ligand and receptor upon binding. The biasing effects of the restraining potentials are rigorously removed. The loss of conformational freedom of the ligand upon binding is determined by introducing a potential of mean force (PMF) as a function of the root-mean-square deviation (rmsd) of the ligand relative to its conformation in the bound state. To reduce the size of the simulated system, the binding pocket of the ribosome is simulated in the framework of the generalized solvent boundary potential (GSBP). The number of solvent molecules in the buried binding site is treated via grand canonical Monte Carlo (GCMC) during the FEP/MD simulations. The correlation coefficient between the calculated and measured binding free energies is 0.96, and the experimentally observed ranking order for the binding affinities of the six ligands is reproduced. However, while the calculated affinities of the strong binders agree well with the experimental values, those for the weak binders are underestimated.

Structural dynamics of the box C/D RNA kink-turn and its complex with proteins: the role of the A-minor 0 interaction, long-residency water bridges, and structural ion-binding sites revealed by molecular simulations

Kink-turns (K-turns) are recurrent elbow-like RNA motifs that participate in protein-assisted RNA folding and contribute to RNA dynamics. We carried out a set of molecular dynamics (MD) simulations using parm99 and parmbsc0 force fields to investigate structural dynamics of the box C/D RNA and its complexes with two proteins: native archaeal L7ae protein and human 15.5 kDa protein, originally bound to very similar structure of U4 snRNA. The box C/D RNA forms K-turn with A-minor 0 tertiary interaction between its canonical (C) and noncanonical (NC) stems. The local K-turn architecture is thus different from the previously studied ribosomal K-turns 38 and 42 having A-minor I interaction. The simulations reveal visible structural dynamics of this tertiary interaction involving altogether six substates which substantially contribute to the elbow-like flexibility of the K-turn. The interaction can even temporarily shift to the A-minor I type pattern; however, this is associated with distortion of the G/A base pair in the NC-stem of the K-turn. The simulations show reduction of the K-turn flexibility upon protein binding. The protein interacts with the apex of the K-turn and with the NC-stem. The protein-RNA interface includes long-residency hydration sites. We have also found long-residency hydration sites and major ion-binding sites associated with the K-turn itself. The overall topology of the K-turn remains stable in all simulations. However, in simulations of free K-turn, we observed instability of the key C16(O2')-A7(N1) H-bond, which is a signature interaction of K-turns and which was visibly more stable in simulations of K-turns possessing A-minor I interaction. It may reflect either some imbalance of the force field or it may be a correct indication of early stages of unfolding since this K-turn requires protein binding for its stabilization. Interestingly, the 16(O2')-7(N1) H- bond is usually not fully lost since it is replaced by a water bridge with a tightly bound water, which is adenine-specific similarly as the original interaction. The 16(O2')-7(N1) H-bond is stabilized by protein binding. The stabilizing effect is more visible with the human 15.5 kDa protein, which is attributed to valine to arginine substitution in the binding site. The behavior of the A-minor interaction is force-field-dependent because the parmbsc0 force field attenuates the A-minor fluctuations compared to parm99 simulations. Behavior of other regions of the box C/D RNA is not sensitive to the force field choice. Simulation with net-neutralizing Na(+) and 0.2 M excess salt conditions appear in all aspects equivalent. The simulations show loss of a hairpin tetraloop, which is not part of the K-turn. This was attributed to force field limitations.

Roles of DEAD-box proteins in RNA and RNP Folding

RNAs and RNA-protein complexes (RNPs) traverse rugged energy landscapes as they fold to their native structures, and many continue to undergo conformational rearrangements as they function. Due to the inherent stability of local RNA structure, proteins are required to assist with RNA conformational transitions during initial folding and in exchange between functional structures. DEAD-box proteins are superfamily 2 RNA helicases that are ubiquitously involved in RNA-mediated processes. Some of these proteins use an ATP-dependent cycle of conformational changes to disrupt RNA structure nonprocessively, accelerating structural transitions of RNAs and RNPs in a manner that bears a strong resemblance to the activities of certain groups of protein chaperones. This review summarizes recent work using model substrates and tractable self-splicing intron RNAs, which has given new insights into how DEAD-box proteins promote RNA folding steps and conformational transitions, and it summarizes recent progress in identifying sites and mechanisms of DEAD-box protein activity within more complex cellular targets.

Yeast 18 S rRNA is directly involved in the ribosomal response to stringent AUG selection during translation initiation

In eukaryotes, the 40 S ribosomal subunit serves as the platform of initiation factor assembly, to place itself precisely on the AUG start codon. Structural arrangement of the 18 S rRNA determines the overall shape of the 40 S subunit. Here, we present genetic evaluation of yeast 18 S rRNA function using 10 point mutations altering the polysome profile. All the mutants reduce the abundance of the mutant 40 S, making it limiting for translation initiation. Two of the isolated mutations, G875A, altering the core of the platform domain that binds eIF1 and eIF2, and A1193U, changing the h31 loop located below the P-site tRNA(i)(Met), show phenotypes indicating defective regulation of AUG selection. Evidence is provided that these mutations reduce the interaction with the components of the preinitiation complex, thereby inhibiting its function at different steps. These results indicate that the 18 S rRNA mutations impair the integrity of scanning-competent preinitiation complex, thereby altering the 40 S subunit response to stringent AUG selection. Interestingly, nine of the mutations alter the body/platform domains of 18 S rRNA, potentially affecting the bridges to the 60 S subunit, but they do not change the level of 18 S rRNA intermediates. Based on these results, we also discuss the mechanism of the selective degradation of the mutant 40 S subunits.

A mutate-and-map strategy for inferring base pairs in structured nucleic acids: proof of concept on a DNA/RNA helix

We propose a rapid chemical strategy for identifying base pairs in structured nucleic acid systems. The approach goes beyond traditional chemical mapping approaches by monitoring perturbations of each residue's chemical accessibility in response to systematic mutagenesis of residues that are distant in sequence but nearby in three dimensions. As a proof of concept, we present high-throughput dimethyl sulfate accessibility data for a chimeric DNA/RNA system in which every possible sequence variation and deletion in a 20 bp region has been synthesized and tested. The data demonstrate that 88% of the system's base pairs can be robustly inferred, with A/A and T/C DNA/RNA mismatches giving the strongest signals. These results point to the feasibility of rapid base pair inference in larger and more complex nucleic acid systems with unknown structure.

Functional elucidation of a key contact between tRNA and the large ribosomal subunit rRNA during decoding

The selection of cognate tRNAs during translation is specified by a kinetic discrimination mechanism driven by distinct structural states of the ribosome. While the biochemical steps that drive the tRNA selection process have been carefully documented, it remains unclear how recognition of matched codon:anticodon helices in the small subunit facilitate global rearrangements in the ribosome complex that efficiently promote tRNA decoding. Here we use an in vitro selection approach to isolate tRNA(Trp) miscoding variants that exhibit a globally perturbed tRNA tertiary structure. Interestingly, the most substantial distortions are positioned in the elbow region of the tRNA that closely approaches helix 69 (H69) of the large ribosomal subunit. The importance of these specific interactions to tRNA selection is underscored by our kinetic analysis of both tRNA and rRNA variants that perturb the integrity of this interaction.

Enhanced functional potential of nucleic acid aptamer libraries patterned to increase secondary structure

The in vitro selection of nucleic acid libraries has driven the discovery of RNA and DNA receptors (aptamers) and catalysts with tailor-made functional properties. Functional nucleic acids emerging from selections have been observed to possess an unusually high degree of secondary structure. In this study, we experimentally examined the relationship between the degree of secondary structure in a nucleic acid library and its ability to yield aptamers that bind protein targets. We designed a patterned nucleic acid library (denoted R*Y*) to enhance the formation of stem-loop structures without imposing any specific sequence or secondary structural requirement. This patterned library was predicted computationally to contain a significantly higher average folding energy compared to a standard, unpatterned N(60) library of the same length. We performed three different iterated selections for protein binding using patterned and unpatterned libraries competing in the same solution. In all three cases, the patterned R*Y* library was enriched relative to the unpatterned library over the course of the 9- to 10-round selection. Characterization of individual aptamer clones emerging from the three selections revealed that the highest affinity aptamer assayed arose from the patterned library for two protein targets, while in the third case, the highest affinity aptamers from the patterned and random libraries exhibited comparable affinity. We identified the binding motif requirements for the most active aptamers generated against two of the targets. The two binding motifs are 3.4- and 27-fold more likely to occur in the R*Y* library than in the N(60) library. Collectively, our findings suggest that researchers performing selections for nucleic acid aptamers and catalysts should consider patterned libraries rather than commonly used N(m) libraries to increase both the likelihood of isolating functional molecules and the potential activities of the resulting molecules.

Comparisons between chemical mapping and binding to isoenergetic oligonucleotide microarrays reveal unexpected patterns of binding to the Bacillus subtilis RNase P RNA specificity domain

Microarrays with isoenergetic pentamer and hexamer 2'-O-methyl oligonucleotide probes with LNA (locked nucleic acid) and 2,6-diaminopurine substitutions were used to probe the binding sites on the RNase P RNA specificity domain of Bacillus subtilis. Unexpected binding patterns were revealed. Because of their enhanced binding free energies, isoenergetic probes can break short duplexes, merge adjacent loops, and/or induce refolding. This suggests new approaches to the rational design of short oligonucleotide therapeutics but limits the utility of microarrays for providing constraints for RNA structure determination. The microarray results are compared to results from chemical mapping experiments, which do provide constraints. Results from both types of experiments indicate that the RNase P RNA folds similarly in 1 M Na(+) and 10 mM Mg(2+).

Disease-associated mutations that alter the RNA structural ensemble

Genome-wide association studies (GWAS) often identify disease-associated mutations in intergenic and non-coding regions of the genome. Given the high percentage of the human genome that is transcribed, we postulate that for some observed associations the disease phenotype is caused by a structural rearrangement in a regulatory region of the RNA transcript. To identify such mutations, we have performed a genome-wide analysis of all known disease-associated Single Nucleotide Polymorphisms (SNPs) from the Human Gene Mutation Database (HGMD) that map to the untranslated regions (UTRs) of a gene. Rather than using minimum free energy approaches (e.g. mFold), we use a partition function calculation that takes into consideration the ensemble of possible RNA conformations for a given sequence. We identified in the human genome disease-associated SNPs that significantly alter the global conformation of the UTR to which they map. For six disease-states (Hyperferritinemia Cataract Syndrome, β-Thalassemia, Cartilage-Hair Hypoplasia, Retinoblastoma, Chronic Obstructive Pulmonary Disease (COPD), and Hypertension), we identified multiple SNPs in UTRs that alter the mRNA structural ensemble of the associated genes. Using a Boltzmann sampling procedure for sub-optimal RNA structures, we are able to characterize and visualize the nature of the conformational changes induced by the disease-associated mutations in the structural ensemble. We observe in several cases (specifically the 5′ UTRs of FTL and RB1) SNP–induced conformational changes analogous to those observed in bacterial regulatory Riboswitches when specific ligands bind. We propose that the UTR and SNP combinations we identify constitute a “RiboSNitch,” that is a regulatory RNA in which a specific SNP has a structural consequence that results in a disease phenotype. Our SNPfold algorithm can help identify RiboSNitches by leveraging GWAS data and an analysis of the mRNA structural ensemble.

Genome-wide association studies identify mutations in the human genome that correlate with a particular disease. It is common to find mutations associated with disease in the non-coding region of the genome. These non-coding mutations are more difficult to interpret at a molecular level, because they do not affect the protein sequence. In this study, we analyze disease-associated mutations in non-coding regions of our genome in the context of their structural effect on the message of genetic information in our cells, Ribonucleic Acid (RNA). We focus in particular on the regulatory parts of our genes known as untranslated regions. We find that certain disease-associated mutations in these regulatory untranslated regions have a significant effect on the structure of the RNA message. We call these elements “RiboSNitches,” because they act like switches turning on and off genes, but are caused by Single Nucleotide Polymorphisms (SNPs), which are single point mutations in our genome. The RiboSNitches we identify are potentially a new class of pharmaceutical targets, as it is possible to change the structure of RNA with small drug-like molecules.

On the role of the cis Hoogsteen:sugar-edge family of base pairs in platforms and triplets-quantum chemical insights into RNA structural biology

Base pairs belonging to the cis Hoogsteen:sugar-edge (H:S) family play important structural roles in folded RNA molecules. Several of these are present in internal loops, where they are involved in interactions leading to planar dinucleotide platforms which stabilize higher order structures such as base triplets and quartets. We report results of analysis of 30 representative examples spanning 16 possible base pair combinations, with several of them showing multimodality of base pairing geometry. The geometries of 23 of these base pairs were modeled directly from coordinates extracted from RNA crystal structures. The other seven were predicted structures which were modeled on the basis of observed isosteric analogues. After appropriate satisfaction of residual valencies, these structures were relaxed using the B3LYP/6-31G(d,p) method and interaction energies were derived at the RIMP2/aug-cc-pVDZ level of theory. The geometries for each of the studied base pairs have been characterized in terms of the number and nature of H-bonds, rmsd values observed on optimization, base pair geometrical parameters, and sugar pucker analysis. In addition to its evaluation, the nature of intermolecular interaction in these complexes was also analyzed using Morokuma decomposition. The gas phase interaction energies range between -5.2 and -20.6 kcal/mol and, in contrast to the H:S trans base pairs, show enhanced relative importance of the electron correlation component, indicative of the greater role of dispersion energy in stabilization of these base pairs. The rich variety of hydrogen bonding pattern, involving the flexible sugar edge, appears to hold the key to several features of structural motifs, such as planarity and propensity to participate in triplets, observed in this family of base pairs. This work explores these aspects by integrating database analysis, and detailed base pairing geometry analysis at the atomistic level, with ab initio computation of interaction energies. The study, involving alternative classification of base pairs and triplets, provides insights into intrinsic properties of these base pairs and their possible structural and functional roles.

Quantifying the relationship between sequence and three-dimensional structure conservation in RNA

Background

In recent years, the number of available RNA structures has rapidly grown reflecting the increased interest on RNA biology. Similarly to the studies carried out two decades ago for proteins, which gave the fundamental grounds for developing comparative protein structure prediction methods, we are now able to quantify the relationship between sequence and structure conservation in RNA.

Results

Here we introduce an all-against-all sequence- and three-dimensional (3D) structure-based comparison of a representative set of RNA structures, which have allowed us to quantitatively confirm that: (i) there is a measurable relationship between sequence and structure conservation that weakens for alignments resulting in below 60% sequence identity, (ii) evolution tends to conserve more RNA structure than sequence, and (iii) there is a twilight zone for RNA homology detection.

Discussion

The computational analysis here presented quantitatively describes the relationship between sequence and structure for RNA molecules and defines a twilight zone region for detecting RNA homology. Our work could represent the theoretical basis and limitations for future developments in comparative RNA 3D structure prediction.

Calculation of the standard binding free energy of sparsomycin to the ribosomal peptidyl-transferase P-site using molecular dynamics simulations with restraining potentials

The standard (absolute) binding free energy of the antibiotic sparsomycin with the 50S bacteria ribosomal subunit is calculated using molecular dynamics (MD) free energy perturbation (FEP) simulations with restraining potentials developed by Wang et al. [Biophys. J. 91:2798-2814 (2006)]. In the simulation protocol, restraining potentials are activated for the orientational and translational movements of the ligand relative to the binding site when it is decoupled from the binding pocket, and then released once the ligand fully interacts with the rest of the system. A reduced system is simulated to decrease the computational cost of the FEP/MD calculations and the effects of the surrounding atoms are incorporated using the generalized solvent boundary potential (GSBP) method. The loss of conformational freedom of the ligand upon binding is characterized using the potential of mean force (PMF) as a function of the root-mean-square deviation (RMSD) relative to the bound conformation. The number of water molecules in the binding pocket is allowed to fluctuate dynamically in response to the ligand during the calculations by combining FEP/MD with grand canonical Monte Carlo (GCMC) simulations. The calculated binding free energy is about -6 kcal/mol, which is in reasonable agreement with the experimental value. The information gleaned from this study provides new insight on the recognition of ribosome by sparsomycin and highlights the challenges in calculations of absolute binding free energies in these systems.

The Ribosome Comes Alive

This essay is a reflection on the ways the X-ray structures of the ribosome are helping in the interpretation of cryogenic electron microscopy (cryo-EM) density maps showing the translating ribosome in motion. Through advances in classification methods, cryo-EM and single-particle reconstruction methods have recently evolved to the point where they can yield an array of structures from a single sample ("story in a sample"), providing snapshots of an entire subprocess of translation, such as translocation or decoding.

On the role of Hoogsteen:Hoogsteen interactions in RNA: ab initio investigations of structures and energies

We use a combination of database analysis and quantum chemical studies to investigate the role of cis and trans Hoogsteen:Hoogsteen (H:H) base pairs and associated higher-order structures in RNA. We add three new examples to the list of previously identified base-pair combinations belonging to these families and, in addition to contextual classification and characterization of their structural and energetic features, we compare their interbase interaction energies and propensities toward participation in triplets and quartets. We find that some base pairs, which are nonplanar in their isolated minimum energy geometries, attain planarity and stability upon triplet formation. A:A H:H trans is the most frequent H:H combination in RNA structures. This base pair occurs at many distinct positions in known rRNA structures, where it helps in the interaction of ribosomal domains in the 50S subunit. It is also present as a part of tertiary interaction in tRNA structures. Although quantum chemical studies suggest an intrinsically nonplanar geometry for this base pair in isolated form, it has the tendency to attain planar geometry in RNA crystal structures by forming higher-order tertiary interactions or in the presence of additional base-phosphate interactions. The tendency of this base pair to form such additional interactions may be helpful in bringing together different segments of RNA, thus making it suitable for the role of facilitator for RNA folding. This also explains the high occurrence frequency of this base pair among all H:H interactions.

The transition from noncoded to coded protein synthesis: did coding mRNAs arise from stability-enhancing binding partners to tRNA?

BACKGROUND

Understanding the origin of protein synthesis has been notoriously difficult. We have taken as a starting premise Wolf and Koonin's view that "evolution of the translation system is envisaged to occur in a compartmentalized ensemble of replicating, co-selected RNA segments, i.e., in an RNA world containing ribozymes with versatile activities".

PRESENTATION OF THE HYPOTHESIS

We propose that coded protein synthesis arose from a noncoded process in an RNA world as a natural consequence of the accumulation of a range of early tRNAs and their serendipitous RNA binding partners. We propose that, initially, RNA molecules with 3' CCA termini that could be aminoacylated by ribozymes, together with an ancestral peptidyl transferase ribozyme, produced small peptides with random or repetitive sequences. Our concept is that the first tRNA arose in this context from the ligation of two RNA hairpins and could be similarly aminoacylated at its 3' end to become a substrate for peptidyl transfer catalyzed by the ancestral ribozyme. Within this RNA world we hypothesize that proto-mRNAs appeared first simply as serendipitous binding partners, forming complementary base pair interactions with the anticodon loops of tRNA pairs. Initially this may have enhanced stability of the paired tRNA molecules so they were held together in close proximity, better positioning the 3' CCA termini for peptidyl transfer and enhancing the rate of peptide synthesis. If there were a selective advantage for the ensemble through the peptide products synthesized, it would provide a natural pathway for the evolution of a coding system with the expansion of a cohort of different tRNAs and their binding partners. The whole process could have occurred quite unremarkably for such a profound acquisition.

TESTING THE HYPOTHESIS

It should be possible to test the different parts of our model using the isolated contemporary 50S ribosomal subunit initially, and then with RNAs transcribed in vitro together with a minimal set of ribosomal proteins that are required today to support protein synthesis.

IMPLICATIONS OF THE HYPOTHESIS

This model proposes that genetic coding arose de novo from complementary base pair interactions between tRNAs and single-stranded RNAs present in the immediate environment.

REVIEWERS

This article was reviewed by Eugene Koonin, Rob Knight and Berthold Kastner (nominated by Laura Landweber).

RNA targeting therapeutics: molecular mechanisms of antisense oligonucleotides as a therapeutic platform

Dramatic advances in understanding of the roles RNA plays in normal health and disease have greatly expanded over the past 10 years and have made it clear that scientists are only beginning to comprehend the biology of RNAs. It is likely that RNA will become an increasingly important target for therapeutic intervention; therefore, it is important to develop strategies for therapeutically modulating RNA function. Antisense oligonucleotides are perhaps the most direct therapeutic strategy to approach RNA. Antisense oligonucleotides are designed to bind to the target RNA by well-characterized Watson-Crick base pairing, and once bound to the target RNA, modulate its function through a variety of postbinding events. This review focuses on the molecular mechanisms by which antisense oligonucleotides can be designed to modulate RNA function in mammalian cells and how synthetic oligonucleotides behave in the body.

The RNA backbone plays a crucial role in mediating the intrinsic stability of the GpU dinucleotide platform and the GpUpA/GpA miniduplex

The side-by-side interactions of nucleobases contribute to the organization of RNA, forming the planar building blocks of helices and mediating chain folding. Dinucleotide platforms, formed by side-by-side pairing of adjacent bases, frequently anchor helices against loops. Surprisingly, GpU steps account for over half of the dinucleotide platforms observed in RNA-containing structures. Why GpU should stand out from other dinucleotides in this respect is not clear from the single well-characterized H-bond found between the guanine N2 and the uracil O4 groups. Here, we describe how an RNA-specific H-bond between O2′(G) and O2P(U) adds to the stability of the GpU platform. Moreover, we show how this pair of oxygen atoms forms an out-of-plane backbone ‘edge’ that is specifically recognized by a non-adjacent guanine in over 90% of the cases, leading to the formation of an asymmetric miniduplex consisting of ‘complementary’ GpUpA and GpA subunits. Together, these five nucleotides constitute the conserved core of the well-known loop-E motif. The backbone-mediated intrinsic stabilities of the GpU dinucleotide platform and the GpUpA/GpA miniduplex plausibly underlie observed evolutionary constraints on base identity. We propose that they may also provide a reason for the extreme conservation of GpU observed at most 5′-splice sites.

Arrangement of 3D structural motifs in ribosomal RNA

Structural 3D motifs in RNA play an important role in the RNA stability and function. Previous studies have focused on the characterization and discovery of 3D motifs in RNA secondary and tertiary structures. However, statistical analyses of the distribution of 3D motifs along the RNA appear to be lacking. Herein, we present a novel strategy for evaluating the distribution of 3D motifs along the RNA chain and those motifs whose distributions are significantly non-random are identified. By applying it to the X-ray structure of the large ribosomal subunit from Haloarcula marismortui, helical motifs were found to cluster together along the chain and in the 3D structure, whereas the known tetraloops tend to be sequentially and spatially dispersed. That the distribution of key structural motifs such as tetraloops differ significantly from a random one suggests that our method could also be used to detect novel 3D motifs of any size in sufficiently long/large RNA structures. The motif distribution type can help in the prediction and design of 3D structures of large RNA molecules.