Using an RNA secondary structure partition function to determine confidence in base pairs predicted by free energy minimization

A partition function calculation for RNA secondary structure is presented that uses a current set of nearest neighbor parameters for conformational free energy at 37 degrees C, including coaxial stacking. For a diverse database of RNA sequences, base pairs in the predicted minimum free energy structure that are predicted by the partition function to have high base pairing probability have a significantly higher positive predictive value for known base pairs. For example, the average positive predictive value, 65.8%, is increased to 91.0% when only base pairs with probability of 0.99 or above are considered. The quality of base pair predictions can also be increased by the addition of experimentally determined constraints, including enzymatic cleavage, flavin mono-nucleotide cleavage, and chemical modification. Predicted secondary structures can be color annotated to demonstrate pairs with high probability that are therefore well determined as compared to base pairs with lower probability of pairing.     

Improved RNA secondary structure prediction by maximizing expected pair accuracy

Free energy minimization has been the most popular method for RNA secondary structure prediction for decades. It is based on a set of empirical free energy change parameters derived from experiments using a nearest-neighbor model. In this study, a program, MaxExpect, that predicts RNA secondary structure by maximizing the expected base-pair accuracy, is reported. This approach was first pioneered in the program CONTRAfold, using pair probabilities predicted with a statistical learning method. Here, a partition function calculation that utilizes the free energy change nearest-neighbor parameters is used to predict base-pair probabilities as well as probabilities of nucleotides being single-stranded. MaxExpect predicts both the optimal structure (having highest expected pair accuracy) and suboptimal structures to serve as alternative hypotheses for the structure. Tested on a large database of different types of RNA, the maximum expected accuracy structures are, on average, of higher accuracy than minimum free energy structures. Accuracy is measured by sensitivity, the percentage of known base pairs correctly predicted, and positive predictive value (PPV), the percentage of predicted pairs that are in the known structure. By favoring double-strandedness or single-strandedness, a higher sensitivity or PPV of prediction can be favored, respectively. Using MaxExpect, the average PPV of optimal structure is improved from 66% to 68% at the same sensitivity level (73%) compared with free energy minimization.     

A set of nearest neighbor parameters for predicting the enthalpy change of RNA secondary structure formation

A complete set of nearest neighbor parameters to predict the enthalpy change of RNA secondary structure formation was derived. These parameters can be used with available free energy nearest neighbor parameters to extend the secondary structure prediction of RNA sequences to temperatures other than 37°C. The parameters were tested by predicting the secondary structures of sequences with known secondary structure that are from organisms with known optimal growth temperatures. Compared with the previous set of enthalpy nearest neighbor parameters, the sensitivity of base pair prediction improved from 65.2 to 68.9% at optimal growth temperatures ranging from 10 to 60°C. Base pair probabilities were predicted with a partition function and the positive predictive value of structure prediction is 90.4% when considering the base pairs in the lowest free energy structure with pairing probability of 0.99 or above. Moreover, a strong correlation is found between the predicted melting temperatures of RNA sequences and the optimal growth temperatures of the host organism. This indicates that organisms that live at higher temperatures have evolved RNA sequences with higher melting temperatures.     

Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure.

An improved dynamic programming algorithm is reported for RNA secondary structure prediction by free energy minimization. Thermodynamic parameters for the stabilities of secondary structure motifs are revised to include expanded sequence dependence as revealed by recent experiments. Additional algorithmic improvements include reduced search time and storage for multibranch loop free energies and improved imposition of folding constraints. An extended database of 151,503 nt in 955 structures? determined by comparative sequence analysis was assembled to allow optimization of parameters not based on experiments and to test the accuracy of the algorithm. On average, the predicted lowest free energy structure contains 73 % of known base-pairs when domains of fewer than 700 nt are folded; this compares with 64 % accuracy for previous versions of the algorithm and parameters. For a given sequence, a set of 750 generated structures contains one structure that, on average, has 86 % of known base-pairs. Experimental constraints, derived from enzymatic and flavin mononucleotide cleavage, improve the accuracy of structure predictions.     

Fast and accurate determination of sites along the FUT2 in vitro transcript that are accessible to antisense oligonucleotides by application of secondary structure predictions and RNase H in combination with MALDI-TOF mass spectrometry

Alteration of gene expression by use of antisense oligonucleotides has considerable potential for therapeutic purposes and scientific studies. Although applied for almost 25 years, this technique is still associated with difficulties in finding antisense-effective regions along the target mRNA. This is mainly due to strong secondary structures preventing binding of antisense oligonucleotides and RNase H, playing a major role in antisense-mediated degradation of the mRNA. These difficulties make empirical testing of a large number of sequences complementary to various sites in the target mRNA a very lengthy and troublesome procedure. To overcome this problem, more recent strategies to find efficient antisense sites are based on secondary structure prediction and RNase H-dependent mechanisms. We were the first who directly combined these two strategies; antisense oligonucleotides complementary to predicted unpaired target mRNA regions were designed and hybridized to the corresponding RNAs. Incubation with RNase H led to cleavage of the RNA at the respective hybridization sites. Analysis of the RNA fragments by matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) mass spectrometry, which has not been used in this context before, allowed exact determination of the cleavage site. Thus the technique described here is very promising when searching for effective antisense sites.     

Thermodynamic characterization of the complete set of sequence symmetric tandem mismatches in RNA and an improved model for predicting the free energy contribution of sequence asymmetric tandem mismatches

Because of the availability of an abundance of RNA sequence information, the ability to rapidly and accurately predict the secondary structure of RNA from sequence is becoming increasingly important. A common method for predicting RNA secondary structure from sequence is free energy minimization. Therefore, accurate free energy contributions for every RNA secondary structure motif are necessary for accurate secondary structure predictions. Tandem mismatches are prevalent in naturally occurring sequences and are biologically important. A common method for predicting the stability of a sequence asymmetric tandem mismatch relies on the stabilities of the two corresponding sequence symmetric tandem mismatches [Mathews, D. H., Sabina, J., Zuker, M., and Turner, D. H. (1999) J. Mol. Biol. 288, 911-940]. To improve the prediction of sequence asymmetric tandem mismatches, the experimental thermodynamic parameters for the 22 previously unmeasured sequence symmetric tandem mismatches are reported. These new data, however, do not improve prediction of the free energy contributions of sequence asymmetric tandem mismatches. Therefore, a new model, independent of sequence symmetric tandem mismatch free energies, is proposed. This model consists of two penalties to account for destabilizing tandem mismatches, two bonuses to account for stabilizing tandem mismatches, and two penalties to account for A-U and G-U adjacent base pairs. This model improves the prediction of asymmetric tandem mismatch free energy contributions and is likely to improve the prediction of RNA secondary structure from sequence.     

Hybridization specificity, enzymatic activity and biological (Ha-ras) activity of oligonucleotides containing 2,4-dideoxy-beta-D-erythro-hexopyranosyl nucleosides.

Antisense oligonucleotides with a 2,4-dideoxyhexopyranosyl nucleoside incorporated at the 3'-end and at a mutation site of the Ha-ras oncogene mRNA were synthesized. Melting temperature studies revealed that an A*-G mismatch is more stable than an A*-T mismatch with these hexopyranosyl nucleosides incorporated at the mutation site. The oligonucleotides are stable against enzymatic degradation. RNase H mediated cleavage studies revealed selective cleavage of mutated Ha-ras mRNA. The oligonucleotide containing two pyranose nucleosides at the penultimate position activates RNase H more strongly than natural oligonucleotides. No correlation, however, was found between DNA - DNA or RNA - DNA melting temperatures and RNase H mediated cleavage capacity. Although the A*-G mismatch gives more stable hybridization than the A*-T base pairing, only the oligonucleotides containing an A*-T base pair are recognized by RNase H. This modification is situated 3 base pairs upstream to the cleavage site. Finally, the double pyranose modified oligonucleotide was able to reduce the growth of T24 cells (bladder carcinoma) while the unmodified antisense oligonucleotide was not.     

Prediction of RNA secondary structure by free energy minimization

RNA secondary structure is often predicted from sequence by free energy minimization. Over the past two years, advances have been made in the estimation of folding free energy change, the mapping of secondary structure and the implementation of computer programs for structure prediction. The trends in computer program development are: efficient use of experimental mapping of structures to constrain structure prediction; use of statistical mechanics to improve the fidelity of structure prediction; inclusion of pseudoknots in secondary structure prediction; and use of two or more homologous sequences to find a common structure.     

Predicting a set of minimal free energy RNA secondary structures common to two sequences

MOTIVATION: Function derives from structure, therefore, there is need for methods to predict functional RNA structures. RESULTS: The Dynalign algorithm, which predicts the lowest free energy secondary structure common to two unaligned RNA sequences, is extended to the prediction of a set of low-energy structures. Dot plots can be drawn to show all base pairs in structures within an energy increment. Dynalign predicts more well-defined structures than structure prediction using a single sequence; in 5S rRNA sequences, the average number of base pairs in structures with energy within 20% of the lowest energy structure is 317 using Dynalign, but 569 using a single sequence. Structure prediction with Dynalign can also be constrained according to experiment or comparative analysis. The accuracy, measured as sensitivity and positive predictive value, of Dynalign is greater than predictions with a single sequence. AVAILABILITY: Dynalign can be downloaded at http://rna.urmc.rochester.edu     

Anti-Influenza Virus Activities of Nicked and Circular Dumbbell RNA/DNA Chimeric Oligonucleotides

Abstract

We have designed a new type of antisense oligonucleotide, containing two hairpin loop structures with RNA/DNA base pairs (sense (RNA) and antisense (DNA)) in the double helical stem (nicked and circular dumbbell DNA/RNA chimeric oligonucleotides). The reaction of the nicked and circular dumbbell DNA/RNA chimeric oligonucleotides with RNase H gave the corresponding anti-DNA together with the sense RNA cleavage products. These oligonucleotides were more resistant to exonuclease attack. We also describe the anti-Fluv activities of nicked and circular dumbbell DNMA chimeric oligonucleotides.     

Identification of c-myc coding region determinant RNA sequences and structures cleaved by an RNase1-like endoribonuclease

The coding region of c-myc mRNA encompassing the coding region determinant (CRD) nucleotides (nts) 1705-1792 is critical in regulating c-myc mRNA stability. This is in part due to the susceptibility of c-myc CRD RNA to attack by an endoribonuclease. We have previously purified and characterized a mammalian endoribonuclease that cleaves c-myc CRD RNA in vitro. This enzyme is tentatively identified as a 35 kDa RNase1-like endonuclease. In an effort to understand the sequence and secondary structure requirements for RNA cleavage by this enzyme, we have determined the secondary structure of the c-myc CRD RNA nts 1705-1792 using RNase probing technique. The secondary structure of c-myc CRD RNA possesses five stems; two of which contain 4 base pairs (stems I and V) and three consisting of 3 base pairs (stems II, III, and IV). Endonucleolytic assays using the c-myc CRD and several c-myc CRD mutants as substrates led to the following conclusions: (i) the enzyme prefers to cleave in between the dinucleotides UA, CA, and UG in single-stranded regions; (ii) the enzyme is more specific towards UA dinucleotides. These properties further distinguish the enzyme from previously described mammalian endonuclease that cleaves c-myc mRNA in vitro.     

Secondary structure models of D2-D3 expansion segments of 28S rRNA for Hoplolaiminae species

The D2-D3 expansion segments of the 28S ribosomal RNA (rRNA) were sequenced and compared to predict secondary structures for Hoplolaiminae species based on free energy minimization and comparative sequence analysis. The free energy based prediction method provides putative stem regions within primary structure and these base pairings in stems were confirmed manually by compensatory base changes among closely and distantly related species. Sequence differences ranged from identical between Hoplolaimus columbus and H. seinhorsti to 20.8% between Scutellonema brachyurum and H. concaudajuvencus. The comparative sequence analysis and energy minimization method yielded 9 stems in the D2 and 6 stems in the D3 which showed complete or partial compensatory base changes. At least 75% of nucleotides in the D2 and 68% of nucleotides in the D3 were related with formation of base pairings to maintain secondary structure. GC contents in stems ranged from 61 to 73% for the D2 and from 64 to 71% for the D3 region. These ranges are higher than G-C contents in loops which ranged from 37 to 48% in the D2 and 33-45% in the D3. In stems, G-C/C-G base pairings were the most common in the D2 and the D3 and also non-canonical base pairs including A•A and U•U, C•U/U•C, and G•A/A•G occurred in stems. The predicted secondary model and new sequence alignment based on predicted secondary structures for the D2 and D3 expansion segments provide useful information to assign positional nucleotide homology and reconstruction of more reliable phylogenetic trees.     

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.     

Properties of cloned and expressed human RNase H1.

We have characterized cloned His-tag human RNase H1. The activity of the enzyme exhibited a bell-shaped response to divalent cations and pH. The optimum conditions for catalysis consisted of 1 mM Mg(2+) and pH 7-8. In the presence of Mg(2+), Mn(2+) was inhibitory. Human RNase H1 shares many enzymatic properties with Escherichia coli RNase H1. The human enzyme cleaves RNA in a DNA-RNA duplex resulting in products with 5'-phosphate and 3'-hydroxy termini, can cleave overhanging single strand RNA adjacent to a DNA-RNA duplex, and is unable to cleave substrates in which either the RNA or DNA strand has 2' modifications at the cleavage site. Human RNase H1 binds selectively to "A-form"-type duplexes with approximately 10-20-fold greater affinity than that observed for E. coli RNase H1. The human enzyme displays a greater initial rate of cleavage of a heteroduplex-containing RNA-phosphorothioate DNA than an RNA-DNA duplex. Unlike the E. coli enzyme, human RNase H1 displays a strong positional preference for cleavage, i.e. it cleaves between 8 and 12 nucleotides from the 5'-RNA-3'-DNA terminus of the duplex. Within the preferred cleavage site, the enzyme displays modest sequence preference with GU being a preferred dinucleotide. The enzyme is inhibited by single-strand phosphorothioate oligonucleotides and displays no evidence of processivity. The minimum RNA-DNA duplex length that supports cleavage is 6 base pairs, and the minimum RNA-DNA "gap size" that supports cleavage is 5 base pairs.     

Predicted structures of apolipoprotein II mRNA constrained by nuclease and dimethyl sulfate reactivity: stable secondary structures occur predominantly in local domains via intraexonic base pairing

Analyses of apolipoprotein II mRNA with chemical and enzymatic probes showed that double- and single-stranded regions were distributed uniformly along the mRNA except for a large (72 nucleotides) single-stranded region containing the translation stop codon. Secondary structure models constrained by the experimental data were made by varying the distance (along the mRNA) over which base pairing was allowed. Four prominent secondary structures were seen with restrictions of 165, 330, or 659 nucleotides suggesting that such structures from via local interactions over distances of 50-120 nucleotides. Predicted long range interactions involve only 2-3 base pairs while local interactions involve helices of 4-10 base pairs. Predicted helices of greater than or equal to 4 base pairs occur primarily within exons, raising the possibility that prominent secondary structures in mRNAs may be largely due to intraexonic base pairing. Tests of single- and double-stranded domains by oligonucleotide-directed RNase H cleavage and primer extension were in accord with the structure model and with nuclease and chemical modification data. The model predicting base pairing between the coding and the 3' noncoding regions was tested by RNase H cleavage followed by oligo(dT)-cellulose chromatography to separate 5' and 3' mRNA fragments. Most (82%) of the 5' fragment remained associated with the 3' noncoding region in a structure with a tm = 50 degrees C in 0.2 M Na+ suggesting that this stem could be stable in vivo. This stem may be stable in the isolated mRNA, but would likely occur transiently in polyribosomal apolipoprotein II mRNA due to ribosome transit through the 5' side of the stem. Alternate structures may occur in this region during ribosome transit and play a role in translation termination or in determining the susceptibility of the mRNA to degradation.     

Targeting Highly Structured RNA by Cooperative Action of siRNAs and Helper Antisense Oligomers in Living Cells

RNA target accessibility is one of the most important factors limiting the efficiency of RNA interference-mediated RNA degradation. However, targeting RNA viruses in their poorly accessible, highly structured regions can be advantageous because these regions are often conserved in sequence and thus less prone to viral escape. We developed an experimental strategy to attack highly structured RNA by means of pairs of specifically designed small interfering RNAs and helper antisense oligonucleotides using the 5’ untranslated region (5’UTR) of coxsackievirus B3 as a model target. In the first step, sites accessible to hybridization of complementary oligonucleotides were identified using two mapping methods with random libraries of short DNA oligomers. Subsequently, the accessibility of the mapped regions for hybridization of longer DNA 16-mers was confirmed by an RNase H assay. Using criteria for the design of efficient small interfering RNAs (siRNA) and a secondary structure model of the viral 5’UTR, several DNA 19-mers were designed against partly double-stranded RNA regions. Target sites for DNA 19-mers were located opposite the sites which had been confirmed as accessible for hybridization. Three pairs of DNA 19-mers and the helper 2’-O-methyl-16-mers were able to effectively induce RNase H cleavage in vitro. For cellular assays, the DNA 19-mers were replaced by siRNAs, and the corresponding three pairs of siRNA-helper oligomer tools were found to target 5’UTR efficiently in a reporter construct in HeLa cells. Addition of the helper oligomer improved silencing capacity of the respective siRNA. We assume that the described procedure will generally be useful for designing of nucleic acid-based tools to silence highly structured RNA targets.     

Structural constraints identified with covariation analysis in ribosomal RNA

Covariation analysis is used to identify those positions with similar patterns of sequence variation in an alignment of RNA sequences. These constraints on the evolution of two positions are usually associated with a base pair in a helix. While mutual information (MI) has been used to accurately predict an RNA secondary structure and a few of its tertiary interactions, early studies revealed that phylogenetic event counting methods are more sensitive and provide extra confidence in the prediction of base pairs. We developed a novel and powerful phylogenetic events counting method (PEC) for quantifying positional covariation with the Gutell lab's new RNA Comparative Analysis Database (rCAD). The PEC and MI-based methods each identify unique base pairs, and jointly identify many other base pairs. In total, both methods in combination with an N-best and helix-extension strategy identify the maximal number of base pairs. While covariation methods have effectively and accurately predicted RNAs secondary structure, only a few tertiary structure base pairs have been identified. Analysis presented herein and at the Gutell lab's Comparative RNA Web (CRW) Site reveal that the majority of these latter base pairs do not covary with one another. However, covariation analysis does reveal a weaker although significant covariation between sets of nucleotides that are in proximity in the three-dimensional RNA structure. This reveals that covariation analysis identifies other types of structural constraints beyond the two nucleotides that form a base pair.     

Facilitating RNA structure prediction with microarrays

Determining RNA secondary structure is important for understanding structure-function relationships and identifying potential drug targets. This paper reports the use of microarrays with heptamer 2'-O-methyl oligoribonucleotides to probe the secondary structure of an RNA and thereby improve the prediction of that secondary structure. When experimental constraints from hybridization results are added to a free-energy minimization algorithm, the prediction of the secondary structure of Escherichia coli 5S rRNA improves from 27 to 92% of the known canonical base pairs. Optimization of buffer conditions for hybridization and application of 2'-O-methyl-2-thiouridine to enhance binding and improve discrimination between AU and GU pairs are also described. The results suggest that probing RNA with oligonucleotide microarrays can facilitate determination of secondary structure.     

Secondary structure of the RNA component of a nuclear/mitochondrial ribonucleoprotein.

RNase mitochondrial RNA processing (MRP) is a site-specific endoribonuclease located in both the nucleus and mitochondria of vertebrate cells. The enzyme is a ribonucleoprotein whose RNA component has been shown to be encoded by a nuclear gene. Because RNase MRP is particular in its substrate requirement, RNA-RNA interaction has been proposed as important for the cleavage reaction. A secondary structure of this RNA from mouse cells has been derived by chemical modification of in vivo MRP RNA in ribonucleoprotein form, as isolated free RNA, and as RNA synthesized in vitro. Full-length MRP RNA appears to adopt a conformation containing a significant number of single-stranded residues and may form a pseudoknot. The data are consistent with both the RNA within the ribonucleoprotein and the free RNA possessing comparable secondary structures and suggest a possible site of interaction between enzyme and substrate. The human MRP RNA can be folded into a conformation very similar to that predicted for the mouse MRP RNA. A more limited analysis of human MRP RNA is consistent with the structure proposed for the mouse species.