Structural features of a six-nucleotide RNA hairpin loop found in ribosomal RNA

Secondary structure prediction for RNA sequences including N6-methyladenosine

There is increasing interest in the roles of covalently modified nucleotides in RNA. There has been, however, an inability to account for modifications in secondary structure prediction because of a lack of software and thermodynamic parameters. We report the solution for these issues for N⁶-methyladenosine (m⁶A), allowing secondary structure prediction for an alphabet of A, C, G, U, and m⁶A. The RNAstructure software now works with user-defined nucleotide alphabets of any size. We also report a set of nearest neighbor parameters for helices and loops containing m⁶A, using experiments. Interestingly, N⁶-methylation decreases folding stability for adenosines in the middle of a helix, has little effect on folding stability for adenosines at the ends of helices, and increases folding stability for unpaired adenosines stacked on a helix. We demonstrate predictions for an N⁶-methylation-activated protein recognition site from MALAT1 and human transcriptome-wide effects of N⁶-methylation on the probability of adenosine being buried in a helix.

RNA folding free energy nearest neighbor parameters were determined for sequences with the nucleotide m⁶A. The RNAstructure software package can accommodate modified nucleotides, enabling secondary structure prediction of sequences with m⁶A.

bpRNA: large-scale automated annotation and analysis of RNA secondary structure

While RNA secondary structure prediction from sequence data has made remarkable progress, there is a need for improved strategies for annotating the features of RNA secondary structures. Here, we present bpRNA, a novel annotation tool capable of parsing RNA structures, including complex pseudoknot-containing RNAs, to yield an objective, precise, compact, unambiguous, easily-interpretable description of all loops, stems, and pseudoknots, along with the positions, sequence, and flanking base pairs of each such structural feature. We also introduce several new informative representations of RNA structure types to improve structure visualization and interpretation. We have further used bpRNA to generate a web-accessible meta-database, ‘bpRNA-1m’, of over 100 000 single-molecule, known secondary structures; this is both more fully and accurately annotated and over 20-times larger than existing databases. We use a subset of the database with highly similar (≥90% identical) sequences filtered out to report on statistical trends in sequence, flanking base pairs, and length. Both the bpRNA method and the bpRNA-1m database will be valuable resources both for specific analysis of individual RNA molecules and large-scale analyses such as are useful for updating RNA energy parameters for computational thermodynamic predictions, improving machine learning models for structure prediction, and for benchmarking structure-prediction algorithms.

Chemically Accurate Relative Folding Stability of RNA Hairpins from Molecular Simulations

To benchmark RNA force fields, we compared the folding stabilities of three 12-nucleotide hairpin stem loops estimated by simulation to stabilities determined by experiment. We used umbrella sampling and a reaction coordinate of end-to-end (5' to 3' hydroxyl oxygen) distance to estimate the free energy change of the transition from the native conformation to a fully extended conformation with no hydrogen bonds between non-neighboring bases. Each simulation was performed four times using the AMBER FF99+bsc0+χ_OL3 force field, and each window, spaced at 1 Å intervals, was sampled for 1 μs, for a total of 552 μs of simulation. We compared differences in the simulated free energy changes to analogous differences in free energies from optical melting experiments using thermodynamic cycles where the free energy change between stretched and random coil sequences is assumed to be sequence-independent. The differences between experimental and simulated ΔΔ G° are, on average, 0.98 ± 0.66 kcal/mol, which is chemically accurate and suggests that analogous simulations could be used predictively. We also report a novel method to identify where replica free energies diverge along a reaction coordinate, thus indicating where additional sampling would most improve convergence. We conclude by discussing methods to more economically perform these simulations.

Nucleobases Undergo Dynamic Rearrangements during RNA Tertiary Folding

The tertiary structure of the GTPase center (GAC) of 23S ribosomal RNA (rRNA) as seen in cocrystals is extremely compact. It is stabilized by long-range hydrogen bonds and nucleobase stacking and by a triloop that forms within its three-way junction. Its folding pathway from secondary structure to tertiary structure has not been previously observed, but it was shown to require Mg²⁺ ions in equilibrium experiments. The fluorescent nucleotide 2-aminopurine was substituted at selected sites within the 60-nt GAC. Fluorescence intensity changes upon addition of MgCl₂ were monitored over a time-course from 1ms to 100s as the RNA folds. The folding pathway is revealed here to be hierarchical through several intermediates. Observation of the nucleobases during folding provides a new perspective on the process and the pathway, revealing the dynamics of nucleobase conformational exchange during the folding transitions.

Formation of Tertiary Interactions during rRNA GTPase Center Folding

The 60-nt GTPase center (GAC) of 23S rRNA has a phylogenetically conserved secondary structure with two hairpin loops and a 3-way junction. It folds into an intricate tertiary structure upon addition of Mg(2+) ions, which is stabilized by the L11 protein in cocrystal structures. Here, we monitor the kinetics of its tertiary folding and Mg(2+)-dependent intermediate states by observing selected nucleobases that contribute specific interactions to the GAC tertiary structure in the cocrystals. The fluorescent nucleobase 2-aminopurine replaced three individual adenines, two of which make long-range stacking interactions and one that also forms hydrogen bonds. Each site reveals a unique response to Mg(2+) addition and temperature, reflecting its environmental change from secondary to tertiary structure. Stopped-flow fluorescence experiments revealed that kinetics of tertiary structure formation upon addition of MgCl2 are also site specific, with local conformational changes occurring from 5 ms to 4s and with global folding from 1 to 5s. Site-specific substitution with (15)N-nucleobases allowed observation of stable hydrogen bond formation by NMR experiments. Equilibrium titration experiments indicate that a stable folding intermediate is present at stoichiometric concentrations of Mg(2+) and suggest that there are two initial sites of Mg(2+) ion association.

NMR characterization of an oligonucleotide model of the miR-21 pre-element

We have used NMR spectroscopy to characterize an oligonucleotide stem loop structure based on the pre-element of an oncogenic microRNA, miR-21. This predicted stem-loop structure is cleaved from the precursor of miR-21 (pre-miR-21) by the nuclease Dicer. It is also a critical feature recognized by the protein complex that converts the primary transcript (pri-miR-21) into the pre-miRNA. The secondary structure of the native sequence is poorly defined by NMR due to rapid exchange of imino protons with solvent; however, replacement of two adjacent putative G•U base pairs with G•C base pairs retains the conformation of the hairpin observed by chemical probing and stabilizes it sufficiently to observe most of the imino proton resonances of the molecule. The observed resonances are consistent with the predicted secondary structure. In addition, a peak due to a loop uridine suggests an interaction between it and a bulged uridine in the stem. Assignment of non-exchangeable proton resonances and characterization of NOEs and coupling constants allows inference of the following features of the structure: extrahelicity of a bulged adenosine, deviation from A-form geometry in a base-paired stem, and consecutive stacking of the adenosines in the 5′ side of the loop, the guanosine of the closing base pair, and a cross-strand adenosine. Modeling of the structure by restrained molecular dynamics suggests a basis for the interaction between the loop uridine, the bulged uridine in the stem, and an A•U base pair in the stem.

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.

The Amber ff99 Force Field Predicts Relative Free Energy Changes for RNA Helix Formation

The ability of the Amber ff99 force field to predict relative free energies of RNA helix formation was investigated. The test systems were three hexaloop RNA hairpins with identical loops and varying stems. The potential of mean force of stretching the hairpins from the native state to an extended conformation was calculated with umbrella sampling. Because the hairpins have identical loop sequence, the differences in free energy changes are only from the stem composition. The Amber ff99 force field was able to correctly predict the order of stabilities of the hairpins, although the magnitude of the free energy change is larger than that determined by optical melting experiments. The two measurements cannot be compared directly because the unfolded state in the optical melting experiments is a random coil, while the end state in the umbrella sampling simulations was an elongated chain. The calculations can be compared to reference data by using a thermodynamic cycle. By applying the thermodynamic cycle to the transitions between the hairpins using simulations and nearest neighbor data, agreement was found to be within the sampling error of simulations, thus demonstrating that ff99 force field is able to accurately predict relative free energies of RNA helix formation.

Nuclear magnetic resonance structure of the prohead RNA E-loop hairpin

The Bacillus subtilis phage phi29 packaging motor requires prohead RNA for genome encapsidation. The nuclear magnetic resonance structure of the prohead RNA E-loop hairpin, r(5'AUUGAGUU), is presented and compared to predictions from MC-SYM. The prohead RNA E-loop hairpins contain sequences similar to rRNA hairpins. Comparison of predicted and experimentally determined prohead and ribosomal hairpin structures reveals that sequence similarity is a stronger determinant of hairpin structural similarity than grouping similar types of RNA. All the hairpins contain a U-turn motif but differ in the first noncanonical pair and backbone orientation. These structures provide benchmarks for further improvements in RNA structure predictions from sequence.

Correlation of RNA secondary structure statistics with thermodynamic stability and applications to folding

The accurate prediction of the secondary and tertiary structure of an RNA with different folding algorithms is dependent on several factors, including the energy functions. However, an RNA higher-order structure cannot be predicted accurately from its sequence based on a limited set of energy parameters. The inter- and intramolecular forces between this RNA and other small molecules and macromolecules, in addition to other factors in the cell such as pH, ionic strength, and temperature, influence the complex dynamics associated with transition of a single stranded RNA to its secondary and tertiary structure. Since all of the factors that affect the formation of an RNAs 3D structure cannot be determined experimentally, statistically derived potential energy has been used in the prediction of protein structure. In the current work, we evaluate the statistical free energy of various secondary structure motifs, including base-pair stacks, hairpin loops, and internal loops, using their statistical frequency obtained from the comparative analysis of more than 50,000 RNA sequences stored in the RNA Comparative Analysis Database (rCAD) at the Comparative RNA Web (CRW) Site. Statistical energy was computed from the structural statistics for several datasets. While the statistical energy for a base-pair stack correlates with experimentally derived free energy values, suggesting a Boltzmann-like distribution, variation is observed between different molecules and their location on the phylogenetic tree of life. Our statistical energy values calculated for several structural elements were utilized in the Mfold RNA-folding algorithm. The combined statistical energy values for base-pair stacks, hairpins and internal loop flanks result in a significant improvement in the accuracy of secondary structure prediction; the hairpin flanks contribute the most.

Bacteriophage P22 antitermination boxB sequence requirements are complex and overlap with those of lambda

Transcription antitermination in phages lambda and P22 uses N proteins that bind to similar boxB RNA hairpins in regulated transcripts. In contrast to the lambda N-boxB interaction, the P22 N-boxB interaction has not been extensively studied. A nuclear magnetic resonance structure of the P22 N peptide boxB(left) complex and limited mutagenesis have been reported but do not reveal a consensus sequence for boxB. We have used a plasmid-based antitermination system to screen boxBs with random loops and to test boxB mutants. We find that P22 N requires boxB to have a GNRA-like loop with no simple requirements on the remaining sequences in the loop or stem. U:A or A:U base pairs are strongly preferred adjacent to the loop and appear to modulate N binding in cooperation with the loop and distal stem. A few GNRA-like hexaloops have moderate activity. Some boxB mutants bind P22 and lambda N, indicating that the requirements imposed on boxB by P22 N overlap those imposed by lambda N. Point mutations can dramatically alter boxB specificity between P22 and lambda N. A boxB specific for P22 N can be mutated to lambda N specificity by a series of single mutations via a bifunctional intermediate, as predicted by neutral theories of evolution.

Base pairing within the psi32,psi39-modified anticodon arm of Escherichia coli tRNA(Phe)

The base-base hydrogen bond interactions of the psi32,psi39-modified anticodon arm of Escherichia coli tRNAPhe have been investigated using heteronuclear NMR spectroscopy. psi32 and psi39 were enzymatically introduced into a [13C,15N]-isotopically enriched RNA sequence corresponding to the tRNAPhe anticodon arm. Both the psi32-A38 and A31-psi39 nucleotide pairs form Watson-Crick base pairing schemes and the anticodon nucleotides adopt a triloop conformation. Similar effects were observed previously with D2-isopentenyl modification of the A37 N6 that also is native to the tRNAPhe anticodon arm. These results demonstrate that the individual modifications are not sufficient to produce the 32-38 bifurcated hydrogen bond or the U-turn motifs that are observed in crystal structures of tRNAs and tRNA-protein complexes. Thus the formation of these conserved structural features in solution likely require the synergistic interaction of multiple modifications.

The naturally occurring N6-threonyl adenine in anticodon loop of Schizosaccharomyces pombe tRNAi causes formation of a unique U-turn motif

Modified nucleosides play an important role in structure and function of tRNA. We have determined the solution structure of the anticodon stem-loop (ASL) of initiator tRNA of Schizosaccharomyces pombe. The incorporation of N6-threonylcarbamoyladenosine at the position 3' to the anticodon triplet (t6A37) results in the formation of a U-turn motif and enhances stacking interactions within the loop and stem regions (i.e. between A35 and t6A37) by bulging out U36. This conformation was not observed in a crystal structure of tRNAi including the same modification in its anticodon loop, nor in the solution structure of the unmodified ASL. A t6A modification also occurs in the well studied anti-stem-loop of lys-tRNA(UUU). A comparison of this stem-loop with our structure demonstrates different effects of the modification depending on the loop sequence.

High-resolution structure of a picornaviral internal cis-acting RNA replication element (cre)

Picornaviruses constitute a medically important family of RNA viruses in which genome replication critically depends on a small RNA element, the cis-acting replication element (cre), that templates 3D(pol) polymerase-catalyzed uridylylation of the protein primer for RNA synthesis, VPg. We report the solution structure of the 33-nt cre of human rhinovirus 14 under solution conditions optimal for uridylylation in vitro. The cre adopts a stem-loop conformation with an extended duplex stem supporting a novel 14-nt loop that derives stability from base-stacking interactions. Base-pair interactions are absent within the loop, and base substitutions within the loop that favor such interactions are detrimental to viral RNA replication. Conserved adenosines in the 5' loop sequence that participate in a slide-back mechanism of VPg-pUpU synthesis are oriented to the inside of the loop but are available for base templating during uridylation. The structure explains why substitutions of the 3' loop nucleotides have little impact on conformation of the critical 5' loop bases and accounts for wide variation in the sequences of cres from different enteroviruses and rhinoviruses.

Structural mimicry in the phage phi21 N peptide-boxB RNA complex

We determined the solution structure of a 22-amino-acid peptide from the amino-terminal domain of the bacteriophage phi21 N protein in complex with its cognate 24-mer boxB RNA hairpin using heteronuclear magnetic resonance spectroscopy. The N peptide binds as an alpha-helix and interacts predominately with the major groove side of the 5' half of the boxB RNA stem-loop. This binding interface is defined by surface complementarity of polar and nonpolar interactions, and little sequence-specific recognition. The phi21 boxB loop (CUAACC) has hydrogen bond and backbone torsions typical of the "U-turn" motif, as well as base stacking of the last 4 nt, and a hydrogen bonded C:C pair closing the loop. The exposed face of the phi21 boxB loop, in complex with the N peptide, is strikingly similar to the GNRA tetraloop-like folds of the related lambda and P22 bacteriophage N peptide-boxB RNA complexes. The N peptide-boxB complexes of the various phage, while individually distinct, provide similar structural features for interactions with the Escherichia coli host factors to enable antitermination.

NMR structure of a ribosomal RNA hairpin containing a conserved CUCAA pentaloop

The structure of a 23 nt RNA sequence, rGGACCCGGGCUCAACCUGGGUCC, was elucidated using homonuclear NMR, distance geometry and restrained molecular dynamics. This RNA is analogous to residues 612-628 of the Escherichia coli 16S rRNA. The structure of the RNA reveals the presence of a pentaloop closed by a duplex stem in typical A-form conformation. The loop does not form a U-turn motif, as previously predicted. A non-planar A.C.A triple base interaction (hydrogen bonds A13 NH6-C10 O2 and C10 N3-A14 NH6) stabilizing the loop structure is inferred from structure calculations. The CUCAA loop structure is asymmetrical, characterized by a reversal of the phosphodiester backbone at the UC step (hydrogen bond C12 NH4-C10 O2') and 3'-stacking within the CAA segment. Loop base U11 is oriented towards the major groove and the consecutive adenosines on the 3'-end of the loop are well stacked, exposing their reactive functional groups in the minor groove defined by the duplex stem. The solution structure of the loop resembles that seen in the 3.3 A X-ray structure of the entire 30S subunit, where the analogous loop interacts with a ribosomal protein and a receptor RNA helix.

Structure and function of the conserved 690 hairpin in Escherichia coli 16 S ribosomal RNA. III. Functional analysis of the 690 loop

An instant-evolution experiment was performed on the eight nucleotides comprising the loop region of the 690 hairpin in Escherichia coli 16 S ribosomal RNA. Positions 690 to 697 were randomly mutated and 101 unique functional mutants were isolated, sequenced and analyzed for function in vivo. Non-random nucleotide distributions were observed at each of the mutated positions except 693 and 694. Nucleotide identity significantly affected ribosome function at positions 690, 695, 696 and 697. Pyrimidines were absent at position 696 in the instant-evolution pool as were C at position 691 and G at position 697. A highly significant covariation was observed between nucleotides 690 and 697. No functional double mutants at positions 691 and 696 were obtained from the instant-evolution pool. In our NMR structure of the 690 loop, both the G690.U697 and G691.A696 form sheared hydrogen-bonded mismatches. To further examine the functional constraints between these paired nucleotides, one set of site-directed mutations was constructed at positions 690:697 and another set was constructed at positions 691:696. Functional analysis of the site-directed mutants is consistent with our instant-evolution findings and revealed constraints on the placement of specific functional groups observed in the NMR structure. Ten instant-evolution mutants were isolated that are more functional than the wild-type. Hyperactivity in these mutants correlates with a single mutation at position 693.

Structure and function of the conserved 690 hairpin in Escherichia coli 16 S ribosomal RNA. II. NMR solution structure

The solution structure of the conserved 690 hairpin from Escherichia coli 16 S rRNA was determined by NMR spectroscopy. The 690 loop is located at the surface of the 30 S subunit in the platform region and has been implicated in interactions with P-site bound tRNA, E-site mRNA, S11 binding, IF3 binding, and in RNA-RNA interactions with the 790 loop of 16 S rRNA and domain IV of 23 S rRNA. The structure reveals a novel sheared type G690.U697 base-pair with a single hydrogen bond from the G690 amino to U697-04. G691 and A696 also form a sheared pair and U692 forms a U-turn with an H-bond to the A695 non-bridging phosphate oxygen. The sheared pairs and U-turn result in the continuous single-stranded stacking of five residues from 6693 to U697 with their Watson-Crick functional groups exposed in the minor groove. The overall fold of the 690 hairpin is similar to the anticodon loop of tRNA. The structure provides an explanation for chemical protection patterns in the loop upon interaction with tRNA, the 50 S subunit, and S11. In vivo genetic studies demonstrate the functional importance of the motifs observed in the solution structure of the 690 hairpin.

A story: unpaired adenosine bases in ribosomal RNAs

In 1985 an analysis of the Escherichia coli 16 S rRNA covariation-based structure model revealed a strong bias for unpaired adenosines. The same analysis revealed that the majority of the G, C, and U bases were paired. These biases are (now) consistent with the high percentage of unpaired adenosine nucleotides in several structure motifs. An analysis of a larger set of bacterial comparative 16 S and 23 S rRNA structure models has substantiated this initial finding and revealed new biases in the distribution of adenosine nucleotides in loop regions. The majority of the adenosine nucleotides are unpaired, while the majority of the G, C, and U bases are paired in the covariation-based structure model. The unpaired adenosine nucleotides predominate in the middle and at the 3' end of loops, and are the second most frequent nucleotide type at the 5' end of loops (G is the most common nucleotide). There are additional biases for unpaired adenosine nucleotides at the 3' end of loops and adjacent to a G at the 5' end of the helix. The most prevalent consecutive nucleotides are GG, GA, AG, and AA. A total of 70 % of the GG sequences are within helices, while more than 70 % of the AA sequences are unpaired. Nearly 50 % of the GA sequences are unpaired, and approximately one-third of the AG sequences are within helices while another third are at the 3' loop.5' helix junction. Unpaired positions with an adenosine nucleotide in more than 50 % of the sequences at the 3' end of 16 S and 23 S rRNA loops were identified and arranged into the A-motif categories XAZ, AAZ, XAG, AAG, and AAG:U, where G or Z is paired, G:U is a base-pair, and X is not an A and Z is not a G in more than 50 % of the sequences. These sequence motifs were associated with several structural motifs, such as adenosine platforms, E and E-like loops, A:A and A:G pairings at the end of helices, G:A tandem base-pairs, GNRA tetraloop hairpins, and U-turns.

Predicting U-turns in ribosomal RNA with comparative sequence analysis

The U-turn is a well-known RNA motif characterized by a sharp reversal of the RNA backbone following a single-stranded uridine base. In experimentally determined U-turn motifs, the nucleotides 3' to the turn are frequently involved in tertiary interactions, rendering this motif particularly attractive in RNA modeling and functional studies. The U-turn signature is composed of an UNR sequence pattern flanked by a Y:Y, Y:A (Y=pyrimidine) or G:A base juxtaposition. We have identified 33 potential UNR-type U-turns and 25 related GNRA-type U-turns in a large set of aligned 16 S and 23 S rRNA sequences. U-turn candidates occur in hairpin loops (34 times) as well as in internal and multi-stem loops (24 times). These are classified into ten families based on loop type, sequence pattern (UNR or GNRA) and the nature of the closing base juxtaposition. In 13 cases, the bases on the 3' side of the turn, or on the immediate 5' side, are involved in tertiary covariations, making these sites strong candidates for tertiary interactions.

Structural motifs in RNA

An RNA motif is a discrete sequence or combination of base juxtapositions found in naturally occurring RNAs in unexpectedly high abundance. Because all the motifs examined so far have three-dimensional structures independent of the context in which they are embedded, they are important components of the "kit" of structural elements from which RNAs are constructed. This review discusses the structures of the motifs that have been identified so far and speculates on the importance of their role in determining RNA conformation and their evolutionary origin.

Aminoglycoside-RNA interactions

The structural and physico-chemical parameters promoting the binding of aminoglycosides to RNAs are becoming clear. The strength of the interaction is dominated by electrostatics, with the positively charged aminoglycosides displacing metal ions. Although aminoglycosides inhibit most known ribozymes, aminoglycosides or polyamines are able to catalyze specific RNA cleavage in the absence of metal ions.

Determining the influence of structure on hybridization using oligonucleotide arrays

We have studied the effects of structure on nucleic acid heteroduplex formation by analyzing hybridization of tRNAphe to a complete set of complementary oligonucleotides, ranging from single nucleotides to dodecanucleotides. The analysis points to features in tRNA that determine heteroduplex yield. All heteroduplexes that give high yield include both double-stranded stems as well as single-stranded regions. Bases in the single-stranded regions are stacked onto the stems, and heteroduplexes terminate at potential interfaces for coaxial stacking. Heteroduplex formation is disfavored by sharp turns or a lack of helical order in single-stranded regions, competition from bases displaced from a stem, and stable tertiary interactions. The study is relevant to duplex formation on oligonucleotide microarrays and to antisense technologies.

Crystal structure of a conserved ribosomal protein-RNA complex

The structure of a highly conserved complex between a 58-nucleotide domain of large subunit ribosomal RNA and the RNA-binding domain of ribosomal protein L11 has been solved at 2.8 angstrom resolution. It reveals a precisely folded RNA structure that is stabilized by extensive tertiary contacts and contains an unusually large core of stacked bases. A bulge loop base from one hairpin of the RNA is intercalated into the distorted major groove of another helix; the protein locks this tertiary interaction into place by binding to the intercalated base from the minor groove side. This direct interaction with a key ribosomal RNA tertiary interaction suggests that part of the role of L11 is to stabilize an unusual RNA fold within the ribosome.

A detailed view of a ribosomal active site: the structure of the L11-RNA complex

We report the crystal structure of a 58 nucleotide fragment of 23S ribosomal RNA bound to ribosomal protein L11. This highly conserved ribonucleoprotein domain is the target for the thiostrepton family of antibiotics that disrupt elongation factor function. The highly compact RNA has both familiar and novel structural motifs. While the C-terminal domain of L11 binds RNA tightly, the N-terminal domain makes only limited contacts with RNA and is proposed to function as a switch that reversibly associates with an adjacent region of RNA. The sites of mutations conferring resistance to thiostrepton and micrococcin line a narrow cleft between the RNA and the N-terminal domain. These antibiotics are proposed to bind in this cleft, locking the putative switch and interfering with the function of elongation factors.

Stabilization of the anticodon stem-loop of tRNALys,3 by an A+-C base-pair and by pseudouridine

NMR spectroscopy was used to determine the solution structures of RNA oligonucleotides comprising the anticodon domain of tRNALys,3. The structural effects of the pseudouridine modification at position 39 were investigated and are well correlated with changes in thermodynamic parameters derived from temperature dependent UV measurements. The pseudouridine-containing hairpin is thermodynamically more stable than the unmodified hairpin by 5 degreesC, and this corresponds with increased base stacking on the 3' side of the tRNA anticodon loop. An A+38-C32 base-pair also forms at the base of the anticodon stem with an approximate pKa of 6 for A38. Formation of the A+-C base-pair increases the Tm of both pseudouridine modified and unmodified RNA hairpins by 5-6 degreesC, and decreases the DeltaG degrees for hairpin formation by 1 kcal/mol. Solution structures were determined for both psi39 and unmodified hairpins under limiting pH conditions at pH 5 and pH 7 to assess the structural effects of both psi modification and the additional A+-C base-pair on tRNALys,3 structure. The A+38-C32 base-pair strengthens the 31-39 base-pair, and induces formation of a dynamic U33-A37 base-pair that effectively reduces the normal seven nucleotide anticodon loop to a three nucleotide UUU loop. These undermodified tRNALys,3 anticodon loops are distinctly different from those seen for other tRNAs exemplified by tRNAPhe. The conformation of the tRNA loop has important implications for the role of nucleoside modification in codon-anticodon recognition and for utilization of tRNALys,3 by HIV-1 as the native reverse transcriptase primer.

Secondary structure dimorphism and interconversion between hairpin and duplex form of oligoribonucleotides

RNA hairpins can alternatively form a dimer with a bulged loop flanked by regularly base paired regions. [1H]NMR spectroscopy and native gel electrophoresis were used to study how the sequence of nucleotides in the loop of the hairpin affect the hairpin-duplex interconversion. As a model system, a hairpin containing 7 nucleotides in the loop and 5 base pairs in the stem was used. The loop size was gradually reduced from 7 to 4 nucleotides, yielding finally the stable UNCG tetraloop. Single nucleotide mutations were performed to investigate the influence of the self-complementarity of the loop sequence on the dimerization. The results demonstrate that (1) the initial fraction of hairpin is determined by concentration of the oligonucleotide, the annealing procedure, and the relative stability of the loop, (2) the degree of self-complementarity of the loop sequence of the hairpin governs the dimerization kinetics, and (3) oligonucleotides complementary to the loop sequence decrease the dimerization rate. We propose a secondary structure-based model for the dimerization reaction of RNA hairpins in which the formation of intermolecular base pairs between self-complementary nucleotides of the loops represents the nucleation step.

The three-dimensional structure of the ribosome and its components

Exciting progress has been made in the last decade by those who use physical methods to study the structure of the ribosome and its components. The structures of 10 ribosomal proteins and three isolated ribosomal protein domains are known, and the conformations of a significant number of rRNA sequences have been determined. Electron microscopists have made major advances in the analysis of images of ribosomes, and microscopically derived ribosome models at resolutions approaching 10A are likely quite soon. Furthermore, ribosome crystallographers are on the verge of phasing the diffraction patterns they have had for several years, and near-atomic resolution models for entire ribosomal subunits could emerge from this source at any time. The literature relevant to these developments is reviewed below.

Structure, recognition and adaptive binding in RNA aptamer complexes

Novel features of RNA structure, recognition and discrimination have been recently elucidated through the solution structural characterization of RNA aptamers that bind cofactors, aminoglycoside antibiotics, amino acids and peptides with high affinity and specificity. This review presents the solution structures of RNA aptamer complexes with adenosine monophosphate, flavin mononucleotide, arginine/citrulline and tobramycin together with an example of hydrogen exchange measurements of the base-pair kinetics for the AMP-RNA aptamer complex. A comparative analysis of the structures of these RNA aptamer complexes yields the principles, patterns and diversity associated with RNA architecture, molecular recognition and adaptive binding associated with complex formation.

The structure of an essential splicing element: stem loop IIa from yeast U2 snRNA

BACKGROUND

Eukaryotic genes are usually transcribed as precursor mRNAs which are then spliced, removing introns to produce functional mRNAs. Splicing is performed by the spliceosome and provides an important level of post-translational control of gene expression. Stem loop IIa from U2 small nuclear (sn)RNA is required for the efficient association of the U2 small nuclear ribonuclear protein (snRNP) with the nascent spliceosome in yeast. Genetic analysis suggests that stem loop IIa is involved in RNA-protein interactions early in splicing, and it may also interact with other RNA sequences in U2. The sequence of loop IIa is well conserved, consistent with the idea that this loop is important for function.

RESULTS

We have solved the structure of U2A, a 20-base analogue of stem loop IIa from Saccharomyces cerevisiae, using NMR and restrained molecular dynamics. In the process, we have demonstrated the efficacy of a new structure calculation protocol, torsion angle molecular dynamics. The structure that has emerged, which is consistent with the in vivo chemical protection data available for stem loop IIa in the context of intact U2 snRNA, contains a sheared GA pair followed by a U-turn in the loop. The U-turn conformation, which resembles the U-turns in tRNA anticodon loops, makes this stretch of U2 snRNA an obvious target for interactions with proteins and/or other RNA sequences.

CONCLUSIONS

The phenotypes of many stem loop IIa mutants can be rationalized assuming that the U-turn conformation in the loop must be preserved for efficient splicing. This observation, combined with the phylogenetic conservation of its sequence, suggests that the conformation of the loop of stem loop IIa is essential for its function in pre-mRNA splicing.

The detailed structure of tandem G.A mismatched base-pair motifs in RNA duplexes is context dependent

The solution structure of the RNA duplex (rGGGCUGAAGCCCU), containing tandem G.A mismatches has been determined by NMR spectroscopy and restrained molecular dynamics. A homonuclear 3D TOCSY-NOESY was used to derive 18 to 30 distance restraints per nucleotide, as well as all gamma torsion angles and sugar puckers for the central UGAA part of the molecule. Using these constraints, together with cross-strand distances, involving exchangeable imino protons, and essentially all other torsion angles that can accurately be determined (i.e. beta, epsilon) otherwise, the structure of the UGAA domain could be determined with high precision (r.m.s.d. 0.62 A), without the aid of isotopically enriched RNA. The G.A base-pairs are of the sheared pairing type, with both nucleotides in the anti conformation, and hydrogen bonds between the guanine 2-amino and the adenine N7 and between the guanine N3 and the adenine 6-amino. Surprisingly the sugar of the guanosine of the G.A. mismatch adopts a 2'-endo sugar pucker conformation. Comparison with other RNA structures, in which two such G.A base-pairs are formed reveals that this detailed structure depends on the identity of the base 5' to the guanosine in the tandem G.A base-pairs. A geometrical model for the incorporation of sheared tandem G.A base-pairs in A-form helices is formulated, which explains the distinct different stacking properties and helical parameters in sequences containing tandem, sheared G.A base-pairs.

The structure of the isolated, central hairpin of the HDV antigenomic ribozyme: novel structural features and similarity of the loop in the ribozyme and free in solution

The structure of an RNA hairpin containing a seven-nucleotide loop that is present in the self-cleaving sequence of hepatitis delta virus antigenomic RNA was determined by high resolution NMR spectroscopy. The loop, which is composed of only one purine and six pyrimidines, has a suprisingly stable structure, mainly supported by sugar hydroxyl hydrogen bonds and base-base and base-phosphate stacking interactions. Compared with the structurally well-determined, seven-membered anticodon loop in tRNA, the sharp turn which affects the required 180 degrees change in direction of the sugar-phosphate backbone in the loop is shifted one nucleotide in the 3' direction. This change in direction can be characterized as a reversed U-turn. It is expected that the reversed U-turn may be found frequently in other molecules as well. There is evidence for a new non-Watson-Crick UC base pair formed between the first and the last residue in the loop, while most of the other bases in the loop are pointing outwards making them accessible to solvent. From chemical modification, mutational and photocrosslinking studies, a similar picture develops for the structure of the hairpin in the active ribozyme indicating that the loop structure in the isolated hairpin and in the ribozyme is very similar.

Structural characterization of three RNA hexanucleotide loops from the internal ribosome entry site of polioviruses

Structural characteristics of three RNA hairpins from the internal ribosome entry site of poliovirus mRNAs have been determined in solution by NMR. Complete proton, phosphorus and carbon resonance assignments were made for the three 16 nt hairpins. The loop sequences, 5'-AAUCCA , AAACCA and GAACCA, have been shown to be essential for viral mRNA translation. NOESY spectra for the three oligomers were very similar indicating a common three dimensional structure. Stems were A-type duplexes with C3'-endo sugar pucker. In the loops, sequential base stacking interactions were detected for all bases except between U8/A8 and C9, indicating a turn in the phosphodiester backbone at this point. Only one nucleotide, U8/A8, had a sugar pucker which deviated appreciably from C3'-endo. The final base in the loop, A11, exhibited an unusual gauche (-) gamma angle. An ensemble of 10 structures calculated for one hairpin using restrained molecular dynamics shows that the first three bases of the loop are turned so as to be exposed to the exterior of the molecule, while the remaining three bases are in an orientation approximating a continuation of the stem helix. Structure calculations and NMR relaxation measurements indicate that the loop apex is subject to considerable local dynamics.

The conformation of ribosomes and rRNA

During the past 18 months, electron microscopists have published two reconstructions of the Escherichia coli ribosome, independently derived from images of unstained particles. The resolutions of their images are 20-25 A-much higher than any previously available. During the same time, NMR spectroscopists have provided an atomic-resolution model for the A-site region of 16S rRNA complexed with paromomycin that explains much of what is known about the interaction of aminoglycoside antibiotics with ribosomes.

Recent solution structures of RNA and its complexes with drugs, peptides and proteins

The past two years have seen remarkable progress in the study of RNA structure: the predicted era of RNA structural biology has arrived. Crystallographic structures of the hammerhead ribozyme and of a large subunit of a group I self-splicing intron have begun to reveal the structural basis of RNA enzymatic activity. A remarkable number of structures of small RNAs and of complexes with drugs, peptides and one protein domain have been determined by NMR.

Solution structure of the conserved 16 S-like ribosomal RNA UGAA tetraloop

The solution structure of the highly conserved UGAA tetraloop found at the 3' end of eukaryotic 16 S-like ribosomal RNA has been solved by nuclear magnetic resonance spectroscopy in the form of the 12 nucleotide hairpin 5'-GGUG[UGAA]CACC. The UGAA tetraloop displays a novel fold. The backbone turn occurs between the G and the third A in the loop, with the U and G in a 5' stack and the As in a 3' stacking arrangement. The loop is closed by a U-A mismatch in which the O2, 2'OH, and O4' groups of the U are within hydrogen bonding distance of the amino group of the A. The tetraloop does not make a uridine-turn, even though its sequence is identical to a U-turn found within the anticodon loop of tRNA(Phe). The hydrogen bonding pattern in the tetraloop provides insight into the function of base modifications found in vivo within this portion of 16 S-like rRNA.

On the conformation of the anticodon loops of initiator and elongator methionine tRNAs

The solution conformations of analogues of initiator and elongator tRNA anticodon stem-loops have been compared by NMR. The data indicate that both have conformations closely similar to those reported for crystalline elongator tRNAs. The two loops differ in their dynamics, however: that of the elongator analogue is more flexible than its initiator counterpart. The anticodon stem-loops of initiator tRNAs are more likely to be distinguished from those of elongator tRNAs during initiation on the basis of their distinctive stem sequences, than they are by differences in the conformations of their anticodon loops.

Remarkable morphological variability of a common RNA folding motif: the GNRA tetraloop-receptor interaction

One of the most common RNA tertiary interactions involves the docking of GNRA hairpin loops into stem-loop structures on other regions of RNA. Domain 5 of the group II intron interacts with Domain 1 through such an interaction, which has been characterized thermodynamically and kinetically for the ai5g intron. Using this system, it was possible to test the morphological tolerances of the GNRA tetraloop involved in tertiary interactions. The data presented herein show that a GNRA tetraloop can still participate in tertiary interaction after being physically cut at any phosphodiester linkage within the loop. The "nicked tetraloop" can be expanded by many nucleotides in either direction and the covalently continuous loop can also be expanded without loss of interaction energy. In the context of the nicked tetraloop, the second nucleotide of the tetraloop sequence can be completely deleted without loss of function. By examining radical alterations in tetraloop structure, this study helps define the minimal sequence and structural requirements of a GNRA motif involved in long-range tertiary interaction. It shows that "tetraloop"-like structures capable of forming tertiary interactions can be imbedded in unexpected contexts, such as internal loops and apparently open structure within RNA. It demonstrates that pentaloops and hexaloops can form the same type of interaction, with almost equal affinity, as a tetraloop. Taken together, these data suggest a more generic term for the GNRA tetraloop-receptor interaction: It is proposed herein that the term "GNRA tetraloop" be replaced by "GNn/RA", where n represents a variable number of nucleotides and / indicates that the loop can be divided and interrupted by other sequences.