Comparative structure analysis of vertebrate ribonuclease P RNA

Chemical reversible crosslinking enables measurement of RNA 3D distances and alternative conformations in cells

Three-dimensional (3D) structures dictate the functions of RNA molecules in a wide variety of biological processes. However, direct determination of RNA 3D structures in vivo is difficult due to their large sizes, conformational heterogeneity, and dynamics. Here we present a method, Spatial 2'-Hydroxyl Acylation Reversible Crosslinking (SHARC), which uses chemical crosslinkers of defined lengths to measure distances between nucleotides in cellular RNA. Integrating crosslinking, exonuclease (exo) trimming, proximity ligation, and high throughput sequencing, SHARC enables transcriptome-wide tertiary structure contact maps at high accuracy and precision, revealing heterogeneous RNA structures and interactions. SHARC data provide constraints that improves Rosetta-based RNA 3D structure modeling at near-nanometer resolution. Integrating SHARC-exo with other crosslinking-based methods, we discover compact folding of the 7SK RNA, a critical regulator of transcriptional elongation. These results establish a strategy for measuring RNA 3D distances and alternative conformations in their native cellular context.

Structural studies of RNase P

Ribonuclease P (RNase P) is one of the first ribozymes discovered and it is found in all phylogenetic groups. It is responsible for processing the 5' end of pre-tRNAs as well as other RNA molecules. RNase P is formed by an RNA molecule responsible for catalysis and one or more proteins. Structural studies of the proteins from different organisms, the bacterial RNA component, and a bacterial RNase P holoenzyme/tRNA complex provide insights into the mechanism of this universal ribozyme. Together with the existing wealth of biochemical information, these studies provide atomic-level information on the mechanism of RNase P and continue to expand our understanding of the structure and architecture of large RNA molecules and ribonucleoprotein complexes, the nature of catalysis by ribozymes, the structural basis of recognition of RNA by RNA molecules, and the evolution of enzymes from the prebiotic, RNA-based world to the modern world.

RNAiFOLD: a constraint programming algorithm for RNA inverse folding and molecular design

Synthetic biology is a rapidly emerging discipline with long-term ramifications that range from single-molecule detection within cells to the creation of synthetic genomes and novel life forms. Truly phenomenal results have been obtained by pioneering groups--for instance, the combinatorial synthesis of genetic networks, genome synthesis using BioBricks, and hybridization chain reaction (HCR), in which stable DNA monomers assemble only upon exposure to a target DNA fragment, biomolecular self-assembly pathways, etc. Such work strongly suggests that nanotechnology and synthetic biology together seem poised to constitute the most transformative development of the 21st century. In this paper, we present a Constraint Programming (CP) approach to solve the RNA inverse folding problem. Given a target RNA secondary structure, we determine an RNA sequence which folds into the target structure; i.e. whose minimum free energy structure is the target structure. Our approach represents a step forward in RNA design--we produce the first complete RNA inverse folding approach which allows for the specification of a wide range of design constraints. We also introduce a Large Neighborhood Search approach which allows us to tackle larger instances at the cost of losing completeness, while retaining the advantages of meeting design constraints (motif, GC-content, etc.). Results demonstrate that our software, RNAiFold, performs as well or better than all state-of-the-art approaches; nevertheless, our approach is unique in terms of completeness, flexibility, and the support of various design constraints. The algorithms presented in this paper are publicly available via the interactive webserver http://bioinformatics.bc.edu/clotelab/RNAiFold; additionally, the source code can be downloaded from that site.

Evolutionary patterns of non-coding RNAs

A plethora of new functions of non-coding RNAs (ncRNAs) have been discovered in past few years. In fact, RNA is emerging as the central player in cellular regulation, taking on active roles in multiple regulatory layers from transcription, RNA maturation, and RNA modification to translational regulation. Nevertheless, very little is known about the evolution of this "Modern RNA World" and its components. In this contribution, we attempt to provide at least a cursory overview of the diversity of ncRNAs and functional RNA motifs in non-translated regions of regular messenger RNAs (mRNAs) with an emphasis on evolutionary questions. This survey is complemented by an in-depth analysis of examples from different classes of RNAs focusing mostly on their evolution in the vertebrate lineage. We present a survey of Y RNA genes in vertebrates and study the molecular evolution of the U7 snRNA, the snoRNAs E1/U17, E2, and E3, the Y RNA family, the let-7 microRNA (miRNA) family, and the mRNA-like evf-1 gene. We furthermore discuss the statistical distribution of miRNAs in metazoans, which suggests an explosive increase in the miRNA repertoire in vertebrates. The analysis of the transcription of ncRNAs suggests that small RNAs in general are genetically mobile in the sense that their association with a hostgene (e.g. when transcribed from introns of a mRNA) can change on evolutionary time scales. The let-7 family demonstrates, that even the mode of transcription (as intron or as exon) can change among paralogous ncRNA.

PNPASE regulates RNA import into mitochondria

RNA import into mammalian mitochondria is considered essential for replication, transcription, and translation of the mitochondrial genome but the pathway(s) and factors that control this import are poorly understood. Previously, we localized polynucleotide phosphorylase (PNPASE), a 3' --> 5' exoribonuclease and poly-A polymerase, in the mitochondrial intermembrane space, a location lacking resident RNAs. Here, we show a new role for PNPASE in regulating the import of nuclear-encoded RNAs into the mitochondrial matrix. PNPASE reduction impaired mitochondrial RNA processing and polycistronic transcripts accumulated. Augmented import of RNase P, 5S rRNA, and MRP RNAs depended on PNPASE expression and PNPASE-imported RNA interactions were identified. PNPASE RNA processing and import activities were separable and a mitochondrial RNA targeting signal was isolated that enabled RNA import in a PNPASE-dependent manner. Combined, these data strongly support an unanticipated role for PNPASE in mediating the translocation of RNAs into mitochondria.

Structure and function of eukaryotic Ribonuclease P RNA

Ribonuclease P (RNase P) is the ribonucleoprotein endonuclease that processes the 5' ends of precursor tRNAs. Bacterial and eukaryal RNase P RNAs had the same primordial ancestor; however, they were molded differently by evolution. RNase P RNAs of eukaryotes, in contrast to bacterial RNAs, are not catalytically active in vitro without proteins. By comparing the bacterial and eukaryal RNAs, we can begin to understand the transitions made between the RNA and protein-dominated worlds. We report, based on crosslinking studies, that eukaryal RNAs, although catalytically inactive alone, fold into functional forms and specifically bind tRNA even in the absence of proteins. Based on the crosslinking results and crystal structures of bacterial RNAs, we develop a tertiary structure model of the eukaryal RNase P RNA. The eukaryal RNA contains a core structure similar to the bacterial RNA but lacks specific features that in bacterial RNAs contribute to catalysis and global stability of tertiary structure.

Ribonuclease P: the evolution of an ancient RNA enzyme

Ribonuclease P (RNase P) is an ancient and essential endonuclease that catalyses the cleavage of the 5' leader sequence from precursor tRNAs (pre-tRNAs). The enzyme is one of only two ribozymes which can be found in all kingdoms of life (Bacteria, Archaea, and Eukarya). Most forms of RNase P are ribonucleoproteins; the bacterial enzyme possesses a single catalytic RNA and one small protein. However, in archaea and eukarya the enzyme has evolved an increasingly more complex protein composition, whilst retaining a structurally related RNA subunit. The reasons for this additional complexity are not currently understood. Furthermore, the eukaryotic RNase P has evolved into several different enzymes including a nuclear activity, organellar activities, and the evolution of a distinct but closely related enzyme, RNase MRP, which has different substrate specificities, primarily involved in ribosomal RNA biogenesis. Here we examine the relationship between the bacterial and archaeal RNase P with the eukaryotic enzyme, and summarize recent progress in characterizing the archaeal enzyme. We review current information regarding the nuclear RNase P and RNase MRP enzymes in the eukaryotes, focusing on the relationship between these enzymes by examining their composition, structure and functions.

INFO-RNA--a fast approach to inverse RNA folding

MOTIVATION

The structure of RNA molecules is often crucial for their function. Therefore, secondary structure prediction has gained much interest. Here, we consider the inverse RNA folding problem, which means designing RNA sequences that fold into a given structure.

RESULTS

We introduce a new algorithm for the inverse folding problem (INFO-RNA) that consists of two parts; a dynamic programming method for good initial sequences and a following improved stochastic local search that uses an effective neighbor selection method. During the initialization, we design a sequence that among all sequences adopts the given structure with the lowest possible energy. For the selection of neighbors during the search, we use a kind of look-ahead of one selection step applying an additional energy-based criterion. Afterwards, the pre-ordered neighbors are tested using the actual optimization criterion of minimizing the structure distance between the target structure and the mfe structure of the considered neighbor. We compared our algorithm to RNAinverse and RNA-SSD for artificial and biological test sets. Using INFO-RNA, we performed better than RNAinverse and in most cases, we gained better results than RNA-SSD, the probably best inverse RNA folding tool on the market.

AVAILABILITY

www.bioinf.uni-freiburg.de?Subpages/software.html.

Interactions between RNase P protein subunits in archaea

A yeast two-hybrid system was used to identify protein-protein interactions between the ribonuclease P (RNase P) protein subunits Mth11p, Mth687p, Mth688p and Mth1618p from the archaeon Methanothermobacter thermoautotrophicus. Clear interactions between Mth688p and Mth687p, and between Mth1618p and Mth11p, were confirmed by HIS3 and LacZ reporter expression. Weaker interactions of Mth687p and Mth688p with Mth 11p, and Mth11p with itself, are also suggested. These interactions resemble, and confirm, those previously seen among the homologs of these proteins in the more complex yeast RNase P holoenzyme.

Structural implications of novel diversity in eucaryal RNase P RNA

Previous eucaryotic RNase P RNA secondary structural models have been based on limited diversity, representing only two of the approximately 30 phylogenetic kingdoms of the domain Eucarya. To elucidate a more generally applicable structure, we used biochemical, bioinformatic, and molecular approaches to obtain RNase P RNA sequences from diverse organisms including representatives of six additional kingdoms of eucaryotes. Novel sequences were from acanthamoeba (Acathamoeba castellanii, Balamuthia mandrillaris, Filamoeba nolandi), animals (Caenorhabditis elegans, Drosophila melanogaster), alveolates (Theileria annulata, Babesia bovis), conosids (Dictyostelium discoideum, Physarum polycephalum), trichomonads (Trichomonas vaginalis), microsporidia (Encephalitozoon cuniculi), and diplomonads (Giardia intestinalis). An improved alignment of eucaryal RNase P RNA sequences was assembled and used for statistical and comparative structural analysis. The analysis identifies a conserved core structure of eucaryal RNase P RNA that has been maintained throughout evolution and indicates that covariation in size occurs between some structural elements of the RNA. Eucaryal RNase P RNA contains regions of highly variable length and structure reminiscent of expansion segments found in rRNA. The eucaryal RNA has been remodeled through evolution as a simplified version of the structure found in bacterial and archaeal RNase P RNAs.

Antisense inhibition of Escherichia coli RNase P RNA: mechanistic aspects

The ribonucleoprotein enzyme RNase P catalyzes endonucleolytic 5'-maturation of tRNA primary transcripts in all domains of life. The indispensability of RNase P for bacterial cell growth and the large differences in structure and function between bacterial and eukaryotic RNase P enzymes comply with the basic requirements for a bacterial enzyme to be suitable as a potential novel drug target. We have identified RNA oligonucleotides that start to show an inhibitory effect on bacterial RNase P RNAs of the structural type A (for example, the Escherichia coli or Klebsiella pneumoniae enzymes) at subnanomolar concentrations in our in vitro precursor tRNA (ptRNA) processing assay. These oligonucleotides are directed against the so-called P15 loop region of RNase P RNA known to interact with the 3'-CCA portion of ptRNA substrates. Lead probing experiments demonstrate that a complementary RNA or DNA 14-mer fully invades the P15 loop region and thereby disrupts local structure in the catalytic core of RNase P RNA. Binding of the RNA 14-mer is essentially irreversible because of a very low dissociation rate. The association rate of this oligonucleotide is on the order of 10(4) M(-1) s(-1) and is thus comparable to those of many other artificial antisense oligonucleotides. The remarkable inhibition efficacy is attributable to the dual effect of direct interference with substrate binding to the RNase P RNA active site and induction of misfolding of the catalytic core of RNase P RNA. Based on our findings, the P15 loop region of bacterial RNase P RNAs of the structural type A can be considered the "Achilles' heel" of the ribozyme and therefore represents a promising target for combatting multiresistant bacterial pathogens.

Isolation and Molecular Evolution of the Selenocysteine tRNA (Cf TRSP) and RNase P RNA (Cf RPPH1) Genes in the Dog Family, Canidae

In an effort to identify rapidly evolving nuclear sequences useful for phylogenetic analyses of closely related species, we isolated two genes transcribed by RNA polymerase III (pol III), the selenocysteine tRNA gene (TRSP) and an RNase P RNA (RPPH1) gene from the domestic dog (Canis familiaris). We focus on genes transcribed by pol III because their coding regions are small (generally 100-300 base pairs [bp]) and their essential promoter elements are located within a couple of hundred bps upstream of the coding region. Therefore, we predicted that regions flanking the coding region and outside of the promoter elements would be free of constraint and would evolve rapidly. We amplified TRSP from 23 canids and RPPH1 from 12 canids and analyzed the molecular evolution of these genes and their utility as phylogenetic markers for resolving relationships among species in Canidae. We compared the rate of evolution of the gene-flanking regions to other noncoding regions of nuclear DNA (introns) and to the mitochondrial encoded COII gene. Alignment of TRSP from 23 canids revealed that regions directly adjacent to the coding region display high sequence variability. We discuss this pattern in terms of functional mechanisms of transcription. Although the flanking regions evolve no faster than introns, both genes were found to be useful phylogenetic markers, in part, because of the synapomorphic indels found in the flanking regions. Gene trees generated from the TRSP and RPPH1 loci were generally in agreement with the published mtDNA phylogeny and are the first phylogeny of Canidae based on nuclear sequences.

A conserved element in the yeast RNase MRP RNA subunit can participate in a long-range base-pairing interaction

RNase MRP is a ribonucleoprotein endoribonuclease involved in eukaryotic pre-rRNA processing. The enzyme possesses a putatively catalytic RNA subunit, structurally related to that of RNase P. A thorough structure analysis of Saccharomyces cerevisiae MRP RNA, entailing enzymatic and chemical probing, mutagenesis and thermal melting, identifies a previously unrecognised stem that occupies a position equivalent to the P7 stem of RNase P. Inclusion of this P7-like stem confers on yeast MRP RNA a greater degree of similarity to the core RNase P RNA structure than that described previously and better delimits domain 2, the proposed specificity domain. The additional stem is created by participation of a conserved sequence element (ymCR-II) in a long-range base-pairing interaction. There is potential for this base-pairing throughout the known yeast MRP RNA sequences. Formation of a P7-like stem is not required, however, for the pre-rRNA processing or essential function of RNase MRP. Mutants that can base-pair are nonetheless detrimental to RNase MRP function, indicating that the stem will form in vivo but that only the wild-type pairing is accommodated. Although the alternative MRP RNA structure described is clearly not part of the active RNase MRP enzyme, it would be the more stable structure in the absence of protein subunits and the probability that it represents a valid intermediate species in the process of yeast RNase MRP assembly is discussed.

A new algorithm for RNA secondary structure design

The function of many RNAs depends crucially on their structure. Therefore, the design of RNA molecules with specific structural properties has many potential applications, e.g. in the context of investigating the function of biological RNAs, of creating new ribozymes, or of designing artificial RNA nanostructures. Here, we present a new algorithm for solving the following RNA secondary structure design problem: given a secondary structure, find an RNA sequence (if any) that is predicted to fold to that structure. Unlike the (pseudoknot-free) secondary structure prediction problem, this problem appears to be hard computationally. Our new algorithm, "RNA Secondary Structure Designer (RNA-SSD)", is based on stochastic local search, a prominent general approach for solving hard combinatorial problems. A thorough empirical evaluation on computationally predicted structures of biological sequences and artificially generated RNA structures as well as on empirically modelled structures from the biological literature shows that RNA-SSD substantially out-performs the best known algorithm for this problem, RNAinverse from the Vienna RNA Package. In particular, the new algorithm is able to solve structures, consistently, for which RNAinverse is unable to find solutions. The RNA-SSD software is publically available under the name of RNA Designer at the RNASoft website (www.rnasoft.ca).

Phylogenetic analysis of the structure of RNase MRP RNA in yeasts

RNase MRP is a ribonucleoprotein enzyme involved in processing precursor rRNA in eukaryotes. To facilitate our structure-function analysis of RNase MRP from Saccharomyces cerevisiae, we have determined the likely secondary structure of the RNA component by a phylogenetic approach in which we sequenced all or part of the RNase MRP RNAs from 17 additional species of the Saccharomycetaceae family. The structure deduced from these sequences contains the helices previously suggested to be common to the RNA subunit of RNase MRP and the related RNA subunit of RNase P, an enzyme cleaving tRNA precursors. However, outside this common region, the structure of RNase MRP RNA determined here differs from a previously proposed universal structure for RNase MRPs. Chemical and enzymatic structure probing analyses were consistent with our revised secondary structure. Comparison of all known RNase MRP RNA sequences revealed three regions with highly conserved nucleotides. Two of these regions are part of a helix implicated in RNA catalysis in RNase P, suggesting that RNase MRP may cleave rRNA using a similar catalytic mechanism.

Eukaryotic ribonuclease P: a plurality of ribonucleoprotein enzymes

Ribonuclease P (RNase P) is an essential endonuclease that acts early in the tRNA biogenesis pathway. This enzyme catalyzes cleavage of the leader sequence of precursor tRNAs (pre-tRNAs), generating the mature 5' end of tRNAs. RNase P activities have been identified in Bacteria, Archaea, and Eucarya, as well as organelles. Most forms of RNase P are ribonucleoproteins, i.e., they consist of an essential RNA subunit and protein subunits, although the composition of the enzyme in mitochondria and chloroplasts is still under debate. The recent purification of the eukaryotic nuclear RNase P has demonstrated a significantly larger protein content compared to the bacterial enzyme. Moreover, emerging evidence suggests that the eukaryotic RNase P has evolved into at least two related nuclear enzymes with distinct functions, RNase P and RNase MRP. Here we review current information on RNase P, with emphasis on the composition, structure, and functions of the eukaryotic nuclear holoenzyme, and its relationship with RNase MRP.

An essential protein-binding domain of nuclear RNase P RNA

Eukaryotic RNase P and RNase MRP are endoribonucleases composed of RNA and protein subunits. The RNA subunits of each enzyme share substantial secondary structural features, and most of the protein subunits are shared between the two. One of the conserved RNA subdomains, designated P3, has previously been shown to be required for nucleolar localization. Phylogenetic sequence analysis suggests that the P3 domain interacts with one of the proteins common to RNase P and RNase MRP, a conclusion strengthened by an earlier observation that the essential domain can be interchanged between the two enzymes. To examine possible functions of the P3 domain, four conserved nucleotides in the P3 domain of Saccharomyces cerevisiae RNase P RNA (RPR1) were randomized to create a library of all possible sequence combinations at those positions. Selection of functional genes in vivo identified permissible variations, and viable clones that caused yeast to exhibit conditional growth phenotypes were tested for defects in RNase P RNA and tRNA biosynthesis. Under nonpermissive conditions, the mutants had reduced maturation of the RPR1 RNA precursor, an expected phenotype in cases where RNase P holoenzyme assembly is defective. This loss of RPR1 RNA maturation coincided, as expected, with a loss of pre-tRNA maturation characteristic of RNase P defects. To test whether mutations at the conserved positions inhibited interactions with a particular protein, specific binding of the individual protein subunits to the RNA subunit was tested in yeast using the three-hybrid system. Pop1p, the largest subunit shared by RNases P and MRP, bound specifically to RPR1 RNA and the isolated P3 domain, and this binding was eliminated by mutations at the conserved P3 residues. These results indicate that Pop1p interacts with the P3 domain common to RNases P and MRP, and that this interaction is critical in the maturation of RNase P holoenzyme.

Eukaryotic ribonuclease P: increased complexity to cope with the nuclear pre-tRNA pathway

Ribonuclease P is an ancient enzyme that cleaves pre-tRNAs to generate mature 5' ends. It contains an essential RNA subunit in Bacteria, Archaea, and Eukarya, but the degree to which the RNA subunit relies on proteins to supplement catalysis is highly variable. The eukaryotic nuclear holoenzyme has recently been found to contain almost twenty times the protein content of the bacterial enzymes, in addition to having split into at least two related enzymes with distinct substrate specificity. In this review, recent progress in understanding the molecular architecture and functions of nuclear forms of RNase P will be considered.

Phylogenetic-comparative analysis of the eukaryal ribonuclease P RNA

Ribonuclease P (RNase P) is the ribonucleoprotein enzyme that cleaves 5'-leader sequences from precursor-tRNAs. Bacterial and eukaryal RNase P RNAs differ fundamentally in that the former, but not the latter, are capable of catalyzing pre-tRNA maturation in vitro in the absence of proteins. An explanation of these functional differences will be assisted by a detailed comparison of bacterial and eukaryal RNase P RNA structures. However, the structures of eukaryal RNase P RNAs remain poorly characterized, compared to their bacterial and archaeal homologs. Hence, we have taken a phylogenetic-comparative approach to refine the secondary structures of eukaryal RNase P RNAs. To this end, 20 new RNase P RNA sequences have been determined from species of ascomycetous fungi representative of the genera Arxiozyma, Clavispora, Kluyveromyces, Pichia, Saccharomyces, Saccharomycopsis, Torulaspora, Wickerhamia, and Zygosaccharomyces. Phylogenetic-comparative analysis of these and other sequences refines previous eukaryal RNase P RNA secondary structure models. Patterns of sequence conservation and length variation refine the minimum-consensus model of the core eukaryal RNA structure. In comparison to bacterial RNase P RNAs, the eukaryal homologs lack RNA structural elements thought to be critical for both substrate binding and catalysis. Nonetheless, the eukaryal RNA retains the main features of the catalytic core of the bacterial RNase P. This indicates that the eukaryal RNA remains intrinsically a ribozyme.

Evidence for an RNA-based catalytic mechanism in eukaryotic nuclear ribonuclease P

Ribonuclease P is the enzyme responsible for removing the 5'-leader segment of precursor transfer RNAs in all organisms. All eukaryotic nuclear RNase Ps are ribonucleoproteins in which multiple protein components and a single RNA species are required for activity in vitro as well as in vivo. It is not known, however, which subunits participate directly in phosphodiester-bond hydrolysis. The RNA subunit of nuclear RNase P is evolutionarily related to its catalytically active bacterial counterpart, prompting speculation that in eukaryotes the RNA may be the catalytic component. In the bacterial RNase P reaction, Mg(II) is required to coordinate the nonbridging phosphodiester oxygen(s) of the scissile bond. As a consequence, bacterial RNase P cannot cleave pre-tRNA in which the pro-Rp nonbridging oxygen of the scissile bond is replaced by sulfur. In contrast, the RNase P reaction in plant chloroplasts is catalyzed by a protein enzyme whose mechanism does not involve Mg(II) coordinated by the pro-Rp oxygen. To determine whether the mechanism of nuclear RNase P resembles more closely an RNA- or a protein-catalyzed reaction, we analyzed the ability of Saccharomyces cerevisiae nuclear RNase P to cleave pre-tRNA containing a sulfur substitution of the pro-Rp oxygen at the cleavage site. Sulfur substitution at this position prohibits correct cleavage of pre-tRNA. Cleavage by eukaryotic RNase P thus depends on the presence of a thio-sensitive ligand to the pro-Rp oxygen of the scissile bond, and is consistent with a common, RNA-based mechanism for the bacterial and eukaryal enzymes.

PCR primers and amplification methods for 12S ribosomal DNA, the control region, cytochrome oxidase I, and cytochrome b in bufonids and other frogs, and an overview of PCR primers which have amplified DNA in amphibians successfully

New primers (N = 24) for the amplification and sequencing of the complete or near complete 12S ribosomal DNA, about 1000 bp of the control region, 390 bp of cytochrome oxidase I, and the near complete cytochrome b are described. The 12S ribosomal DNA primers successfully amplify DNA in tetrapods; other primers successfully amplify DNA in bufonoids and other anurans. An overview of published literature and sequence data banks identified 170 mitochondrial and 96 nuclear DNA primers that have been used or are highly likely to be useful in amphibians. Primer sequences, their locations within genes, and sequence location and identity in Xenopus and human and/or mouse are presented for each primer. The utility of each primer was estimated by identifying the smallest, yet most inclusive, taxonomic category within which each primer has been successful. Primers from all published sources are mapped together. We hope that these new primers, as well as the list of primers that have been useful in amphibians, will encourage further systematic and population genetic studies of amphibians.