Computer-aided prediction of RNA secondary structures

General RNA-binding proteins have a function in poly(A)-binding protein-dependent translation

The interaction between the poly(A)-binding protein (PABP) and eukaryotic translational initiation factor 4G (eIF4G), which brings about circularization of the mRNA, stimulates translation. General RNA-binding proteins affect translation, but their role in mRNA circularization has not been studied before. Here, we demonstrate that the major mRNA ribonucleoprotein YB-1 has a pivotal function in the regulation of eIF4F activity by PABP. In cell extracts, the addition of YB-1 exacerbated the inhibition of 80S ribosome initiation complex formation by PABP depletion. Rabbit reticulocyte lysate in which PABP weakly stimulates translation is rendered PABP-dependent after the addition of YB-1. In this system, eIF4E binding to the cap structure is inhibited by YB-1 and stimulated by a nonspecific RNA. Significantly, adding PABP back to the depleted lysate stimulated eIF4E binding to the cap structure more potently if this binding had been downregulated by YB-1. Conversely, adding nonspecific RNA abrogated PABP stimulation of eIF4E binding. These data strongly suggest that competition between YB-1 and eIF4G for mRNA binding is required for efficient stimulation of eIF4F activity by PABP.

Capped mRNAs with reduced secondary structure can function in extracts from poliovirus-infected cells

Extracts from poliovirus-infected HeLa cells were used to study ribosome binding of native and denatured reovirus mRNAs and translation of capped mRNAs with different degrees of secondary structure. Here, we demonstrate that ribosomes in extracts from poliovirus-infected cells could form initiation complexes with denatured reovirus mRNA, in contrast to their inability to bind native reovirus mRNA. Furthermore, the capped alfalfa mosaic virus 4 RNA, which is most probably devoid of stable secondary structure at its 5' end, could be translated at much higher efficiency than could other capped mRNAs in extracts from poliovirus-infected cells.

Capped mRNAs with reduced secondary structure can function in extracts from poliovirus-infected cells

Extracts from poliovirus-infected HeLa cells were used to study ribosome binding of native and denatured reovirus mRNAs and translation of capped mRNAs with different degrees of secondary structure. Here, we demonstrate that ribosomes in extracts from poliovirus-infected cells could form initiation complexes with denatured reovirus mRNA, in contrast to their inability to bind native reovirus mRNA. Furthermore, the capped alfalfa mosaic virus 4 RNA, which is most probably devoid of stable secondary structure at its 5' end, could be translated at much higher efficiency than could other capped mRNAs in extracts from poliovirus-infected cells.

The Therapeutic Potential of Ribozymes

Ribozymes are catalytic RNA molecules that recognize their target RNA in a highly sequence-specific manner. They can therefore be used to inhibit deleterious gene expression (by cleavage of the target mRNA) or even repair mutant cellular RNAs. Targets such as the mRNAs of oncogenes (resulting from base mutations or chromosome translocations, eg, ras or bcr-abl) and viral genomes and transcripts (human immunodeficiency virus–type 1 [HIV-1]) are ideal targets for such sequence-specific agents. The aim of this review is therefore to introduce the different classes of ribozymes, highlighting some of the chemistry of the reactions they catalyze, to address the specific inhibition of genes by ribozymes, the problems yet to be resolved, and how new developments in the field give hope to the future for ribozymes in the therapeutic field.

The Therapeutic Potential of Ribozymes

Ribozymes are catalytic RNA molecules that recognize their target RNA in a highly sequence-specific manner. They can therefore be used to inhibit deleterious gene expression (by cleavage of the target mRNA) or even repair mutant cellular RNAs. Targets such as the mRNAs of oncogenes (resulting from base mutations or chromosome translocations, eg, ras or bcr-abl) and viral genomes and transcripts (human immunodeficiency virus–type 1 [HIV-1]) are ideal targets for such sequence-specific agents. The aim of this review is therefore to introduce the different classes of ribozymes, highlighting some of the chemistry of the reactions they catalyze, to address the specific inhibition of genes by ribozymes, the problems yet to be resolved, and how new developments in the field give hope to the future for ribozymes in the therapeutic field.

The equilibrium partition function and base pair binding probabilities for RNA secondary structure

A novel application of dynamic programming to the folding problem for RNA enables one to calculate the full equilibrium partition function for secondary structure and the probabilities of various substructures. In particular, both the partition function and the probabilities of all base pairs are computed by a recursive scheme of polynomial order N3 in the sequence length N. The temperature dependence of the partition function gives information about melting behavior for the secondary structure. The pair binding probabilities, the computation of which depends on the partition function, are visually summarized in a "box matrix" display and this provides a useful tool for examining the full ensemble of probable alternative equilibrium structures. The calculation of this ensemble representation allows a proper application and assessment of the predictive power of the secondary structure method, and yields important information on alternatives and intermediates in addition to local information about base pair opening and slippage. The results are illustrated for representative tRNA, 5S RNA, and self-replicating and self-splicing RNA molecules, and allow a direct comparison with enzymatic structure probes. The effect of changes in the thermodynamic parameters on the equilibrium ensemble provides a further sensitivity check to the predictions.

In vitro secondary structure analysis of mRNA from lacZ translation initiation mutants

mRNA secondary structure can be an important determinant of the efficiency of translation initiation. To study the effect of secondary structure on translation initiation, in vitro secondary structure analysis was performed on 32 lacZ RNA transcripts that differ in their in vivo translation initiation efficiencies because of mutations. We have shown that well-translated RNA has a relatively unstructured translation initiation region in vitro. In contrast, the translation initiation region of many of the poorly translated RNA transcripts is involved in a stem-loop structure. Mutations that decrease the in vitro stability of the stem-loop increase the frequency of translation initiation. The sequences responsible for forming this stem-loop structure were localized to a small region of RNA. The results confirm some of the previous predictions of the RNA secondary structure of the mutant RNAs based on computer modeling, but they disagree with some of the predicted long-range interactions.

Human promyelocytic leukemia HL-60 cell proliferation and c-myc protein expression are inhibited by an antisense pentadecadeoxynucleotide targeted against c-myc mRNA

The human promyelocytic leukemia cell line HL-60 overexpresses the c-myc protooncogene. A calculated secondary structure for c-myc mRNA placed the initiation codon in a bulge of a weakly base-paired region. Treatment of HL-60 cells with 5' d(AACGTTGAGGGGCAT) 3', complementary to the initiation codon and the next four codons of c-myc mRNA, inhibited c-myc protein expression in a dose-dependent manner. However, treatment of HL-60 cells with 5' d(TTGGGATAACACTTA) 3', complementary to nucleotides 17-31 of vesicular stomatitis virus matrix protein mRNA, displayed no such effects. These results agree with analogous studies of normal human T lymphocytes [Heikkila, R., Schwab, G., Wickstrom, E., Loke, S. L., Pluznik, D. H., Watt, R. & Neckers, L. M. (1987) Nature (London) 328, 445-449], except that only one-third as much oligomer was needed for a comparable effect. Proliferation of HL-60 cells in culture was inhibited in a sequence-specific, dose-dependent manner by the c-myc-complementary oligomer, but neither the oligomer complementary to vesicular stomatitis virus matrix protein mRNA nor 5' d(CATTTCTTGCTCTCC) 3', complementary to nucleotides 5399-5413 of human immunodeficiency virus tat gene mRNA, inhibited proliferation. It thus appears that antisense oligodeoxynucleotides added to myc-transformed cells via culture medium are capable of eliciting sequence-specific, dose-dependent inhibition of c-myc protein expression and cell proliferation.

Enzymatic amplification of translation inhibition of rabbit beta-globin mRNA mediated by anti-messenger oligodeoxynucleotides covalently linked to intercalating agents

The effects of anti-messenger oligodeoxynucleotides, covalently linked to an intercalating agent, on translation of rabbit beta-globin mRNA, were investigated both in wheat germ extract and in microinjected Xenopus oocytes. A specific inhibition of beta-globin synthesis was observed in both expression systems with a modified 11-mer covalently linked to an acridine derivative. In injected oocytes a more efficient block was observed with this modified oligonucleotide than with its unsubstituted homolog. This was ascribed to stacking interactions of the intercalating agent with base pairs which provide an additional stabilization of the [mRNA/DNA] hybrid. We demonstrated that in wheat germ extract, the modified and unmodified oligonucleotides behaved similarly due to the presence of a high RNaseH activity. RNaseH was also present, although to a lesser extent, in the oocyte cytoplasm. This anti-messenger DNA-induced degradation of target mRNA resulted in amplified efficiency of hybrid-arrested translation. This additional mechanism might provide anti-sense DNAs with an advantage over anti-sense RNAs.

A computer model of evolutionary optimization

Molecular evolution is viewed as a typical combinatorial optimization problem. We analyse a chemical reaction model which considers RNA replication including correct copying and point mutations together with hydrolytic degradation and the dilution flux of a flow reactor. The corresponding stochastic reaction network is implemented on a computer in order to investigate some basic features of evolutionary optimization dynamics. Characteristic features of real molecular systems are mimicked by folding binary sequences into unknotted two-dimensional structures. Selective values are derived from these molecular 'phenotypes' by an evaluation procedure which assigns numerical values to different elements of the secondary structure. The fitness function obtained thereby contains nontrivial long-range interactions which are typical for real systems. The fitness landscape also reveals quite involved and bizarre local topologies which we consider also representative of polynucleotide replication in actually occurring systems. Optimization operates on an ensemble of sequences via mutation and natural selection. The strategy observed in the simulation experiments is fairly general and resembles closely a heuristic widely applied in operations research areas. Despite the relative smallness of the system--we study 2000 molecules of chain length v = 70 in a typical simulation experiment--features typical for the evolution of real populations are observed as there are error thresholds for replication, evolutionary steps and quasistationary sequence distributions. The relative importance of selectively neutral or almost neutral variants is discussed quantitatively. Four characteristic ensemble properties, entropy of the distribution, ensemble correlation, mean Hamming distance and diversity of the population, are computed and checked for their sensitivity in recording major optimization events during the simulation.

The secondary structure of human 28S rRNA: the structure and evolution of a mosaic rRNA gene

We have determined the secondary structure of the human 28S rRNA molecule based on comparative analysis of available eukaryotic cytoplasmic and prokaryotic large-rRNA gene sequences. Examination of large-rRNA sequences of both distantly and closely related species has enabled us to derive a structure that accounts both for highly conserved sequence tracts and for previously unanalyzed variable-sequence tracts that account for the evolutionary differences in size among the large rRNAs. Human 28S rRNA is composed of two different types of sequence tracts: conserved and variable. They differ in composition, degree of conservation, and evolution. The conserved regions demonstrate a striking constancy of size and sequence. We have confirmed that the conserved regions of large-rRNA molecules are capable of forming structures that are superimposable on one another. The variable regions contain the sequences responsible for the 83% increase in size of the human large-rRNA molecule over that of Escherichia coli. Their locations in the gene are maintained during evolution. They are G + C rich and largely nonhomologous, contain simple repetitive sequences, appear to evolve by frequent recombinational events, and are capable of forming large, stable hairpins. The secondary-structure model presented here is in close agreement with existing prokaryotic 23S rRNA secondary-structure models. The introduction of this model helps resolve differences between previously proposed prokaryotic and eukaryotic large-rRNA secondary-structure models.

Improved free-energy parameters for predictions of RNA duplex stability

Thermodynamic parameters for prediction of RNA duplex stability are reported. One parameter for duplex initiation and 10 parameters for helix propagation are derived from enthalpy and free-energy changes for helix formation by 45 RNA oligonucleotide duplexes. The oligomer sequences were chosen to maximize reliability of secondary structure predictions. Each of the 10 nearest-neighbor sequences is well-represented among the 45 oligonucleotides, and the sequences were chosen to minimize experimental errors in delta GO at 37 degrees C. These parameters predict melting temperatures of most oligonucleotide duplexes within 5 degrees C. This is about as good as can be expected from the nearest-neighbor model. Free-energy changes for helix propagation at dangling ends, terminal mismatches, and internal G X U mismatches, and free-energy changes for helix initiation at hairpin loops, internal loops, or internal bulges are also tabulated.

Secondary structure model for mouse beta Maj globin mRNA derived from enzymatic digestion data, comparative sequence and computer analysis

A model for the secondary structure of mouse beta Maj globin messenger RNA is presented based on enzymatic digestion data, comparative sequence and computer analysis. Using 5'-32P-end-labeled beta globin mRNA as a substrate, single-stranded regions were determined with S1 and T1 nucleases and double-stranded regions with V1 ribonuclease from cobra venom. The structure data obtained for ca. 75% of the molecule was introduced into a computer algorithm which predicts secondary structures of minimum free energy consistent with the enzymatic data. Two prominent base paired regions independently derived by phylogenetic analysis were also present in the computer generated structure lending support for the model. An interesting feature of the model is the presence of long-range base pairing interactions which permit the beta globin mRNA to fold back on itself, thereby bringing the 5'- and 3'-noncoding regions within close proximity. This feature is consistent with data from other laboratories suggesting an interaction of the 5'- and 3'-domains in the mammalian globin mRNAs.

Improved free energies for G.C base-pairs

Thermodynamic parameters of helix formation are reported for seven oligoribonucleotides containing only G.C pairs. These data are used with the nearest-neighbor model to calculate enthalpies and free energies of base-pair formation for G.C pairs. For helix initiation, the free energy change at 37 degrees C, delta G(0)37, is +3.9 kcal/mol; for helix propagation, the delta G(0)37 values are -2.3, -3.2 and -3.3 kcal/mol for C-G, G-G and G-C neighbors, respectively.

U1 snRNA: the evolution of its primary and secondary structure

In this paper we first show that the primary structure of U1 snRNA is homologous to that of tandem repeated pre-tRNA. Two sets of polymerase III promoter sites (the a and b boxes) are clearly recognisable at the appropriate positions in U1, although neither is functional; these sites occur in a degenerate form and their transcription is initiated by polymerase II. Moreover, several of the conserved subsequences of tRNAs that are not associated with transcription initiation (and supposedly are conserved because of their role in translation) are conserved in U1 as well, one of them being the pattern Py-Py-anticodon-Pu-Pu (for both anticodons of tandem tRNA). Second, we show that the secondary structure of U1 is apparently formed after fixation of the "B-hairpin loop' by one of the associated proteins. If and only if this hairpin loop is fixed, a consensus secondary structure is produced by the minimisation-of-free-energy technique. Moreover, we show that this B-hairpin loop has been destabilised relatively recently in evolutionary time by deletions (e.g., in the polymerase III box). If we reinsert the deleted bases, the so constructed hypothetical "ancestral" molecule folds into the consensus secondary structure by unconstrained energy minimisation (i.e., without fixation of the B-loop). Some features of the secondary structure of tandem repeated pre-tRNA are conserved in U1, but the overall structure has changed dramatically. Like tRNA, U1 has a cloverleaf-like structure, but its overall size has doubled. By comparing their secondary structures and by alignment of the sequences, we trace the local events associated with the global change in secondary structure (and apparently in the function of the molecule). Finally, we discuss our results from the perspective of informatic prerequisites for heterarchical multilevel evolution.

Solvent effects on the stability of A7U7p

The thermodynamics of double-helix formation were measured spectrophotometrically for A7U7 in water at 1 M NaCl and for A7U7p in a variety of solvent mixtures and salt. Comparison of the A7U7 results with calorimetric measurements indicates duplex formation involves intermediate states. For A7U7p between 0.06 and 0.55 M Na+, dTm/d(log [Na+]) = 17.4 degrees C, similar to the value of 19.6 degrees C for poly-(A).poly(U) [Krakauer, H., & Sturtevant, J. M. (1968) Biopolymers 6, 491-512]. At 1 M NaCl, the A7U7p duplex is most stable in 100% water. For 10 mol % solutions, the order for A7U7p duplex stability is ethylene glycol greater than glycerol greater than ethanol greater than 2-propanol greater than dimethyl sulfoxide greater than 1-propanol greater than formamide greater than N,N-dimethylformamide greater than urea greater than dioxane. Comparison of changes in stability and thermodynamic parameters with literature results for proteins suggests proteins and A7U7p interact differently with solvent. The results suggest hydrophobic bonding is not a major contributor to the stability of the A7U7p duplex. Comparisons with bulk solvent surface tension suggest the energy of cavity formation is also not a major contributor to duplex stability.

A new principle of RNA folding based on pseudoknotting

Tertiary interactions involving hairpin or interior loops of RNA can lead to extended quasi-continuous double helical stem regions, consisting of coaxially stacked segments of duplex RNA, bridged by single-stranded connections. This type of compact folding plays a role in various strategic regions of RNA molecules. Their role in ribosome functioning, RNA splicing and recognition of tRNA-like structures is discussed.

The secondary structure of oocyte and somatic 5S ribosomal RNAs of the fish Misgurnus fossilis L. from nuclease hydrolyses and chemical modification data

We have studied the accessibility of 5'- 32P labeled oocyte and somatic 5S rRNAs from the fish Misgurnus fossilis L. to S1, T1 and cobra venom nucleases and have found that the cleavage sites of 5S rRNAs closely related in primary structures differ in these molecules. The data of nuclease hydrolyses revealed the existence of two conformers corresponding to renatured and partially denatured somatic 5S rRNA and capable of mutual interconversions. The exposed cytosine residues were located in oocyte and somatic 5S rRNAs converted into uridine ones by sodium bisulfite treatment. The data have been used to construct the secondary structure models of somatic and oocyte 5S rRNAs by means of specially devised computer program. These models differ in their 5'-halves which contain all the nucleotide substitutions in the primary structure, all differences in location of the exposed cytosine residues, and finally, in the cleavage pattern by the nucleases used.

Thermodynamic studies of RNA stability

Enthalpies and entropies of helix stabilization due to addition of 3' terminal unpaired nucleotides to a CCGG or GGCC core double helix are derived from UV melting studies. The results suggest stacking provides a significant fraction of the free energy of a terminal base pair. The effects of temperature, aggregation, and ionic strength on the determination of thermodynamic parameters are considered. Helix propagation parameters are revised and extended based on recent additions to the data set.

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.

Generating non-overlapping displays of nucleic acid secondary structure

A new algorithm is presented which permits the display of nucleic acid secondary structure by computer. This algorithm circumvents the problem of overlapping portions of the molecule which is inherent in some other drawing programs. The results from this algorithm may also be used as input to the drawing algorithm previously reported in this journal [1] to untangle most of a drawing. The algorithm also represents the molecule in a form which makes visual comparisons for similarity quite easy since it guarantees that comparable features will reside in the same relative position in the drawings when the drawings are normalized.

Computer-aided nucleic acid secondary structure modeling incorporating enzymatic digestion data

We present a computer-aided method for determining nucleic acid secondary structure. The method utilizes a program which has the capability to filter matrix diagonal data on the basis of diagonal length, stabilization energy, and chemical and enzymatic data. The program also allows the user to assign selected regions of the structure as uniquely single-stranded or paired, and to filter out "trade-off" structures on the basis of such pairing. In order to demonstrate the utility of the program we present a preliminary secondary structure for the 3' end of alfalfa mosaic virus RNA 4 (AMV-4 RNA). This structure is based on an analysis which includes the use of in vitro partial enzymatic digestion of the RNA.

An accelerated algorithm for calculating the secondary structure of single stranded RNAs

We describe a code designed for secondary structure computation of single stranded RNA molecules. While it incorporates the same principles as the original algorithm of Nussinov et al (1978), its restructuring improves the logic and the approach of the codes based on it. For long sequences the code is at least an order of magnitude faster. For a chain n nucleotides long, references to computer disk memory are reduced from n3 to less than n2. For n much greater than 100, disk references behave like n3/6000.

The combinatorial distance geometry method for the calculation of molecular conformation. I. A new approach to an old problem

A new approach to the long-standing local minimum problem of molecular energy minimization is proposed. The approach relies upon a field of computer mathematics known as combinatorial optimization, together with methods of conformational analysis derived from distance geometry. The advantages over the usual numerical techniques of optimization are, first, that the algorithms derived are globally convergent, and second, that the mathematical problems involved are well-posed and suitable for study within the modern theory of computational complexity. In this paper we introduce the approach, and describe a computer program based on it.

Characterization of the yeast mitochondrial locus necessary for tRNA biosynthesis: DNA sequence analysis and identification of a new transcript

Most components necessary for the biosynthesis of mitochondrial tRNAs are coded by nuclear genes, but one mitochondrial locus other than the tRNA genes themselves is required to make functional tRNAs in the yeast Saccharomyces cerevisiae. DNA sequence analysis of this yeast mitochondrial tRNA synthesis locus is reported here. This region of mitochondrial DNA is almost exclusively A+T-rich DNA with one G+C-rich element. Despite the unusual structure of the DNA in this region, we have demonstrated that it codes for a heretofore unidentified mitochondrial transcript about 450 bases in length. Since this RNA is the only RNA encoded by the tRNA synthesis locus, it must be the active agent of the locus. This RNA could either act autonomously through RNA-RNA interactions or as part of an RNA-protein complex to effect tRNA biosynthesis.

Chemical modification of guanine residues of mouse 5 S ribosomal RNA with kethoxal. (Nucleosides and nucleotides 46)

Chemical modification of mouse 5 S rRNA with kethoxal was carried out to examine the secondary structure. The guanine residues located at positions 37, 41, 56, 66, 75 and 89 were modified. The relative rates of reaction are in the order G37, G56, G89, G66, G41, G75 at 28 degrees C and G37, G41, G56, G89, G75, G66 at 35 degrees C. These results support a secondary structure model containing 5 helices and 5 loops and indicate that the region around position 37 is the most exposed in higher-order structure.

Chemical modification of cytosine residues of mouse 5 S ribosomal RNA with hydrogen sulfide. (Nucleosides and nucleotides 43)

Cytosine residues of nucleic acids were converted to 4-thiouracil residues with hydrogen sulfide in pyridine and water to examine the secondary and tertiary structures of mouse 5 S rRNA. The cytosine residues at positions 10, 24, 34 (or 36), 39, 44 (or 46) and 63 were converted preferentially when the treatment was carried out at 28 degrees C. This result supports the model of the secondary structure of 5 S rRNA of Nishikawa, K. and Takemura, S. ((1974) FEBS Lett. 40, 106-109) consisting of five helices and five loops. As the temperature was increased to 35 degrees C, additional cytosine residues in positions 26, 52 and 78 were modified to moderate extents.

An algorithm for the display of nucleic acid secondary structure

A simple algorithm is presented for the graphic display of nucleic acid secondary structure. Examples of secondary structure displays are given for tRNA, 5S RNA and part of the 16S RNA. Due to its speed, this algorithm could easily be used in conjunction with secondary structure programs which calculate various alternate structures.

An interactive technique for the display of nucleic acid secondary structure

The ability to visualize nucleic acid secondary structure has become quite important since the advent of computer prediction and biochemical techniques that depict such structures. Manually drawing the conformations can be quite time consuming and tedious. Thus, the ability to draw with the aid of a computer the secondary structure of nucleic acid molecules is quite advantageous. This paper describes an interactive algorithm that permits one to generate such drawings which may then be used for further analysis and/or publications.