Consecutive G·U base pairs inside RNA helices can be destabilizing, while those at the ends of helices are thermodynamically stabilizing. To determine if this paradox could be explained by differences in base stacking, we determined the high-resolution (1.32 Å) crystal structure of (5′-GGUGGCUGUU-3′)2 and studied three sequences with four consecutive terminal G·U pairs by NMR spectroscopy. In the crystal structure of (5′-GGUGGCUGUU-3′)2, the helix is overwound but retains the overall features of A-form RNA. The penultimate base steps at each end of the helix have high base overlap and contribute to the unexpectedly favorable energetic contribution for the 5′-GU-3′/3′-UG-5′ motif in this helix position. The balance of base stacking and helical twist contributes to the positional dependence of G·U pair stabilities. The energetic stabilities and similarity to A-form RNA helices suggest that consecutive G·U pairs would be recognized by RNA helix binding proteins, such as Dicer and Ago. Thus, these results will aid future searches for target sites of small RNAs in gene regulation.

Viral RNA structure prediction is a valuable tool for development of drugs against viral disease. This work discusses different approaches to predicting encapsidated viral RNA and highlights satellite tobacco mosaic virus (STMV) RNA as a model system with excellent crystallography data. Fundamentally important issues for debate include thermodynamic versus kinetic control of virus assembly and the possible consequences of quasi-species in the primary structure on RNA secondary structure prediction of a single structure or an ensemble of structures. Multiple computational tools and chemical reagents are now available for improved viral RNA structure prediction. Two different predicted structures for encapsidated STMV RNA result from differences in three main areas: a different approach and philosophy to studying encapsidated viral RNA, an emphasis on different RNA motifs, and technical differences in computational methods and chemical reagents. The experiments with traditional chemical probing and SHAPE reagents are compared in terms of chemistry, results, and interpretation for STMV RNA as well as other RNA protein assemblies, such as the 5′UTR of HIV and the ribosome. This discussion of the challenges of viral RNA structure prediction will lead to new experiments and improved future predictions for viral RNA.

Consecutive terminal mismatches add thermodynamic stability to RNA duplexes and occur frequently in microRNA-mRNA interactions. Accurate thermodynamic stabilities of consecutive terminal mismatches contribute to the development of specific, high-affinity siRNA therapeutics. Consecutive terminal adenosine mismatches (TAMS) are studied at different salt concentrations, with polyethylene glycol cosolutes, and with locked nucleic acid (LNA) substitutions. These measurements provide benchmarks for the application of thermodynamic predictions to different physiological or therapeutic conditions. The salt dependence for RNA duplex stability is similar for TAMS, internal loops, and Watson–Crick duplexes. A unified model for predicting the free energy of an RNA duplex with or without loops and mismatches at lower sodium concentrations is presented. The destabilizing effects of PEG 200 are larger for TAMS than internal loops or Watson–Crick duplexes, which may result from different base stacking conformations, dynamics, and water hydration. In contrast, LNA substitutions stabilize internal loops much more than TAMS. Surprisingly, the average per adenosine increase in stability for LNA substitutions in internal loops is −1.82 kcal/mol and only −0.20 kcal/mol for TAMS. The stabilities of TAMS and internal loops with LNA substitutions have similar favorable free energies. Thus, the unfavorable free energy of adenosine internal loops is largely an entropic effect. The favorable stabilities of TAMS result mainly from base stacking. The ability of RNA duplexes to form extended terminal mismatches in the absence of proteins such as argonaute and identifying the enthalpic contributions to terminal mismatch stabilities provide insight into the physical basis of microRNA–mRNA molecular recognition and specificity.

Resistance mutations to antibiotics targeting rRNA can be far from the drug-binding site. Crystallography studies revealed that the antibiotic resistance mutation G2482A (G2447A in E. coli) in Haloarcula marismortui 50S ribosomes does not directly contact the drug or introduce changes to the ribosomal structure except for losing a potassium ion coordinated to a base triple at the drug-binding site. Using molecular dynamics simulations, we tested hypotheses regarding the effects of the G2482A mutation and ion coordination on the conformational dynamics of the 50S ribosome. Simulations show that the mutation enhances conformational fluctuation at the antibiotic binding site, weakens the hydrogen-bonding network, and increases flexibility at the 50S peptidyl transferase center (PTC). Our data supports the view that distant mutations can perturb the dynamic network in the ribosomal PTC, thereby raising the entropic cost of antibiotic binding. These results underscore the importance of considering conformational dynamics in rational drug design.Keywords: antibiotic resistance; molecular dynamics simulations; ribosome dynamics;

The Bacillus subtilis phage ϕ29 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.

Consecutive GU pairs at the ends of RNA helices provide significant thermodynamic stability between −1.0 and −3.8 kcal/mol at 37 °C, which is equivalent to approximately 2 orders of magnitude in the value of a binding constant. The thermodynamic stabilities of GU pairs depend on the sequence, stacking orientation, and position in the helix. In contrast to GU pairs in the middle of a helix that may be destabilizing, all consecutive terminal GU pairs contribute favorable thermodynamic stability. This work presents measured thermodynamic stabilities for 30 duplexes containing two, three, or four consecutive GU pairs at the ends of RNA helices and a model to predict the thermodynamic stabilities of terminal GU pairs. Imino proton NMR spectra show that the terminal GU nucleotides form hydrogen-bonded pairs. Different orientations of terminal GU pairs can have different conformations with equivalent thermodynamic stabilities. These new data and prediction model will help improve RNA secondary structure prediction, identification of miRNA target sequences with GU pairs, and efforts to understand the fundamental physical forces directing RNA structure and energetics.

The thermodynamic stabilities of consecutive mismatches at the ends of RNA helices are determined by the 3′ terminal nucleotides. More than 40 RNA duplexes containing terminal motifs of 3 or more nucleotides were studied by optical melting experiments. Up to three noncanonical pairs of nucleotides at the end of RNA helices provide additional thermodynamic stability. 3′ nucleotides contribute more stability than 5′ nucleotides, and purines contribute more stability than pyrimidines. The additional stability of a second or third 3′ nucleotide stacking on a purine is the same for both dangling ends and consecutive terminal mismatches. Current predictions underestimate RNA duplex stabilities with terminal motifs by 1.4 kcal/mol on average, which is an order of magnitude in a binding constant at 37 °C. Accurate thermodynamic parameters for these terminal motifs will contribute to improvements in RNA secondary structure predictions, identification of microRNA targets, and design of siRNA therapeutics with fewer off-target effects.

Viral genomic RNA adopts many conformations during its life cycle as the genome is replicated, translated, and encapsidated. The high-resolution crystallographic structure of the satellite tobacco mosaic virus (STMV) particle reveals 30 helices of well-ordered RNA. The crystallographic data provide global constraints on the possible secondary structures for the encapsidated RNA. Traditional free energy minimization methods of RNA secondary structure prediction do not generate structures consistent with the crystallographic data, and to date no complete STMV RNA basepaired secondary structure has been generated. RNA-protein interactions and tertiary interactions may contribute a significant degree of stability, and the kinetics of viral assembly may dominate the folding process. The computational tools, Helix Find & Combine, Crumple, and Sliding Windows and Assembly, evaluate and explore the possible secondary structures for encapsidated STMV RNA. All possible hairpins consistent with the experimental data and a cotranscriptional folding and assembly hypothesis were generated, and the combination of hairpins that was most consistent with experimental data is presented as the best representative structure of the ensemble. Multiple solutions to the genome packaging problem could be an evolutionary advantage for viruses. In such cases, an ensemble of structures that share favorable global features best represents the RNA fold.