MD simulation studies to investigate iso-energetic conformational behaviour of modified nucleosides m(2)G and m(2) 2G present in tRNA

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop

The N6-methyladenosine modification is one of the most abundant post-transcriptional modifications in ribonucleic acid (RNA) molecules. Using molecular dynamics simulations and alchemical free-energy calculations, we studied the structural and energetic implications of incorporating this modification in an adenine mononucleotide and an RNA hairpin structure. At the mononucleotide level, we found that the syn configuration is more favorable than the anti configuration by 2.05 ± 0.15 kcal/mol. The unfavorable effect of methylation was due to the steric overlap between the methyl group and a nitrogen atom in the purine ring. We then probed the effect of methylation in an RNA hairpin structure containing an AUCG tetraloop, which is recognized by a "reader" protein (YTHDC1) to promote transcriptional silencing of long noncoding RNAs. While methylation had no significant conformational effect on the hairpin stem, the methylated tetraloop showed enhanced conformational flexibility compared to the unmethylated tetraloop. The increased flexibility was associated with the outward flipping of two bases (A6 and U7) which formed stacking interactions with each other and with the C8 and G9 bases in the tetraloop, leading to a conformation similar to that in the RNA/reader protein complex. Therefore, methylation-induced conformational flexibility likely facilitates RNA recognition by the reader protein.

Human TRMT1 catalyzes m2G or m22G formation on tRNAs in a substrate-dependent manner

TRMT1 is an N2-methylguanosine (m2G) and N2,N2-methylguanosine (m22G) methyltransferase that targets G26 of both cytoplasmic and mitochondrial tRNAs. In higher eukaryotes, most cytoplasmic tRNAs with G26 carry m22G26, although the majority of mitochondrial G26-containing tRNAs carry m2G26 or G26, suggesting differences in the mechanisms by which TRMT1 catalyzes modification of these tRNAs. Loss-of-function mutations of human TRMT1 result in neurological disorders and completely abrogate tRNA:m22G26 formation. However, the mechanism underlying the independent catalytic activity of human TRMT1 and identity of its specific substrate remain elusive, hindering a comprehensive understanding of the pathogenesis of neurological disorders caused by TRMT1 mutations. Here, we showed that human TRMT1 independently catalyzes formation of the tRNA:m2G26 or m22G26 modification in a substrate-dependent manner, which explains the distinct distribution of m2G26 and m22G26 on cytoplasmic and mitochondrial tRNAs. For human TRMT1-mediated tRNA:m22G26 formation, the semi-conserved C11:G24 serves as the determinant, and the U10:A25 or G10:C25 base pair is also required, while the size of the variable loop has no effect. We defined the requirements of this recognition mechanism as the "m22G26 criteria". We found that the m22G26 modification occurred in almost all the higher eukaryotic tRNAs conforming to these criteria, suggesting the "m22G26 criteria" are applicable to other higher eukaryotic tRNAs.

Epitranscriptome: Review of Top 25 Most-Studied RNA Modifications

The alphabet of building blocks for RNA molecules is much larger than the standard four nucleotides. The diversity is achieved by the post-transcriptional biochemical modification of these nucleotides into distinct chemical entities that are structurally and functionally different from their unmodified counterparts. Some of these modifications are constituent and critical for RNA functions, while others serve as dynamic markings to regulate the fate of specific RNA molecules. Together, these modifications form the epitranscriptome, an essential layer of cellular biochemistry. As of the time of writing this review, more than 300 distinct RNA modifications from all three life domains have been identified. However, only a few of the most well-established modifications are included in most reviews on this topic. To provide a complete overview of the current state of research on the epitranscriptome, we analyzed the extent of the available information for all known RNA modifications. We selected 25 modifications to describe in detail. Summarizing our findings, we describe the current status of research on most RNA modifications and identify further developments in this field.

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

RFhy-m2G: Identification of RNA N2-methylguanosine modification sites based on random forest and hybrid features

N2-methylguanosine is a post-transcriptional modification of RNA that is found in eukaryotes and archaea. The biological function of m2G modification discovered so far is to control and stabilize the three-dimensional structure of tRNA and the dynamic barrier of reverse transcription. To discover additional biological functions of m2G, it is necessary to develop time-saving and labor-saving calculation tools to identify m2G. In this paper, based on hybrid features and a random forest, a novel predictor, RFhy-m2G, was developed to identify the m2G modification sites for three species. The hybrid feature used by the predictor is used to fuse the three features of ENAC, PseDNC, and NPPS. These three features include primary sequence derivation properties, physicochemical properties, and position-specific properties. Since there are redundant features in hybrid features, MRMD2.0 is used for optimal feature selection. Through feature analysis, it is found that the optimal hybrid features obtained still contain three kinds of properties, and the hybrid features can more accurately identify m2G modification sites and improve prediction performance. Based on five-fold cross-validation and independent testing to evaluate the prediction model, the accuracies obtained were 0.9982 and 0.9417, respectively. The robustness of the predictor is demonstrated by comparisons with other predictors.

A Virtual Screening Platform Identifies Chloroethylagelastatin A as a Potential Ribosomal Inhibitor

Chloroethylagelastatin A (CEAA) is an analogue of agelastatin A (AA), a natural alkaloid derived from a marine sponge. It is under development for therapeutic use against brain tumors as it has excellent central nervous system (CNS) penetration and pre-clinical therapeutic activity against brain tumors. Recently, AA was shown to inhibit protein synthesis by binding to the ribosomal A-site. In this study, we developed a novel virtual screening platform to perform a comprehensive screening of various AA analogues showing that AA analogues with proven therapeutic activity including CEAA have significant ribosomal binding capacity whereas therapeutically inactive analogues show poor ribosomal binding and revealing structural fingerprint features essential for drug-ribosome interactions. In particular, CEAA was found to have greater ribosomal binding capacity than AA. Biological tests showed that CEAA binds the ribosome and contributes to protein synthesis inhibition. Our findings suggest that CEAA may possess ribosomal inhibitor activity and that our virtual screening platform may be a useful tool in discovery and development of novel ribosomal inhibitors.

Direct Sequencing of tRNA by 2D-HELS-AA MS Seq Reveals Its Different Isoforms and Dynamic Base Modifications

Post-transcriptional modifications are intrinsic to RNA structure and function. However, methods to sequence RNA typically require a cDNA intermediate and are either not able to sequence these modifications or are tailored to sequence one specific nucleotide modification only. Interestingly, some of these modifications occur with <100% frequency at their particular sites, and site-specific quantification of their stoichiometries is another challenge. Here, we report a direct method for sequencing tRNA^Phe without cDNA by integrating a two-dimensional hydrophobic RNA end-labeling strategy with an anchor-based algorithm in mass spectrometry-based sequencing (2D-HELS-AA MS Seq). The entire tRNA^Phe was sequenced and the identity, location, and stoichiometry of all eleven different RNA modifications was determined, five of which were not 100% modified, including a 2'-O-methylated G (Gm) in the wobble anticodon position as well as an N², N²-dimethylguanosine (m²₂G), a 7-methylguanosine (m⁷G), a 1-methyladenosine (m¹A), and a wybutosine (Y), suggesting numerous post-transcriptional regulations in tRNA. Two truncated isoforms at the 3'-CCA tail of the tRNA^Phe (75 nt with a 3'-CC tail (80% abundance) and 74 nt with a 3'-C tail (3% abundance)) were identified in addition to the full-length 3'-CCA-tailed tRNA^Phe (76 nt, 17% abundance). We discovered a new isoform with A-G transitions/editing at the 44 and 45 positions in the tRNA^Phe variable loop, and discuss possible mechanisms related to the emergence and functions of the isoforms with these base transitions or editing. Our method revealed new isoforms, base modifications, and RNA editing as well as their stoichiometries in the tRNA that cannot be determined by current cDNA-based methods, opening new opportunities in the field of epitranscriptomics.

An Emerging Role for isomiRs and the microRNA Epitranscriptome in Neovascularization

Therapeutic neovascularization can facilitate blood flow recovery in patients with ischemic cardiovascular disease, the leading cause of death worldwide. Neovascularization encompasses both angiogenesis, the sprouting of new capillaries from existing vessels, and arteriogenesis, the maturation of preexisting collateral arterioles into fully functional arteries. Both angiogenesis and arteriogenesis are highly multifactorial processes that require a multifactorial regulator to be stimulated simultaneously. MicroRNAs can regulate both angiogenesis and arteriogenesis due to their ability to modulate expression of many genes simultaneously. Recent studies have revealed that many microRNAs have variants with altered terminal sequences, known as isomiRs. Additionally, endogenous microRNAs have been identified that carry biochemically modified nucleotides, revealing a dynamic microRNA epitranscriptome. Both types of microRNA alterations were shown to be dynamically regulated in response to ischemia and are able to influence neovascularization by affecting the microRNA’s biogenesis, or even its silencing activity. Therefore, these novel regulatory layers influence microRNA functioning and could provide new opportunities to stimulate neovascularization. In this review we will highlight the formation and function of isomiRs and various forms of microRNA modifications, and discuss recent findings that demonstrate that both isomiRs and microRNA modifications directly affect neovascularization and vascular remodeling.

iRNA-m2G: Identifying N2-methylguanosine Sites Based on Sequence-Derived Information

RNA N²-methylguanosine (m2G) is one kind of posttranscriptional modification and plays crucial roles in the control and stabilization of tRNA. However, our knowledge about the biological functions of m2G is still limited. The key step of revealing its new function is to recognize the m2G sites in the transcriptome. Since there is no effective method for detecting m2G sites, it is desirable to develop new methods to identify m2G sites. In this study, a computational predictor called iRNA-m2G was proposed to identify m2G sites in eukaryotic transcriptomes. In iRNA-m2G, the RNA sequences were encoded by using nucleotide chemical property and accumulated nucleotide frequency. iRNA-m2G was not only validated by the rigorous jackknife test on the benchmark dataset but also examined by performing cross-species validations. In addition, iRNA-m2G was also tested on an independent dataset. It was found that the accuracies obtained by iRNA-m2G were all quite promising in these tests, indicating that the proposed method could become a powerful tool for identifying m2G sites. Finally, a user-friendly web server for iRNA-m2G is freely accessible at http://lin-group.cn/server/iRNA-m2G.php.

The emerging impact of tRNA modifications in the brain and nervous system

A remarkable number of neurodevelopmental disorders have been linked to defects in tRNA modifications. These discoveries place tRNA modifications in the spotlight as critical modulators of gene expression pathways that are required for proper organismal growth and development. Here, we discuss the emerging molecular and cellular functions of the diverse tRNA modifications linked to cognitive and neurological disorders. In particular, we describe how the structure and location of a tRNA modification influences tRNA folding, stability, and function. We then highlight how modifications in tRNA can impact multiple aspects of protein translation that are instrumental for maintaining proper cellular proteostasis. Importantly, we describe how perturbations in tRNA modification lead to a spectrum of deleterious biological outcomes that can disturb neurodevelopment and neurological function. Finally, we summarize the biological themes shared by the different tRNA modifications linked to cognitive disorders and offer insight into the future questions that remain to decipher the role of tRNA modifications. This article is part of a Special Issue entitled: mRNA modifications in gene expression control edited by Dr. Soller Matthias and Dr. Fray Rupert.

TRMT1-Catalyzed tRNA Modifications Are Required for Redox Homeostasis To Ensure Proper Cellular Proliferation and Oxidative Stress Survival

Mutations in the tRNA methyltransferase 1 (TRMT1) gene have been identified as the cause of certain forms of autosomal-recessive intellectual disability (ID). However, the molecular pathology underlying ID-associated TRMT1 mutations is unknown, since the biological role of the encoded TRMT1 protein remains to be determined. Here, we have elucidated the molecular targets and function of TRMT1 to uncover the cellular effects of ID-causing TRMT1 mutations. Using human cells that have been rendered deficient in TRMT1, we show that TRMT1 is responsible for catalyzing the dimethylguanosine (m2,2G) base modification in both nucleus- and mitochondrion-encoded tRNAs. TRMT1-deficient cells exhibit decreased proliferation rates, alterations in global protein synthesis, and perturbations in redox homeostasis, including increased endogenous ROS levels and hypersensitivity to oxidizing agents. Notably, ID-causing TRMT1 variants are unable to catalyze the formation of m2,2G due to defects in RNA binding and cannot rescue oxidative stress sensitivity. Our results uncover a biological role for TRMT1-catalyzed tRNA modification in redox metabolism and show that individuals with TRMT1-associated ID are likely to have major perturbations in cellular homeostasis due to the lack of m2,2G modifications.

Conformational Preferences of Modified Nucleoside 5-Taurinomethyluridine, τm(5)U Occur at 'wobble' 34th Position in the Anticodon Loop of tRNA

Conformational preferences of hypermodified nucleoside 5-taurinomethyluridine 5'-monophoshate 'p-τm(5)U' (-CH2-NH2(+)-CH2-CH2-SO3(-)) have been investigated using semi-empirical RM1 method. Automated geometry optimization using ab initio molecular orbital HF-SCF (6-31G**) and DFT (B3LYP/6-31G**) calculations have also been made to compare the salient features. The RM1 preferred most stable conformation of 'p-τm(5)U' has been stabilized by hydrogen bonding interactions between O(11a)…HN(8), O1P(34)…HN(8), and O1P(34)…HC(10). Another conformational study of 5-taurinomethyluridine side chain has also been performed in context of anticodon loop bases of E. coli tRNA(Leu). The atom O(11a) of τm(5)U(34) side chain interacts with adenosine (A35) as well as ribose-phosphate backbone which might provide structural stability to the anticodon loop. The glycosyl torsion angle of τm(5)U retains 'anti'-conformation. The solvent accessible surface area calculations revealed the role of τm(5)U in tRNA(Leu) anticodon loop. MD simulation results are found in agreement with RM1 preferred stable structure. The MEPs calculations of τm(5)U(34):G3 model show unique potential tunnels between the hydrogen bond donor and acceptor atoms as compared to τm(5)U(34):A3 model. Thus, these results could pave the way to understand the role of τm(5)U(34) to recognize UUG/UUA codons at atomic level in the mitochondrial disease, MELAS.

The influence of hypermodified nucleosides lysidine and t(6)A to recognize the AUA codon instead of AUG: a molecular dynamics simulation study

Hypermodified nucleosides lysidine (L) and N(6)-threonylcarbamoyladenosine (t(6)A) influence codon-anticodon interactions during the protein biosynthesis process. Lysidine prevents the misrecognition of the AUG codon as isoleucine and that of AUA as methionine. The structural significance of these modified bases has not been studied in detail at the atomic level. Hence, in the present study we performed multiple molecular dynamics (MD) simulations of anticodon stem loop (ASL) of tRNA(Ile) in the presence and absence of modified bases 'L' and 't(6)A' at the 34th and 37th positions respectively along with trinucleotide 'AUA' and 'AUG' codons. Hydrogen bonding interactions formed by the tautomeric form of lysidine may assist in reading the third base adenine of the 'AUA' codon, unlike the guanine of the 'AUG' codon. Such interactions might be useful to restrict codon specificity to recognize isoleucine tRNA instead of methionine tRNA. The t(6)A side chain interacts with the purine ring of the first codon nucleotide adenine, which might provide base stacking interactions and could be responsible for restricting extended codon-anticodon recognition. We found that ASL tRNA(Ile) in the absence of modifications at the 34th and 37th positions cannot establish proper hydrogen bonding interactions to recognize the isoleucine codon 'AUA' and subsequently disturbs the anticodon loop structure. The binding free energy calculations revealed that tRNA(Ile) ASL with modified nucleosides prefers the codon AUA over AUG. Thus, these findings might be useful to understand the role of modified bases L and t(6)A to recognize the AUA codon instead of AUG.

Comparative Structural Dynamics of tRNA(Phe) with Respect to Hinge Region Methylated Guanosine: A Computational Approach

Transfer RNAs (tRNAs) contain various uniquely modified nucleosides thought to be useful for maintaining the structural stability of tRNAs. However, their significance for upholding the tRNA structure has not been investigated in detail at the atomic level. In this study, molecular dynamic simulations have been performed to assess the effects of methylated nucleic acid bases, N (2)-methylguanosine (m(2)G) and N (2)-N (2)-dimethylguanosine (m 2 (2) G) at position 26, i.e., the hinge region of E. coli tRNA(Phe) on its structure and dynamics. The results revealed that tRNA(Phe) having unmodified guanosine in the hinge region (G26) shows structural rearrangement in the core of the molecule, resulting in lack of base stacking interactions, U-turn feature of the anticodon loop, and TΨC loop. We show that in the presence of the unmodified guanosine, the overall fold of tRNA(Phe) is essentially not the same as that of m(2)G26 and m 2 (2) G26 containing tRNA(Phe). This structural rearrangement arises due to intrinsic factors associated with the weak hydrogen-bonding patterns observed in the base triples of the tRNA(Phe) molecule. The m(2)G26 and m 2 (2) G26 containing tRNA(Phe) retain proper three-dimensional fold through tertiary interactions. Single-point energy and molecular electrostatics potential calculation studies confirmed the structural significance of tRNAs containing m(2)G26 and m 2 (2) G26 compared to tRNA with normal G26, showing that the mono-methylated (m(2)G26) and dimethylated (m 2 (2) G26) modifications are required to provide structural stability not only in the hinge region but also in the other parts of tRNA(Phe). Thus, the present study allows us to better understand the effects of modified nucleosides and ionic environment on tRNA folding.