Major reorientation of tRNA substrates defines specificity of dihydrouridine synthases

Exploring a unique class of flavoenzymes: Identification and biochemical characterization of ribosomal RNA dihydrouridine synthase

Dihydrouridine (D), a prevalent and evolutionarily conserved base in the transcriptome, primarily resides in tRNAs and, to a lesser extent, in mRNAs. Notably, this modification is found at position 2449 in the Escherichia coli 23S rRNA, strategically positioned near the ribosome's peptidyl transferase site. Despite the prior identification, in E. coli genome, of three dihydrouridine synthases (DUS), a set of NADPH and FMN-dependent enzymes known for introducing D in tRNAs and mRNAs, characterization of the enzyme responsible for D2449 deposition has remained elusive. This study introduces a rapid method for detecting D in rRNA, involving reverse transcriptase-blockage at the rhodamine-labeled D2449 site, followed by PCR amplification (RhoRT-PCR). Through analysis of rRNA from diverse E. coli strains, harboring chromosomal or single-gene deletions, we pinpoint the yhiN gene as the ribosomal dihydrouridine synthase, now designated as RdsA. Biochemical characterizations uncovered RdsA as a unique class of flavoenzymes, dependent on FAD and NADH, with a complex structural topology. In vitro assays demonstrated that RdsA dihydrouridylates a short rRNA transcript mimicking the local structure of the peptidyl transferase site. This suggests an early introduction of this modification before ribosome assembly. Phylogenetic studies unveiled the widespread distribution of the yhiN gene in the bacterial kingdom, emphasizing the conservation of rRNA dihydrouridylation. In a broader context, these findings underscore nature's preference for utilizing reduced flavin in the reduction of uridines and their derivatives.

RNA modifying enzymes shape tRNA biogenesis and function

Transfer RNAs (tRNAs) are the most highly modified cellular RNAs, both with respect to the proportion of nucleotides that are modified within the tRNA sequence and with respect to the extraordinary diversity in tRNA modification chemistry. However, the functions of many different tRNA modifications are only beginning to emerge. tRNAs have two general clusters of modifications. The first cluster is within the anticodon stem-loop including several modifications essential for protein translation. The second cluster of modifications is within the tRNA elbow, and roles for these modifications are less clear. In general, tRNA elbow modifications are typically not essential for cell growth, but nonetheless several tRNA elbow modifications have been highly conserved throughout all domains of life. In addition to forming modifications, many tRNA modifying enzymes have been demonstrated or hypothesized to also play an important role in folding tRNA acting as tRNA chaperones. In this review, we summarize the known functions of tRNA modifying enzymes throughout the lifecycle of a tRNA molecule, from transcription to degradation. Thereby, we describe how tRNA modification and folding by tRNA modifying enzymes enhance tRNA maturation, tRNA aminoacylation, and tRNA function during protein synthesis, ultimately impacting cellular phenotypes and disease.

Functional redundancy in tRNA dihydrouridylation

Dihydrouridine (D) is a common modified base found predominantly in transfer RNA (tRNA). Despite its prevalence, the mechanisms underlying dihydrouridine biosynthesis, particularly in prokaryotes, have remained elusive. Here, we conducted a comprehensive investigation into D biosynthesis in Bacillus subtilis through a combination of genetic, biochemical, and epitranscriptomic approaches. Our findings reveal that B. subtilis relies on two FMN-dependent Dus-like flavoprotein homologs, namely DusB1 and DusB2, to introduce all D residues into its tRNAs. Notably, DusB1 exhibits multisite enzyme activity, enabling D formation at positions 17, 20, 20a and 47, while DusB2 specifically catalyzes D biosynthesis at positions 20 and 20a, showcasing a functional redundancy among modification enzymes. Extensive tRNA-wide D-mapping demonstrates that this functional redundancy impacts the majority of tRNAs, with DusB2 displaying a higher dihydrouridylation efficiency compared to DusB1. Interestingly, we found that BsDusB2 can function like a BsDusB1 when overexpressed in vivo and under increasing enzyme concentration in vitro. Furthermore, we establish the importance of the D modification for B. subtilis growth at suboptimal temperatures. Our study expands the understanding of D modifications in prokaryotes, highlighting the significance of functional redundancy in this process and its impact on bacterial growth and adaptation.

Sequence- and Structure-Specific tRNA Dihydrouridylation by hDUS2

The post-transcriptional reduction of uridine to dihydrouridine (D) by dihydrouridine synthase (DUS) enzymes is among the most ubiquitous transformations in RNA biology. D is found at multiple sites in tRNAs, and studies in yeast have proposed that each of the four eukaryotic DUS enzymes modifies a different site; however, the molecular basis for this exquisite selectivity is unknown, and human DUS enzymes have remained largely uncharacterized. Here we investigate the substrate specificity of human dihydrouridine synthase 2 (hDUS2) using mechanism-based cross-linking with 5-bromouridine (5-BrUrd)-modified oligonucleotide probes and in vitro dihydrouridylation assays. We find that hDUS2 exclusively modifies U20 across diverse tRNA substrates and identify a minimal GU sequence within the tRNA D loop that underlies selective substrate modification. Further, we use our mechanism-based platform to screen small molecule inhibitors of hDUS2, a potential anticancer target. Our work elucidates the principles of substrate modification by a conserved DUS and provides a general platform for studying RNA modifying enzymes with sequence-defined activity-based probes.

Diversity, Distribution, and Chromosomal Rearrangements of TRIP1 Repeat Sequences in Escherichia coli

The bacterial genome contains numerous repeated sequences that greatly affect its genomic plasticity. The Escherichia coli K-12 genome contains three copies of the TRIP1 repeat sequence (TRIP1a, TRIP1b, and TRIP1c). However, the diversity, distribution, and role of the TRIP1 repeat sequence in the E. coli genome are still unclear. In this study, after screening 6725 E. coli genomes, the TRIP1 repeat was found in the majority of E. coli strains (96%: 6454/6725). The copy number and direction of the TRIP1 repeat sequence varied in each genome. Overall, 2449 genomes (36%: 2449/6725) had three copies of TRIP1 (TRIP1a, TRIP1b, and TRIP1c), which is the same as E. coli K-12. Five types of TRIP1 repeats, including two new types (TRIP1d and TRIP1e), are identified in E. coli genomes, located in 4703, 3529, 5741, 1565, and 232 genomes, respectively. Each type of TRIP1 repeat is localized to a specific locus on the chromosome. TRIP1 repeats can cause intra-chromosomal rearrangements. A total of 156 rearrangement events were identified, of which 88% (137/156) were between TRIP1a and TRIP1c. These findings have important implications for future research on TRIP1 repeats.

A minimal sequence motif drives selective tRNA dihydrouridylation by hDUS2

The post-transcriptional reduction of uridine to dihydrouridine (D) by dihydrouridine synthase (DUS) enzymes is among the most ubiquitous transformations in RNA biology. D is found at multiple sites in tRNAs and studies in yeast have proposed that each of the four eukaryotic DUS enzymes modifies a different site, however the molecular basis for this exquisite selectivity is unknown and human DUS enzymes have remained largely uncharacterized. Here we investigate the substrate specificity of human dihydrouridine synthase 2 (hDUS2) using mechanism-based crosslinking with 5-bromouridine (5-BrUrd)-modified oligonucleotide probes and in vitro dihydrouridylation assays. We find that hDUS2 modifies U20 in the D loop of diverse tRNA substrates and identify a minimal GU motif within the tRNA tertiary fold required for directing its activity. Further, we use our mechanism-based platform to screen small molecule inhibitors of hDUS2, a potential anti-cancer target. Our work elucidates the principles of substrate modification by a conserved DUS and provides a general platform to studying RNA modifying enzymes with sequence-defined activity-based probes.

The diverse structural modes of tRNA binding and recognition

tRNAs are short noncoding RNAs responsible for decoding mRNA codon triplets, delivering correct amino acids to the ribosome, and mediating polypeptide chain formation. Due to their key roles during translation, tRNAs have a highly conserved shape and large sets of tRNAs are present in all living organisms. Regardless of sequence variability, all tRNAs fold into a relatively rigid three-dimensional L-shaped structure. The conserved tertiary organization of canonical tRNA arises through the formation of two orthogonal helices, consisting of the acceptor and anticodon domains. Both elements fold independently to stabilize the overall structure of tRNAs through intramolecular interactions between the D- and T-arm. During tRNA maturation, different modifying enzymes posttranscriptionally attach chemical groups to specific nucleotides, which not only affect translation elongation rates but also restrict local folding processes and confer local flexibility when required. The characteristic structural features of tRNAs are also employed by various maturation factors and modification enzymes to assure the selection, recognition, and positioning of specific sites within the substrate tRNAs. The cellular functional repertoire of tRNAs continues to extend well beyond their role in translation, partly, due to the expanding pool of tRNA-derived fragments. Here, we aim to summarize the most recent developments in the field to understand how three-dimensional structure affects the canonical and noncanonical functions of tRNA.

Self-attention enabled deep learning of dihydrouridine (D) modification on mRNAs unveiled a distinct sequence signature from tRNAs

Dihydrouridine (D) is a modified pyrimidine nucleotide universally found in viral, prokaryotic, and eukaryotic species. It serves as a metabolic modulator for various pathological conditions, and its elevated levels in tumors are associated with a series of cancers. Precise identification of D sites on RNA is vital for understanding its biological function. A number of computational approaches have been developed for predicting D sites on tRNAs; however, none have considered mRNAs. We present here DPred, the first computational tool for predicting D on mRNAs in yeast from the primary RNA sequences. Built on a local self-attention layer and a convolutional neural network (CNN) layer, the proposed deep learning model outperformed classic machine learning approaches (random forest, support vector machines, etc.) and achieved reasonable accuracy and reliability with areas under the curve of 0.9166 and 0.9027 in jackknife cross-validation and on an independent testing dataset, respectively. Importantly, we showed that distinct sequence signatures are associated with the D sites on mRNAs and tRNAs, implying potentially different formation mechanisms and putative divergent functionality of this modification on the two types of RNA. DPred is available as a user-friendly Web server.

Evolutionary Diversity of Dus2 Enzymes Reveals Novel Structural and Functional Features among Members of the RNA Dihydrouridine Synthases Family

Dihydrouridine (D) is an abundant modified base found in the tRNAs of most living organisms and was recently detected in eukaryotic mRNAs. This base confers significant conformational plasticity to RNA molecules. The dihydrouridine biosynthetic reaction is catalyzed by a large family of flavoenzymes, the dihydrouridine synthases (Dus). So far, only bacterial Dus enzymes and their complexes with tRNAs have been structurally characterized. Understanding the structure-function relationships of eukaryotic Dus proteins has been hampered by the paucity of structural data. Here, we combined extensive phylogenetic analysis with high-precision 3D molecular modeling of more than 30 Dus2 enzymes selected along the tree of life to determine the evolutionary molecular basis of D biosynthesis by these enzymes. Dus2 is the eukaryotic enzyme responsible for the synthesis of D20 in tRNAs and is involved in some human cancers and in the detoxification of β-amyloid peptides in Alzheimer's disease. In addition to the domains forming the canonical structure of all Dus, i.e., the catalytic TIM-barrel domain and the helical domain, both participating in RNA recognition in the bacterial Dus, a majority of Dus2 proteins harbor extensions at both ends. While these are mainly unstructured extensions on the N-terminal side, the C-terminal side extensions can adopt well-defined structures such as helices and beta-sheets or even form additional domains such as zinc finger domains. 3D models of Dus2/tRNA complexes were also generated. This study suggests that eukaryotic Dus2 proteins may have an advantage in tRNA recognition over their bacterial counterparts due to their modularity.

Dihydrouridine in the Transcriptome: New Life for This Ancient RNA Chemical Modification

Until recently, post-transcriptional modifications of RNA were largely restricted to noncoding RNA species. However, this belief seems to have quickly dissipated with the growing number of new modifications found in mRNA that were originally thought to be primarily tRNA-specific, such as dihydrouridine. Recently, transcriptomic profiling, metabolic labeling, and proteomics have identified unexpected dihydrouridylation of mRNAs, greatly expanding the catalog of novel mRNA modifications. These data also implicated dihydrouridylation in meiotic chromosome segregation, protein translation rates, and cell proliferation. Dihydrouridylation of tRNAs and mRNAs are introduced by flavin-dependent dihydrouridine synthases. In this review, we will briefly outline the current knowledge on the distribution of dihydrouridines in the transcriptome, their chemical labeling, and highlight structural and mechanistic aspects regarding the dihydrouridine synthases enzyme family. A special emphasis on important research directions to be addressed will also be discussed. This new entry of dihydrouridine into mRNA modifications has definitely added a new layer of information that controls protein synthesis.

The Dihydrouridine landscape from tRNA to mRNA: a perspective on synthesis, structural impact and function

The universal dihydrouridine (D) epitranscriptomic mark results from a reduction of uridine by the Dus family of NADPH-dependent reductases and is typically found within the eponym D-loop of tRNAs. Despite its apparent simplicity, D is structurally unique, with the potential to deeply affect the RNA backbone and many, if not all, RNA-connected processes. The first landscape of its occupancy within the tRNAome was reported 20 years ago. Its potential biological significance was highlighted by observations ranging from a strong bias in its ecological distribution to the predictive nature of Dus enzymes overexpression for worse cancer patient outcomes. The exquisite specificity of the Dus enzymes revealed by a structure-function analyses and accumulating clues that the D distribution may expand beyond tRNAs recently led to the development of new high-resolution mapping methods, including Rho-seq that established the presence of D within mRNAs and led to the demonstration of its critical physiological relevance.

Structural and functional characterization of TrmM in m6 A modification of bacterial tRNA

N⁶ -methyladenosine (m⁶ A), widely distributed in both coding and noncoding RNAs, regulates the epigenetic signals and RNA metabolism in eukaryotes. Although this posttranscriptional modification is frequently observed in messenger and ribosomal RNA, it is relatively rare in transfer RNA. In Escherichia coli, TrmM encoded by yfiC is the tRNA-specific N⁶ methyltransferase, which modifies the A37 residue of tRNA^Val (cmo⁵ UAC) using S-adenosyl-l-methionine as a methyl donor. However, the structure-function relationship of this enzyme is not completely understood. In this report, we determined two x-ray crystal structures of Mycoplasma capricolum TrmM with and without S-adenosyl-l-homocysteine, which is a reaction product. We also demonstrated the cellular and in vitro activities of this enzyme in the m⁶ A modification of tRNA and the requirement of a divalent metal ion for its function, which is unprecedented in other RNA N⁶ methyltransferases, including the E. coli TrmM. Our results reveal that the dimeric form of M. capricolum TrmM is important for efficient tRNA binding and catalysis, thereby offering insights into the distinct substrate specificity of the monomeric E. coli homolog.

Unique anticodon loop conformation with the flipped-out wobble nucleotide in the crystal structure of unbound tRNAVal

During protein synthesis on ribosome, tRNA recognizes its cognate codon of mRNA through base-pairing with the anticodon. The 5'-end nucleotide of the anticodon is capable of wobble base-pairing, offering a molecular basis for codon degeneracy. The wobble nucleotide is often targeted for post-transcriptional modification, which affects the specificity and fidelity of the decoding process. Flipping-out of a wobble nucleotide in the anticodon loop has been proposed to be necessary for modifying enzymes to access the target nucleotide, which has been captured in selective structures of protein-bound complexes. Meanwhile, all other structures of free or ribosome-bound tRNA display anticodon bases arranged in stacked conformation. We report the X-ray crystal structure of unbound tRNA^Val1 to a 2.04 Å resolution showing two different conformational states of wobble uridine in the anticodon loop, one stacked on the neighboring base and the other swiveled out toward solvent. In addition, the structure reveals a rare magnesium ion coordination to the nitrogen atom of a nucleobase, which has been sampled very rarely among known structures of nucleic acids.

De novo crystal structure determination of double stranded RNA binding domain using only the sulfur anomalous diffraction in SAD phasing

Single-wavelength anomalous dispersion (SAD)-phasing using sulfur as the unique anomalous scatterer is a powerful method to solve the phase problem in protein crystallography. However, it is not yet widely used by non-expert crystallographers. We report here the structure determination of the double stranded RNA binding domain of human dihydrouridine synthase using the sulfur-SAD method and highly redundant data collected at 1.8 ​Å ("off-edge"), at which the estimated overall anomalous signal was 1.08%. High multiplicity data were collected on a single crystal rotated along the ϕ or ω axis at different κ angles, with the primary beam intensity being attenuated from 50% to 95%, compared to data collection at 0.98 ​Å, to reduce radiation damage. SHELXD succeeded to locate 14 out 15 sulfur sites only using the data sets recorded with highest beam attenuation, which provided phases sufficient for structure solving. In an attempt to stimulate the use of sulfur-SAD phasing by a broader community of crystallographers, we describe our experimental strategy together with a compilation of previous successful cases, suggesting that sulfur-SAD phasing should be attempted for determining the de novo structure of any protein with average sulfur content diffracting better than 3 ​Å resolution.

Dihydrouridine synthesis in tRNAs is under reductive evolution in Mollicutes

Dihydrouridine (D) is a tRNA-modified base conserved throughout all kingdoms of life and assuming an important structural role. The conserved dihydrouridine synthases (Dus) carries out D-synthesis. DusA, DusB and DusC are bacterial members, and their substrate specificity has been determined in Escherichia coli. DusA synthesizes D20/D20a while DusB and DusC are responsible for the synthesis of D17 and D16, respectively. Here, we characterize the function of the unique dus gene encoding a DusB detected in Mollicutes, which are bacteria that evolved from a common Firmicute ancestor via massive genome reduction. Using in vitro activity tests as well as in vivo E. coli complementation assays with the enzyme from Mycoplasma capricolum (DusB_MCap), a model organism for the study of these parasitic bacteria, we show that, as expected for a DusB homolog, DusB_MCap modifies U17 to D17 but also synthetizes D20/D20a combining therefore both E. coli DusA and DusB activities. Hence, this is the first case of a Dus enzyme able to modify up to three different sites as well as the first example of a tRNA-modifying enzyme that can modify bases present on the two opposite sides of an RNA-loop structure. Comparative analysis of the distribution of DusB homologs in Firmicutes revealed the existence of three DusB subgroups namely DusB1, DusB2 and DusB3. The first two subgroups were likely present in the Firmicute ancestor, and Mollicutes have retained DusB1 and lost DusB2. Altogether, our results suggest that the multisite specificity of the M. capricolum DusB enzyme could be an ancestral property.

Posttranscriptional modifications at the 37th position in the anticodon stem-loop of tRNA: structural insights from MD simulations

Transfer RNA (tRNA) is the most diversely modified RNA. Although the strictly conserved purine position 37 in the anticodon stem-loop undergoes modifications that are phylogenetically distributed, we do not yet fully understand the roles of these modifications. Therefore, molecular dynamics simulations are used to provide molecular-level details for how such modifications impact the structure and function of tRNA. A focus is placed on three hypermodified base families that include the parent i6A, t6A, and yW modifications, as well as derivatives. Our data reveal that the hypermodifications exhibit significant conformational flexibility in tRNA, which can be modulated by additional chemical functionalization. Although the overall structure of the tRNA anticodon stem remains intact regardless of the modification considered, the anticodon loop must rearrange to accommodate the bulky, dynamic hypermodifications, which includes changes in the nucleotide glycosidic and backbone conformations, and enhanced or completely new nucleobase-nucleobase interactions compared to unmodified tRNA or tRNA containing smaller (m1G) modifications at the 37th position. Importantly, the extent of the changes in the anticodon loop is influenced by the addition of small functional groups to parent modifications, implying each substituent can further fine-tune tRNA structure. Although the dominant conformation of the ASL is achieved in different ways for each modification, the molecular features of all modified tRNA drive the ASL domain to adopt the functional open-loop conformation. Importantly, the impact of the hypermodifications is preserved in different sequence contexts. These findings highlight the likely role of regulating mRNA structure and translation.

Cyclodipeptide Synthases of the NYH Subfamily Recognize tRNA Using an α-Helix Enriched with Positive Residues

Cyclodipeptide synthases (CDPSs) perform nonribosomal protein synthesis using two aminoacyl-tRNA substrates to produce cyclodipeptides. At present, there are no structural details of the CDPS:tRNA interaction available. Using AlbC, a CDPS that produces cyclo(l-Phe-l-Phe), the interaction between AlbC and its Phe-tRNA substrate was investigated. Simulations of models of the AlbC:tRNA complex, proposed by rigid-body docking or homology modeling, demonstrated that interactions with residues of an AlbC α-helix, α4, significantly contribute to the free energy of binding of AlbC to tRNA. Individual residue contributions to the tRNA binding free energy of the discovered binding mode explain well the available biochemical data, and the results of in vivo assay experiments performed in this work and guided by simulations. In molecular dynamics simulations, the phenylalanyl group predominantly occupied the two positions observed in the experimental structure of AlbC in the dipeptide intermediate state, suggesting that tRNAs of the first and second substrates interact with AlbC in a similar manner. Overall, given the high degree of sequence and structural similarity among the members of the CDPS NYH protein subfamily, the mechanism of the protein:tRNA interaction is expected to be pertinent to a wide range of proteins interacting with tRNA.

Structural basis of the interaction between cyclodipeptide synthases and aminoacylated tRNA substrates

Cyclodipeptide synthases (CDPSs) catalyze the synthesis of various cyclodipeptides by using two aminoacyl-tRNA (aa-tRNA) substrates in a sequential mechanism. Here, we studied binding of phenylalanyl-tRNAPhe to the CDPS from Candidatus Glomeribacter gigasporarum (Cglo-CDPS) by gel filtration and electrophoretic mobility shift assay. We determined the crystal structure of the Cglo-CDPS:Phe-tRNAPhe complex to 5 Å resolution and further studied it in solution using small-angle X-ray scattering (SAXS). The data show that the major groove of the acceptor stem of the aa-tRNA interacts with the enzyme through the basic β2 and β7 strands of CDPSs belonging to the XYP subfamily. A bending of the CCA extremity enables the amino acid moiety to be positioned in the P1 pocket while the terminal A76 adenosine occupies the P2 pocket. Such a positioning indicates that the present structure illustrates the binding of the first aa-tRNA. In cells, CDPSs and the elongation factor EF-Tu share aminoacylated tRNAs as substrates. The present study shows that CDPSs and EF-Tu interact with opposite sides of tRNA. This may explain how CDPSs hijack aa-tRNAs from canonical ribosomal protein synthesis.

The Thermus thermophilus DEAD-box protein Hera is a general RNA binding protein and plays a key role in tRNA metabolism

RNA helicases catalyze the ATP-dependent destabilization of RNA duplexes. DEAD-box helicases share a helicase core that mediates ATP binding and hydrolysis, RNA binding and unwinding. Most members of this family contain domains flanking the core that can confer RNA substrate specificity and guide the helicase to a specific RNA. However, the in vivo RNA substrates of most helicases are currently not defined. The DEAD-box helicase Hera from Thermus thermophilus contains a helicase core, followed by a dimerization domain and an RNA binding domain that folds into an RNA recognition motif (RRM). The RRM mediates high affinity binding to an RNA hairpin, and an adjacent duplex is then unwound by the helicase core. Hera is a cold-shock protein, and has been suggested to act as an RNA chaperone under cold-shock conditions. Using crosslinking immunoprecipitation of Hera/RNA complexes and sequencing, we show that Hera binds to a large fraction of T. thermophilus RNAs under normal-growth and cold-shock conditions without a strong sequence preference, in agreement with a structure-specific recognition of RNAs and a general function in RNA metabolism. Under cold-shock conditions, Hera is recruited to RNAs with high propensities to form stable secondary structures. We show that selected RNAs identified, including a set of tRNAs, bind to Hera in vitro, and activate the Hera helicase core. Gene ontology analysis reveals an enrichment of genes related to translation, including mRNAs of ribosomal proteins, tRNAs, tRNA ligases, and tRNA-modifying enzymes, consistent with a key role of Hera in ribosome and tRNA metabolism.

iRNAD: a computational tool for identifying D modification sites in RNA sequence

MOTIVATION

Dihydrouridine (D) is a common RNA post-transcriptional modification found in eukaryotes, bacteria and a few archaea. The modification can promote the conformational flexibility of individual nucleotide bases. And its levels are increased in cancerous tissues. Therefore, it is necessary to detect D in RNA for further understanding its functional roles. Since wet-experimental techniques for the aim are time-consuming and laborious, it is urgent to develop computational models to identify D modification sites in RNA.

RESULTS

We constructed a predictor, called iRNAD, for identifying D modification sites in RNA sequence. In this predictor, the RNA samples derived from five species were encoded by nucleotide chemical property and nucleotide density. Support vector machine was utilized to perform the classification. The final model could produce the overall accuracy of 96.18% with the area under the receiver operating characteristic curve of 0.9839 in jackknife cross-validation test. Furthermore, we performed a series of validations from several aspects and demonstrated the robustness and reliability of the proposed model.

AVAILABILITY AND IMPLEMENTATION

A user-friendly web-server called iRNAD can be freely accessible at http://lin-group.cn/server/iRNAD, which will provide convenience and guide to users for further studying D modification.

Prospective serum metabolomic profiling of lethal prostate cancer

Impaired metabolism may play an important role in the pathogenesis of lethal prostate cancer, yet there is a paucity of evidence regarding the association. We conducted a large prospective serum metabolomic analysis of lethal prostate cancer in 523 cases and 523 matched controls nested within the Alpha-Tocopherol, Beta-Carotene Cancer Prevention (ATBC) Study. Median time from baseline fasting serum collection to prostate cancer death was 18 years (maximum 30 years). We identified 860 known biochemicals through an ultrahigh-performance LC-MS/MS platform. Conditional logistic regression models estimated odds ratios (OR) and 95% confidence intervals of risk associated with 1-standard deviation (s.d.) increases in log-metabolite signals. We identified 34 metabolites associated with lethal prostate cancer with a false discovery rate (FDR) < 0.15. Notably, higher serum thioproline, and thioproline combined with two other cysteine-related amino acids and redox metabolites, cystine and cysteine, were associated with reduced risk (1-s.d. OR = 0.75 and 0.71, respectively; p ≤ 8.2 × 10-5 ). By contrast, the dipeptide leucylglycine (OR = 1.36, p = 8.2 × 10-5 ), and three gamma-glutamyl amino acids (OR = 1.28-1.30, p ≤ 4.6 × 10-4 ) were associated with increased risk of lethal prostate cancer. Cases with metastatic disease at diagnosis (n = 179) showed elevated risk for several lipids, including especially the ketone body 3-hydroxybutyrate (BHBA), acyl carnitines, and dicarboxylic fatty acids (1.37 ≤ OR ≤ 1.49, FDR < 0.15). These findings provide a prospective metabolomic profile of lethal prostate cancer characterized by altered biochemicals in the redox, dipeptide, pyrimidine, and gamma-glutamyl amino acid pathways, whereas ketone bodies and fatty acids were associated specifically with metastatic disease.

Molecular basis for transfer RNA recognition by the double-stranded RNA-binding domain of human dihydrouridine synthase 2

Double stranded RNA-binding domain (dsRBD) is a ubiquitous domain specialized in the recognition of double-stranded RNAs (dsRNAs). Present in many proteins and enzymes involved in various functional roles of RNA metabolism, including RNA splicing, editing, and transport, dsRBD generally binds to RNAs that lack complex structures. However, this belief has recently been challenged by the discovery of a dsRBD serving as a major tRNA binding module for human dihydrouridine synthase 2 (hDus2), a flavoenzyme that catalyzes synthesis of dihydrouridine within the complex elbow structure of tRNA. We here unveil the molecular mechanism by which hDus2 dsRBD recognizes a tRNA ligand. By solving the crystal structure of this dsRBD in complex with a dsRNA together with extensive characterizations of its interaction with tRNA using mutagenesis, NMR and SAXS, we establish that while hDus2 dsRBD retains a conventional dsRNA recognition capability, the presence of an N-terminal extension appended to the canonical domain provides additional residues for binding tRNA in a structure-specific mode of action. Our results support that this extension represents a feature by which the dsRBD specializes in tRNA biology and more broadly highlight the importance of structural appendages to canonical domains in promoting the emergence of functional diversity.

Conformational Stability Adaptation of a Double-Stranded RNA-Binding Domain to Transfer RNA Ligand

The double-stranded RNA-binding domain (dsRBD) is a broadly distributed domain among RNA-maturing enzymes. Although this domain recognizes dsRNA's structures via a conserved canonical structure adopting an α₁-β₁β₂β₃-α₂ topology, several dsRBDs can accommodate discrete structural extensions expanding further their functional repertoire. How these structural elements engage cooperative communications with the canonical structure and how they contribute to the dsRBD's overall folding are poorly understood. Here, we addressed these issues using the dsRBD of human dihydrouridine synthase-2 (hDus2) (hDus2-dsRBD) as a model. This dsRBD harbors N- and C-terminal extensions, the former being directly involved in the recognition of tRNA substrate of hDus2. These extensions engage residues that form a long-range hydrophobic network (LHN) outside the RNA-binding interface. We show by coarse-grain Brownian dynamics that the Nt-extension and its residues F359 and Y364 rigidify the major folding nucleus of the canonical structure via an indirect effect. hDus2-dsRBD unfolds following a two-state cooperative model, whereas both F359A and Y364A mutants, designed to destabilize this LHN, unfold irreversibly. Structural and computational analyses show that these mutants are unstable due to an increase in the dynamics of the two extensions favoring solvent exposure of α2-helix and weakening the main folding nucleus rigidity. This LHN appears essential for maintaining a thermodynamic stability of the overall system and eventually a functional conformation for tRNA recognition. Altogether, our findings suggest that functional adaptability of extended dsRBDs is promoted by a cooperative hydrophobic coupling between the extensions acting as effectors and the folding nucleus of the canonical structure.

Structural studies on dihydrouridine synthase A (DusA) from Pseudomonas aeruginosa

Dihydrouridination is one of the abundant modifications in tRNA editing. The presence of dihydrouridine is attributed to tRNA stability desired for the efficient gene translation process. The conversion of uridine to dihydrouridine is catalyzed by flavine containing enzyme called dihydrouridine synthase (Dus). We report first ever information about DusA enzyme from Pseudomonas aeruginosa in form of structural and functional studies. The gene coding for DusA from P. aeruginosa (PADusA) was cloned, expressed and purified, using recombinant DNA technology methods. Thermal and chemical stability of PADusA was determined with respect to temperature and urea-induced equilibrium unfolding experiments, with monitoring the change of ellipticity at 200-260 nm by Circular Dichroism (CD) spectroscopy. Unfolding studies revealed that PADusA has acquired a stable tertiary structure fold with a T_m value of 46.2 °C and C_m of 2.7 M for urea. The enzyme contains 43% α-helices and 16% β-strands. The three dimensional structure of PADusA was modeled using insilico methods. In order to understand the mechanism of substrate recognition and catalysis, tRNA and puromycin were docked on PADusA structure and their binding was analyzed. The structural features suggested that PADusA may also form a novel target for structure based drug design of antimicrobial agents.

Transfer RNA Modification Enzymes from Thermophiles and Their Modified Nucleosides in tRNA

To date, numerous modified nucleosides in tRNA as well as tRNA modification enzymes have been identified not only in thermophiles but also in mesophiles. Because most modified nucleosides in tRNA from thermophiles are common to those in tRNA from mesophiles, they are considered to work essentially in steps of protein synthesis at high temperatures. At high temperatures, the structure of unmodified tRNA will be disrupted. Therefore, thermophiles must possess strategies to stabilize tRNA structures. To this end, several thermophile-specific modified nucleosides in tRNA have been identified. Other factors such as RNA-binding proteins and polyamines contribute to the stability of tRNA at high temperatures. Thermus thermophilus, which is an extreme-thermophilic eubacterium, can adapt its protein synthesis system in response to temperature changes via the network of modified nucleosides in tRNA and tRNA modification enzymes. Notably, tRNA modification enzymes from thermophiles are very stable. Therefore, they have been utilized for biochemical and structural studies. In the future, thermostable tRNA modification enzymes may be useful as biotechnology tools and may be utilized for medical science.

A rationale for tRNA modification circuits in the anticodon loop

The numerous post-transcriptional modifications of tRNA play a crucial role in tRNA function. While most modifications are introduced to tRNA independently, several sets of modifications are found to be interconnected such that the presence of one set of modifications drives the formation of another modification. The vast majority of these modification circuits are found in the anticodon loop (ACL) region where the largest variety and highest density of modifications occur compared to the other parts of the tRNA and where there is relatively limited sequence and structural information. We speculate here that the modification circuits in the ACL region arise to enhance enzyme modification specificity by direct or indirect use of the first modification in the circuit as an additional recognition element for the second modification. We also describe the five well-studied modification circuits in the ACL, and outline possible mechanisms by which they may act. The prevalence of these modification circuits in the ACL and the phylogenetic conservation of some of them suggest that a number of other modification circuits will be found in this region in different organisms.

Unveiling structural and functional divergences of bacterial tRNA dihydrouridine synthases: perspectives on the evolution scenario

Post-transcriptional base modifications are important to the maturation process of transfer RNAs (tRNAs). Certain modifications are abundant and present at several positions in tRNA as for example the dihydrouridine, a modified base found in the three domains of life. Even though the function of dihydrourine is not well understood, its high content in tRNAs from psychrophilic bacteria or cancer cells obviously emphasizes a central role in cell adaptation. The reduction of uridine to dihydrouridine is catalyzed by a large family of flavoenzymes named dihydrouridine synthases (Dus). Prokaryotes have three Dus (A, B and C) wherein DusB is considered as an ancestral protein from which the two others derived via gene duplications. Here, we unequivocally established the complete substrate specificities of the three Escherichia coli Dus and solved the crystal structure of DusB, enabling for the first time an exhaustive structural comparison between these bacterial flavoenzymes. Based on our results, we propose an evolutionary scenario explaining how substrate specificities has been diversified from a single structural fold.

tRNA Modifications: Impact on Structure and Thermal Adaptation

Transfer RNAs (tRNAs) are central players in translation, functioning as adapter molecules between the informational level of nucleic acids and the functional level of proteins. They show a highly conserved secondary and tertiary structure and the highest density of post-transcriptional modifications among all RNAs. These modifications concentrate in two hotspots-the anticodon loop and the tRNA core region, where the D- and T-loop interact with each other, stabilizing the overall structure of the molecule. These modifications can cause large rearrangements as well as local fine-tuning in the 3D structure of a tRNA. The highly conserved tRNA shape is crucial for the interaction with a variety of proteins and other RNA molecules, but also needs a certain flexibility for a correct interplay. In this context, it was shown that tRNA modifications are important for temperature adaptation in thermophilic as well as psychrophilic organisms, as they modulate rigidity and flexibility of the transcripts, respectively. Here, we give an overview on the impact of modifications on tRNA structure and their importance in thermal adaptation.

Structural effects of modified ribonucleotides and magnesium in transfer RNAs

Modified nucleotides are ubiquitous and important to tRNA structure and function. To understand their effect on tRNA conformation, we performed a series of molecular dynamics simulations on yeast tRNA^Phe and tRNA_init, Escherichia coli tRNA_init and HIV tRNA^Lys. Simulations were performed with the wild type modified nucleotides, using the recently developed CHARMM compatible force field parameter set for modified nucleotides (J. Comput. Chem.2016, 37, 896), or with the corresponding unmodified nucleotides, and in the presence or absence of Mg²⁺. Results showed a stabilizing effect associated with the presence of the modifications and Mg²⁺ for some important positions, such as modified guanosine in position 37 and dihydrouridines in 16/17 including both structural properties and base interactions. Some other modifications were also found to make subtle contributions to the structural properties of local domains. While we were not able to investigate the effect of adenosine 37 in tRNA_init and limitations were observed in the conformation of E. coli tRNA_init, the presence of the modified nucleotides and of Mg²⁺ better maintained the structural features and base interactions of the tRNA systems than in their absence indicating the utility of incorporating the modified nucleotides in simulations of tRNA and other RNAs.

From Prebiotics to Probiotics: The Evolution and Functions of tRNA Modifications

All nucleic acids in cells are subject to post-transcriptional chemical modifications. These are catalyzed by a myriad of enzymes with exquisite specificity and that utilize an often-exotic array of chemical substrates. In no molecule are modifications more prevalent than in transfer RNAs. In the present document, we will attempt to take a chemical rollercoaster ride from prebiotic times to the present, with nucleoside modifications as key players and tRNA as the centerpiece that drove the evolution of biological systems to where we are today. These ideas will be put forth while touching on several examples of tRNA modification enzymes and their modus operandi in cells. In passing, we submit that the choice of tRNA is not a whimsical one but rather highlights its critical function as an essential invention for the evolution of protein enzymes.

The tRNA Elbow in Structure, Recognition and Evolution

Prominent in the L-shaped three-dimensional structure of tRNAs is the "elbow" where their two orthogonal helical stacks meet. It has a conserved structure arising from the interaction of the terminal loops of the D- and T-stem-loops, and presents to solution a flat face of a tertiary base pair between the D- and T-loops. In addition to the ribosome, which interacts with the elbow in all three of its tRNA binding sites, several cellular RNAs and many proteins are known to recognize the elbow. At least three classes of non-coding RNAs, namely 23S rRNA, ribonuclease P, and the T-box riboswitches, recognize the tRNA elbow employing an identical structural motif consisting of two interdigitated T-loops. In contrast, structural solutions to tRNA-elbow recognition by proteins are varied. Some enzymes responsible for post-transcriptional tRNA modification even disrupt the elbow structure in order to access their substrate nucleotides. The evolutionary origin of the elbow is mysterious, but, because it does not explicitly participate in the flow of genetic information, it has been proposed to be a late innovation. Regardless, it is biologically essential. Even some viruses that hijack the cellular machinery using tRNA decoys have convergently evolved near-perfect mimics of the tRNA elbow.

Correlation between the stability of tRNA tertiary structure and the catalytic efficiency of a tRNA-modifying enzyme, archaeal tRNA-guanine transglycosylase

In many archaeal tRNAs, archaeosine is found at position 15. During archaeosine biosynthesis, archaeal tRNA-guanine transglycosylase (ArcTGT) first replaces the guanine base at position 15 with 7-cyano-7-deazaguanine (preQ0). In this study, we investigated whether modified nucleosides in tRNA substrates would affect ArcTGT incorporation of preQ0. We prepared a series of hypomodified tRNAs(Ser)(GGA) from Escherichia coli strains lacking each tRNA-modifying enzyme. Measurement of ArcTGT kinetic parameters with the various tRNAs(Ser)(GGA) as substrates showed that the Km decreased due to the lack of modified nucleosides. The tRNAs(Ser)(GGA) melting profiles resulted in experimental evidence showing that each modified nucleoside in tRNA(Ser)(GGA) enhanced tRNA stability. Furthermore, the ArcTGT K(m) strongly correlated with the melting temperature (T(m)), suggesting that the unstable tRNA containing fewer modified nucleosides served as a better ArcTGT substrate. These results show that preQ0 incorporation into tRNA by ArcTGT takes place early in the archaeal tRNA modification process.

An extended dsRBD is required for post-transcriptional modification in human tRNAs

In tRNA, dihydrouridine is a conserved modified base generated by the post-transcriptional reduction of uridine. Formation of dihydrouridine 20, located in the D-loop, is catalyzed by dihydrouridine synthase 2 (Dus2). Human Dus2 (HsDus2) expression is upregulated in lung cancers, offering a growth advantage throughout its ability to interact with components of the translation apparatus and inhibit apoptosis. Here, we report the crystal structure of the individual domains of HsDus2 and their functional characterization. HsDus2 is organized into three major modules. The N-terminal catalytic domain contains the flavin cofactor involved in the reduction of uridine. The second module is the conserved α-helical domain known as the tRNA binding domain in HsDus2 homologues. It is connected via a flexible linker to an unusual extended version of a dsRNA binding domain (dsRBD). Enzymatic assays and yeast complementation showed that the catalytic domain binds selectively NADPH but cannot reduce uridine in the absence of the dsRBD. While in Dus enzymes from bacteria, plants and fungi, tRNA binding is essentially achieved by the α-helical domain, we showed that in HsDus2 this function is carried out by the dsRBD. This is the first reported case of a tRNA-modifying enzyme carrying a dsRBD used to bind tRNAs.

From bacterial to human dihydrouridine synthase: automated structure determination

The reduction of uridine to dihydrouridine at specific positions in tRNA is catalysed by dihydrouridine synthase (Dus) enzymes. Increased expression of human dihydrouridine synthase 2 (hDus2) has been linked to pulmonary carcinogenesis, while its knockdown decreased cancer cell line viability, suggesting that it may serve as a valuable target for therapeutic intervention. Here, the X-ray crystal structure of a construct of hDus2 encompassing the catalytic and tRNA-recognition domains (residues 1-340) determined at 1.9 Å resolution is presented. It is shown that the structure can be determined automatically by phenix.mr_rosetta starting from a bacterial Dus enzyme with only 18% sequence identity and a significantly divergent structure. The overall fold of the human Dus2 is similar to that of bacterial enzymes, but has a larger recognition domain and a unique three-stranded antiparallel β-sheet insertion into the catalytic domain that packs next to the recognition domain, contributing to domain-domain interactions. The structure may inform the development of novel therapeutic approaches in the fight against lung cancer.