Prediction of pairing schemes in RNA molecules-loop contributions and energy of wobble and non-wobble pairs

Folding and finding RNA secondary structure

Optimal exploitation of the expanding database of sequences requires rapid finding and folding of RNAs. Methods are reviewed that automate folding and discovery of RNAs with algorithms that couple thermodynamics with chemical mapping, NMR, and/or sequence comparison. New functional noncoding RNAs in genome sequences can be found by combining sequence comparison with the assumption that functional noncoding RNAs will have more favorable folding free energies than other RNAs. When a new RNA is discovered, experiments and sequence comparison can restrict folding space so that secondary structure can be rapidly determined with the help of predicted free energies. In turn, secondary structure restricts folding in three dimensions, which allows modeling of three-dimensional structure. An example from a domain of a retrotransposon is described. Discovery of new RNAs and their structures will provide insights into evolution, biology, and design of therapeutics. Applications to studies of evolution are also reviewed.

Diversity of 16S rRNA genes within individual prokaryotic genomes

Analysis of intragenomic variation of 16S rRNA genes is a unique approach to examining the concept of ribosomal constraints on rRNA genes; the degree of variation is an important parameter to consider for estimation of the diversity of a complex microbiome in the recently initiated Human Microbiome Project (http://nihroadmap.nih.gov/hmp). The current GenBank database has a collection of 883 prokaryotic genomes representing 568 unique species, of which 425 species contained 2 to 15 copies of 16S rRNA genes per genome (2.22 +/- 0.81). Sequence diversity among the 16S rRNA genes in a genome was found in 235 species (from 0.06% to 20.38%; 0.55% +/- 1.46%). Compared with the 16S rRNA-based threshold for operational definition of species (1 to 1.3% diversity), the diversity was borderline (between 1% and 1.3%) in 10 species and >1.3% in 14 species. The diversified 16S rRNA genes in Haloarcula marismortui (diversity, 5.63%) and Thermoanaerobacter tengcongensis (6.70%) were highly conserved at the 2 degrees structure level, while the diversified gene in B. afzelii (20.38%) appears to be a pseudogene. The diversified genes in the remaining 21 species were also conserved, except for a truncated 16S rRNA gene in "Candidatus Protochlamydia amoebophila." Thus, this survey of intragenomic diversity of 16S rRNA genes provides strong evidence supporting the theory of ribosomal constraint. Taxonomic classification using the 16S rRNA-based operational threshold could misclassify a number of species into more than one species, leading to an overestimation of the diversity of a complex microbiome. This phenomenon is especially seen in 7 bacterial species associated with the human microbiome or diseases.

Complete suboptimal folding of RNA and the stability of secondary structures

An algorithm is presented for generating rigorously all suboptimal secondary structures between the minimum free energy and an arbitrary upper limit. The algorithm is particularly fast in the vicinity of the minimum free energy. This enables the efficient approximation of statistical quantities, such as the partition function or measures for structural diversity. The density of states at low energies and its associated structures are crucial in assessing from a thermodynamic point of view how well-defined the ground state is. We demonstrate this by exploring the role of base modification in tRNA secondary structures, both at the level of individual sequences from Escherichia coli and by comparing artificially generated ensembles of modified and unmodified sequences with the same tRNA structure. The two major conclusions are that (1) base modification considerably sharpens the definition of the ground state structure by constraining energetically adjacent structures to be similar to the ground state, and (2) sequences whose ground state structure is thermodynamically well defined show a significant tendency to buffer single point mutations. This can have evolutionary implications, since selection pressure to improve the definition of ground states with biological function may result in increased neutrality.

Prediction of sequentially optimal RNA secondary structures

A rigorous mathematical modeling of the RNA sequential folding process during transcription is proposed. It is based, at each transcription step, on a homogeneous markovian jump process, the state space of which is the set of structures constructible on the part of the RNA already transcribed. A theoretical formula permitting the computation of the structures probabilities at the end of the RNA transcription is derived. Successive approximations, aimed at reducing the size of the state space, permit the design of a prediction algorithm. The algorithm is tested on some structural RNAs (tRNA, 5S, 16S, hammerhead, ...), results are discussed and possible improvements are proposed.

Direct 5S rRNA Assay for Monitoring Mixed-Culture Bioprocesses

This study demonstrates the efficacy of a direct 5S rRNA assay for the characterization of mixed microbial populations by using as an example the bacteria associated with acidic mining environments. The direct 5S rRNA assay described herein represents a nonselective, direct molecular method for monitoring and characterizing the predominant, metabolically active members of a microbial population. The foundation of the assay is high-resolution denaturing gradient gel electrophoresis (DGGE), which is used to separate 5S rRNA species extracted from collected biomass. Separation is based on the unique migration behavior of each 5S rRNA species during electrophoresis in denaturing gradient gels. With mixtures of RNA extracted from laboratory cultures, the upper practical limit for detection in the current experimental system has been estimated to be greater than 15 different species. With this method, the resolution was demonstrated to be effective at least to the species level. The strength of this approach was demonstrated by the ability to discriminate between Thiobacillus ferrooxidans ATCC 19859 and Thiobacillus thiooxidans ATCC 8085, two very closely related species. Migration patterns for the 5S rRNA from members of the genus Thiobacillus were readily distinguishable from those of the genera Acidiphilium and Leptospirillum. In conclusion, the 5S rRNA assay represents a powerful method by which the structure of a microbial population within acidic environments can be assessed.

Chemical and computer probing of RNA structure

Ribonucleic acids (RNAs) are one of the most important types of biopolymers. RNAs play key roles in the storage and multiplication of genetic information. They are important in catalysis and RNA splicing and are the most important steps of translation. This chapter describes experimental methods for probing RNA structure and theoretical methods allowing the prediction of thermodynamically favorable RNA folding. These methods are complementary and together they provide a powerful approach to determine the structure of RNAs. The three-dimensional (tertiary) structure of RNA is formed by hydrogen-bonding among functional groups of nucleosides in different regions of the molecule, by coordination of polyvalent cations, and by stacking between the double-stranded regions present in the RNA. The tertiary structures of only some small RNAs have been determined by high-resolution X-ray crystallographic analysis and nuclear magnetic resonance analysis. The most widely used approach for the investigation of RNA structure is chemical and enzymatic probing, in combination with theoretical methods and phylogenetic studies allowing the prediction of variants of RNA folding. Investigations of RNA structures with different enzymatic and chemical probes can provide detailed data allowing the identification of double-stranded regions of the molecules and nucleotides involved in tertiary interactions.

A periodic table of symmetric tandem mismatches in RNA

The stabilities and structures of a series of RNA octamers containing symmetric tandem mismatches were studied by UV melting and imino proton NMR. The free energy increments for tandem mismatch formation are found to depend upon both mismatch sequence and adjacent base pairs. The observed sequence dependence of tandem mismatch stability is UGGU > GUUG > GAAG > or = AGGA > UUUU > CAAC > or = CUUC approximately UCCU approximately CCCC approximately ACCA approximately AAAA, and the closing base pair dependence is 5'G3'C > 5'C3'G > 5'U3'A approximately 5'A3'U. These results differ from expectations based on models used in RNA folding algorithms and from the sequence dependence observed for folding of RNA hairpins. Imino proton NMR results indicate the sequence dependence is partially due to hydrogen bonding within mismatches.

Potential secondary structure at the translational start domain of eukaryotic and prokaryotic mRNAs

In order to identify conserved potential secondary structures within translational start sites, mRNA sequences derived from different species were studied with programs able to depict such features. The potential secondary structure of 71 bases around the initiator AUG or AUGs in the coding sequences of 290 eukaryotic mRNAs was first examined and compared to 290 similarly analyzed regions derived from prokaryotic mRNA sequences (Nucleic Acids Res (1987) 15, 345-360). In both sets of sequences the initiator codon was often found to be in an open potential structure whereas a denser region characterized by nearly-periodic spacings defined the coding regions. Randomization of the sequences obliterated the observed patterns suggesting that the structure of the mRNA may determine these differences. Three sets of eukaryotic and prokaryotic mRNAs of approximately equal length were analyzed and found to preserve an open unpaired non-coding region 5' to the start codon. The start codon was found free of potential secondary structure in over 80% of all the sequences analyzed. These data, and study of mutants that restrict the accessibility of the start codon to the ribosomal initiation complex, suggest that both the prokaryotic and eukaryotic mRNA start sites must occur free of potential secondary structure for efficient initiation. A striking difference of the eukaryotic mRNA sequences analyzed was the high propensity of the coding region vicinal to the start codon to form secondary structures. Certain translation-defective mutants exhibit impaired formation of these secondary structures suggesting that the structure of the coding regions adjacent to the start codons of eukaryotic mRNAs may be an important, thus far unexamined, determinant of initiation. We propose that, for all genes studied, the transition in secondary structure between the coding and non-coding regions may be an important determinant of initiation.

Effect of deletions 5' to the translation initiation sequence on the expression of an mRNA in animal cells

To learn if an mRNA.18S rRNA interaction or a special secondary structure in the mRNA start region is essential for translation in eukaryotic cells, we constructed recombinant plasmids with the SV40 early promoter 5' to part of the Escherichia coli tufB-lacZ gene. Deletion of bases potentially complementary to the 18S rRNA highly increased the transient beta-galactosidase expressed in transfected CHO cells. Deletion of bases that fostered formation of potential hairpins with the mRNA 5'-terminus or altered the structure of the coding region reduced beta-galactosidase activity suggesting that these features of the mRNA secondary structure may be essential for initiation of translation. Computer aided analysis of the potential structure of 290 mRNAs suggests these are conserved features of the initiation region.

An apparent pause site in the transcription unit of the rabbit alpha-globin gene

Transcription of the rabbit alpha-globin gene begins primarily at the cap site, although some upstream start sites are also observed. Analysis by RNA polymerase run-on assays in nuclei shows that transcription continues at a high level past the polyadenylation site, after which the polymerase density actually increases in a region of about 400 nucleotides, followed by a gradual decline over the 700 nucleotides. These features are also observed in the transcription unit of the rabbit beta-globin gene. The region with the unexpectedly high nascent RNA hybridization signal in the 3' flank contains a conserved sequence, KGCAGCWGGR (K = G or T, W = A or T, R = A or G), followed by an inverted repeat. The inverted repeat (perhaps with the conserved sequence) may be a pause site for RNA polymerase II, thus accounting for the increase in polymerase density. This sequence and inverted repeat are found in the 3' flank of several globin genes and the simian virus 40 (SV40) early genes, as well as in the regions implicated in pausing or termination of transcription of eight different genes. Deletion of the conserved sequence and inverted repeat from the 3' flank of the SV40 early region causes a small increase in the levels of transcription downstream from this site. Replacement with the conserved sequence and inverted repeat from the rabbit alpha-globin gene causes an accumulation of polymerases, supporting the hypothesis that polymerases pause at this site. This proposed pause site may affect the efficiency of termination at some sites further downstream, perhaps by loss of a processivity factor.

A proposal for the secondary structure of a variable area of eukaryotic small ribosomal subunit RNA involving the existence of a pseudoknot

Eukaryotic small ribosomal subunit RNAs contain an area of variable structure, V4, which comprises about 250 nucleotides in most species, whereas the corresponding area in bacterial small ribosomal subunit RNAs consists of about 64 nucleotides folded into a single hairpin. There is no consensus on the secondary structure of area V4 in eukaryotes, about 10 different models having been proposed. The prediction of a model on a comparative basis poses special problems because, due to the variability of the area in length as well as sequence, a dependable alignment is very difficult to achieve. A new model was derived by systematic examination of all combinations of helices that have been hitherto proposed, plus some new ones. The following properties of the helices were examined: transposability to all presently known sequences, presence of compensating substitutions, and thermodynamic stability. A model was selected by ranking all possible combinations of transposable helices according to the number of compensating substitutions scored. The optimal model comprises a pseudoknot and four hairpin structures. Certain species contain additional hairpins inserted between these structural elements, while in others the structure is partially or entirely deleted.

Speeding up the dynamic algorithm for planar RNA folding

The simplest dynamic algorithm for planar RNA folding searches for the maximum number of base pairs. The algorithm uses O(n3) steps. The more general case, where different weights (energies) are assigned to stacked base pairs and to the various types of single-stranded region topologies, requires a considerably longer computation time because of the partial backtracking involved. Limiting the loop size reduces the running time back to O(n3). Reduction in the number of steps in the calculations of the various RNA topologies has recently been suggested, thereby improving the time behavior. Here we show how a "jumping" procedure can be used to speed up the computation, not only for the maximal number of base pairs algorithm, but for the minimal energy algorithm as well.

GC balance in the internal transcribed spacers ITS 1 and ITS 2 of nuclear ribosomal RNA genes

The internal transcribed spacer (ITS) 1 and 2, the 5.8S rRNA gene, and adjacent 18S rRNA and 25S rRNA coding regions of two Cucurbitaceae (Cucurbita pepo, zucchini, ITS 1: 187 bp, and ITS 2: 252 bp in length, and Cucumis sativus, cucumber, ITS 1: 229 bp, and ITS 2: 245 bp in length) have been sequenced. The evolutionary pattern shown by the ITSs of these plants is different from that found in vertebrates. Deletions, insertions, and base substitutions have occurred in both spacers; however, it is obvious that some selection pressure is responsible for the preservation of stem-loop structures. The dissimilarity of the 5' region of ITS 2 found in higher plants has consequences for proposed models on U3 snRNA-ITS 2 interaction in higher eukaryotes. The two investigated Cucurbitaceae species show a G + C content of ITS 1 that nearly equals that of ITS 2. An analysis of the ITS sequences reveals that in 19 out of 20 organisms published, the G + C content of ITS 1 nearly equals that of ITS 2, although it ranges from 20% to 90% in different organisms (GC balance). Moreover, the balanced G + C content of the ITSs in a given species seems to be similar to that of so-called expansion segments (ESs) in the 25/28S rRNA coding region. Thus, ITSs show a phenomenon called molecular coevolution with respect to each other and to the ESs. In the ITSs of Cucurbitaceae the balanced G + C composition is at least partly achieved by C to T transitions, via deamination of 5-methylcytosine. Other mutational events must be taken into account. The appearance of this phenomenon is discussed in terms of functional constraints linked to the structures of these spacers.

String analysis and energy minimization in the partition of DNA sequences

Two approaches to the understanding of biological sequences are confronted. While the recognition of particular signals in sequences relies on complex physical interactions, the problem is often analysed in terms of the presence or absence of literal motifs (strings) in the sequence. We present here a test-case for evaluating the potential of this approach. We classify DNA sequences as positive or negative depending on whether they contain a single melted domain in the middle of the sequence, which is a global physical property. Two sets of positive "biological" sequences were generated by a computer simulation of evolutionary divergence along the branches of a phylogenetic tree, under the constraint that each intermediate sequence be positive. These two sets and a set of random positive sequences were subjected to pattern analysis. The observed local patterns were used to construct expert systems to discriminate positive from negative sequences. The experts achieved 79% to 90% success on random positive sequences and up to 99% on the biological sets, while making less than 2% errors on negative sequences. Thus, the global constraints imposed on sequences by a physical process may generate local patterns that are sufficient to predict, with a reasonable probability, the behaviour of the sequences. However, rather large sets of biological sequences are required to generate patterns free of illegitimate constraints. Furthermore, depending upon the initial sequence, the sets of sequences generated on a phylogenetic tree may be amenable or refractory to string analysis, while obeying identical physical constraints. Our study clarifies the relationship between experts' errors on positive and negative sequences, and the contributions of legitimate and illegitimate patterns to these errors. The test-case appears suitable both for further investigations of problems in the theory of sequence evolution and for further testing of pattern analysis techniques.

Electrostatic potential of macromolecules measured by pKa shift of a fluorophore. 2. Transfer RNA

The procedures developed earlier (Friedrich and Woolley, preceding paper in this journal) for probing electrostatic potential with the fluorescein label were applied to transfer RNA. By using tRNA species that contain chemically reactive bases we were able to label these bases with fluorescein derivatives and thus to 'map' the electrostatic potential around the molecule. Both the electrostatic potential and the fluorescence emission anisotropy data that were obtained at the same time could be understood in terms of the well-known, paradigmatic crystal structure of tRNA(Phe). However, within the distribution of the various tRNA species, tRNA(Met)f appeared to occupy an extreme position, which suggests a relation between the conformation in solution and the initiation function of this molecule. Comparison with theoretical predictions by others of the electrostatic potential map of tRNA showed agreement in respect of trends, but the values of the potentials measured were orders of magnitude lower than predicted. This we attribute primarily to solvation.

Variable base pairing in a helix of eubacterial 5 S ribosomal RNA points to the existence of a conformational switch

A survey of 160 published sequences of eubacterial 5 S rRNAs shows that there exists structural variability in one of the helices of the generally accepted secondary structure model. Four structural variants are found, which differ with respect to the position and the number of bulges present. Most eubacterial 5 S RNAs fit into at least two of these conformations. A reaction scheme connecting the four observed conformations by changes in the base pairing scheme is proposed. Each of the known 5 S RNA sequences fits into conformations interconvertible by the proposed reactions.

Substrate recognition and identification of splice sites by the tRNA-splicing endonuclease and ligase from Saccharomyces cerevisiae

We have examined the substrate requirements for efficient and accurate splicing of tRNA precursors in Saccharomyces cerevisiae. The effects of Schizosaccharomyces pombe tRNASer gene mutations on the two steps in splicing, intron excision and joining of tRNA halves, were determined independently by using partially purified splicing endonuclease and tRNA ligase from S. cerevisiae. Two mutations (G14 and A46) reduced the efficiency of excision and joining in parallel, whereas two others (U47:7 and C33) produced differential effects on these two steps; U47:7 affected primarily the excision reaction, and C33 had a greater impact on ligation. These data indicate that endonuclease and ligase recognize both common and unique features of their substrates. Another two mutations (Ai26 and A37:13) induced miscutting, although with converse effects on the two splice sites. Thus, the two cutting events appear to be independent. Finally, we suggest that splice sites may be determined largely through their position relative to sites within the tRNA-like domain of the precursors. Several of these important sites were identified, and others are proposed based on the data described here.

Potential secondary structure at translation-initiation sites

Since translational start codons also occur internally, more-complex features within mRNA must determine initiation. We compare the potential secondary structure of 123 prokaryotic mRNA start regions to that of regions coding for internal methionines. The latter display an unexpectedly-uniform, almost-periodic pattern of pairing potential. In contrast, sequences 5' to start codons have little self-pairing, and do not pair extensively with the proximal coding region. Pairing potential surrounding start codons was found to be less than half of that found near internal AUGs. In groups of random sequences where the distribution of nucleotides at each position, or of trinucleotides at each in-frame codon position, matched the observed natural distribution, there was no periodicity in the pairing potential of the internal sequences. Randomized internal sequences had less pairing: the ratio of pairing intensity between internals and starts was reduced from 2.0 to 1.6 by randomization. We propose that the transition from the relatively-unstructured start domains to the highly-structured internal sequences may be an important determinant of translational start-site recognition.

Substrate recognition and identification of splice sites by the tRNA-splicing endonuclease and ligase from Saccharomyces cerevisiae

We have examined the substrate requirements for efficient and accurate splicing of tRNA precursors in Saccharomyces cerevisiae. The effects of Schizosaccharomyces pombe tRNASer gene mutations on the two steps in splicing, intron excision and joining of tRNA halves, were determined independently by using partially purified splicing endonuclease and tRNA ligase from S. cerevisiae. Two mutations (G14 and A46) reduced the efficiency of excision and joining in parallel, whereas two others (U47:7 and C33) produced differential effects on these two steps; U47:7 affected primarily the excision reaction, and C33 had a greater impact on ligation. These data indicate that endonuclease and ligase recognize both common and unique features of their substrates. Another two mutations (Ai26 and A37:13) induced miscutting, although with converse effects on the two splice sites. Thus, the two cutting events appear to be independent. Finally, we suggest that splice sites may be determined largely through their position relative to sites within the tRNA-like domain of the precursors. Several of these important sites were identified, and others are proposed based on the data described here.

Nuclease S1 analysis of eubacterial 5S rRNA secondary structure

Single-strand-specific nuclease S1 was employed as a structural probe to confirm locations of unpaired nucleotide bases in 5S rRNAs purified from prokaryotic species of rRNA superfamily I. Limited nuclease S1 digests of 3'- and 5'-end-labeled [32P]5S rRNAs were electrophoresed in parallel with reference endoribonuclease digests on thin sequencing gels. Nuclease S1 primary hydrolysis patterns were comparable for 5S rRNAs prepared from all 11 species examined in this study. The locations of base-paired regions determined by enzymatic analysis corroborate the general features of the proposed universal five-helix model for prokaryotic 5S rRNA, although the results of this study suggest a significant difference between prokaryotic and eukaryotic 5S rRNAs in the evolution of helix IV. Furthermore, the extent of base-pairing predicted by helix IV needs to be reevaluated for eubacterial species. Clipping patterns in helices II and IV appear to be consistent with a secondary structural model that undergoes a conformational rearrangement between two (or more) structures. Primary clipping patterns in the helix II region, obtained by S1 analysis, may provide useful information concerning the tertiary structure of the 5S rRNA molecule.

Temperature jump relaxation studies on the interactions between transfer RNAs with complementary anticodons. The effect of modified bases adjacent to the anticodon triplet

We have used the temperature-jump relaxation technique to determine the kinetic and thermodynamic parameters for the association between the following tRNAs pairs having complementary anticodons: tRNA(Ser) with tRNA(Gly), tRNA(Cys) with tRNA(Ala) and tRNA(Trp) with tRNA(Pro). The anticodon sequence of E. coli tRNA(Ser), GGA, is complementary to the U*CC anticodon of E. coli tRNA(Gly(2] (where U* is a still unknown modified uridine base) and A37 is not modified in none of these two tRNAs. E. coli tRNA(Ala) has a VGC anticodon (V is 5-oxyacetic acid uridine) while tRNA(Cys) has the complementary GCA anticodon with a modified adenine on the 3' side, namely 2-methylthio N6-isopentenyl adenine (mS2i6A37) in E. Coli tRNA(Cys) and N6-isopentenyl adenine (i6A37) in yeast tRNA(Cys). The brewer yeast tRNA(Trp) (anticodon CmCA) differs from the wild type E. coli tRNA(Trp) (anticodon CCA) in several positions of the nucleotide sequence. Nevertheless, in the anticodon loop, only two interesting differences are present: A37 is not modified while C34 at the first anticodon position is modified into a ribose 2'-O methyl derivative (Cm). The corresponding complementary tRNA is E.coli tRNA(Pro) with the VGG anticodon. Our results indicate a dominant effect of the nature and sequence of the anticodon bases and their nearest neighbor in the anticodon loop (particularly at position 37 on the 3' side); no detectable influence of modifications in the other tRNA stems has been detected. We found a strong stabilizing effect of the methylthio group on i6A37 as compared to isopentenyl modification of the same residue. We have not been able so far to assess the effect of isopentenyl modification alone in comparison to unmodified A37. The results obtained with the complex yeast tRNA(Trp)-E.coli tRNA(Pro) also suggest that a modification of C34 to Cm34 does not significantly increase the stability of tRNA(Trp) association with its complementary anticodon in tRNA(Pro). The observations are discussed in the light of inter- and intra-strand stacking interactions among the anticodon triplets and with the purine base adjacent to them, and of possible biological implications.

Base-base mismatches. Thermodynamics of double helix formation for dCA3XA3G + dCT3YT3G (X, Y = A,C,G,T)

Thermodynamic parameters for double strand formation have been measured for the sixteen double helices of the sequence dCA3XA3G.dCT3YT3G, with each of the bases A, C, G and T at the positions labelled X and Y. The results are analyzed in terms of nearest-neighbors and are compared with thermodynamic parameters for RNA secondary structure. At room temperature the sequence (Formula: see text) is more stable than (Formula: see text) and is similar in stability to (Formula: see text) and (Formula: see text) are least stable. At higher temperatures the sequences containing a G.C base pair become more stable than those containing only A.T. All molecules containing mismatches are destabilized with respect to those with only Watson-Crick pairing, but there is a wide range of destabilization. At room temperature the most stable mismatches are those containing guanine (G.T, G.G, G.A); the least stable contain cytosine (C.A, C.C). At higher temperatures pyrimidine-pyrimidine mismatches become the least stable.

Nucleotide sequence, secondary structure and evolution of the 5S ribosomal RNA from five bacterial species

The nucleotide sequences of the 5S ribosomal RNAs of the bacteria Agrobacterium tumefaciens, Alcaligenes faecalis, Pseudomonas cepacia, Aquaspirillum serpens and Acinetobacter calcoaceticus have been determined. The sequences fit in a generally accepted model for 5S RNA secondary structure. However, a closer comparative examination of these and other bacterial 5S RNA primary structures reveals the potential of additional base pairing and of multiple equilibria between a set of slightly different alternative secondary structures in one area of the molecule. The phylogenetic position of the examined bacteria is derived from a 5S RNA sequence alignment by a clustering method and compared with the position derived on the basis of 16S ribosomal RNA oligonucleotide catalogs.

[Prediction of secondary structures of nucleic acids: algorithmic and physical aspects]

Prediction of secondary structures in nucleic acids requires both an adequate physical model and powerful calculation algorithms. In our approach, we cut the molecules in sections of which the contributions to the global energy are context-dependent but roughly additive. The structure of minimum energy is obtained by a tree search under constraints of binary incompatibilities. Our algorithm of the "incompatibility islets" is shown to be more powerful than the "bit parallel forward checking" algorithm, well known in Artificial Intelligence. Recurrent algorithms, proposed by other authors are even more rapid, but often miss the correct structures, for they demand a strict additivity of the energetic contributions, physically unjustified. New strategies, required to deal with molecules of more than 200 nucleotides are discussed. Our physical model has been improved by considering the special case of internal loops beginning with a G-A opposition. A bonus of 1.5 kcal. is attributed to such a feature, at each side of an internal loop. To illustrate our programs, we give the computed schemes for the 3' termini of the small subunit ribosomal RNA.

A universal model for the secondary structure of 5.8S ribosomal RNA molecules, their contact sites with 28S ribosomal RNAs, and their prokaryotic equivalent

The phylogenetic approach (ref. 1) has been utilized in construction of a universal 5.8S rRNA secondary structure model, in which about 65% of the residues exist in paired structures. Conserved nucleotides primarily occupy unpaired regions. Multiple compensating base changes are demonstrated to be present in each of the five postulated helices, thereby forming a major basis for their proof. The results of chemical and enzymatic probing of 5.8S rRNAs (ref. 13, 32) are fully consistent with, and support, our model. This model differs in several ways from recently proposed 5.8S rRNA models (ref. 3, 4), which are discussed. Each of the helices in our model has been extended to the corresponding bacterial, chloroplast and mitochondrial sequences, which are demonstrated to be positionally conserved by alignment with their eukaryotic counterparts. This extension is also made for the base paired 5.8S/28S contact points, and their prokaryotic and organelle counterparts. The demonstrated identity of secondary structure in these diverse molecules strongly suggests that they perform equivalent functions in prokaryotic and eukaryotic ribosomes.

Equilibria in 5-S ribosomal RNA secondary structure. Bulges and interior loops in 5-S RNA secondary structure may serve as articulations for a flexible molecule

The basic assumption in this paper is that the secondary structure of a 5-S ribosomal RNA cannot be represented by a single model. We propose that the molecule can adopt, at least within the ribosome, a series of slightly different structures of nearly equal stability. The different structures arise from the existence of ambiguous base-pairing opportunities in bulged helices and the adjacent interior loops. In eubacterial 5-S RNAs there is one such an area, in eukaryotic 5-S RNAs two such areas that can give rise to structural switches. We explain how a change in secondary structure in these areas may influence the relative orientation of the surrounding helices, in other words how bulges and interior loops may serve as articulations and give rise to a flexible tertiary structure.

The nucleotide sequence of the cytoplasmic 5S rRNA from the horsetail, Equisetum arvense

Using 3'- and 5'-end labelling sequencing techniques, the following sequence for the cytoplasmic 5S rRNA of the horsetail Equisetum arvense could be determined: (sequence in text). This sequence exhibits all features expected for higher plant cytoplasmic 5S rRNAs, and can be fitted to the secondary structure model for 5S rRNA proposed by De Wachter et al. (15).

An energy model that predicts the correct folding of both the tRNA and the 5S RNA molecules

A new set of energy values to predict the secondary structures in RNA molecules has been derived through a multiple-step refinement procedure. It achieves more than 80% success in predicting the cloverleaf pattern in tRNA (200 sequences tested) and more than 60% success in predicting the consensus folding of 5S RNA (100 sequences). Improvements in our initial program for predicting secondary structures, based on the principle of the "incompatibility islets" made possible the work on 5S RNA. The program was speeded up by introducing a dynamic grouping of the islets into three disjoint blocks. The novel features in the energy model include i) an evaluation of the contribution of odd pairs according to their position within a segment ii) a penalty for internal loops related to their dissymmetry iii) a bonus for bulge loops when the two terminal paired bases at the junction point are both pyrimidines.

Detailed analysis of the higher-order structure of 16S-like ribosomal ribonucleic acids

Images

Secondary structure of the Dictyostelium discoideum small subunit ribosomal RNA

We have used comparative analyses of prokaryotic and eukaryotic small subunit ribosomal RNAs to deduce a secondary structure for the Dictyostelium discoideum 18S rRNA. Most of the duplex regions are evolutionarily conserved in all organisms. We have taken advantage of the variation to the D. discoideum sequence (relative to the yeast and frog 19S rRNAs) to identify additional helical regions which are common to the eukaryotic 18S rRNAs.

Nucleotide sequences of the 5.8S rRNAs of a mollusc and a porifer, and considerations regarding the secondary structure of 5.8S rRNA and its interaction with 28S rRNA

We report the primary structures of the 5.8 S ribosomal RNAs isolated from the sponge Hymeniacidon sanguinea and the snail Arion rufus. We had previously proposed (Ursi et al., Nucl. Acids Res. 10, 3517-3530 (1982)) a secondary structure model on the basis of a comparison of twelve 5.8 S RNA sequences then known, and a matching model for the interaction of 5.8 S RNA with 26 S RNA in yeast. Here we show that the secondary structure model can be extended to the 25 sequences presently available, and that the interaction model can be extended to the binding of 5.8 S RNA to the 5'-terminal domain of 28 S (26 S) RNA in three species.

Secondary structure of the Tetrahymena ribosomal RNA intervening sequence: structural homology with fungal mitochondrial intervening sequences

Splicing of the ribosomal RNA precursor of Tetrahymena is an autocatalytic reaction, requiring no enzyme or other protein in vitro. The structure of the intervening sequence (IVS) appears to direct the cleavage/ligation reactions involved in pre-rRNA splicing and IVS cyclization. We have probed this structure by treating the linear excised IVS RNA under nondenaturing conditions with various single- and double-strand-specific nucleases and then mapping the cleavage sites by using sequencing gel electrophoresis. A computer program was then used to predict the lowest-free-energy secondary structure consistent with the nuclease cleavage data. The resulting structure is appealing in that the ends of the IVS are in proximity; thus, the IVS can help align the adjacent coding regions (exons) for ligation, and IVS cyclization can occur. The Tetrahymena IVS has several sequences in common with those of fungal mitochondrial mRNA and rRNA IVSs, sequences that by genetic analysis are known to be important cis-acting elements for splicing of the mitochondrial RNAs. In the predicted structure of the Tetrahymena IVS, these sequences interact in a pairwise manner similar to that postulated for the mitochondrial IVSs. These findings suggest a common origin of some nuclear and mitochondrial introns and common elements in the mechanism of their splicing.

The nucleotide sequences of the 5S rRNAs of four mushrooms and their use in studying the phylogenetic position of basidiomycetes among the eukaryotes

The nucleotide sequences of the 5 S ribosomal RNAs of the mushrooms Russula cyanoxantha, Pleurotus ostreatus, Agaricus edulis, and Auricularia auricula-judae were determined. The sequences fit in a universal five-helix secondary structure model for 5 S RNA. As in most other 5 S RNAs, some helical areas contain non-standard base pairs. A clustering method was used to reconstruct an evolutionary tree from 82 eukaryotic 5 S RNA sequences. It allows to make a choice between alternative systematic classifications for basidiomycetes and reveals that the fungal kingdom is highly polyphyletic.

Computer building and folding of fictitious transfer-RNA sequences

In order to evaluate the common occurrence with which polynucleotides may adopt the cloverleaf configuration, 1150 random sequences were computer built and folded into their most stable secondary structure. Various constraints modulated the generation of the sequences: i) the base-pairing pattern, ii) the nucleotide composition, iii) the presence of assigned bases (modified or not) at certain sites, and iv) the chain length. In many cases, artificial tRNAs appear to require a more complex organization than a cloverleaf pairing scheme to achieve, as do natural molecules, the corresponding secondary structure. Moreover, the preferred foldings of sequences from 50 to 90 nucleotide long without an imposed pairing pattern usually contain two rather than three hairpin-loops. Implications concerning the emergence and the evolution of the protein-synthesis apparatus are discussed.

Sequences of the 5S rRNAs of Azotobacter vinelandii, Pseudomonas aeruginosa and Pseudomonas fluorescens with some notes on 5S RNA secondary structure

Recently published alignments of available 5 S rRNA sequences have shown that a rigid base pairing pattern, pointing to the existence of a universal five-helix secondary structure for all 5 S RNAs, can be superimposed on such alignments. For a few species, the alignment and the base pairing pattern show distortions with respect to the large majority of sequences. Their 5 S RNAs may form exceptional secondary structures, or there may just be errors in the published sequences. We have examined such a case, Pseudomonas fluorescens, and found the sequence to be in error. The corrected sequence, as well as those of the related species Azotobacter vinelandii and Pseudomonas aeruginosa, fit perfectly in the 5 S RNA sequence alignment and in the five-helix secondary structure model. There exists comparative evidence for the frequent presence of non-standard base pairs at several points of the 5 S RNA secondary structure.

Comparison of fungal mitochondrial introns reveals extensive homologies in RNA secondary structure

The complete sequences of nine Saccharomyces cerevisiae mitochondrial introns, six of which carry long open reading frames, have already been published. We have recently determined the sequence of an intron in the large ribosomal mitochondrial RNA of Kluyveromyces thermotolerans (Jacquier et al., in preparation), which we found to be closely related to its S. cerevisiae counterpart. This latter result prompted us to undertake a systematic search for possible homologous elements in the other, available sequences with the help of an original computer program. A previously unsuspected wealth of evolutionarily conserved sequences and secondary structures was thus uncovered. Seven at least of the available sequences may be folded up into elaborate secondary structure models, the cores of which are nearly identical. These models result in bringing together the exon-intron junctions into relatively close spatial proximity and looping out either all or most of the sequences in open reading frame, when present. These results and their possible implications with respect to the mechanism of splicing are discussed in the light of available genetic and biochemical data.