MODOMICS: a database of RNA modification pathways. 2008 update

The architecture of eukaryotic translation

Translation in baker’s yeast involves the coordinated interaction of 200 000 ribosomes, 3 000 000 tRNAs and between 15 000 and 60 000 mRNAs. It is currently unknown whether this specific constellation of components has particular relevance for the requirements of the yeast proteome, or whether this is simply a frozen accident. Our study uses a computational simulation model of the genome-wide translational apparatus of yeast to explore quantitatively which combinations of mRNAs, ribosomes and tRNAs can produce viable proteomes. Surprisingly, we find that if we only consider total translational activity over time without regard to composition of the proteome, then there are many and widely differing combinations that can generate equivalent synthesis yields. In contrast, translational activity required for generating specific proteomes can only be achieved within a much more constrained parameter space. Furthermore, we find that strongly ribosome limited regimes are optimal for cells in that they are resource efficient and simplify the dynamics of the system.

Structural significance of hypermodified nucleic acid base hydroxywybutine (OHyW) which occur at 37th position in the anticodon loop of yeast tRNA(Phe)

Conformational preferences of hypermodified nucleic acid base hydroxywybutine (OHyW) have been studied using quantum chemical single point semi-empirical PM3 method. Automated geometry optimization using semi-empirical RM1, molecular mechanics force field (MMFF) along with ab-initio HF-SCF (6-31G** basis set) and DFT (B3LYP/6-31G** basis set) calculations have also been made to compare the salient features. Molecular electrostatic potentials (MEPs) depict the polarities of hydroxywybutine (OHyW) side chain. Another conformational study showed that hydroxywybutosine side chain interacts with adjacent bases within the anticodon loop of tRNA(Phe). The solvent accessible surface area (SASA) calculations revealed the structural role of hydroxywybutine in anticodon loop. Explicit molecular dynamics (MD) simulation has been done over the PM3 most stable structure of OHyW. The hydroxywybutine side chain prefers 'distal' conformation i.e. spreads away from the cyclic five membered imidazole moiety of modified tricyclic guanine base. The predicted preferred conformation of hydroxywybutine may prevent extended Watson-Crick base pairing during protein biosynthesis process. This conformation of OHyW stabilized by intramolecular interactions between O(6)⋯HO(16), O(6)⋯HC(15) and O(20)⋯HC(17). Further stabilization is also expected from interactions between O(22)⋯HC(16) and O(23)⋯HC(15). Explicit molecular dynamics (MD) simulation over the PM3 most stable structure of OHyW support the preferred geometry by preserving the 'distal' orientation of hydroxywybutine side chain and intramolecular hydrogen bonding interactions. MD simulation study revealed the role of hydroxyl group of OHyW to avoid fluctuations and prevent multiple iso-energetic conformations of hydroxywybutine side chain as compared to wybutine (yW).

Reprogramming of tRNA modifications controls the oxidative stress response by codon-biased translation of proteins

Selective translation of survival proteins is an important facet of the cellular stress response. We recently demonstrated that this translational control involves a stress-specific reprogramming of modified ribonucleosides in tRNA. Here we report the discovery of a step-wise translational control mechanism responsible for survival following oxidative stress. In yeast exposed to hydrogen peroxide, there is a Trm4 methyltransferase-dependent increase in the proportion of tRNA^LEU(CAA) containing m⁵C at the wobble position, which causes selective translation of mRNA from genes enriched in the TTG codon. Of these genes, oxidative stress increases protein expression from the TTG-enriched ribosomal protein gene RPL22A, but not its unenriched paralog. Loss of either TRM4 or RPL22A confers hypersensitivity to oxidative stress. Proteomic analysis reveals that oxidative stress causes a significant translational bias toward proteins coded by TTG-enriched genes. These results point to stress-induced reprogramming of tRNA modifications and consequential reprogramming of ribosomes in translational control of cell survival.

Structure and function of noncanonical nucleobases

DNA and RNA contain, next to the four canonical nucleobases, a number of modified nucleosides that extend their chemical information content. RNA is particularly rich in modifications, which is obviously an adaptation to their highly complex and variable functions. In fact, the modified nucleosides and their chemical structures establish a second layer of information which is of central importance to the function of the RNA molecules. Also the chemical diversity of DNA is greater than originally thought. Next to the four canonical bases, the DNA of higher organisms contains a total of four epigenetic bases: m(5) dC, hm(5) dC, f(5) dC und ca(5) dC. While all cells of an organism contain the same genetic material, their vastly different function and properties inside complex higher organisms require the controlled silencing and activation of cell-type specific genes. The regulation of the underlying silencing and activation process requires an additional layer of epigenetic information, which is clearly linked to increased chemical diversity. This diversity is provided by the modified non-canonical nucleosides in both DNA and RNA.

Molecular evolution of dihydrouridine synthases

BACKGROUND

Dihydrouridine (D) is a modified base found in conserved positions in the D-loop of tRNA in Bacteria, Eukaryota, and some Archaea. Despite the abundant occurrence of D, little is known about its biochemical roles in mediating tRNA function. It is assumed that D may destabilize the structure of tRNA and thus enhance its conformational flexibility. D is generated post-transcriptionally by the reduction of the 5,6-double bond of a uridine residue in RNA transcripts. The reaction is carried out by dihydrouridine synthases (DUS). DUS constitute a conserved family of enzymes encoded by the orthologous gene family COG0042. In protein sequence databases, members of COG0042 are typically annotated as "predicted TIM-barrel enzymes, possibly dehydrogenases, nifR3 family".

RESULTS

To elucidate sequence-structure-function relationships in the DUS family, a comprehensive bioinformatic analysis was carried out. We performed extensive database searches to identify all members of the currently known DUS family, followed by clustering analysis to subdivide it into subfamilies of closely related sequences. We analyzed phylogenetic distributions of all members of the DUS family and inferred the evolutionary tree, which suggested a scenario for the evolutionary origin of dihydrouridine-forming enzymes. For a human representative of the DUS family, the hDus2 protein suggested as a potential drug target in cancer, we generated a homology model. While this article was under review, a crystal structure of a DUS representative has been published, giving us an opportunity to validate the model.

CONCLUSIONS

We compared sequences and phylogenetic distributions of all members of the DUS family and inferred the phylogenetic tree, which provides a framework to study the functional differences among these proteins and suggests a scenario for the evolutionary origin of dihydrouridine formation. Our evolutionary and structural classification of the DUS family provides a background to study functional differences among these proteins that will guide experimental analyses.

RNA mediated Toll-like receptor stimulation in health and disease

Besides their well known functions in storage and translation of information nucleic acids have emerged as a target of pattern recognition receptors that drive activation of innate immunity. Due to the paucity of building block monomers used in nucleic acids, discrimination of host and microbial nucleic acids as a means of self/foreign discrimination is a complicated task. Pattern recognition receptors rely on discrimination by sequence, structural features and spatial compartmentalization to differentiate microbial derived nucleic acids from host ones. Microbial nucleic acid detection is important for the sensing of infectious danger and initiating an immune response to microbial attack. Failures in the underlying recognitions systems can have severe consequences: thus, inefficient recognition of microbial nucleic acids may increase susceptibility to infectious diseases. On the other hand, excessive immune responses as a result of failed self/foreign discrimination are associated with autoimmune diseases. This review gives a general overview over the underlying concepts of nucleic acid sensing by Toll-like receptors. Within this general framework, we focus on bacterial RNA and synthetic RNA oligomers.

Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq

An extensive repertoire of modifications is known to underlie the versatile coding, structural and catalytic functions of RNA, but it remains largely uncharted territory. Although biochemical studies indicate that N(6)-methyladenosine (m(6)A) is the most prevalent internal modification in messenger RNA, an in-depth study of its distribution and functions has been impeded by a lack of robust analytical methods. Here we present the human and mouse m(6)A modification landscape in a transcriptome-wide manner, using a novel approach, m(6)A-seq, based on antibody-mediated capture and massively parallel sequencing. We identify over 12,000 m(6)A sites characterized by a typical consensus in the transcripts of more than 7,000 human genes. Sites preferentially appear in two distinct landmarks--around stop codons and within long internal exons--and are highly conserved between human and mouse. Although most sites are well preserved across normal and cancerous tissues and in response to various stimuli, a subset of stimulus-dependent, dynamically modulated sites is identified. Silencing the m(6)A methyltransferase significantly affects gene expression and alternative splicing patterns, resulting in modulation of the p53 (also known as TP53) signalling pathway and apoptosis. Our findings therefore suggest that RNA decoration by m(6)A has a fundamental role in regulation of gene expression.

The archaeal COG1901/DUF358 SPOUT-methyltransferase members, together with pseudouridine synthase Pus10, catalyze the formation of 1-methylpseudouridine at position 54 of tRNA

The methylation of pseudouridine (Ψ) at position 54 of tRNA, producing m(1)Ψ, is a hallmark of many archaeal species, but the specific methylase involved in the formation of this modification had yet to be characterized. A comparative genomics analysis had previously identified COG1901 (DUF358), part of the SPOUT superfamily, as a candidate for this missing methylase family. To test this prediction, the COG1901 encoding gene, HVO_1989, was deleted from the Haloferax volcanii genome. Analyses of modified base contents indicated that while m(1)Ψ was present in tRNA extracted from the wild-type strain, it was absent from tRNA extracted from the mutant strain. Expression of the gene encoding COG1901 from Halobacterium sp. NRC-1, VNG1980C, complemented the m(1)Ψ minus phenotype of the ΔHVO_1989 strain. This in vivo validation was extended with in vitro tests. Using the COG1901 recombinant enzyme from Methanocaldococcus jannaschii (Mj1640), purified enzyme Pus10 from M. jannaschii and full-size tRNA transcripts or TΨ-arm (17-mer) fragments as substrates, the sequential pathway of m(1)Ψ54 formation in Archaea was reconstituted. The methylation reaction is AdoMet dependent. The efficiency of the methylase reaction depended on the identity of the residue at position 55 of the TΨ-loop. The presence of Ψ55 allowed the efficient conversion of Ψ54 to m(1)Ψ54, whereas in the presence of C55, the reaction was rather inefficient and no methylation reaction occurred if a purine was present at this position. These results led to renaming the Archaeal COG1901 members as TrmY proteins.

Mass spectrometry in the biology of RNA and its modifications

Many powerful analytical techniques for investigation of nucleic acids exist in the average modern molecular biology lab. The current review will focus on questions in RNA biology that have been answered by the use of mass spectrometry, which means that new biological information is the purpose and outcome of most of the studies we refer to. The review begins with a brief account of the subject "MS in the biology of RNA" and an overview of the prevalent RNA modifications identified to date. Fundamental considerations about mass spectrometric analysis of RNA are presented with the aim of detailing the analytical possibilities and challenges relating to the unique chemical nature of nucleic acids. The main biological topics covered are RNA modifications and the enzymes that perform the modifications. Modifications of RNA are essential in biology, and it is a field where mass spectrometry clearly adds knowledge of biological importance compared to traditional methods used in nucleic acid research. The biological applications are divided into analyses exclusively performed at the building block (mainly nucleoside) level and investigations involving mass spectrometry at the oligonucleotide level. We conclude the review discussing aspects of RNA identification and quantifications, which are upcoming fields for MS in RNA research. This article is part of a Special Section entitled: Understanding genome regulation and genetic diversity by mass spectrometry.

Identification of modifications in microbial, native tRNA that suppress immunostimulatory activity

2′-O-methylation of guanosine 18 is a naturally occurring tRNA modification that can suppress immune TLR7 responses.

Naturally occurring nucleotide modifications within RNA have been proposed to be structural determinants for innate immune recognition. We tested this hypothesis in the context of native nonself-RNAs. Isolated, fully modified native bacterial transfer RNAs (tRNAs) induced significant secretion of IFN-α from human peripheral blood mononuclear cells in a manner dependent on TLR7 and plasmacytoid dendritic cells. As a notable exception, tRNA^Tyr from Escherichia coli was not immunostimulatory, as were all tested eukaryotic tRNAs. However, the unmodified, 5′-unphosphorylated in vitro transcript of tRNA^Tyr induced IFN-α, thus revealing posttranscriptional modifications as a factor suppressing immunostimulation. Using a molecular surgery approach based on catalytic DNA, a panel of tRNA^Tyr variants featuring differential modification patterns was examined. Out of seven modifications present in this tRNA, 2′-O-methylated G_m18 was identified as necessary and sufficient to suppress immunostimulation. Transplantation of this modification into the scaffold of yeast tRNA^Phe also resulted in blocked immunostimulation. Moreover, an RNA preparation of an E. colitrmH mutant that lacks G_m18 2′-O-methyltransferase activity was significantly more stimulatory than the wild-type sample. The experiments identify the single methyl group on the 2′-oxygen of G_m18 as a natural modification in native tRNA that, beyond its primary structural role, has acquired a secondary function as an antagonist of TLR7.

RNA folding with soft constraints: reconciliation of probing data and thermodynamic secondary structure prediction

Thermodynamic folding algorithms and structure probing experiments are commonly used to determine the secondary structure of RNAs. Here we propose a formal framework to reconcile information from both prediction algorithms and probing experiments. The thermodynamic energy parameters are adjusted using ‘pseudo-energies’ to minimize the discrepancy between prediction and experiment. Our framework differs from related approaches that used pseudo-energies in several key aspects. (i) The energy model is only changed when necessary and no adjustments are made if prediction and experiment are consistent. (ii) Pseudo-energies remain biophysically interpretable and hold positional information where experiment and model disagree. (iii) The whole thermodynamic ensemble of structures is considered thus allowing to reconstruct mixtures of suboptimal structures from seemingly contradicting data. (iv) The noise of the energy model and the experimental data is explicitly modeled leading to an intuitive weighting factor through which the problem can be seen as folding with ‘soft’ constraints of different strength. We present an efficient algorithm to iteratively calculate pseudo-energies within this framework and demonstrate how this approach can be used in combination with SHAPE chemical probing data to improve secondary structure prediction. We further demonstrate that the pseudo-energies correlate with biophysical effects that are known to affect RNA folding such as chemical nucleotide modifications and protein binding.

Recent advances in developing small molecules targeting RNA

RNAs are underexploited targets for small molecule drugs or chemical probes of function. This may be due, in part, to a fundamental lack of understanding of the types of small molecules that bind RNA specifically and the types of RNA motifs that specifically bind small molecules. In this review, we describe recent advances in the development and design of small molecules that bind to RNA and modulate function that aim to fill this void.

YhiQ is RsmJ, the methyltransferase responsible for methylation of G1516 in 16S rRNA of E. coli

Ten methyltransferases and one pseudouridine synthase are required for complete modification of the small ribosomal subunit in Escherichia coli. Nine methyltransferases, as well as the pseudouridine synthase, are already known. Here, we identify RsmJ, the last unknown methyltransferase required for methylation of m(2)G1516 in 16S ribosomal RNA (rRNA), as the protein encoded by yhiQ. Reverse transcription primer extension analysis reveals that rRNA extracted from a yhiQ deletion strain is not methylated at G1516. Moreover, methylation is restored upon gene complementation. Also, purified recombinant YhiQ specifically methylates 30S subunits extracted from the deletion strain. The absence of the yhiQ gene leads to a cold-sensitive phenotype. Based on these data, we propose that the yhiQ gene be renamed rsmJ.

Decoding in Candidatus Riesia pediculicola, close to a minimal tRNA modification set?

A comparative genomic analysis of the recently sequenced human body louse unicellular endosymbiont Candidatus Riesia pediculicola with a reduced genome (582 Kb), revealed that it is the only known organism that might have lost all post-transcriptional base and ribose modifications of the tRNA body, retaining only modifications of the anticodon-stem-loop essential for mRNA decoding. Such a minimal tRNA modification set was not observed in other insect symbionts or in parasitic unicellular bacteria, such as Mycoplasma genitalium (580 Kb), that have also evolved by considerably reducing their genomes. This could be an example of a minimal tRNA modification set required for life, a question that has been at the center of the field for many years, especially for understanding the emergence and evolution of the genetic code.

Human Mitochondrial tRNAs: Biogenesis, Function, Structural Aspects, and Diseases

Mitochondria are eukaryotic organelles that generate most of the energy in the cell by oxidative phosphorylation (OXPHOS). Each mitochondrion contains multiple copies of a closed circular double-stranded DNA genome (mtDNA). Human (mammalian) mtDNA encodes 13 essential subunits of the inner membrane complex responsible for OXPHOS. These mRNAs are translated by the mitochondrial protein synthesis machinery, which uses the 22 species of mitochondrial tRNAs (mt tRNAs) encoded by mtDNA. The unique structural features of mt tRNAs distinguish them from cytoplasmic tRNAs bearing the canonical cloverleaf structure. The genes encoding mt tRNAs are highly susceptible to point mutations, which are a primary cause of mitochondrial dysfunction and are associated with a wide range of pathologies. A large number of nuclear factors involved in the biogenesis and function of mt tRNAs have been identified and characterized, including processing endonucleases, tRNA-modifying enzymes, and aminoacyl-tRNA synthetases. These nuclear factors are also targets of pathogenic mutations linked to various diseases, indicating the functional importance of mt tRNAs for mitochondrial activity.

Of P and Z: mitochondrial tRNA processing enzymes

Mitochondrial tRNAs are generally synthesized as part of polycistronic transcripts. Release of tRNAs from these precursors is thus not only required to produce functional adaptors for translation, but also responsible for the maturation of other mitochondrial RNA species. Cleavage of mitochondrial tRNAs appears to be exclusively accomplished by endonucleases. 5'-end maturation in the mitochondria of different Eukarya is achieved by various kinds of RNase P, representing the full range of diversity found in this enzyme family. While ribonucleoprotein enzymes with RNA components of bacterial-like appearance are found in a few unrelated protists, algae, and fungi, highly degenerate RNAs of dramatic size variability are found in the mitochondria of many fungi. The majority of mitochondrial RNase P enzymes, however, appear to be pure protein enzymes. Human mitochondrial RNase P, the first to be identified and possibly the prototype of all animal mitochondrial RNases P, is composed of three proteins. Homologs of its nuclease subunit MRPP3/PRORP, are also found in plants, algae and several protists, where they are apparently responsible for RNase P activity in mitochondria (and beyond) without the help of extra subunits. The diversity of RNase P enzymes is contrasted by the uniformity of mitochondrial RNases Z, which are responsible for 3'-end processing. Only the long form of RNase Z, which is restricted to eukarya, is found in mitochondria, even when an additional short form is present in the same organism. Mitochondrial tRNA processing thus appears dominated by new, eukaryal inventions rather than bacterial heritage. This article is part of a Special Issue entitled: Mitochondrial Gene Expression.

Comprehensive survey and geometric classification of base triples in RNA structures

Base triples are recurrent clusters of three RNA nucleobases interacting edge-to-edge by hydrogen bonding. We find that the central base in almost all triples forms base pairs with the other two bases of the triple, providing a natural way to geometrically classify base triples. Given 12 geometric base pair families defined by the Leontis-Westhof nomenclature, combinatoric enumeration predicts 108 potential geometric base triple families. We searched representative atomic-resolution RNA 3D structures and found instances of 68 of the 108 predicted base triple families. Model building suggests that some of the remaining 40 families may be unlikely to form for steric reasons. We developed an on-line resource that provides exemplars of all base triples observed in the structure database and models for unobserved, predicted triples, grouped by triple family, as well as by three-base combination (http://rna.bgsu.edu/Triples). The classification helps to identify recurrent triple motifs that can substitute for each other while conserving RNA 3D structure, with applications in RNA 3D structure prediction and analysis of RNA sequence evolution.

Informatics for mass spectrometry-based RNA analysis

Mass spectrometry (MS) allows the sensitive and direct characterization of biological macromolecules and therefore has the potential to complement the more conventional genetic and biochemical methods used for RNA characterization. Although MS has been used much less frequently for RNA research than it has been for protein research, recent technical improvements in both instrumentation and software make MS a powerful tool for RNA analysis because it can now be used to sequence, quantify, and chemically analyze RNAs. Mass spectrometry is particularly well suited for the characterization of RNAs associated with ribonucleoprotein complexes. This review focuses on the software and databases that can be used for MS-based RNA studies. Software for the processing of raw mass spectra, the identification and characterization of RNAs by mass mapping, de novo sequencing, and tandem MS-based database searching are available.

RNAcentral: A vision for an international database of RNA sequences

During the last decade there has been a great increase in the number of noncoding RNA genes identified, including new classes such as microRNAs and piRNAs. There is also a large growth in the amount of experimental characterization of these RNA components. Despite this growth in information, it is still difficult for researchers to access RNA data, because key data resources for noncoding RNAs have not yet been created. The most pressing omission is the lack of a comprehensive RNA sequence database, much like UniProt, which provides a comprehensive set of protein knowledge. In this article we propose the creation of a new open public resource that we term RNAcentral, which will contain a comprehensive collection of RNA sequences and fill an important gap in the provision of biomedical databases. We envision RNA researchers from all over the world joining a federated RNAcentral network, contributing specialized knowledge and databases. RNAcentral would centralize key data that are currently held across a variety of databases, allowing researchers instant access to a single, unified resource. This resource would facilitate the next generation of RNA research and help drive further discoveries, including those that improve food production and human and animal health. We encourage additional RNA database resources and research groups to join this effort. We aim to obtain international network funding to further this endeavor.

Identification and modeling of a phosphatase-like domain in a tRNA 2'-O-ribosyl phosphate transferase Rit1p

Cytoplasmic initiator tRNAs from plants and fungi are excluded from participating in translational elongation by the presence of a unique 2'-phosphoribosyl modification of purine 64, introduced posttranscriptionally by the enzyme Rit1p. Members of the Rit1p family show no obvious similarity to other proteins or domains, there is no structural information available to guide experimental analyses, and the mechanism of action of this enzyme remains a mystery. Using protein fold recognition, we identified a phosphatase-like domain in the C-terminal part of Rit1p. A comparative model of the C-terminal domain was constructed and used to predict the function of conserved residues and to propose the mechanism of action of Rit1p. The model will facilitate experimental analyses of Rit1p and its interactions with the initiator tRNA substrate.

Potential for interdependent development of tRNA determinants for aminoacylation and ribosome decoding

Although the nucleotides in tRNA required for aminoacylation are conserved in evolution, bacterial aminoacyl-tRNA synthetases (aaRSs) are unable to acylate eukaryotic tRNA. The cross-species barrier may be due to the absence of eukaryote-specific domains from bacterial aaRSs. Here we show that while E. coli CysRS cannot acylate human tRNA^Cys, the fusion of a eukaryote-specific domain of human CysRS overcomes the cross-species barrier in human tRNA^Cys. In addition to enabling recognition of the sequence differences in the tertiary core of tRNA^Cys, the fused eukaryotic domain redirects the specificity of E. coli CysRS from the A37 present in bacterial tRNA^Cys to the G37 in mammals. Further experiments show that the accuracy of codon recognition on the ribosome was also highly sensitive to the A37-to-G37 transition in tRNA^Cys. These results raise the possibility of the development of tRNA nucleotide determinants for aminoacylation being interdependent with those for ribosome decoding.

RNA tertiary structure prediction with ModeRNA

Noncoding RNAs perform important roles in the cell. As their function is tightly connected with structure, and as experimental methods are time-consuming and expensive, the field of RNA structure prediction is developing rapidly. Here, we present a detailed study on using the ModeRNA software. The tool uses the comparative modeling approach and can be applied when a structural template is available and an alignment of reasonable quality can be performed. We guide the reader through the entire process of modeling Escherichia coli tRNA(Thr) in a conformation corresponding to the complex with an aminoacyl-tRNA synthetase (aaRS). We describe the choice of a template structure, preparation of input files, and explore three possible modeling strategies. In the end, we evaluate the resulting models using six alternative benchmarks. The ModeRNA software can be freely downloaded from http://iimcb.genesilico.pl/moderna/ under the conditions of the General Public License. It runs under LINUX, Windows and Mac OS. It is also available as a server at http://iimcb.genesilico.pl/modernaserver/. The models and the script to reproduce the study from this article are available at http://www.genesilico.pl/moderna/examples/.

Evolution of the deaminase fold and multiple origins of eukaryotic editing and mutagenic nucleic acid deaminases from bacterial toxin systems

The deaminase-like fold includes, in addition to nucleic acid/nucleotide deaminases, several catalytic domains such as the JAB domain, and others involved in nucleotide and ADP-ribose metabolism. Using sensitive sequence and structural comparison methods, we develop a comprehensive natural classification of the deaminase-like fold and show that its ancestral version was likely to operate on nucleotides or nucleic acids. Consequently, we present evidence that a specific group of JAB domains are likely to possess a DNA repair function, distinct from the previously known deubiquitinating peptidase activity. We also identified numerous previously unknown clades of nucleic acid deaminases. Using inference based on contextual information, we suggest that most of these clades are toxin domains of two distinct classes of bacterial toxin systems, namely polymorphic toxins implicated in bacterial interstrain competition and those that target distantly related cells. Genome context information suggests that these toxins might be delivered via diverse secretory systems, such as Type V, Type VI, PVC and a novel PrsW-like intramembrane peptidase-dependent mechanism. We propose that certain deaminase toxins might be deployed by diverse extracellular and intracellular pathogens as also endosymbionts as effectors targeting nucleic acids of host cells. Our analysis suggests that these toxin deaminases have been acquired by eukaryotes on several independent occasions and recruited as organellar or nucleo-cytoplasmic RNA modifiers, operating on tRNAs, mRNAs and short non-coding RNAs, and also as mutators of hyper-variable genes, viruses and selfish elements. This scenario potentially explains the origin of mutagenic AID/APOBEC-like deaminases, including novel versions from Caenorhabditis, Nematostella and diverse algae and a large class of fast-evolving fungal deaminases. These observations greatly expand the distribution of possible unidentified mutagenic processes catalyzed by nucleic acid deaminases.

The 2-thiouridylase function of the human MTU1 (TRMU) enzyme is dispensable for mitochondrial translation

MTU1 (TRMU) is a mitochondrial enzyme responsible for the 2-thiolation of the wobble U in tRNA(Lys), tRNA(Glu) and tRNA(Gln), a post-transcriptional modification believed to be important for accurate and efficient synthesis of the 13 respiratory chain subunits encoded by mtDNA. Mutations in MTU1 are associated with acute infantile liver failure, and this has been ascribed to a transient lack of cysteine, the sulfur donor for the thiouridylation reaction, resulting in a mitochondrial translation defect during early development. A mutation in tRNA(Lys) that causes myoclonic epilepsy with ragged-red fibers (MERRF) is also reported to prevent modification of the wobble U. Here we show that mitochondrial translation is unaffected in fibroblasts from an MTU1 patient, in which MTU1 is undetectable by immunoblotting, despite the severe reduction in the 2-thiolation of mitochondrial tRNA(Lys), tRNA(Glu) and tRNA(Gln). The only respiratory chain abnormality that we could observe in these cells was an accumulation of a Complex II assembly intermediate, which, however, did not affect the level of the fully assembled enzyme. The identical phenotype was observed by siRNA-mediated knockdown of MTU1 in HEK 293 cells. Further, the mitochondrial translation deficiencies present in myoblasts from mitochondrial encephalomyopathy, lactic acidosis and stroke-like episode and MERRF patients, which are associated with defects in post-transcriptional modification of mitochondrial tRNAs, did not worsen following knockdown of MTU1 in these cells. This study demonstrates that MTU1 is not required for mitochondrial translation at normal steady-state levels of tRNAs, and that it may possess an as yet uncharacterized function in another sulfur-trafficking pathway.

A toolbox for developing bioinformatics software

Creating useful software is a major activity of many scientists, including bioinformaticians. Nevertheless, software development in an academic setting is often unsystematic, which can lead to problems associated with maintenance and long-term availibility. Unfortunately, well-documented software development methodology is difficult to adopt, and technical measures that directly improve bioinformatic programming have not been described comprehensively. We have examined 22 software projects and have identified a set of practices for software development in an academic environment. We found them useful to plan a project, support the involvement of experts (e.g. experimentalists), and to promote higher quality and maintainability of the resulting programs. This article describes 12 techniques that facilitate a quick start into software engineering. We describe 3 of the 22 projects in detail and give many examples to illustrate the usage of particular techniques. We expect this toolbox to be useful for many bioinformatics programming projects and to the training of scientific programmers.

Identification of N6,N6-dimethyladenosine in transfer RNA from Mycobacterium bovis Bacille Calmette-Guérin

There are more than 100 different ribonucleoside structures incorporated as post-transcriptional modifications, mainly in tRNA and rRNA of both prokaryotes and eukaryotes, and emerging evidence suggests that these modifications function as a system in the translational control of cellular responses. However, our understanding of this system is hampered by the paucity of information about the complete set of RNA modifications present in individual organisms. To this end, we have employed a chromatography-coupled mass spectrometric approach to define the spectrum of modified ribonucleosides in microbial species, starting with Mycobacterium bovis BCG. This approach revealed a variety of ribonucleoside candidates in tRNA from BCG, of which 12 were definitively identified based on comparisons to synthetic standards and 5 were tentatively identified by exact mass comparisons to RNA modification databases. Among the ribonucleosides observed in BCG tRNA was one not previously described in tRNA, which we have now characterized as N⁶,N⁶-dimethyladenosine.

Functional consequences of T-stem mutations in E. coli tRNAThrUGU in vitro and in vivo

The binding affinities between Escherichia coli EF-Tu and 34 single and double base-pair changes in the T stem of E. coli tRNA(Thr)(UGU) were compared with similar data obtained previously for several aa-tRNAs binding to Thermus thermophilus EF-Tu. With a single exception, the two proteins bound to mutations in three T-stem base pairs in a quantitatively identical manner. However, tRNA(Thr) differs from other tRNAs by also using its rare A52-C62 pair as a negative specificity determinant. Using a plasmid-based tRNA gene replacement strategy, we show that many of the tRNA(Thr)(UGU) T-stem changes are either unable to support growth of E. coli or are less effective than the wild-type sequence. Since the inviable T-stem sequences are often present in other E. coli tRNAs, it appears that T-stem sequences in each tRNA body have evolved to optimize function in a different way. Although mutations of tRNA(Thr) can substantially increase or decrease its affinity to EF-Tu, the observed affinities do not correlate with the growth phenotype of the mutations in any simple way. This may either reflect the different conditions used in the two assays or indicate that the T-stem mutants affect another step in the translation mechanism.

Biosynthesis and function of tRNA modifications in Archaea

tRNA modifications are important for decoding, translation accuracy, and structural integrity of tRNAs. Archaeal tRNAs contain at least 47 different tRNA modifications, some of them, including archaeosine, agmatidine, and mimG, are specific to the archaeal domain. The biosynthetic pathways for these complex signature modifications have recently been elucidated and are extensively described in this review. Archaeal organisms still lag Escherichia coli and Saccharomyces cerevisiae in terms of genetic characterization and in vivo function of tRNA modifications. However, recent advances in the model Haloferax volcanii, described here, should allow closing this gap soon. Consequently, an update on experimental characterizations of archaeal tRNA modification genes and proteins is given to set the stage for future work in this field.

RNA nucleotide methylation

Methylation of RNA occurs at a variety of atoms, nucleotides, sequences and tertiary structures. Strongly related to other posttranscriptional modifications, methylation of different RNA species includes tRNA, rRNA, mRNA, tmRNA, snRNA, snoRNA, miRNA, and viral RNA. Different catalytic strategies are employed for RNA methylation by a variety of RNA-methyltransferases which fall into four superfamilies. This review outlines the different functions of methyl groups in RNA, including biophysical, biochemical and metabolic stabilization of RNA, quality control, resistance to antibiotics, mRNA reading frame maintenance, deciphering of normal and altered genetic code, selenocysteine incorporation, tRNA aminoacylation, ribotoxins, splicing, intracellular trafficking, immune response, and others. Connections to other fields including gene regulation, DNA repair, stress response, and possibly histone acetylation and exocytosis are pointed out. WIREs RNA 2011 2 611-631 DOI: 10.1002/wrna.79 For further resources related to this article, please visit the WIREs website.

ModeRNA: a tool for comparative modeling of RNA 3D structure

RNA is a large group of functionally important biomacromolecules. In striking analogy to proteins, the function of RNA depends on its structure and dynamics, which in turn is encoded in the linear sequence. However, while there are numerous methods for computational prediction of protein three-dimensional (3D) structure from sequence, with comparative modeling being the most reliable approach, there are very few such methods for RNA. Here, we present ModeRNA, a software tool for comparative modeling of RNA 3D structures. As an input, ModeRNA requires a 3D structure of a template RNA molecule, and a sequence alignment between the target to be modeled and the template. It must be emphasized that a good alignment is required for successful modeling, and for large and complex RNA molecules the development of a good alignment usually requires manual adjustments of the input data based on previous expertise of the respective RNA family. ModeRNA can model post-transcriptional modifications, a functionally important feature analogous to post-translational modifications in proteins. ModeRNA can also model DNA structures or use them as templates. It is equipped with many functions for merging fragments of different nucleic acid structures into a single model and analyzing their geometry. Windows and UNIX implementations of ModeRNA with comprehensive documentation and a tutorial are freely available.

Natural history of eukaryotic DNA methylation systems

Methylation of cytosines and adenines in DNA is a widespread epigenetic mark in both prokaryotes and eukaryotes. In eukaryotes, it has a profound influence on chromatin structure and dynamics. Recent advances in genomics and biochemistry have considerably elucidated the functions and provenance of these DNA modifications. DNA methylases appear to have emerged first in bacterial restriction-modification (R-M) systems from ancient RNA-modifying enzymes, in transitions that involved acquisition of novel catalytic residues and DNA-recognition features. DNA adenine methylases appear to have been acquired by ciliates, heterolobosean amoeboflagellates, and certain chlorophyte algae. Six distinct clades of cytosine methylases, including the DNMT1, DNMT2, and DNMT3 clades, were acquired by eukaryotes through independent lateral transfer of their precursors from bacteria or bacteriophages. In addition to these, multiple adenine and cytosine methylases were acquired by several families of eukaryotic transposons. In eukaryotes, the DNA-methylase module was often combined with distinct modified and unmodified peptide recognition domains and other modules mediating specialized interactions, for example, the RFD module of DNMT1 which contains a permuted Sm domain linked to a helix-turn-helix domain. In eukaryotes, the evolution of DNA methylases appears to have proceeded in parallel to the elaboration of histone-modifying enzymes and the RNAi system, with functions related to counter-viral and counter-transposon defense, and regulation of DNA repair and differential gene expression being their primary ancestral functions. Diverse DNA demethylation systems that utilize base-excision repair via DNA glycosylases and cytosine deaminases appear to have emerged in multiple eukaryotic lineages. Comparative genomics suggests that the link between cytosine methylation and DNA glycosylases probably emerged first in a novel R-M system in bacteria. Recent studies suggest that the 5mC is not a terminal DNA modification, with enzymes of the Tet/JBP family of 2-oxoglutarate- and iron-dependent dioxygenases further hydroxylating it to form 5-hydroxymethylcytosine (5hmC). These enzymes emerged first in bacteriophages and appear to have been transferred to eukaryotes on one or more occasions. Eukaryotes appear to have recruited three major types of DNA-binding domains (SRA/SAD, TAM/MBD, and CXXC) in discriminating DNA with methylated or unmethylated cytosines. Analysis of the domain architectures of these domains and the DNA methylases suggests that early in eukaryotic evolution they developed a close functional link with SET-domain methylases and Jumonji-related demethylases that operate on peptides in chromatin proteins. In several eukaryotes, other functional connections were elaborated in the form of various combinations between domains related to DNA methylation and those involved in ATP-dependent chromatin remodeling and RNAi. In certain eukaryotes, such as mammals and angiosperms, novel dependencies on the DNA methylation system emerged, which resulted in it affecting unexpected aspects of the biology of these organisms such as parent-offspring interactions. In genomic terms, this was reflected in the emergence of new proteins related to methylation, such as Stella. The well-developed methylation systems of certain heteroloboseans, stramenopiles, chlorophytes, and haptophyte indicate that these might be new model systems to explore the relevance of DNA modifications in eukaryotes.

A quantitative systems approach reveals dynamic control of tRNA modifications during cellular stress

Decades of study have revealed more than 100 ribonucleoside structures incorporated as post-transcriptional modifications mainly in tRNA and rRNA, yet the larger functional dynamics of this conserved system are unclear. To this end, we developed a highly precise mass spectrometric method to quantify tRNA modifications in Saccharomyces cerevisiae. Our approach revealed several novel biosynthetic pathways for RNA modifications and led to the discovery of signature changes in the spectrum of tRNA modifications in the damage response to mechanistically different toxicants. This is illustrated with the RNA modifications Cm, m(5)C, and m(2) (2)G, which increase following hydrogen peroxide exposure but decrease or are unaffected by exposure to methylmethane sulfonate, arsenite, and hypochlorite. Cytotoxic hypersensitivity to hydrogen peroxide is conferred by loss of enzymes catalyzing the formation of Cm, m(5)C, and m(2) (2)G, which demonstrates that tRNA modifications are critical features of the cellular stress response. The results of our study support a general model of dynamic control of tRNA modifications in cellular response pathways and add to the growing repertoire of mechanisms controlling translational responses in cells.

Mechanism of N-methylation by the tRNA m1G37 methyltransferase Trm5

Trm5 is a eukaryal and archaeal tRNA methyltransferase that catalyzes methyl transfer from S-adenosylmethionine (AdoMet) to the N(1) position of G37 directly 3' to the anticodon. While the biological role of m(1)G37 in enhancing translational fidelity is well established, the catalytic mechanism of Trm5 has remained obscure. To address the mechanism of Trm5 and more broadly the mechanism of N-methylation to nucleobases, we examined the pH-activity profile of an archaeal Trm5 enzyme, and performed structure-guided mutational analysis. The data reveal a marked dependence of enzyme-catalyzed methyl transfer on hydrogen ion equilibria: the single-turnover rate constant for methylation increases by one order of magnitude from pH 6.0 to reach a plateau at pH 7.0. This suggests a mechanism involving proton transfer from G37 as the key element in catalysis. Consideration of the kinetic data in light of the Trm5-tRNA-AdoMet ternary cocrystal structure, determined in a precatalytic conformation, suggests that proton transfer is associated with an induced fit rearrangement of the complex that precedes formation of the reactive configuration in the active site. Key roles for the conserved R145 side chain in stabilizing a proposed oxyanion at G37-O(6), and for E185 as a general base to accept the proton from G37-N(1), are suggested based on the mutational analysis.

The RNA Modification Database, RNAMDB: 2011 update

Since its inception in 1994, The RNA Modification Database (RNAMDB, http://rna-mdb.cas.albany.edu/RNAmods/) has served as a focal point for information pertaining to naturally occurring RNA modifications. In its current state, the database employs an easy-to-use, searchable interface for obtaining detailed data on the 109 currently known RNA modifications. Each entry provides the chemical structure, common name and symbol, elemental composition and mass, CA registry numbers and index name, phylogenetic source, type of RNA species in which it is found, and references to the first reported structure determination and synthesis. Though newly transferred in its entirety to The RNA Institute, the RNAMDB continues to grow with two notable additions, agmatidine and 8-methyladenosine, appended in the last year. The RNA Modification Database is staying up-to-date with significant improvements being prepared for inclusion within the next year and the following year. The expanded future role of The RNA Modification Database will be to serve as a primary information portal for researchers across the entire spectrum of RNA-related research.

REPAIRtoire--a database of DNA repair pathways

REPAIRtoire is the first comprehensive database resource for systems biology of DNA damage and repair. The database collects and organizes the following types of information: (i) DNA damage linked to environmental mutagenic and cytotoxic agents, (ii) pathways comprising individual processes and enzymatic reactions involved in the removal of damage, (iii) proteins participating in DNA repair and (iv) diseases correlated with mutations in genes encoding DNA repair proteins. REPAIRtoire provides also links to publications and external databases. REPAIRtoire contains information about eight main DNA damage checkpoint, repair and tolerance pathways: DNA damage signaling, direct reversal repair, base excision repair, nucleotide excision repair, mismatch repair, homologous recombination repair, nonhomologous end-joining and translesion synthesis. The pathway/protein dataset is currently limited to three model organisms: Escherichia coli, Saccharomyces cerevisiae and Homo sapiens. The DNA repair and tolerance pathways are represented as graphs and in tabular form with descriptions of each repair step and corresponding proteins, and individual entries are cross-referenced to supporting literature and primary databases. REPAIRtoire can be queried by the name of pathway, protein, enzymatic complex, damage and disease. In addition, a tool for drawing custom DNA-protein complexes is available online. REPAIRtoire is freely available and can be accessed at http://repairtoire.genesilico.pl/.

Crystal structure of a novel JmjC-domain-containing protein, TYW5, involved in tRNA modification

Wybutosine (yW) is a hypermodified nucleoside found in position 37 of tRNA(Phe), and is essential for correct phenylalanine codon translation. yW derivatives widely exist in eukaryotes and archaea, and their chemical structures have many species-specific variations. Among them, its hydroxylated derivative, hydroxywybutosine (OHyW), is found in eukaryotes including human, but the modification mechanism remains unknown. Recently, we identified a novel Jumonji C (JmjC)-domain-containing protein, TYW5 (tRNA yW-synthesizing enzyme 5), which forms the OHyW nucleoside by carbon hydroxylation, using Fe(II) ion and 2-oxoglutarate (2-OG) as cofactors. In this work, we present the crystal structures of human TYW5 (hTYW5) in the free and complex forms with 2-OG and Ni(II) ion at 2.5 and 2.8 Å resolutions, respectively. The structure revealed that the catalytic domain consists of a β-jellyroll fold, a hallmark of the JmjC domains and other Fe(II)/2-OG oxygenases. hTYW5 forms a homodimer through C-terminal helix bundle formation, thereby presenting a large, positively charged patch involved in tRNA binding. A comparison with the structures of other JmjC-domain-containing proteins suggested a mechanism for substrate nucleotide recognition. Functional analyses of structure-based mutants revealed the essential Arg residues participating in tRNA recognition by TYW5. These findings extend the repertoire of the tRNA modification enzyme into the Fe(II)/2-OG oxygenase superfamily.

Fidelity in archaeal information processing

A key element during the flow of genetic information in living systems is fidelity. The accuracy of DNA replication influences the genome size as well as the rate of genome evolution. The large amount of energy invested in gene expression implies that fidelity plays a major role in fitness. On the other hand, an increase in fidelity generally coincides with a decrease in velocity. Hence, an important determinant of the evolution of life has been the establishment of a delicate balance between fidelity and variability. This paper reviews the current knowledge on quality control in archaeal information processing. While the majority of these processes are homologous in Archaea, Bacteria, and Eukaryotes, examples are provided of nonorthologous factors and processes operating in the archaeal domain. In some instances, evidence for the existence of certain fidelity mechanisms has been provided, but the factors involved still remain to be identified.

The critical role of RNA processing and degradation in the control of gene expression

The continuous degradation and synthesis of prokaryotic mRNAs not only give rise to the metabolic changes that are required as cells grow and divide but also rapid adaptation to new environmental conditions. In bacteria, RNAs can be degraded by mechanisms that act independently, but in parallel, and that target different sites with different efficiencies. The accessibility of sites for degradation depends on several factors, including RNA higher-order structure, protection by translating ribosomes and polyadenylation status. Furthermore, RNA degradation mechanisms have shown to be determinant for the post-transcriptional control of gene expression. RNases mediate the processing, decay and quality control of RNA. RNases can be divided into endonucleases that cleave the RNA internally or exonucleases that cleave the RNA from one of the extremities. Just in Escherichia coli there are >20 different RNases. RNase E is a single-strand-specific endonuclease critical for mRNA decay in E. coli. The enzyme interacts with the exonuclease polynucleotide phosphorylase (PNPase), enolase and RNA helicase B (RhlB) to form the degradosome. However, in Bacillus subtilis, this enzyme is absent, but it has other main endonucleases such as RNase J1 and RNase III. RNase III cleaves double-stranded RNA and family members are involved in RNA interference in eukaryotes. RNase II family members are ubiquitous exonucleases, and in eukaryotes, they can act as the catalytic subunit of the exosome. RNases act in different pathways to execute the maturation of rRNAs and tRNAs, and intervene in the decay of many different mRNAs and small noncoding RNAs. In general, RNases act as a global regulatory network extremely important for the regulation of RNA levels.

New archaeal methyltransferases forming 1-methyladenosine or 1-methyladenosine and 1-methylguanosine at position 9 of tRNA

Two archaeal tRNA methyltransferases belonging to the SPOUT superfamily and displaying unexpected activities are identified. These enzymes are orthologous to the yeast Trm10p methyltransferase, which catalyses the formation of 1-methylguanosine at position 9 of tRNA. In contrast, the Trm10p orthologue from the crenarchaeon Sulfolobus acidocaldarius forms 1-methyladenosine at the same position. Even more surprisingly, the Trm10p orthologue from the euryarchaeon Thermococcus kodakaraensis methylates the N¹-atom of either adenosine or guanosine at position 9 in different tRNAs. This is to our knowledge the first example of a tRNA methyltransferase with a broadened nucleoside recognition capability. The evolution of tRNA methyltransferases methylating the N¹ atom of a purine residue is discussed.

Functional characterization of the YmcB and YqeV tRNA methylthiotransferases of Bacillus subtilis

Methylthiotransferases (MTTases) are a closely related family of proteins that perform both radical-S-adenosylmethionine (SAM) mediated sulfur insertion and SAM-dependent methylation to modify nucleic acid or protein targets with a methyl thioether group (-SCH(3)). Members of two of the four known subgroups of MTTases have been characterized, typified by MiaB, which modifies N(6)-isopentenyladenosine (i(6)A) to 2-methylthio-N(6)-isopentenyladenosine (ms(2)i(6)A) in tRNA, and RimO, which modifies a specific aspartate residue in ribosomal protein S12. In this work, we have characterized the two MTTases encoded by Bacillus subtilis 168 and find that, consistent with bioinformatic predictions, ymcB is required for ms(2)i(6)A formation (MiaB activity), and yqeV is required for modification of N(6)-threonylcarbamoyladenosine (t(6)A) to 2-methylthio-N(6)-threonylcarbamoyladenosine (ms(2)t(6)A) in tRNA. The enzyme responsible for the latter activity belongs to a third MTTase subgroup, no member of which has previously been characterized. We performed domain-swapping experiments between YmcB and YqeV to narrow down the protein domain(s) responsible for distinguishing i(6)A from t(6)A and found that the C-terminal TRAM domain, putatively involved with RNA binding, is likely not involved with this discrimination. Finally, we performed a computational analysis to identify candidate residues outside the TRAM domain that may be involved with substrate recognition. These residues represent interesting targets for further analysis.

RNA FRABASE 2.0: an advanced web-accessible database with the capacity to search the three-dimensional fragments within RNA structures

BACKGROUND

Recent discoveries concerning novel functions of RNA, such as RNA interference, have contributed towards the growing importance of the field. In this respect, a deeper knowledge of complex three-dimensional RNA structures is essential to understand their new biological functions. A number of bioinformatic tools have been proposed to explore two major structural databases (PDB, NDB) in order to analyze various aspects of RNA tertiary structures. One of these tools is RNA FRABASE 1.0, the first web-accessible database with an engine for automatic search of 3D fragments within PDB-derived RNA structures. This search is based upon the user-defined RNA secondary structure pattern. In this paper, we present and discuss RNA FRABASE 2.0. This second version of the system represents a major extension of this tool in terms of providing new data and a wide spectrum of novel functionalities. An intuitionally operated web server platform enables very fast user-tailored search of three-dimensional RNA fragments, their multi-parameter conformational analysis and visualization.

DESCRIPTION

RNA FRABASE 2.0 has stored information on 1565 PDB-deposited RNA structures, including all NMR models. The RNA FRABASE 2.0 search engine algorithms operate on the database of the RNA sequences and the new library of RNA secondary structures, coded in the dot-bracket format extended to hold multi-stranded structures and to cover residues whose coordinates are missing in the PDB files. The library of RNA secondary structures (and their graphics) is made available. A high level of efficiency of the 3D search has been achieved by introducing novel tools to formulate advanced searching patterns and to screen highly populated tertiary structure elements. RNA FRABASE 2.0 also stores data and conformational parameters in order to provide "on the spot" structural filters to explore the three-dimensional RNA structures. An instant visualization of the 3D RNA structures is provided. RNA FRABASE 2.0 is freely available at http://rnafrabase.cs.put.poznan.pl.

CONCLUSIONS

RNA FRABASE 2.0 provides a novel database and powerful search engine which is equipped with new data and functionalities that are unavailable elsewhere. Our intention is that this advanced version of the RNA FRABASE will be of interest to all researchers working in the RNA field.

Origin and evolution of peptide-modifying dioxygenases and identification of the wybutosine hydroxylase/hydroperoxidase

Unlike classical 2-oxoglutarate and iron-dependent dioxygenases, which include several nucleic acid modifiers, the structurally similar jumonji-related dioxygenase superfamily was only known to catalyze peptide modifications. Using comparative genomics methods, we predict that a family of jumonji-related enzymes catalyzes wybutosine hydroxylation/peroxidation at position 37 of eukaryotic tRNAPhe. Identification of this enzyme raised questions regarding the emergence of protein- and nucleic acid-modifying activities among jumonji-related domains. We addressed these with a natural classification of DSBH domains and reconstructed the precursor of the dioxygenases as a sugar-binding domain. This precursor gave rise to sugar epimerases and metal-binding sugar isomerases. The sugar isomerase active site was exapted for catalysis of oxygenation, with a radiation of these enzymes in bacteria, probably due to impetus from the primary oxygenation event in Earth’s history. 2-Oxoglutarate-dependent versions appear to have further expanded with rise of the tricarboxylic acid cycle. We identify previously under-appreciated aspects of their active site and multiple independent innovations of 2-oxoacid-binding basic residues among these superfamilies. We show that double-stranded β-helix dioxygenases diversified extensively in biosynthesis and modification of halogenated siderophores, antibiotics, peptide secondary metabolites and glycine-rich collagen-like proteins in bacteria. Jumonji-related domains diversified into three distinct lineages in bacterial secondary metabolism systems and these were precursors of the three major clades of eukaryotic enzymes. The specificity of wybutosine hydroxylase/peroxidase probably relates to the structural similarity of the modified moiety to the ancestral amino acid substrate of this superfamily.

Human AlkB homolog ABH8 Is a tRNA methyltransferase required for wobble uridine modification and DNA damage survival

tRNA nucleosides are extensively modified to ensure their proper function in translation. However, many of the enzymes responsible for tRNA modifications in mammals await identification. Here, we show that human AlkB homolog 8 (ABH8) catalyzes tRNA methylation to generate 5-methylcarboxymethyl uridine (mcm(5)U) at the wobble position of certain tRNAs, a critical anticodon loop modification linked to DNA damage survival. We find that ABH8 interacts specifically with tRNAs containing mcm(5)U and that purified ABH8 complexes methylate RNA in vitro. Significantly, ABH8 depletion in human cells reduces endogenous levels of mcm(5)U in RNA and increases cellular sensitivity to DNA-damaging agents. Moreover, DNA-damaging agents induce ABH8 expression in an ATM-dependent manner. These results expand the role of mammalian AlkB proteins beyond that of direct DNA repair and support a regulatory mechanism in the DNA damage response pathway involving modulation of tRNA modification.

Combination of the loss of cmnm5U34 with the lack of s2U34 modifications of tRNALys, tRNAGlu, and tRNAGln altered mitochondrial biogenesis and respiration

Yeast Saccharomyces cerevisiae MTO2, MTO1, and MSS1 genes encoded highly conserved tRNA modifying enzymes for the biosynthesis of carboxymethylaminomethyl (cmnm)(5)s(2)U(34) in mitochondrial tRNA(Lys), tRNA(Glu), and tRNA(Gln). In fact, Mto1p and Mss1p are involved in the biosynthesis of the cmnm(5) group (cmnm(5)U(34)), while Mto2p is responsible for the 2-thiouridylation (s(2)U(34)) of these tRNAs. Previous studies showed that partial modifications at U(34) in mitochondrial tRNA enabled mto1, mto2, and mss1 strains to respire. In this report, we investigated the functional interaction between MTO2, MTO1, and MSS1 genes by using the mto2, mto1, and mss1 single, double, and triple mutants. Strikingly, the deletion of MTO2 was synthetically lethal with a mutation of MSS1 or deletion of MTO1 on medium containing glycerol but not on medium containing glucose. Interestingly, there were no detectable levels of nine tRNAs including tRNA(Lys), tRNA(Glu), and tRNA(Gln) in mto2/mss1, mto2/mto1, and mto2/mto1/mss1 strains. Furthermore, mto2/mss1, mto2/mto1, and mto2/mto1/mss1 mutants exhibited extremely low levels of COX1 and CYTB mRNA and 15S and 21S rRNA as well as the complete loss of mitochondrial protein synthesis. The synthetic enhancement combinations likely resulted from the completely abolished modification at U(34) of tRNA(Lys), tRNA(Glu), and tRNA(Gln), caused by the combination of eliminating the 2-thiouridylation by the mto2 mutation with the absence of the cmnm(5)U(34) by the mto1 or mss1 mutation. The complete loss of modifications at U(34) of tRNAs altered mitochondrial RNA metabolisms, causing a degradation of mitochondrial tRNA, mRNA, and rRNAs. As a result, failures in mitochondrial RNA metabolisms were responsible for the complete loss of mitochondrial translation. Consequently, defects in mitochondrial protein synthesis caused the instability of their mitochondrial genomes, thus producing the respiratory-deficient phenotypes. Therefore, our findings demonstrated a critical role of modifications at U(34) of tRNA(Lys), tRNA(Glu), and tRNA(Gln) in maintenance of mitochondrial genome, mitochondrial RNA stability, translation, and respiratory function.