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.     

Automated classification of RNA 3D motifs and the RNA 3D Motif Atlas

The analysis of atomic-resolution RNA three-dimensional (3D) structures reveals that many internal and hairpin loops are modular, recurrent, and structured by conserved non-Watson–Crick base pairs. Structurally similar loops define RNA 3D motifs that are conserved in homologous RNA molecules, but can also occur at nonhomologous sites in diverse RNAs, and which often vary in sequence. To further our understanding of RNA motif structure and sequence variability and to provide a useful resource for structure modeling and prediction, we present a new method for automated classification of internal and hairpin loop RNA 3D motifs and a new online database called the RNA 3D Motif Atlas. To classify the motif instances, a representative set of internal and hairpin loops is automatically extracted from a nonredundant list of RNA-containing PDB files. Their structures are compared geometrically, all-against-all, using the FR3D program suite. The loops are clustered into motif groups, taking into account geometric similarity and structural annotations and making allowance for a variable number of bulged bases. The automated procedure that we have implemented identifies all hairpin and internal loop motifs previously described in the literature. All motif instances and motif groups are assigned unique and stable identifiers and are made available in the RNA 3D Motif Atlas (http://rna.bgsu.edu/motifs), which is automatically updated every four weeks. The RNA 3D Motif Atlas provides an interactive user interface for exploring motif diversity and tools for programmatic data access.     

The 3′ end of Turnip crinkle virus contains a highly interactive structure including a translational enhancer that is disrupted by binding to the RNA-dependent RNA polymerase

Precise temporal control is needed for RNA viral genomes to translate sufficient replication-required products before clearing ribosomes and initiating replication. A 3′ translational enhancer in Turnip crinkle virus (TCV) overlaps an internal T-shaped structure (TSS) that binds to 60S ribosomal subunits. The higher-order structure in the region was examined through alteration of critical sequences revealing novel interactions between an H-type pseudoknot and upstream residues, and between the TSS and internal and terminal loops of an upstream hairpin. Our results suggest that the TSS forms a stable scaffold that allows for simultaneous interactions with external sequences through base pairings on both sides of its large internal symmetrical loop. Binding of TCV RNA-dependent RNA polymerase (RdRp) to the region potentiates a widespread conformational shift with substantial rearrangement of the TSS region, including the element required for efficient ribosome binding. Degrading the RdRp caused the RNA to resume its original conformation, suggesting that the initial conformation is thermodynamically favored. These results suggest that the 3′ end of TCV folds into a compact, highly interactive structure allowing RdRp access to multiple elements including the 3′ end, which causes structural changes that potentiate the shift between translation and replication.     

Structural motifs in ribosomal RNAs: implications for RNA design and genomics

The various motifs of RNA molecules are closely related to their structural and functional properties. To better understand the nature and distributions of such structural motifs (i.e., paired and unpaired bases in stems, junctions, hairpin loops, bulges, and internal loops) and uncover characteristic features, we analyze the large 16S and 23S ribosomal RNAs of Escherichia coli. We find that the paired and unpaired bases in structural motifs have characteristic distribution shapes and ranges; for example, the frequency distribution of paired bases in stems declines linearly with the number of bases, whereas that for unpaired bases in junctions has a pronounced peak. Significantly, our survey reveals that the ratio of total (over the entire molecule) unpaired to paired bases (0.75) and the fraction of bases in stems (0.6), junctions (0.16), hairpin loops (0.12), and bulges/internal loops (0.12) are shared by 16S and 23S ribosomal RNAs, suggesting that natural RNAs may maintain certain proportions of bases in various motifs to ensure structural integrity. These findings may help in the design of novel RNAs and in the search (via constraints) for RNA-coding motifs in genomes, problems of intense current focus.     

RNA Structure and Stability

RNA molecules fold into specific base-paired conformations that contain single-stranded regions, A-form double helices, hairpin loops, internal loops, bulges, junctions, pseudoknots, kissing hairpins, and so forth. These structural motifs are recognized by proteins, other RNAs, and other parts of the same RNA. The interactions of these structural elements are crucial to the biological functions of the RNA molecules. We describe the different motifs and discuss their thermodynamic stabilities relative to single strands of RNA. The stabilities determine under what conditions they occur and whether they change when interacting with proteins or other ligands.     

Analysis of Four-Way Junctions in RNA Structures

RNA secondary structures can be divided into helical regions composed of canonical Watson-Crick and related base pairs, as well as single-stranded regions such as hairpin loops, internal loops, and junctions. These elements function as building blocks in the design of diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze existing RNA four-way junctions in terms of base-pair interactions and 3D configurations. Specifically, we identify nine broad junction families according to coaxial stacking patterns and helical configurations. We find that helices within junctions tend to arrange in roughly parallel and perpendicular patterns and stabilize their conformations using common tertiary motifs such as coaxial stacking, loop-helix interaction, and helix packing interaction. Our analysis also reveals a number of highly conserved base-pair interaction patterns and novel tertiary motifs such as A-minor-coaxial stacking combinations and sarcin/ricin motif variants. Such analyses of RNA building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

Extracting stacking interaction parameters for RNA from the data set of native structures

A crucial step in the determination of the three-dimensional native structures of RNA is the prediction of their secondary structures, which are stable independent of the tertiary fold. Accurate prediction of the secondary structure requires context-dependent estimates of the interaction parameters. We have exploited the growing database of natively folded RNA structures in the Protein Data Bank (PDB) to obtain stacking interaction parameters using a knowledge-based approach. Remarkably, the calculated values of the resulting statistical potentials (SPs) are in excellent agreement with the parameters determined using measurements in small oligonucleotides. We validate the SPs by predicting 74% of the base-pairs in a dataset of structures using the ViennaRNA package. Interestingly, this number is similar to that obtained using the measured thermodynamic parameters. We also tested the efficacy of the SP in predicting secondary structure by using gapless threading, which we advocate as an alternative method for rapidly predicting RNA structures. For RNA molecules with less than 700 nucleotides, about 70% of the native base-pairs are correctly predicted. As a further validation of the SPs we calculated Z-scores, which measure the relative stability of the native state with respect to a manifold of higher free energy states. The computed Z-scores agree with estimates made using calorimetric measurements for a few RNA molecules. Structural analysis was used to rationalize the success and failures of SP and experimentally determined parameters. First, from the near perfect linear relationship between the number of native base-pairs and sequence length, we show that nearly 46% of nucleotides are not in stacks. Second, by analyzing the suboptimal structures that are generated in gapless threading we show that the SPs and experimentally determined parameters are most successful in predicting stacks that end in hairpins. These results show that further improvement in secondary structure prediction requires reliable estimates of interaction parameters for loops, bulges, and stacks that do not end in hairpins.     

Solution Structure of NPSL2, A Regulatory Element in the oncomiR-1 RNA

The miR-17 ～ 92a polycistron, also known as oncomiR-1, is commonly overexpressed in multiple cancers and has several oncogenic properties. OncomiR-1 encodes six constituent microRNAs (miRs), each enzymatically processed with different efficiencies. However, the structural mechanism that regulates this differential processing remains unclear. Chemical probing of oncomiR-1 revealed that the Drosha cleavage sites of pri-miR-92a are sequestered in a four-way junction. NPSL2, an independent stem loop element, is positioned just upstream of pri-miR-92a and sequesters a crucial part of the sequence that constitutes the basal helix of pri-miR-92a. Disruption of the NPSL2 hairpin structure could promote the formation of a pri-miR-92a structure that is primed for processing by Drosha. Thus, NPSL2 is predicted to function as a structural switch, regulating pri-miR-92a processing. Here, we determined the solution structure of NPSL2 using solution NMR spectroscopy. This is the first high-resolution structure of an oncomiR-1 element. NPSL2 adopts a hairpin structure with a large, but highly structured, apical and internal loops. The 10-bp apical loop contains a pH-sensitive A⁺·C mismatch. Additionally, several adenosines within the apical and internal loops have elevated pK_a values. The protonation of these adenosines can stabilize the NPSL2 structure through electrostatic interactions. Our study provides fundamental insights into the secondary and tertiary structure of an important RNA hairpin proposed to regulate miR biogenesis.     

An Approach to Identifying the Functionally Important RNA Sites by Complementary Addressed Modification

An approach based on complementary addressed modification of nucleic acids by oligodeoxyribonucleotide derivatives was proposed for changing the spatial structure of particular RNA sites in order to study their role in the biological activity of the total RNA molecule. Hepatitis C virus (HCV) IRES was used as a model. Oligodeoxyribonucleotide derivatives contained a 4-[N-(2-chloroethyl)-N-methylamino]benzylamino group at the 5-phosphate and were complementary to various RNA sites located in regions of hairpins II, IIId, or IIIe. Covalent adducts resulting from RNA alkylation with the derivatives were isolated by denaturing PAGE and tested for binding with the 40S subunit of human ribosomes. Structural alteration of hairpin II had no effect on the binding, whereas that of hairpin IIIe substantially reduced it. The RNA with modified hairpin IIId showed virtually no binding with the 40S subunit. Hairpin IIId was assumed to play the critical role in the binding of HCV IRES with the 40S subunit.     

Fully differentiable coarse-grained and all-atom knowledge-based potentials for RNA structure evaluation

RNA molecules play integral roles in gene regulation, and understanding their structures gives us important insights into their biological functions. Despite recent developments in template-based and parameterized energy functions, the structure of RNA--in particular the nonhelical regions--is still difficult to predict. Knowledge-based potentials have proven efficient in protein structure prediction. In this work, we describe two differentiable knowledge-based potentials derived from a curated data set of RNA structures, with all-atom or coarse-grained representation, respectively. We focus on one aspect of the prediction problem: the identification of native-like RNA conformations from a set of near-native models. Using a variety of near-native RNA models generated from three independent methods, we show that our potential is able to distinguish the native structure and identify native-like conformations, even at the coarse-grained level. The all-atom version of our knowledge-based potential performs better and appears to be more effective at discriminating near-native RNA conformations than one of the most highly regarded parameterized potential. The fully differentiable form of our potentials will additionally likely be useful for structure refinement and/or molecular dynamics simulations.     

RGRNA: prediction of RNA secondary structure based on replacement and growth of stems

Owing to their structural diversity, RNAs perform many diverse biological functions in the cell. RNA secondary structure is thus important for predicting RNA function. Here, we propose a new combinatorial optimization algorithm, named RGRNA, to improve the accuracy of predicting RNA secondary structure. Following the establishment of a stempool, the stems are sorted by length, and chosen from largest to smallest. If the stem selected is the true stem, the secondary structure of this stem when combined with another stem selected at random will have low free energy, and the free energy will tend to gradually diminish. The free energy is considered as a parameter and the structure is converted into binary numbers to determine stem compatibility, for step-by-step prediction of the secondary structure for all combinations of stems. The RNA secondary structure can be predicted by the RGRNA method. Our experimental results show that the proposed algorithm outperforms RNAfold in terms of sensitivity, specificity, and Matthews correlation coefficient value.     

Mg2+ Dependence of the Structure and Thermodynamics of Wheat Germ and Lupin Seeds 5S rRNA

Abstract

The formation and stability of structural elements in two 5S rRNA molecules from wheat germ (WG) and lupin seeds (LS) as a function of Mg²⁺ concentration in solution was determined using the adiabatic differential scanning microcalorimetry (DSC). The experimentally determined thermodynamic parameters are compared with calculations using thermodynamic databases used for prediction of RNA structure. The 5S rRNA molecules which show minor differences in the nucleotide sequence display very different thermal unfolding profiles (DSC profiles). Numerical deconvolution of DSC profiles provided information about structural transformations that take place in both 5S rRNA molecules. A comparative analysis of DSC data and the theoretical thermodynamic models of the structure was used to establish a relationship between the constituting transitions found in the melting profiles and the unfolding of structural domains of the 5S rRNA and stability of its particular helical elements.

Increased concentration of Mg²⁺ ions induces additional internal interactions stabilising 5S rRNA structures found at low Na⁺ concentrations. Observed conformational transitions suggest a structural model in which the extension of helical region E dominates over the postulated tertiary interaction between hairpin loops. We propose that helix E is stabilised by a sequence of non-standard pairings extending this helix by the formation of tetra loop e and an almost total reduction of loop d between helices E and D. Two hairpin structures in both 5S rRNA molecules: the extended C-C' and the extended E-E'-E” hairpins appear as the most stable elements of the structure. The cooperativity of the unfolding of helixes in these 5S rRNA molecules changes already at 2 mM Mg²⁺.     

Substrate-specific kinetics of Dicer-catalyzed RNA processing

The specialized ribonuclease Dicer plays a central role in eukaryotic gene expression by producing small regulatory RNAs—microRNAs (miRNAs) and short interfering RNAs (siRNAs)—from larger double-stranded RNA (dsRNA) substrates. Although Dicer will cleave both imperfectly base-paired hairpin structures (pre-miRNAs) and perfect duplexes (pre-siRNAs) in vitro, it has not been clear whether these are mechanistically equivalent substrates and how dsRNA binding proteins such as trans-activation response (TAR) RNA binding protein (TRBP) influence substrate selection and RNA processing efficiency. We show here that human Dicer is much faster at processing a pre-miRNA substrate compared to a pre-siRNA substrate under both single and multiple turnover conditions. Maximal cleavage rates (V_max) calculated by Michaelis-Menten analysis differed by more than 100-fold under multiple turnover conditions. TRBP was found to enhance dicing of both substrates to similar extents, and this stimulation required the two N-terminal dsRNA binding domains of TRBP. These results demonstrate that multiple factors influence dicing kinetics. While TRBP stimulates dicing by enhancing the stability of Dicer-substrate complexes, Dicer itself generates product RNAs at rates determined at least in part by the structural properties of the substrate.     

Crystal Structure of the Human Rad9-Hus1-Rad1 Clamp

Three evolutionarily conserved proteins, Rad9, Hus1, and Rad1, form a heterotrimeric 9-1-1 complex that plays critical roles in cellular responses to DNA damage by activating checkpoints and by recruiting DNA repair enzymes to DNA lesions. We have determined the crystal structure of the human Rad9 (residues 1-272)-Hus1-Rad1 complex at 2.5 Å resolution. The 9^1-272-1-1 complex forms a closed ring, with each subunit having a similar structure. Despite its high level of similarity to proliferating cell nucleus antigen in terms of overall structure, the 9^1-272-1-1 complex exhibits notable differences in local structures, including interdomain connecting loops, H2 and H3 helices, and loops in the vicinity of the helices of each subunit. These local structural variations provide several unique features to the 9-1-1 heterotrimeric complex—including structures of intermolecular interfaces and the inner surface around the central hole, and different electrostatic potentials at and near the interdomain connecting loops of each 9-1-1 subunit—compared to the proliferating cell nucleus antigen trimer. We propose that these structural features allow the 9-1-1 complex to bind to a damaged DNA during checkpoint control and to serve as a platform for base excision repair. We also show that the 9^1-272-1-1 complex, but not the full-length 9-1-1 complex, forms a stable complex with the 5′ recessed DNA, suggesting that the C-terminal tail of Rad9 is involved in the regulation of the 9-1-1 complex in DNA binding.     

Salt dependence of nucleic acid hairpin stability

Single-stranded junctions/loops are frequently occurring structural motifs in nucleic acid structures. Due to the polyanionic nature of the nucleic acid backbone, metal ions play a crucial role in the loop stability. Here we use the tightly bound ion theory, which can account for the possible ion correlation and ensemble (fluctuation) effects, to predict the ion-dependence of loop and stem-loop (hairpin) free energies. The predicted loop free energy is a function of the loop length, the loop end-to-end distance, and the ion (Na⁺ and Mg²⁺ in this study) concentrations. Based on the statistical mechanical calculations, we derive a set of empirical formulas for the loop thermodynamic parameters as functions of Na⁺ and Mg²⁺ concentrations. For three specific types of loops, namely, hairpin, bulge, and internal loops, the predicted free energies agree with the experimental data. Further applications of these empirical formulas to RNA and DNA hairpin stability lead to good agreements with the available experimental data. Our results indicate that the ion-dependent loop stability makes significant contribution to the overall ion-dependence of the hairpin stability.     

High-throughput genetic identification of functionally important regions of the yeast DEAD-box protein Mss116p

The Saccharomyces cerevisiae DEAD-box protein Mss116p is a general RNA chaperone that functions in splicing mitochondrial group I and group II introns. Recent X-ray crystal structures of Mss116p in complex with ATP analogs and single-stranded RNA show that the helicase core induces a bend in the bound RNA, as in other DEAD-box proteins, while a C-terminal extension (CTE) induces a second bend, resulting in RNA crimping. Here, we illuminate these structures by using high-throughput genetic selections, unigenic evolution, and analyses of in vivo splicing activity to comprehensively identify functionally important regions and permissible amino acid substitutions throughout Mss116p. The functionally important regions include those containing conserved sequence motifs involved in ATP and RNA binding or interdomain interactions, as well as previously unidentified regions, including surface loops that may function in protein-protein interactions. The genetic selections recapitulate major features of the conserved helicase motifs seen in other DEAD-box proteins but also show surprising variations, including multiple novel variants of motif III (SAT). Patterns of amino acid substitutions indicate that the RNA bend induced by the helicase core depends on ionic and hydrogen-bonding interactions with the bound RNA; identify a subset of critically interacting residues; and indicate that the bend induced by the CTE results primarily from a steric block. Finally, we identified two conserved regions—one the previously noted post II region in the helicase core and the other in the CTE—that may help displace or sequester the opposite RNA strand during RNA unwinding.     

Calculation of pKas in RNA: on the structural origins and functional roles of protonated nucleotides

pK(a) calculations based on the Poisson-Boltzmann equation have been widely used to study proteins and, more recently, DNA. However, much less attention has been paid to the calculation of pK(a) shifts in RNA. There is accumulating evidence that protonated nucleotides can stabilize RNA structure and participate in enzyme catalysis within ribozymes. Here, we calculate the pK(a) shifts of nucleotides in RNA structures using numerical solutions to the Poisson-Boltzmann equation. We find that significant shifts are predicted for several nucleotides in two catalytic RNAs, the hairpin ribozyme and the hepatitis delta virus ribozyme, and that the shifts are likely to be related to their functions. We explore how different structural environments shift the pK(a)s of nucleotides from their solution values. RNA structures appear to use two basic strategies to shift pK(a)s: (a) the formation of compact structural motifs with structurally-conserved, electrostatic interactions; and (b) the arrangement of the phosphodiester backbone to focus negative electrostatic potential in specific regions.     

Structure of the RNA Binding Domain of a DEAD-Box Helicase Bound to Its Ribosomal RNA Target Reveals a Novel Mode of Recognition by an RNA Recognition Motif

DEAD-box RNA helicases of the bacterial DbpA subfamily are localized to their biological substrate when a carboxy-terminal RNA recognition motif domain binds tightly and specifically to a segment of 23S ribosomal RNA (rRNA) that includes hairpin 92 of the peptidyl transferase center. A complex between a fragment of 23S rRNA and the RNA binding domain (RBD) of the Bacillus subtilis DbpA protein YxiN was crystallized and its structure was determined to 2.9 Å resolution, revealing an RNA recognition mode that differs from those observed with other RNA recognition motifs. The RBD is bound between two RNA strands at a three-way junction. Multiple phosphates of the RNA backbone interact with an electropositive band generated by lysines of the RBD. Nucleotides of the single-stranded loop of hairpin 92 interact with the RBD, including the guanosine base of G2553, which forms three hydrogen bonds with the peptide backbone. A G2553U mutation reduces the RNA binding affinity by 2 orders of magnitude, confirming that G2553 is a sequence specificity determinant in RNA binding. Binding of the RBD to 23S rRNA in the late stages of ribosome subunit maturation would position the ATP-binding duplex destabilization fragment of the protein for interaction with rRNA in the peptidyl transferase cleft of the subunit, allowing it to “melt out” unstable secondary structures and allow proper folding.