Comparative structure analysis of vertebrate ribonuclease P RNA

Ribonuclease P cleaves 5'-precursor sequences from pre-tRNAs. All cellular RNase P holoenzymes contain homologous RNA elements; the eucaryal RNase P RNA, in contrast to the bacterial RNA, is catalytically inactive in the absence of the protein component(s). To understand the function of eucaryal RNase P RNA, knowledge of its structure is needed. Considerable effort has been devoted to comparative studies of the structure of this RNA from diverse organisms, including eucaryotes, primarily fungi, but also a limited set of vertebrates. The substantial differences in the sequences and structures of the vertebrate RNAs from those of other organisms have made it difficult to align the vertebrate sequences, thus limiting comparative studies. To expand our understanding of the structure of diverse RNase P RNAs, we have isolated by PCR and sequenced 13 partial RNase P RNA genes from 11 additional vertebrate taxa representing most extant major vertebrate lineages. Based on a recently proposed structure of the core elements of RNase P RNA, we aligned the sequences and propose a minimum consensus secondary structure for the vertebrate RNase P RNA.     

Mitochondrial RNase P: the RNA family grows

Recent results of biochemical approaches and genome sequencing approaches has extended the members of the family of mitochondrial RNase P RNA genes. So far all of them are AU rich, and most of their secondary structures are easier to predict than was the structure of the first mitochondrial RNase P RNA from S. cerevisiae. The recently sequenced protozoan R. americana mitochondrial gene displays many of the evolutionarily conserved primary sequence and secondary structure attributes of the well characterized bacterial RNase P RNAs. Continued addition of RNAs to this data base should allow increasingly informative alignments and an understanding of what structural elements are dispensable in the smallest mitochondrial RNAs. The only protein subunits identified to date are Rpm2p from S.cerevisiae and the homologous protein from S. douglasii.     

Derivation of the three-dimensional architecture of bacterial ribonuclease P RNAs from comparative sequence analysis

The secondary structure of bacterial RNase P RNA, a ribozyme responsible for the maturation of the 5' end of tRNAs, is well established on the basis of sequence comparison analysis. RNase P RNA secondary structures fall into two types, A and B, which share a common core formed by the assembly of two main folding domains, but differ in their peripheral elements.A revised alignment of 137 available sequences reveals new covariations allowing for the refinement of both types of secondary structures. Phylogenetic evidence is thus provided for the extension of stems P11, P14, P19, P10.1 and P15.1 through further canonical base-pairs or GAellipsisGA mismatches. These refinements led in turn to a new organization of the catalytic core, with coaxial stackings of helices P2 and P19 as well as P1 and P4. New inter-domain tertiary interactions involve loop L9 and helix P1 and loop L8 with helix P4. These features were incorporated into atomic-scale 3D models of RNase P RNA for representatives of each structural type, namely Escherichia coli and Bacillus subtilis. In each model, the juxtaposition of the core helices creates a cradle onto which the pre-tRNA substrate binds with most evolutionarily conserved residues converging towards the cleavage site. The inner cores of both types are stabilized similarly, albeit by different peripheral elements, emphasizing the modular and hierarchical organisation of the architecture of RNase P RNAs. Similarities are thus apparent between the type A modules, P16/P17/P6 and P13/P14, and their type B analogs, P5.1/P15.1 and P10. 1/P10.1a, respectively. Other noteworthy features of these models include compactness and good agreement with published crosslinking data.     

Differentiation of Chlamydia spp. by sequence determination and restriction endonuclease cleavage of RNase P RNA genes

The amplification of DNA from Chlamydia trachomatis by PCR with degenerated primers yielded a 345-bp fragment of the putative RNase P RNA gene. From the deduced DNA sequence of this gene in C. trachomatis, a modified primer pair was designed. The primer pair was subsequently used to obtain the corresponding gene products from Chlamydia pneumoniae and Chlamydia psittaci. Sequence comparisons revealed similarities of 76.6% between C. trachomatis and C. pneumoniae, 79.5% between C. trachomatis and C. psittaci, and 84.7% between C. pneumoniae and C. psittaci. Furthermore, the three species were differentiated by fragment length polymorphism analysis after restriction enzyme cleavage of the PCR products. Sequence variations among 14 serotypes of C. trachomatis were confined to one purine base substitution in the putative RNase P RNA gene of lymphogranuloma venereum strains L1 to L3. Complete sequence similarity was found for nine strains of C. pneumoniae of different geographic origins. Taken together, our results indicate a possibility of the general application of this method in clinical bacteriology. Analysis of the secondary structures of the putative RNase P RNA genes from the different Chlamydia species suggested that a novel structural element in the domain of RNase P RNA is involved in base pairing with the 3'-terminal CCA motif of a tRNA precursor. This structure has not previously been found among RNase P RNAs of members of the division Bacteria.     

An active domain of the nuclear RNase P RNA

The P10/11-P12 RNA domain of yeast nuclear RNase P RNA has been characterized using genetic and biochemical analysis. This RNA domain contains some of the most conserved nucleotides throughout yeast species and shares considerable homology with the P10-P11-P12 bacterial RNase P RNA domain. Viable yeast variants generated by sequence randomization of the conserved internal loop nucleotides have demonstrated magnesium-sensitive growth defects. Partial purification and characterization of the RNase P holoenzyme from these variants reveals that the mutations affect the catalytic rate of the enzyme and increased magnesium concentrations are required to achieve maximal activity compared to wild type enzyme. Biochemical structure probing has been employed to address the interaction of the RNA domain with magnesium. Several nucleotides within the loop portion of the domain show magnesium-induced changes in reagent accessibility. These include the highly conserved nucleotides shared between yeast and bacteria, which become less accessible in the presence of magnesium. Conversely, accessibility of other regions of the RNA increases. The genetic and biochemical data suggest that the P10/11-P12 RNA domain, and the conserved nucleotides in particular, interacts with magnesium in a manner that affects catalysis by RNase P.     

Automatic detection of conserved RNA structure elements in complete RNA virus genomes

We propose a new method for detecting conserved RNA secondary structures in a family of related RNA sequences. Our method is based on a combination of thermodynamic structure prediction and phylogenetic comparison. In contrast to purely phylogenetic methods, our algorithm can be used for small data sets of approximately 10 sequences, efficiently exploiting the information contained in the sequence variability. The procedure constructs a prediction only for those parts of sequences that are consistent with a single conserved structure. Our implementation produces reasonable consensus structures without user interference. As an example we have analysed the complete HIV-1 and hepatitis C virus (HCV) genomes as well as the small segment of hantavirus. Our method confirms the known structures in HIV-1 and predicts previously unknown conserved RNA secondary structures in HCV.     

Prediction of RNA secondary structure based on helical regions distribution

MOTIVATION: RNAs play an important role in many biological processes and knowing their structure is important in understanding their function. Due to difficulties in the experimental determination of RNA secondary structure, the methods of theoretical prediction for known sequences are often used. Although many different algorithms for such predictions have been developed, this problem has not yet been solved. It is thus necessary to develop new methods for predicting RNA secondary structure. The most-used at present is Zuker's algorithm which can be used to determine the minimum free energy secondary structure. However many RNA secondary structures verified by experiments are not consistent with the minimum free energy secondary structures. In order to solve this problem, a method used to search a group of secondary structures whose free energy is close to the global minimum free energy was developed by Zuker in 1989. When considering a group of secondary structures, if there is no experimental data, we cannot tell which one is better than the others. This case also occurs in combinatorial and heuristic methods. These two kinds of methods have several weaknesses. Here we show how the central limit theorem can be used to solve these problems. RESULTS: An algorithm for predicting RNA secondary structure based on helical regions distribution is presented, which can be used to find the most probable secondary structure for a given RNA sequence. It consists of three steps. First, list all possible helical regions. Second, according to central limit theorem, estimate the occurrence probability of every helical region based on the Monte Carlo simulation. Third, add the helical region with the biggest probability to the current structure and eliminate the helical regions incompatible with the current structure. The above processes can be repeated until no more helical regions can be added. Take the current structure as the final RNA secondary structure. In order to demonstrate the confidence of the program, a test on three RNA sequences: tRNAPhe, Pre-tRNATyr, and Tetrahymena ribosomal RNA intervening sequence, is performed. AVAILABILITY: The program is written in Turbo Pascal 7.0. The source code is available upon request. CONTACT: Wujj@nic.bmi.ac.cn or Liwj@mail.bmi.ac.cn     

Statistics of RNA secondary structures

A statistical reference for RNA secondary structures with minimum free energies is computed by folding large ensembles of random RNA sequences. Four nucleotide alphabets are used: two binary alphabets, AU and GC, the biophysical AUGC and the synthetic GCXK alphabet. RNA secondary structures are made of structural elements, such as stacks, loops, joints, and free ends. Statistical properties of these elements are computed for small RNA molecules of chain lengths up to 100. The results of RNA structure statistics depend strongly on the particular alphabet chosen. The statistical reference is compared with the data derived from natural RNA molecules with similar base frequencies. Secondary structures are represented as trees. Tree editing provides a quantitative measure for the distance dt, between two structures. We compute a structure density surface as the conditional probability of two structures having distance t given that their sequences have distance h. This surface indicates that the vast majority of possible minimum free energy secondary structures occur within a fairly small neighborhood of any typical (random) sequence. Correlation lengths for secondary structures in their tree representations are computed from probability densities. They are appropriate measures for the complexity of the sequence-structure relation. The correlation length also provides a quantitative estimate for the mean sensitivity of structures to point mutations.     

Structure-function analysis in nuclear RNase P RNA

Eukaryotic ribonuclease P (RNase P) enzymes require both RNA and protein subunits for activity in vivo and in vitro. We have undertaken an analysis of the complex RNA subunit of the nuclear holoenzyme in an effort to understand its structure and its similarities to and differences from the bacterial ribozymes. Phylogenetic analysis, structure-sensitive RNA footprinting, and directed mutagenesis reveal conserved secondary and tertiary structures with both strong similarities to the bacterial consensus and distinctive features. The effects of mutations in the most highly conserved positions are being used to dissect the functions of individual subdomains.     

Multiple protein structure alignment

A method was developed to compare protein structures and to combine them into a multiple structure consensus. Previous methods of multiple structure comparison have only concatenated pairwise alignments or produced a consensus structure by averaging coordinate sets. The current method is a fusion of the fast structure comparison program SSAP and the multiple sequence alignment program MULTAL. As in MULTAL, structures are progressively combined, producing intermediate consensus structures that are compared directly to each other and all remaining single structures. This leads to a hierarchic "condensation," continually evaluated in the light of the emerging conserved core regions. Following the SSAP approach, all interatomic vectors were retained with well-conserved regions distinguished by coherent vector bundles (the structural equivalent of a conserved sequence position). Each bundle of vectors is summarized by a resultant, whereas vector coherence is captured in an error term, which is the only distinction between conserved and variable positions. Resultant vectors are used directly in the comparison, which is weighted by their error values, giving greater importance to the matching of conserved positions. The resultant vectors and their errors can also be used directly in molecular modeling. Applications of the method were assessed by the quality of the resulting sequence alignments, phylogenetic tree construction, and databank scanning with the consensus. Visual assessment of the structural superpositions and consensus structure for various well-characterized families confirmed that the consensus had identified a reasonable core.     

Changes in the rabbit temporomandibular joint associated with posterior displacement of the mandible

We report a detailed evolutionary study of the RNase P- and RNase MRP- associated RNAs. The analyses were performed on all the available complete sequences of RNase MRP (vertebrates, yeast, plant), nuclear RNase P (vertebrates, yeast), and mitochondrial RNase P (yeast) RNAs. For the first time the phylogenetic distance between these sequences and the nucleotide substitution rates have been quantitatively measured.The analyses were performed by considering the optimal multiple alignments obtained mostly by maximizing similarity between primary sequences. RNase P RNA and MRP RNA display evolutionary dynamics following the molecular clock. Both have similar rates and evolve about one order of magnitude faster than the corresponding small rRNA sequences which have been, so far, the most common gene markers used for phylogeny. However, small rRNAs evolve too slowly to solve close phylogenetic relationships such as those between mammals. The quicker rate of RNase P and MRP RNA allowed us to assess phylogenetic relationships between mammals and other vertebrate species and yeast strains. The phylogenetic data obtained with yeasts perfectly agree with those obtained by functional assays, thus demonstrating the potential offered by this approach for laboratory experiments.     

Assessing the reliability of RNA folding using statistical mechanics

We have analyzed the base-pairing probability distributions of 16 S and 16 S-like, and 23 S and 23 S-like ribosomal RNAs of Archaea, Bacteria, chloroplasts, mitochondria and Eukarya, as predicted by the partition function approach for RNA folding introduced by McCaskill. A quantitative analysis of the reliability of RNA folding is done by comparing the base-pairing probability distributions with the structures predicted by comparative sequence analysis (comparative structures). We distinguish two factors that show a relationship to the reliability of RNA minimum free energy structure. The first factor is the dominance of one particular base-pair or the absence of base-pairing for a given base within the base-pairing probability distribution (BPPD). We characterize the BPPD per base, including the probability of not base-pairing, by its Shannon entropy (S). The S value indicates the uncertainty about the base-pairing of a base: low S values result from BPPDs that are strongly dominated by a single base-pair or by the absence of base-pairing. We show that bases with low S values have a relatively high probability that their minimum free energy (MFE) structure corresponds to the comparative structure. The BPPDs of prokaryotes that live at high temperatures (thermophilic Archaea and Bacteria) have, calculated at 37 degrees C, lower S values than the BPPDs of prokaryotes that live at lower temperatures (mesophilic and psychrophilic Archaea and Bacteria). This reflects an adaptation of the ribosomal RNAs to the environmental temperature. A second factor that is important to consider with regard to the reliability of MFE structure folding is a variable degree of applicability of the thermodynamic model of RNA folding for different groups of RNAs. Here we show that among the bases that show low S values, the Archaea and Bacteria have similar, high probabilities (0.96 and 0.94 in 16 S and 0.93 and 0.91 in 23 S, respectively) that the MFE structure corresponds to the comparative structure. These probabilities are lower in the chloroplasts (16 S 0.91, 23 S 0.79), mitochondria (16 S-like 0.89, 23 S-like 0.69) and Eukarya (18 S 0.81, 28 S 0.86).     

Enhancement of antitumor immunity by expression of CD70 (CD27 ligand) or CD154 (CD40 ligand) costimulatory molecules in tumor cells

BACKGROUND: Most catalytic RNAs depend on divalent metal ions for folding and catalysis. A thorough structure-function analysis of catalytic RNA therefore requires the identification of the metal-ion-binding sites. Here, we probed the binding sites using Fenton chemistry, which makes use of the ability of Fe2+ to functionally or structurally replace Mg2+ at ion-binding sites and to generate short-lived and highly reactive hydroxyl radicals that can cleave nucleic acid and protein backbones in spatial proximity of these ion-binding sites. RESULTS: Incubation of group I intron RNA with Fe2+, sodium ascorbate and hydrogen peroxide yields distinctly cleaved regions that occur only in the correctly folded RNA in the presence of Mg2+ and can be competed by additional Mg2+, suggesting that Fe2+ and Mg2+ interact with the same sites. Cleaved regions in the catalytic core are conserved for three different group I introns, and there is good correlation between metal-ion-binding sites determined using our method and those determined using other techniques. In a model of the T4 phage-derived td intron, cleaved regions separated in the secondary structure come together in three-dimensional space to form several metal-ion-binding pockets. CONCLUSIONS: In contrast to structural probing with Fe2+/EDTA, cleavage with Fe2+ detects metal-ion-binding sites located primarily in the inside of the RNA. Essentially all metal-ion-binding pockets detected are formed by tertiary structure elements. Using this method, we confirmed proposed metal-ion-binding sites and identified new ones in group I intron RNAs. This approach should allow the localization of metal-ion-binding sites in RNAs of interest.     

A common 40 amino acid motif in eukaryotic RNases H1 and caulimovirus ORF VI proteins binds to duplex RNAs

Eukaryotic RNases H from Saccharomyces cerevisiae , Schizosaccharomyces pombe and Crithidia fasciculata , unlike the related Escherichia coli RNase HI, contain a non-RNase H domain with a common motif. Previously we showed that S.cerevisiae RNase H1 binds to duplex RNAs (either RNA-DNA hybrids or double-stranded RNA) through a region related to the double-stranded RNA binding motif. A very similar amino acid sequence is present in caulimovirus ORF VI proteins. The hallmark of the RNase H/caulimovirus nucleic acid binding motif is a stretch of 40 amino acids with 11 highly conserved residues, seven of which are aromatic. Point mutations, insertions and deletions indicated that integrity of the motif is important for binding. However, additional amino acids are required because a minimal peptide containing the motif was disordered in solution and failed to bind to duplex RNAs, whereas a longer protein bound well. Schizosaccharomyces pombe RNase H1 also bound to duplex RNAs, as did proteins in which the S.cerevisiae RNase H1 binding motif was replaced by either the C.fasciculata or by the cauliflower mosaic virus ORF VI sequence. The similarity between the RNase H and the caulimovirus domain suggest a common interaction with duplex RNAs of these two different groups of proteins.     

Implications of RNA structure on the annealing of a potent antisense RNA directed against the human immunodeficiency virus type 1

Antisense RNA-mediated regulation in bacterial systems is related to the kinetics of RNA-RNA annealing in vitro. Here, we investigated the secondary structure of alphaY69, an effective HIV-directed antisense RNA in human cells. Purified RNA preparations contain a single conformer. The global structure was identified by a cleavage experiment under native conditions using a short complementary oligonucleotide and RNase H. Structural analyses indicate a three-domain structure of alphaY69 consisting of two stem-loop elements connected by a seven-nucleotide single-stranded hinge region. Kinetic data suggest that the formation of base pairs between a CGC triplet of alphaY69 and its target RNA is essential for fast annealing. The complementary sequence stretch of the target folds into a high-energy secondary structure. The relationship between modifications in structural elements of alphaY69 and the annealing kinetics suggested that rate-limiting steps of the annealing involve a single site of alphaY69 and do not involve its 5' or 3'-end. Further, the data indicate that both initial base-specific interactions and duplex formation are dependent on the CGC triplet of the central region of alphaY69. This mechanism represents a specific and efficient way of RNA-RNA annealing that is initiated by the interaction of unstructured RNA regions.     

The catalytic core of RNase P

A deletion mutant of the catalytic RNA component of Escherichia coli RNase P missing residues 87-241 retains the ability to interact with the protein component to form a functional catalyst. The deletion of this phylogenetically conserved region significantly increases the Km, indicating that the deleted structures may be important for binding to the precursor tRNA substrate but not for the cleavage reaction. Under some reaction conditions, this RNase P deletion mutant can become a relatively non-specific nuclease, indicating that this RNA's catalytic center may be more exposed. The catalytic core of the RNase P is formed by less than one third of the 377 residues of the RNase P RNA.