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.

Fluorescent RNA cytochemistry: tracking gene transcripts in living cells

The advent of jellyfish green fluorescent protein and its spectral variants, together with promising new fluorescent proteins from other classes of the Cnidarian phylum (coral and anemones), has greatly enhanced and promises to further boost the detection and localization of proteins in cell biology. It has been less widely appreciated that highly sensitive methods have also recently been developed for detecting the movement and localization in living cells of the very molecules that precede proteins in the gene expression pathway, i.e. RNAs. These approaches include the microinjection of fluorescent RNAs into living cells, the in vivo hybridization of fluorescent oligonucleotides to endogenous RNAs and the expression in cells of fluorescent RNA-binding proteins. This new field of ‘fluorescent RNA cytochemistry’ is summarized in this article, with emphasis on the biological insights it has already provided. These new techniques are likely to soon collaborate with other emerging approaches to advance the investigation of RNA birth, RNA–protein assembly and ribonucleoprotein particle transport in systems such as oocytes, embryos, neurons and other somatic cells, and may even permit the observation of viral replication and transcription pathways as they proceed in living cells, ushering in a new era of nucleic acids research in vivo.

RegA proteins from phage T4 and RB69 have conserved helix–loop groove RNA binding motifs but different RNA binding specificities

The RegA proteins from the bacteriophage T4 and RB69 are translational repressors that control the expression of multiple phage mRNAs. RegA proteins from the two phages share 78% sequence identity; however, in vivo expression studies have suggested that the RB69 RegA protein binds target RNAs with a higher affinity than T4 RegA protein. To study the RNA binding properties of T4 and RB69 RegA proteins more directly, the binding sites of RB69 RegA protein on synthetic RNAs corresponding to the translation initiation region of two RB69 target genes were mapped by RNase protection assays. These assays revealed that RB69 RegA protein protects nucleotides –9 to –3 (relative to the start codon) on RB69 gene 44, which contains the sequence GAAAAUU. On RB69 gene 45, the protected site (nucleotides –8 to –3) contains a similar purine-rich sequence: GAAAUA. Interestingly, T4 RegA protein protected the same nucleotides on these RNAs. To examine the specificity of RNA binding, quantitative RNA gel shift assays were performed with synthetic RNAs corresponding to recognition elements (REs) in three T4 and three RB69 mRNAs. Comparative gel shift assays demonstrated that RB69 RegA protein has an ∼7-fold higher affinity for T4 gene 44 RE RNA than T4 RegA protein. RB69 RegA protein also binds RB69 gene 44 RE RNA with a 4-fold higher affinity than T4 RegA protein. On the other hand, T4 RegA exhibited a higher affinity than RB69 RegA protein for RB69 gene 45 RE RNA. With respect to their affinities for cognate RNAs, both RegA proteins exhibited the following hierarchy of affinities: gene 44 > gene 45 > regA. Interestingly, T4 RegA exhibited the highest affinity towards RB69 gene 45 RE RNA, whereas RB69 RegA protein had the highest affinity for T4 gene 44 RE RNA. The helix–loop groove RNA binding motif of T4 RegA protein is fully conserved in RB69 RegA protein. However, homology modeling of the structure of RB69 RegA protein reveals that the divergent residues are clustered in two areas of the surface, and that there are two large areas of high conservation near the helix–loop groove, which may also play a role in RNA binding.

Kinetics of formation of hypoxanthine containing base pairs by HIV-RT: RNA template effects on the base substitution frequencies

Hypoxanthine (H), the deamination product of adenine, has been implicated in the high frequency of A to G transitions observed in retroviral and other RNA genomes. Although H·C base pairs are thermodynamically more stable than other H·N pairs, polymerase selection may be determined in part by kinetic factors. Therefore, the hypoxanthine induced substitution pattern resulting from replication by viral polymerases may be more complex than that predicted from thermodynamics. We have examined the steady-state kinetics of formation of base pairs opposite template H in RNA by HIV-RT, and for the incorporation of dITP during first- and second-strand synthesis. Hypoxanthine in an RNA template enhances the k2app for pairing with standard dNTPs by factors of 10–1000 relative to adenine at the same sequence position. The order of base pairing preferences for H in RNA was observed to be H·C >> H·T > H·A > H·G. Steady-state kinetics of insertion for all possible mispairs formed with dITP were examined on RNA and DNA templates of identical sequence. Insertion of dITP opposite all bases occurs 2–20 times more frequently on RNA templates. This bias for higher insertion frequencies on RNA relative to DNA templates is also observed for formation of mispairs at template A. This kinetic advantage afforded by RNA templates for mismatches and pairing involving H suggests a higher induction of mutations at adenines during first-strand synthesis by HIV-RT.

RNA: a method to specifically inhibit PCR amplification of known members of a multigene family by degenerate primers

The polymerase chain reaction (PCR) is a versatile method to amplify specific DNA with oligonucleotide primers. By designing degenerate PCR primers based on amino acid sequences that are highly conserved among all known gene family members, new members of a multigene family can be identified. The inherent weakness of this approach is that the degenerate primers will amplify previously identified, in addition to new, family members. To specifically address this problem, we synthesized a specific RNA for each known family member so that it hybridized to one strand of the template, adjacent to the 3′-end of the primer, allowing the degenerate primer to bind yet preventing extension by DNA polymerase. To test our strategy, we used known members of the soluble, nitric oxide-sensitive guanylyl cyclase family as our templates and degenerate primers that discriminate this family from other guanylyl cyclases. We demonstrate that amplification of known members of this family is effectively and specifically inhibited by the corresponding RNAs, alone or in combination. This robust method can be adapted to any application where multiple PCR products are amplified, as long as the sequence of the desired and the undesired PCR product(s) is sufficiently distinct between the primers.

Fluorescent quenching-based quantitative detection of specific DNA/RNA using a BODIPY® FL-labeled probe or primer

We have developed a simple method for the quantitative detection of specific DNA or RNA molecules based on the finding that BODIPY® FL fluorescence was quenched by its interaction with a uniquely positioned guanine. This approach makes use of an oligonucleotide probe or primer containing a BODIPY® FL-modified cytosine at its 5′-end. When such a probe was hybridized with a target DNA, its fluorescence was quenched by the guanine in the target, complementary to the modified cytosine, and the quench rate was proportional to the amount of target DNA. This widely applicable technique will be used directly with larger samples or in conjunction with the polymerase chain reaction to quantify small DNA samples.

RISSC: a novel database for ribosomal 16S–23S RNA genes spacer regions

A novel database, under the acronym RISSC (Ribosomal Intergenic Spacer Sequence Collection), has been created. It compiles more than 1600 entries of edited DNA sequence data from the 16S–23S ribosomal spacers present in most prokaryotes and organelles (e.g. mitochondria and chloroplasts) and is accessible through the Internet (http://ulises.umh.es/RISSC), where systematic searches for specific words can be conducted, as well as BLAST-type sequence searches. Additionally, a characteristic feature of this region, the presence/absence and nature of tRNA genes within the spacer, is included in all the entries, even when not previously indicated in the original database. All these combined features could provide a useful documen­tation tool for studies on evolution, identification, typing and strain characterization, among others.

PseudoBase: structural information on RNA pseudoknots

PseudoBase is a database containing structural, functional and sequence data related to RNA pseudo­knots. It can be reached at http://wwwbio.LeidenUniv.nl/∼Batenburg/PKB.html. For each pseudoknot, thirteen items are stored, for example the relevant sequence, the stem positions of the pseudoknot, the EMBL accession number of the sequence and the support that can be given regarding the reliability of the pseudo­knot. Since the last publication, information on sizes of the stems and the loops in the pseudoknots has been added. Also added are alternative entries that produce surveys of where the pseudoknots are, sorted according to stem size or loop size.

The European Large Subunit Ribosomal RNA Database

The European Large Subunit Ribosomal RNA Database compiles all complete or nearly complete large subunit ribosomal RNA sequences available from public sequence databases. These are provided in aligned format and the secondary structure, as derived by comparative sequence analysis, is included. Additional information about the sequences such as literature references and taxonomic information is also included. The database is available from our WWW server at http://rrna.uia.ac.be/lsu/.

rrndb: the Ribosomal RNA Operon Copy Number Database

The Ribosomal RNA Operon Copy Number Database (rrndb) is an Internet-accessible database containing annotated information on rRNA operon copy number among prokaryotes. Gene redundancy is uncommon in prokaryotic genomes, yet the rRNA genes can vary from one to as many as 15 copies. Despite the widespread use of 16S rRNA gene sequences for identification of prokaryotes, information on the number and sequence of individual rRNA genes in a genome is not readily accessible. In an attempt to understand the evolutionary implications of rRNA operon redundancy, we have created a phylogenetically arranged report on rRNA gene copy number for a diverse collection of prokaryotic microorganisms. Each entry (organism) in the rrndb contains detailed information linked directly to external websites including the Ribosomal Database Project, GenBank, PubMed and several culture collections. Data contained in the rrndb will be valuable to researchers investigating microbial ecology and evolution using 16S rRNA gene sequences. The rrndb web site is directly accessible on the WWW at http://rrndb.cme.msu.edu.

ACTIVITY: a database on DNA/RNA sites activity adapted to apply sequence-activity relationships from one system to another

ACTIVITY is a database on DNA/RNA site sequences with known activity magnitudes, measurement systems, sequence-activity relationships under fixed experimental conditions and procedures to adapt these relationships from one measurement system to another. This database deposits information on DNA/RNA affinities to proteins and cell nuclear extracts, cutting efficiencies, gene transcription activity, mRNA translation efficiencies, mutability and other biological activities of natural sites occurring within promoters, mRNA leaders, and other regulatory regions in pro- and eukaryotic genomes, their mutant forms and synthetic analogues. Since activity magnitudes are heavily system-dependent, the current version of ACTIVITY is supplemented by three novel sub-databases: (i) SYSTEM, measurement systems; (ii) KNOWLEDGE, sequence-activity relationships under fixed experimental conditions; and (iii) CROSS_TEST, procedures adapting a relationship from one measurement system to another. These databases are useful in molecular biology, pharmacogenetics, metabolic engineering, drug design and biotechnology. The databases can be queried using SRS and are available through the Web, http://wwwmgs.bionet.nsc.ru/systems/Activity/.

HyPaLib: a database of RNAs and RNA structural elements defined by hybrid patterns

The database, called HyPaLib (for Hybrid Pattern Library), contains annotated structural elements characteristic for certain classes of structural and/or functional RNAs. These elements are described in a language specifically designed for this purpose. The language allows convenient specification of hybrid patterns, i.e. motifs consisting of sequence features and structural elements together with sequence similarity and thermodynamic constraints. We are currently developing software tools that allow a user to search sequence databases for any pattern in HyPaLib, thus providing functionality which is similar to PROSITE, but dedicated to the more complex patterns in RNA sequences. HyPaLib is available at http://bibiserv.techfak.uni-bielefeld.de/HyPa/.

Interaction of a transcriptional repressor with the RNA polymerase II holoenzyme plays a crucial role in repression

The yeast transcriptional repressor Tup1, tethered to DNA, represses to strikingly different degrees transcription elicited by members of two classes of activators. Repression in both cases is virtually eliminated by mutation of either member of the cyclin-kinase pair Srb10/11. In contrast, telomeric chromatin affects both classes of activators equally, and in neither case is that repression affected by mutation of Srb10/11. In vitro, Tup1 interacts with RNA polymerase II holoenzyme bearing Srb10 as well as with the separated Srb10. These and other findings indicate that at least one aspect of Tup1's action involves interaction with the RNA polymerase II holoenzyme.

RNA editing in Trypanosoma brucei: characterization of gRNA U-tail interactions with partially edited mRNA substrates

Guide RNAs (gRNAs), key components of the RNA editing reaction in Trypanosoma brucei, direct the insertion and deletion of uridylate (U) residues. Analyses of gRNAs reveal three functional elements. The 5′-end of the gRNA contains the anchor, which is responsible for selection and binding to the pre-edited mRNA. The second element (the guiding region) provides the information required for editing. At the 3′-end of the gRNA is a non-encoded U-tail, whose function remains unclear. However, the cleavage–ligation model for editing proposes that the U-tail binds to purine-rich regions upstream of editing sites, thereby strengthening the interaction and holding onto the 5′ cleavage product. Our previous studies demonstrated that the U-tail interacts with upstream sequences and may play roles in both stabilization and tethering. These studies also indicated that the U-tail interactions involved mRNA regions that were to be subsequently edited. This raised the question of what happens to the mRNA–U-tail interaction as editing proceeds in the 3′→5′ direction. We examined gCYb-558 and its U-tail interaction with 5′CYbUT and two partially edited 5′CYb substrates. Our results indicate that the 3′-end of the U-tail interacts with the same sequence in all three mRNAs. Predicted secondary structures using crosslinking data suggest that a similar structure is maintained as editing proceeds. These results indicate that the role of the U-tail may also involve maintenance of important secondary structure motifs.

Ribosomal protein S7 from Escherichia coli uses the same determinants to bind 16S ribosomal RNA and its messenger RNA

Ribosomal protein S7 from Escherichia coli binds to the lower half of the 3′ major domain of 16S rRNA and initiates its folding. It also binds to its own mRNA, the str mRNA, and represses its translation. Using filter binding assays, we show in this study that the same mutations that interfere with S7 binding to 16S rRNA also weaken its affinity for its mRNA. This suggests that the same protein regions are responsible for mRNA and rRNA binding affinities, and that S7 recognizes identical sequence elements within the two RNA targets, although they have dissimilar secondary structures. Overexpression of S7 is known to inhibit bacterial growth. This phenotypic growth defect was relieved in cells overexpressing S7 mutants that bind poorly the str mRNA, confirming that growth impairment is controlled by the binding of S7 to its mRNA. Interestingly, a mutant with a short deletion at the C-terminus of S7 was more detrimental to cell growth than wild-type S7. This suggests that the C-terminal portion of S7 plays an important role in ribosome function, which is perturbed by the deletion.