RNA-modifying proteins as anticancer drug targets

Epitranscriptome in action: RNA modifications in the DNA damage response

Complex pathways involving the DNA damage response (DDR) contend with cell-intrinsic and -extrinsic sources of DNA damage. DDR mis-regulation results in genome instability that can contribute to aging and diseases including cancer and neurodegeneration. Recent studies have highlighted key roles for several RNA species in the DDR, including short RNAs and RNA/DNA hybrids (R-loops) at DNA break sites, all contributing to efficient DNA repair. RNAs can undergo more than 170 distinct chemical modifications. These RNA modifications have emerged as key orchestrators of the DDR. Here, we highlight the function of enzyme- and non-enzyme-induced RNA modifications in the DDR, with particular emphasis on m6A, m5C, and RNA editing. We also discuss stress-induced RNA damage, including RNA alkylation/oxidation, RNA-protein crosslinks, and UV-induced RNA damage. Uncovering molecular mechanisms that underpin the contribution of RNA modifications to DDR and genome stability will have direct application to disease and approaches for therapeutic intervention.

Reserpine, a novel N6-methyladenosine regulator, reverses Lenvatinib resistance in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is an aggressive malignancy and a growing global health problem. Reserpine (Res), a plant-derived hypertension drug, has been reported to possess anti-tumor efficacy. However, the role and function of Res in N6-methyladenosine (m6A) regulation and Lenvatinib (Len) resistance in HCC have not been clarified.

PURPOSE

To verify whether Res can be used as a natural small-molecule regulator of m6A to reverse Len resistance in HCC.

METHODS

Dot blotting, Western blotting and m6A quantification were used to compare and analyze the differential expression of m6A and its methyltransferase METTL3. Western blotting, Real-Time PCR (RT-PCR), cellular thermal shift assay (CETSA) and molecular docking were used to explore the mechanism of interaction between Res and m6A. The effects of Res on the biological characteristics of Lenvatinib-resistant HCC cells were investigated through CCK-8, clone formation, and Transwell assays. Cell line-derived xenograft (CDX) and patient-derived xenograft (PDX) mouse models were used to assess the ability of Res to reverse Len resistance in vivo. MeRIP m6A sequencing, PATHWAY analysis and Western blotting were used to analyze the downstream signaling pathways and genes involved in Res-mediated reversal of Len resistance.

RESULTS

Len resistance in HCC is related to the increased m6A level and the high expression of METTL3. Res affects the activity of METTL3 protein by binding to it, thereby downregulating the level of m6A. In vitro study showed that Res can sensitize HCC cells to the anti-tumor effects of Len treatment, including blocking proliferation, inhibiting migration, and inducing apoptosis. Len-resistant CDX and PDX models revealed that Res can reverse the resistant phenotype, with the tumor inhibition rates of 77.46 % and 62.1 %, respectively, when combined with Len treatment. Analysis of xenograft tissues showed that the combination of Res and Len down-regulates the m6A level, reduces proliferation biomarkers, and induces apoptosis, which is consistent with the in vitro data. Mechanistically, our preliminary results indicate that Res can up-regulate the SMAD3 level by down-regulating m6A in Len-resistant cells.

CONCLUSIONS

Reserpine, a small-molecule regulator of m6A, reverses Lenvatinib-resistant phenotypes, including proliferation, migration and anti-apoptosis, in vitro and in vivo by targeting SMAD3 and down-regulating the m6A level in HCC.

Writers, readers, and erasers RNA modifications and drug resistance in cancer

Drug resistance in cancer cells significantly diminishes treatment efficacy, leading to recurrence and metastasis. A critical factor contributing to this resistance is the epigenetic alteration of gene expression via RNA modifications, such as N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C), 7-methylguanosine (m7G), pseudouridine (Ψ), and adenosine-to-inosine (A-to-I) editing. These modifications are pivotal in regulating RNA splicing, translation, transport, degradation, and stability. Governed by "writers," "readers," and "erasers," RNA modifications impact numerous biological processes and cancer progression, including cell proliferation, stemness, autophagy, invasion, and apoptosis. Aberrant RNA modifications can lead to drug resistance and adverse outcomes in various cancers. Thus, targeting RNA modification regulators offers a promising strategy for overcoming drug resistance and enhancing treatment efficacy. This review consolidates recent research on the role of prevalent RNA modifications in cancer drug resistance, with a focus on m6A, m1A, m5C, m7G, Ψ, and A-to-I editing. Additionally, it examines the regulatory mechanisms of RNA modifications linked to drug resistance in cancer and underscores the existing limitations in this field.

Electrochemical detection of FTO with N3-kethoxal labeling and MazF cleavage

N6-Methyladenosine (m6A) is a prevalent modification in eukaryotic mRNAs and is linked to various human cancers. The fat mass and obesity-associated protein (FTO), a key m6A demethylase, is crucial in m6A regulation, affecting many biological processes and diseases. Detecting FTO is vital for clinical and research applications. Our study leverages the specific cleavage properties of the MazF endoribonuclease to design an electrochemical method with signal amplification guided by streptavidin-horseradish peroxidase (SA-HRP), intended for FTO detection. Initially, the compound N3-kethoxal is employed for its reversible tagging ability, selectively attaching to guanine (G) bases. Subsequently, dibenzocyclooctyne polyethylene glycol biotin (DBCO-PEG4-Biotin), is introduced through a reaction with N3-kethoxal. HRP is then employed to catalyze the redox system to enhance the current response further. A promising linear correlation between the peak current and the FTO concentration was observed within the range of 7.90 × 10-8 to 3.50 × 10-7 M, with a detection limit of 5.80 × 10-8 M. Moreover, this method assessed the FTO inhibitor FB23's inhibitory effect, revealing a final IC50 value of 54.73 nM. This result aligns with the IC50 value of 60 nM obtained through alternative methods and is very close to the values reported in the literature. The study provides reference value for research into obesity, diabetes, cancer, and other FTO-related diseases, as well as for the screening of potential therapeutic drugs.

Regulations of m6A and other RNA modifications and their roles in cancer

RNA modification is an essential component of the epitranscriptome, regulating RNA metabolism and cellular functions. Several types of RNA modifications have been identified to date; they include N6-methyladenosine (m6A), N1-methyladenosine (m1A), 5-methylcytosine (m5C), N7-methylguanosine (m7G), N6,2'-O-dimethyladenosine (m6Am), N4-acetylcytidine (ac4C), etc. RNA modifications, mediated by regulators including writers, erasers, and readers, are associated with carcinogenesis, tumor microenvironment, metabolic reprogramming, immunosuppression, immunotherapy, chemotherapy, etc. A novel perspective indicates that regulatory subunits and post-translational modifications (PTMs) are involved in the regulation of writer, eraser, and reader functions in mediating RNA modifications, tumorigenesis, and anticancer therapy. In this review, we summarize the advances made in the knowledge of different RNA modifications (especially m6A) and focus on RNA modification regulators with functions modulated by a series of factors in cancer, including regulatory subunits (proteins, noncoding RNA or peptides encoded by long noncoding RNA) and PTMs (acetylation, SUMOylation, lactylation, phosphorylation, etc.). We also delineate the relationship between RNA modification regulator functions and carcinogenesis or cancer progression. Additionally, inhibitors that target RNA modification regulators for anticancer therapy and their synergistic effect combined with immunotherapy or chemotherapy are discussed.

Dysregulation of tRNA methylation in cancer: Mechanisms and targeting therapeutic strategies

tRNA is the RNA type that undergoes the most modifications among known RNA, and in recent years, tRNA methylation has emerged as a crucial process in regulating gene translation. Dysregulation of tRNA abundance occurs in cancer cells, along with increased expression and activity of tRNA methyltransferases to raise the level of tRNA modification and stability. This leads to hijacking of translation and synthesis of multiple proteins associated with tumor proliferation, metastasis, invasion, autophagy, chemotherapy resistance, and metabolic reprogramming. In this review, we provide an overview of current research on tRNA methylation in cancer to clarify its involvement in human malignancies and establish a theoretical framework for future therapeutic interventions targeting tRNA methylation processes.

Emerging role of RNA modification and long noncoding RNA interaction in cancer

RNA modification, especially N6-methyladenosine, 5-methylcytosine, and N7-methylguanosine methylation, participates in the occurrence and progression of cancer through multiple pathways. The function and expression of these epigenetic regulators have gradually become a hot topic in cancer research. Mutation and regulation of noncoding RNA, especially lncRNA, play a major role in cancer. Generally, lncRNAs exert tumor-suppressive or oncogenic functions and its dysregulation can promote tumor occurrence and metastasis. In this review, we summarize N6-methyladenosine, 5-methylcytosine, and N7-methylguanosine modifications in lncRNAs. Furthermore, we discuss the relationship between epigenetic RNA modification and lncRNA interaction and cancer progression in various cancers. Therefore, this review gives a comprehensive understanding of the mechanisms by which RNA modification affects the progression of various cancers by regulating lncRNAs, which may shed new light on cancer research and provide new insights into cancer therapy.

Histone lactylation-boosted ALKBH3 potentiates tumor progression and diminished promyelocytic leukemia protein nuclear condensates by m1A demethylation of SP100A

Albeit N1-Methyladenosine (m1A) RNA modification represents an important regulator of RNA metabolism, the role of m1A modification in carcinogenesis remains enigmatic. Herein, we found that histone lactylation enhances ALKBH3 expression and simultaneously attenuates the formation of tumor-suppressive promyelocytic leukemia protein (PML) condensates by removing the m1A methylation of SP100A, promoting the malignant transformation of cancers. First, ALKBH3 is specifically upregulated in high-risk ocular melanoma due to excessive histone lactylation levels, referring to m1A hypomethylation status. Moreover, the multiomics analysis subsequently identified that SP100A, a core component for PML bodies, serves as a downstream candidate target for ALKBH3. Therapeutically, the silencing of ALKBH3 exhibits efficient therapeutic efficacy in melanoma both in vitro and in vivo, which could be reversed by the depletion of SP100A. Mechanistically, we found that YTHDF1 is responsible for recognition of the m1A methylated SP100A transcript, which increases its RNA stability and translational efficacy. Conclusively, we initially demonstrated that m1A modification is necessary for tumor suppressor gene expression, expanding the current understandings of dynamic m1A function during tumor progression. In addition, our results indicate that lactylation-driven ALKBH3 is essential for the formation of PML nuclear condensates, which bridges our knowledge of m1A modification, metabolic reprogramming, and phase-separation events.

NAT10-mediated ac4C-modified ANKZF1 promotes tumor progression and lymphangiogenesis in clear-cell renal cell carcinoma by attenuating YWHAE-driven cytoplasmic retention of YAP1

BACKGROUND

Lymphatic metastasis is one of the most common metastatic routes and indicates a poor prognosis in clear-cell renal cell carcinoma (ccRCC). N-acetyltransferase 10 (NAT10) is known to catalyze N4-acetylcytidine (ac4C) modification of mRNA and participate in many cellular processes. However, its role in the lymphangiogenic process of ccRCC has not been reported. This study aimed to elucidate the role of NAT10 in ccRCC lymphangiogenesis, providing valuable insights into potential therapeutic targets for intervention.

METHODS

ac4C modification and NAT10 expression levels in ccRCC were assessed using public databases and clinical samples. Functional investigations involved manipulating NAT10 expression in cellular and mouse models to study its role in ccRCC. Mechanistic insights were gained through a combination of RNA sequencing, mass spectrometry, co-immunoprecipitation, RNA immunoprecipitation, immunofluorescence, and site-specific mutation analyses.

RESULTS

We found that ac4C modification and NAT10 expression levels increased in ccRCC. NAT10 promoted tumor progression and lymphangiogenesis of ccRCC by enhancing the nuclear import of Yes1-associated transcriptional regulator (YAP1). Subsequently, we identified ankyrin repeat and zinc finger peptidyl tRNA hydrolase 1 (ANKZF1) as the functional target of NAT10, and its upregulation in ccRCC was caused by NAT10-mediated ac4C modification. Mechanistic analyses demonstrated that ANKZF1 interacted with tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein epsilon (YWHAE) to competitively inhibit cytoplasmic retention of YAP1, leading to transcriptional activation of pro-lymphangiogenic factors.

CONCLUSIONS

These results suggested a pro-cancer role of NAT10-mediated acetylation in ccRCC and identified the NAT10/ANKZF1/YAP1 axis as an under-reported pathway involving tumor progression and lymphangiogenesis in ccRCC.

The impact of tRNA modifications on translation in cancer: identifying novel therapeutic avenues

Recent advancements have illuminated the critical role of RNA modifications in post-transcriptional regulation, shaping the landscape of gene expression. This review explores how tRNA modifications emerge as critical players, fine-tuning functionalities that not only maintain the fidelity of protein synthesis but also dictate gene expression and translation profiles. Highlighting their dysregulation as a common denominator in various cancers, we systematically investigate the intersection of both cytosolic and mitochondrial tRNA modifications with cancer biology. These modifications impact key processes such as cell proliferation, tumorigenesis, migration, metastasis, bioenergetics and the modulation of the tumor immune microenvironment. The recurrence of altered tRNA modification patterns across different cancer types underscores their significance in cancer development, proposing them as potential biomarkers and as actionable targets to disrupt tumorigenic processes, offering new avenues for precision medicine in the battle against cancer.

Molecular basis of the reaction mechanism of the methyltransferase HENMT1

PIWI-interacting RNAs (piRNAs) are important for ensuring the integrity of the germline. 3'-terminal 2'-O-methylation is essential for piRNA maturation and to protect them from degradation. HENMT1 (HEN Methyltransferase 1) carries out the 2'-O-methylation, which is of key importance for piRNA stability and functionality. However, neither the structure nor the catalytic mechanism of mammalian HENMT1 have been studied. We have constructed a catalytic-competent HENMT1 complex using computational approaches, in which Mg2+ is primarily coordinated by four evolutionary conserved residues, and is further auxiliary coordinated by the 3'-O and 2'-O on the 3'-terminal nucleotide of the piRNA. Our study suggests that metal has limited effects on substrate and cofactor binding but is essential for catalysis. The reaction consists of deprotonation of the 2'-OH to 2'-O and a methyl transfer from SAM to the 2'-O. The methyl transfer is spontaneous and fast. Our in-depth analysis suggests that the 2'-OH may be deprotonated before entering the active site or it may be partially deprotonated at the active site by His800 and Asp859, which are in a special alignment that facilitates the proton transfer out of the active site. Furthermore, we have developed a detailed potential reaction scenario indicating that HENMT1 is Mg2+ utilizing but is not a Mg2+ dependent enzyme.

RMBase v3.0: decode the landscape, mechanisms and functions of RNA modifications

Although over 170 chemical modifications have been identified, their prevalence, mechanism and function remain largely unknown. To enable integrated analysis of diverse RNA modification profiles, we have developed RMBase v3.0 (http://bioinformaticsscience.cn/rmbase/), a comprehensive platform consisting of eight modules. These modules facilitate the exploration of transcriptome-wide landscape, biogenesis, interactome and functions of RNA modifications. By mining thousands of epitranscriptome datasets with novel pipelines, the 'RNA Modifications' module reveals the map of 73 RNA modifications of 62 species. the 'Genes' module allows to retrieve RNA modification profiles and clusters by gene and transcript. The 'Mechanisms' module explores 23 382 enzyme-catalyzed or snoRNA-guided modified sites to elucidate their biogenesis mechanisms. The 'Co-localization' module systematically formulates potential correlations between 14 histone modifications and 6 RNA modifications in various cell-lines. The 'RMP' module investigates the differential expression profiles of 146 RNA-modifying proteins (RMPs) in 18 types of cancers. The 'Interactome' integrates the interactional relationships between 73 RNA modifications with RBP binding events, miRNA targets and SNPs. The 'Motif' illuminates the enriched motifs for 11 types of RNA modifications identified from epitranscriptome datasets. The 'Tools' introduces a novel web-based 'modGeneTool' for annotating modifications. Overall, RMBase v3.0 provides various resources and tools for studying RNA modifications.

The role of the methyltransferase METTL3 in prostate cancer: a potential therapeutic target

The incidence of prostate cancer (PCa), the most prevalent malignancy, is currently at the forefront. RNA modification is a subfield of the booming field of epigenetics. To date, more than 170 types of RNA modifications have been described, and N6-methyladenosine (m6A) is the most abundant and well-characterized internal modification of mRNAs involved in various aspects of cancer progression. METTL3, the first identified key methyltransferase, regulates human mRNA and non-coding RNA expression in an m6A-dependent manner. This review elucidates the biological function and role of METTL3 in PCa and discusses the implications of METTL3 as a potential therapeutic target for future research directions and clinical applications.

The catalytic mechanism of the RNA methyltransferase METTL3

The complex of methyltransferase-like proteins 3 and 14 (METTL3-14) is the major enzyme that deposits N6-methyladenosine (m6A) modifications on mRNA in humans. METTL3-14 plays key roles in various biological processes through its methyltransferase (MTase) activity. However, little is known about its substrate recognition and methyl transfer mechanism from its cofactor and methyl donor S-adenosylmethionine (SAM). Here, we study the MTase mechanism of METTL3-14 by a combined experimental and multiscale simulation approach using bisubstrate analogues (BAs), conjugates of a SAM-like moiety connected to the N6-atom of adenosine. Molecular dynamics simulations based on crystal structures of METTL3-14 with BAs suggest that the Y406 side chain of METTL3 is involved in the recruitment of adenosine and release of m6A. A crystal structure with a bisubstrate analogue representing the transition state of methyl transfer shows a direct involvement of the METTL3 side chains E481 and K513 in adenosine binding which is supported by mutational analysis. Quantum mechanics/molecular mechanics (QM/MM) free energy calculations indicate that methyl transfer occurs without prior deprotonation of adenosine-N6. Furthermore, the QM/MM calculations provide further support for the role of electrostatic contributions of E481 and K513 to catalysis. The multidisciplinary approach used here sheds light on the (co)substrate binding mechanism, catalytic step, and (co)product release catalysed by METTL3, and suggests that the latter step is rate-limiting. The atomistic information on the substrate binding and methyl transfer reaction of METTL3 can be useful for understanding the mechanisms of other RNA MTases and for the design of transition state analogues as their inhibitors.

The role of RNA modification in urological cancers: mechanisms and clinical potential

RNA modification is a post-transcriptional level of regulation that is widely distributed in all types of RNAs, including mRNA, tRNA, rRNA, miRNA, and lncRNA, where N6-methyladenine (m6A) is the most abundant mRNA methylation modification. Significant evidence has depicted that m6A modifications are closely related to human diseases, especially cancer, and play pivotal roles in RNA transcription, splicing, stabilization, and translation processes. The most common urological cancers include prostate, bladder, kidney, and testicular cancers, accounting for a certain proportion of human cancers, with an ever-increasing incidence and mortality. The recurrence, systemic metastasis, poor prognosis, and drug resistance of urologic tumors have prompted the identification of new therapeutic targets and mechanisms. Research on m6A modifications may provide new solutions to the current puzzles. In this review, we provide a comprehensive overview of the key roles played by RNA modifications, especially m6A modifications, in urologic cancers, as well as recent research advances in diagnostics and molecularly targeted therapies.

Correlation between long non-coding RNA MAFG-AS1 and cancer prognosis: a meta-analysis

Background

MAF transcription factor G antisense RNA 1 (MAFG-AS1), a novel long non-coding RNA discovered recently, was proved to be useful in predicting malignancy prognosis. Nevertheless, its association with cancer prognosis has been inconsistent. Therefore, this meta-analysis aimed to explore the clinicopathological and prognostic significance of MAFG-AS1 in diverse carcinomas.

Methods

Studies focused on MAFG-AS1 expression as a prognostic role in cancers were thoroughly searched in six electronic databases. The value of MAFG-AS1 in malignancies was assessed by hazard ratios (HRs) or odds ratios (ORs). Additionally, the GEPIA database was utilized to further strengthen our conclusion.

Results

A total of 15 studies involving 1187 cases and nine types of cancers were recruited into this meta-analysis. High MAFG-AS1 expression was significantly related to advanced tumor stage (OR = 0.52, 95%CI [0.39, 0.69], P < 0.00001), earlier lymph node metastasis (OR = 3.62, 95%CI [2.19, 5.99], P < 0.00001), worse tumor differentiation (OR = 0.64, 95%CI [0.43, 0.95], P = 0.03), and poor overall survival (HR = 1.94, 95%CI [1.72, 2.19], P < 0.00001). No significant heterogeneity and publication bias was detected across studies. Meanwhile, MAFG-AS1 was significantly elevated in ten kinds of cancers based on the validation of the GEPIA database.

Conclusion

The results of this meta-analysis indicated that high MAFG-AS1 expression is dramatically correlated with unfavorable prognosis in cancers. MAFG-AS1 may be served as a promising biomarker for malignancies.

A novel m7G regulator-based methylation patterns in head and neck squamous cell carcinoma

Abnormal RNA N7-methylguanosine (m7G) modification is known to contribute to effects on tumor occurrence and development. Nevertheless, the mechanisms of its function in immunoregulation, tumor microenvironment (TME) modulation, and tumor promotion remain largely unknown. A series of computer-aided bioinformatic analyses were conducted based on transcriptomic, single-cell sequence, and spatial transcriptomic data to determine the m7G modification patterns in head and neck squamous cell carcinoma (HNSCC). Consensus clustering approach was employed according to the expressions of 33 m7G regulators. ESTIMATE, CIBERSORT, and single sample gene set enrichment analysis algorithms were adopted to investigate the immune cell infiltration features. A prognostic model named m7Gscore was established. Seurat, SingleR, and Monocle2 were used to analyze the single-cell sequence profiling. STUtility was used to integrate multiple spatial transcriptomic datasets. Quantitative reverse transcription polymerase chain reaction, transwell, and wound-healing assay were performed to verify the oncogenes. Here, three different m7G modification patterns were highlighted in HNSCC patients, which were also related to various clinical manifestations and three representative immunophenotypes: immune-excluded, immune-desert, and inflamed, separately. Patients with lower m7Gscore were highlighted by higher immune cell infiltrations, better overall survival rates, lesser tumor mutation burden (TMB), lower sensitivities to target inhibitors therapies, and better immunotherapeutic response. Moreover, DCPS, EIF4E, EIF4E2, LSM1, NCBP2, NUDT1, and NUDT5 were identified to play critical roles in T-cell differentiation. Knockdown of LSM1/NUDT5 could restrain the malignancy of HNSCC cells. Collectively, quantitative assessment of m7G modification patterns in individual HNSCC patients could contribute to identifying more efficient immunotherapeutic approaches and improve the clinical outcome of HNSCC.

Integrative Approach to Probe Alternative Redox Mechanisms in RNA Modifications

RNA modifications found in most RNAs, particularly in tRNAs and rRNAs, reveal an abundance of chemical alterations of nucleotides. Over 150 distinct RNA modifications are known, emphasizing a remarkable diversity of chemical moieties in RNA molecules. These modifications play pivotal roles in RNA maturation, structural integrity, and the fidelity and efficiency of translation processes. The catalysts responsible for these modifications are RNA-modifying enzymes that use a striking array of chemistries to directly influence the chemical landscape of RNA. This diversity is further underscored by instances where the same modification is introduced by distinct enzymes that use unique catalytic mechanisms and cofactors across different domains of life. This phenomenon of convergent evolution highlights the biological importance of RNA modification and the vast potential within the chemical repertoire for nucleotide alteration. While shared RNA modifications can hint at conserved enzymatic pathways, a major bottleneck is to identify alternative routes within species that possess a modified RNA but are devoid of known RNA-modifying enzymes. To address this challenge, a combination of bioinformatic and experimental strategies proves invaluable in pinpointing new genes responsible for RNA modifications. This integrative approach not only unveils new chemical insights but also serves as a wellspring of inspiration for biocatalytic applications and drug design. In this Account, we present how comparative genomics and genome mining, combined with biomimetic synthetic chemistry, biochemistry, and anaerobic crystallography, can be judiciously implemented to address unprecedented and alternative chemical mechanisms in the world of RNA modification. We illustrate these integrative methodologies through the study of tRNA and rRNA modifications, dihydrouridine, 5-methyluridine, queuosine, 8-methyladenosine, 5-carboxymethylamino-methyluridine, or 5-taurinomethyluridine, each dependent on a diverse array of redox chemistries, often involving organic compounds, organometallic complexes, and metal coenzymes. We explore how vast genome and tRNA databases empower comparative genomic analyses and enable the identification of novel genes that govern RNA modification. Subsequently, we describe how the isolation of a stable reaction intermediate can guide the synthesis of a biomimetic to unveil new enzymatic pathways. We then discuss the usefulness of a biochemical "shunt" strategy to study catalytic mechanisms and to directly visualize reactive intermediates bound within active sites. While we primarily focus on various RNA-modifying enzymes studied in our laboratory, with a particular emphasis on the discovery of a SAM-independent methylation mechanism, the strategies and rationale presented herein are broadly applicable for the identification of new enzymes and the elucidation of their intricate chemistries. This Account offers a comprehensive glimpse into the evolving landscape of RNA modification research and highlights the pivotal role of integrated approaches to identify novel enzymatic pathways.

m6A regulator-mediated methylation modification patterns and tumor microenvironment immune infiltration with prognostic analysis in esophageal cancer

Esophageal cancer is a highly malignant disease with poor prognosis. Despite recent advances in the study of esophageal cancer, there has been only limited improvement in the treatment and prognosis. N6-methyladenosine (m6A), a type of RNA modification, has been extensively investigated and is involved in many biological behaviors, including tumorigenesis and progression. Thus, more research on m6A modification may increase our understanding of esophageal cancer pathogenesis and provide potential targets. In our study, we integrated the public data of esophageal cancer from The Cancer Genome Atlas (TCGA) and Gene-Expression Omnibus (GEO) databases. Unsupervised clustering analysis was used to classify patients into different groups. Gene set variation analysis (GSVA) was performed in a nonparametric and unsupervised mode. We evaluated immune cell infiltration by single sample gene set enrichment analysis (ssGSEA). Differentially expressed genes (DEGs) among m6A clusters were identified using Empirical Bayesian approach. Both multivariate and univariate Cox regression models were used for prognostic analysis. We provided an overview of gene variation and expression of 23 m6A regulators in esophageal cancer, as well as their effects on survival. Based on the overall expression level of m6A regulators, patients were classified into three m6A clusters (A-C) with different immune cell infiltration abundance, gene expression signatures and prognosis. Among m6A clusters, we identified 206 DEGs, according to which patients were classified into 4 gene clusters (A-D). Quantitative m6A score was calculated for each patient based on those DEGs with significant impact on survival. The infiltration of all types of immune cells except type 2 T helper (Th2) cells were negatively correlated with m6A score. M6Acluster C exhibited the lowest m6A score, the most abundant immune cell infiltration, and the worst prognosis, suggesting an immune excluded phenotype. Consistently, gene cluster D with the lowest m6A score showed the worst prognosis. In short, patients with esophageal cancer showed different m6A modification patterns. Quantitative scoring indicated that patients with the lowest m6A score exhibited the most abundant immune cell infiltration and the poorest prognosis. This m6A scoring system is promising to assess m6A modification pattern, characterize immune infiltration and guide personalized treatment and prognostic prediction.

Integrative analysis of multi-omics data reveals a pseudouridine-related lncRNA signature for prediction of glioma prognosis and chemoradiotherapy sensitivity

BACKGROUND

Glioblastoma is the most common type of glioma with a high incidence and poor prognosis, and effective medical treatment remains challenging. Pseudouridine (Ψ) is the first post-transcriptional modification discovered and one of the most abundant modifications to RNA. However, the prognostic value of Ψ-related lncRNAs (ΨrLs) for glioma patients has never been systematically evaluated. This study aims to construct a risk model based on ΨrLs signature and to validate the predictive efficiency of the model.

METHOD

Transcriptomic data, genomic data, and relevant clinical data of glioma patients were extracted from the Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). ΨrLs with significant correlation with Ψ-related genes were identified, and univariate Cox regression, least absolute shrinkage and selection operator (LASSO) regression, and multivariate Cox regression were used to further select biomarkers and construct a ΨrLs signature risk model. Then, the expression of lncRNAs of ΨrLs signature in multiple glioma cell lines was detected by qPCR. Further, ROC analysis, stratification analysis, correlation analysis, survival analysis, nomogram, enrichment analysis, immune infiltration analysis, chemoradiotherapy sensitivity analysis, somatic mutation, and recurrent copy number variation (CNV) analysis were used to validate the predictive efficiency of ΨrLs signature in TCGA and CGGA datasets.

RESULTS

A four-lncRNA ΨrLs signature (DNAJC27-AS1, GDNF-AS1, ZBTB20-AS4, and DNMBP-AS1) risk model was constructed. By ROC analysis, stratified analysis, correlation analysis, survival analysis, and nomogram, the signature showed satisfactory predictive efficiency. Functional enrichment analysis revealed the differences in immune-related biological processes between high- and low-risk groups. Immune infiltration analysis showed that the high-risk group had lower tumor purity and higher stromal, immune and ESTIMATE scores. Mitoxantrone was identified as effective drug for low-risk group of glioma patients. Key genes in glioma development, including IDH1, EGFR, PTEN, etc., were differentially mutated between risk groups. The main recurrent CNVs in low-risk groups were 19q13.42 deletion and 7q34 amplification; 10q23.31 deletion and 12q14.1 in the high-risk group.

CONCLUSIONS

Our study identified a four-lncRNA ΨrLs signature that effectively predicts the prognosis of glioma patients and may serve as a diagnostic tool. Risk scores of glioma patients generated by the signature is associated with immune-related biological processes and chemoradiotherapy sensitivity. These findings may inform the development of more targeted and effective therapies for glioma patients.

m6A mRNA modification potentiates Th17 functions to inflame autoimmunity

N6-methyladenosine (m6A), the most common and abundant epigenetic RNA modification, governs mRNA metabolism to determine cell differentiation, proliferation and response to stimulation. m6A methyltransferase METTL3 has been reported to control T cell homeostasis and sustain the suppressive function of regulatory T cells (Tregs). However, the role of m6A methyltransferase in other subtypes of T cells remains unknown. T helper cells 17 (Th17) play a pivotal role in host defense and autoimmunity. Here, we found that the loss of METTL3 in T cells caused serious defect of Th17 cell differentiation, and impeded the development of experimental autoimmune encephalomyelitis (EAE). We generated Mettl3f/fIl17aCre mice and observed that METTL3 deficiency in Th17 cells significantly suppressed the development of EAE and displayed less Th17 cell infiltration into central nervous system (CNS). Importantly, we demonstrated that depletion of METTL3 attenuated IL-17A and CCR5 expression by facilitating SOCS3 mRNA stability in Th17 cells, leading to disrupted Th17 cell differentiation and infiltration, and eventually attenuating the process of EAE. Collectively, our results highlight that m6A modification sustains Th17 cell function, which provides new insights into the regulatory network of Th17 cells, and also implies a potential therapeutic target for Th17 cell mediated autoimmune disease.

N 6-methyladenosine Modification of FZR1 mRNA Promotes Gemcitabine Resistance in Pancreatic Cancer

The therapeutic options for treating pancreatic ductal adenocarcinoma (PDAC) are limited, and resistance to gemcitabine, a cornerstone of PDAC chemotherapy regimens, remains a major challenge. N6-methyladenosine (m6A) is a prevalent modification in mRNA that has been linked to diverse biological processes in human diseases. Herein, by characterizing the global m6A profile in a panel of gemcitabine-sensitive and gemcitabine-insensitive PDAC cells, we identified a key role for elevated m6A modification of the master G0-G1 regulator FZR1 in regulating gemcitabine sensitivity. Targeting FZR1 m6A modification augmented the response to gemcitabine treatment in gemcitabine-resistant PDAC cells both in vitro and in vivo. Mechanistically, GEMIN5 was identified as a novel m6A mediator that specifically bound to m6A-modified FZR1 and recruited the eIF3 translation initiation complex to accelerate FZR1 translation. FZR1 upregulation maintained the G0-G1 quiescent state and suppressed gemcitabine sensitivity in PDAC cells. Clinical analysis further demonstrated that both high levels of FZR1 m6A modification and FZR1 protein corresponded to poor response to gemcitabine. These findings reveal the critical function of m6A modification in regulating gemcitabine sensitivity in PDAC and identify the FZR1-GEMIN5 axis as a potential target to enhance gemcitabine response.

SIGNIFICANCE

Increased FZR1 translation induced by m6A modification engenders a gemcitabine-resistant phenotype by inducing a quiescent state and confers a targetable vulnerability to improve treatment response in PDAC.

RNA modification: mechanisms and therapeutic targets

RNA modifications are dynamic and reversible chemical modifications on substrate RNA that are regulated by specific modifying enzymes. They play important roles in the regulation of many biological processes in various diseases, such as the development of cancer and other diseases. With the help of advanced sequencing technologies, the role of RNA modifications has caught increasing attention in human diseases in scientific research. In this review, we briefly summarized the basic mechanisms of several common RNA modifications, including m6A, m5C, m1A, m7G, Ψ, A-to-I editing and ac4C. Importantly, we discussed their potential functions in human diseases, including cancer, neurological disorders, cardiovascular diseases, metabolic diseases, genetic and developmental diseases, as well as immune disorders. Through the "writing-erasing-reading" mechanisms, RNA modifications regulate the stability, translation, and localization of pivotal disease-related mRNAs to manipulate disease development. Moreover, we also highlighted in this review all currently available RNA-modifier-targeting small molecular inhibitors or activators, most of which are designed against m6A-related enzymes, such as METTL3, FTO and ALKBH5. This review provides clues for potential clinical therapy as well as future study directions in the RNA modification field. More in-depth studies on RNA modifications, their roles in human diseases and further development of their inhibitors or activators are needed for a thorough understanding of epitranscriptomics as well as diagnosis, treatment, and prognosis of human diseases.

Chinese Ecliptae herba (Eclipta prostrata (L.) L.) extract and its component wedelolactone enhances osteoblastogenesis of bone marrow mesenchymal stem cells via targeting METTL3-mediated m6A RNA methylation

ETHNOPHARMACOLOGICAL RELEVANCE

Chinese Ecliptae herba (Eclipta prostrata (L.) L.) is an ethnomedicinal herb, which is used mainly to nourish kidney and thus strengthen bones according to traditional Chinese medicine theory. Pharmacological studies have supported the ethnomedicine use, showing that Ecliptae herba extract has an anti-osteoporotic effect in vivo and promoted osteoblast proliferation and activity in vitro. However, the molecular mechanism of Ecliptae herba on osteoblast differentiation from bone marrow mesenchymal stem cells (BMSC), the progenitors of osteoblasts, is still unclear.

AIM OF THE STUDY

N6-methyladenosine (m6A) mRNA epigenetic modification may play a key role in promoting osteoblastic differentiation, and thus treating osteoporosis. This study sought to assess the mechanism through which Eclipate herba and its component wedelolactone influence m6A modification during the process of osteoblastogenesis from BMSC.

MATERIAL AND METHODS

The alkaline phosphatase (ALP) and Alizarin red S (ARS) staining were applied to determine osteoblastogenesis from BMSC. Western blot and quantitative real-time PCR were performed. RNA sequencing analysis was used to determine the characteristics of m6A methylation. Stable knocking down of METTL3 using lentiviral-based shRNA was performed.

RESULTS

Upon 9 d treatment of BMSC with ethyl acetate extract of Ecliptae herba (MHL), ALP activity and ossification level increased in comparison with osteogenic medium (OS)-treated control. The expression of methyltransferase METTL3 and METTL14 was significantly increased, but WTAP expression had no change in response to MHL treatment. Knocking down of METTL3 resulted in a decrease in MHL-induced ALP activity, ossification level as well as mRNA expression of Osterix and Osteocalcin, two bone formation-related markers. The level of m6A increased when BMSC was treated with MHL for 9 d. RNA sequencing analysis indicated that MHL treatment altered mRNA m6A modification of genes associated with osteoblastogenesis. By kyoto encyclopedia of genes and genomes (KEGG) pathway analysis, HIF-1α, PI3K/Akt, and Hippo signaling pathways were enriched and associated with m6A modification. The expression of m6A-modified genes including HIF-1α, VEGF-A, and RASSF1, was upregulated by MHL, but the upregulation was reversed after METTL3 knockdown. Additionally, the enhanced expression of METTL3 was also observed after treatment with wedelolactone, a component from MHL.

CONCLUSIONS

These results suggested a previously uncharacterized mechanism of MHL and wedelolactone on osteoblastogenesis, by which METTL3-mediated m6A methylation is involved and thus contributes to the enhancement of osteoblastogenesis.

Role of N6-methyladenosine RNA modification in gastric cancer

N6-methyladenosine (m6A) RNA methylation is the most prevalent internal modification of mammalian messenger RNA. The m6A modification affects multiple aspects of RNA metabolism, including processing, splicing, export, stability, and translation through the reversible regulation of methyltransferases (Writers), demethylases (Erasers), and recognition binding proteins (Readers). Accumulating evidence indicates that altered m6A levels are associated with a variety of human cancers. Recently, dysregulation of m6A methylation was shown to be involved in the occurrence and development of gastric cancer (GC) through various pathways. Thus, elucidating the relationship between m6A and the pathogenesis of GC has important clinical implications for the diagnosis, treatment, and prognosis of GC patients. In this review, we evaluate the potential role and clinical significance of m6A-related proteins which function in GC in an m6A-dependent manner. We discuss current issues regarding m6A-targeted inhibition of GC, explore new methods for GC diagnosis and prognosis, consider new targets for GC treatment, and provide a reasonable outlook for the future of GC research.

Advances in targeting RNA modifications for anticancer therapy

Numerous strategies are employed by cancer cells to control gene expression and facilitate tumorigenesis. In the study of epitranscriptomics, a diverse set of modifications to RNA represent a new player of gene regulation in disease and in development. N6-methyladenosine (m6A) is the most common modification on mammalian messenger RNA and tends to be aberrantly placed in cancer. Recognized by a series of reader proteins that dictate the fate of the RNA, m6A-modified RNA could promote tumorigenesis by driving protumor gene expression signatures and altering the immunologic response to tumors. Preclinical evidence suggests m6A writer, reader, and eraser proteins are attractive therapeutic targets. First-in-human studies are currently testing small molecule inhibition against the methyltransferase-like 3 (METTL3)/methyltransferase-like 14 (METTL14) methyltransferase complex. Additional modifications to RNA are adopted by cancers to drive tumor development and are under investigation.

Long noncoding RNA (lncRNA) H19: An essential developmental regulator with expanding roles in cancer, stem cell differentiation, and metabolic diseases

Recent advances in deep sequencing technologies have revealed that, while less than 2% of the human genome is transcribed into mRNA for protein synthesis, over 80% of the genome is transcribed, leading to the production of large amounts of noncoding RNAs (ncRNAs). It has been shown that ncRNAs, especially long non-coding RNAs (lncRNAs), may play crucial regulatory roles in gene expression. As one of the first isolated and reported lncRNAs, H19 has gained much attention due to its essential roles in regulating many physiological and/or pathological processes including embryogenesis, development, tumorigenesis, osteogenesis, and metabolism. Mechanistically, H19 mediates diverse regulatory functions by serving as competing endogenous RNAs (CeRNAs), Igf2/H19 imprinted tandem gene, modular scaffold, cooperating with H19 antisense, and acting directly with other mRNAs or lncRNAs. Here, we summarized the current understanding of H19 in embryogenesis and development, cancer development and progression, mesenchymal stem cell lineage-specific differentiation, and metabolic diseases. We discussed the potential regulatory mechanisms underlying H19's functions in those processes although more in-depth studies are warranted to delineate the exact molecular, cellular, epigenetic, and genomic regulatory mechanisms underlying the physiological and pathological roles of H19. Ultimately, these lines of investigation may lead to the development of novel therapeutics for human diseases by exploiting H19 functions.

Targeted Profiling of Epitranscriptomic Reader, Writer, and Eraser Proteins Regulated by H3K36me3

Trimethylation of lysine 36 on histone H3 (H3K36me3), an epigenetic mark associated with actively transcribed genes, plays an important role in multiple cellular processes, including transcription elongation, DNA methylation, DNA repair, etc. Aberrant expression and mutations of the main methyltransferase for H3K36me3, i.e., SET domain-containing 2 (SETD2), were shown to be associated with various cancers. Here, we performed targeted profiling of 154 epitranscriptomic reader, writer, and eraser (RWE) proteins using a scheduled liquid chromatography-parallel-reaction monitoring (LC-PRM) method coupled with the use of stable isotope-labeled (SIL) peptides as internal standards to investigate how H3K36me3 modulates the chromatin occupancies of epitranscriptomic RWE proteins. Our results showed consistent changes in chromatin occupancies of RWE proteins upon losses of H3K36me3 and H4K16ac and a role of H3K36me3 in recruiting METTL3 to chromatin following induction of DNA double-strand breaks. In addition, protein-protein interaction network and Kaplan-Meier survival analyses revealed the importance of METTL14 and TRMT11 in kidney cancer. Taken together, our work unveiled cross-talks between histone epigenetic marks (i.e., H3K36me3 and H4K16ac) and epitranscriptomic RWE proteins and uncovered the potential roles of these RWE proteins in H3K36me3-mediated biological processes.

The Regulation of m6A Modification in Glioblastoma: Functional Mechanisms and Therapeutic Approaches

Glioblastoma is the most prevalent primary brain tumour and invariably confers a poor prognosis. The immense intra-tumoral heterogeneity of glioblastoma and its ability to rapidly develop treatment resistance are key barriers to successful therapy. As such, there is an urgent need for the greater understanding of the tumour biology in order to guide the development of novel therapeutics in this field. N6-methyladenosine (m6A) is the most abundant of the RNA modifications in eukaryotes. Studies have demonstrated that the regulation of this RNA modification is altered in glioblastoma and may serve to regulate diverse mechanisms including glioma stem-cell self-renewal, tumorigenesis, invasion and treatment evasion. However, the precise mechanisms by which m6A modifications exert their functional effects are poorly understood. This review summarises the evidence for the disordered regulation of m6A in glioblastoma and discusses the downstream functional effects of m6A modification on RNA fate. The wide-ranging biological consequences of m6A modification raises the hope that novel cancer therapies can be targeted against this mechanism.

Over-expression of ADAR1 in mice does not initiate or accelerate cancer formation in vivo

Adenosine to inosine editing (A-to-I) in regions of double stranded RNA (dsRNA) is mediated by adenosine deaminase acting on RNA 1 (ADAR1) or ADAR2. ADAR1 and A-to-I editing levels are increased in many human cancers. Inhibition of ADAR1 has emerged as a high priority oncology target, however, whether ADAR1 overexpression enables cancer initiation or progression has not been directly tested. We established a series of in vivo models to allow overexpression of full-length ADAR1, or its individual isoforms, to test if increased ADAR1 expression was oncogenic. Widespread over-expression of ADAR1 or the p110 or p150 isoforms individually as sole lesions was well tolerated and did not result in cancer initiation. Therefore, ADAR1 overexpression alone is not sufficient to initiate cancer. We demonstrate that endogenous ADAR1 and A-to-I editing increased upon immortalization in murine cells, consistent with the observations from human cancers. We tested if ADAR1 over-expression could co-operate with cancer initiated by loss of tumour suppressors using a model of osteosarcoma. We did not see a disease potentiating or modifying effect of overexpressing ADAR1 or its isoforms in the models assessed. We conclude that increased ADAR1 expression and A-to-I editing in cancers is most likely a consequence of tumor formation.

RNA Folding Based on 5 Beads Model and Multiscale Simulation

RNA folding prediction is very meaningful and challenging. The molecular dynamics simulation (MDS) of all atoms (AA) is limited to the folding of small RNA molecules. At present, most of the practical models are coarse grained (CG) model, and the coarse-grained force field (CGFF) parameters usually depend on known RNA structures. However, the limitation of the CGFF is obvious that it is difficult to study the modified RNA. Based on the 3 beads model (AIMS_RNA_B3), we proposed the AIMS_RNA_B5 model with three beads representing a base and two beads representing the main chain (sugar group and phosphate group). We first run the all atom molecular dynamic simulation (AAMDS), and fit the CGFF parameter with the AA trajectory. Then perform the coarse-grained molecular dynamic simulation (CGMDS). AAMDS is the foundation of CGMDS. CGMDS is mainly to carry out the conformation sampling based on the current AAMDS state and improve the folding speed. We simulated the folding of three RNAs, which belong to hairpin, pseudoknot and tRNA respectively. Compared to the AIMS_RNA_B3 model, the AIMS_RNA_B5 model is more reasonable and performs better.

Construction of an RNA modification-related gene predictive model associated with prognosis and immunity in gastric cancer

BACKGROUND

Gastric cancer (GC) is one of the most common causes of cancer-related fatalities worldwide, and its progression is associated with RNA modifications. Here, using RNA modification-related genes (RNAMRGs), we aimed to construct a prognostic model for patients with GC.

METHODS

Based on RNAMRGs, RNA modification scores (RNAMSs) were obtained for GC samples from The Cancer Genome Atlas and were divided into high- and low-RNAMS groups. Differential analysis and weighted correlation network analysis were performed for the differential expressed genes (DEGs) to obtain the key genes. Next, univariate Cox regression, least absolute shrinkage and selection operator, and multivariate Cox regression analyses were performed to obtain the model. According to the model risk score, samples were divided into high- and low-risk groups. Enrichment analysis and immunoassays were performed for the DEGs in these groups. Four external datasets from Gene Expression Omnibus data base were used to test the accuracy of the predictive model.

RESULTS

We identified SELP and CST2 as key DEGs, which were used to generate the predictive model. The high-risk group had a worse prognosis compared to the low-risk group (p < 0.05). Enrichment analysis and immunoassays revealed that 144 DEGs related to immune cell infiltration were associated with the Wnt signaling pathway and included hub genes such as ELN. Overall mutation levels, tumor mutation burden, and microsatellite instability were lower, but tumor immune dysfunction and exclusion scores were greater (p < 0.05) in the high-risk group than in the low-risk group. The validation results showed that the prediction model score can accurately predict the prognosis of GC patients. Finally, a nomogram was constructed using the risk score combined with the clinicopathological characteristics of patients with GC.

CONCLUSION

This risk score from the prediction model related to the tumor microenvironment and immunotherapy could accurately predict the overall survival of GC patients.

The role of RNA methyltransferase METTL3 in gynecologic cancers: Results and mechanisms

N6-methyladenosine (m6A) methylation is the most prevalent mRNA modification in eukaryotes, and it is defined as the methylation of nitrogen atoms on the six adenine (A) bases of RNA in the presence of methyltransferases. Methyltransferase-like 3 (Mettl3), one of the components of m6A methyltransferase, plays a decisive catalytic role in m6A methylation. Recent studies have confirmed that m6A is associated with a wide spectrum of biological processes and it significantly affects disease progression and prognosis of patients with gynecologic tumors, in which the role of Mettl3 cannot be ignored. Mettl3 is involved in numerous pathophysiological functions, such as embryonic development, fat accumulation, and tumor progression. Moreover, Mettl3 may serve as a potential target for treating gynecologic malignancies, thus, it may benefit the patients and prolong survival. However, there is a need to further study the role and mechanism of Mettl3 in gynecologic malignancies. This paper reviews the recent progression on Mettl3 in gynecologic malignancies, hoping to provide a reference for further research.

Contemporary Progress and Opportunities in RNA-Targeted Drug Discovery

The surprising discovery that RNAs are the predominant gene products to emerge from the human genome catalyzed a renaissance in RNA biology. It is now well-understood that RNAs act as more than just a messenger and comprise a large and diverse family of ribonucleic acids of differing sizes, structures, and functions. RNAs play expansive roles in the cell, contributing to the regulation and fine-tuning of nearly all aspects of gene expression and genome architecture. In line with the significance of these functions, we have witnessed an explosion in discoveries connecting RNAs with a variety of human diseases. Consequently, the targeting of RNAs, and more broadly RNA biology, has emerged as an untapped area of drug discovery, making the search for RNA-targeted therapeutics of great interest. In this Microperspective, I highlight contemporary learnings in the field and present my views on how to catapult us toward the systematic discovery of RNA-targeted medicines.

Epitranscriptomics in metabolic disease

While epigenetic modifications of DNA and histones play main roles in gene transcription regulation, recently discovered post-transcriptional RNA modifications, known as epitranscriptomic modifications, have been found to have a profound impact on gene expression by regulating RNA stability, localization and decoding efficiency. Importantly, genetic variations or environmental perturbations of epitranscriptome modifiers (that is, writers, erasers and readers) are associated with obesity and metabolic diseases, such as type 2 diabetes. The epitranscriptome is closely coupled to epigenetic signalling, adding complexity to our understanding of gene expression in both health and disease. Moreover, the epitranscriptome in the parental generation can affect organismal phenotypes in the next generation. In this Review, we discuss the relationship between epitranscriptomic modifications and metabolic diseases, their relationship with the epigenome and possible therapeutic strategies.

A five-pseudouridylation-associated-LncRNA classifier for primary prostate cancer prognosis prediction

Background: Prostate cancer (PCa) is one of the most common cancers in males around the globe, and about one-third of patients with localized PCa will experience biochemical recurrence (BCR) after radical prostatectomy or radiation therapy. Reportedly, a proportion of patients with BCR had a poor prognosis. Cumulative studies have shown that RNA modifications participate in the cancer-related transcriptome, but the role of pseudouridylation occurring in lncRNAs in PCa remains opaque. Methods: Spearman correlation analysis and univariate Cox regression were utilized to determine pseudouridylation-related lncRNAs with prognostic value in PCa. Prognostic pseudouridylation-related lncRNAs were included in the LASSO (least absolute shrinkage and selection operator) regression algorithm to develop a predictive model. KM (Kaplan-Meier) survival analysis and ROC (receiver operating characteristic) curves were applied to validate the constructed model. A battery of biological cell assays was conducted to confirm the cancer-promoting effects of RP11-468E2.5 in the model. Results: A classifier containing five pseudouridine-related lncRNAs was developed to stratify PCa patients on BCR and named the "ψ-lnc score." KM survival analysis showed patients in the high ψ-lnc score group experienced BCR more than those in the low ψ-lnc score group. ROC curves demonstrated that ψ-lnc score outperformed other clinical indicators in BCR prediction. An external dataset, GSE54460, was utilized to validate the predictive model's efficacy and authenticity. A ceRNA (competitive endogenous RNA) network was constructed to explore the model's potential molecular functions and was annotated through GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses. RP11-468E2.5 was picked for further investigation, including pan-cancer analysis and experimental validation. Preliminarily, RP11-468E2.5 was confirmed as a tumor promoter. Conclusion: We provide some evidence that pseudouridylation in lncRNA played a role in the development of PCa and propose a novel prognostic classifier for clinical practice.

The roles of lncRNA functions and regulatory mechanisms in the diagnosis and treatment of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most frequent and deadly type of liver cancer. While the underlying molecular mechanisms are poorly understood, it is documented that lncRNAs may play key roles. Many HCC-associated lncRNAs have been linked to HBV and HCV infection, mediating gene expression, cell growth, development, and death. Studying the regulatory mechanisms and biological functions of HCC-related lncRNAs will assist our understanding of HCC pathogenesis as well as its diagnosis and management. Here, we address the potential of dysregulated lncRNAs in HCC as diagnostic and therapeutic biomarkers, and we evaluate the oncogenic or tumor-suppressive properties of these lncRNAs.

Dysfunctional tRNA reprogramming and codon-biased translation in cancer

Many cancers hijack translation to increase the synthesis of tumor-driving proteins, the messenger mRNAs of which have specific codon usage patterns. Termed 'codon-biased translation' and originally identified in stress response regulation, this mechanism is supported by diverse studies demonstrating how the 50 RNA modifications of the epitranscriptome, specific tRNAs, and codon-biased mRNAs are used by oncogenic programs to promote proliferation and chemoresistance. The epitranscriptome writers METTL1-WDR4, Elongator complex protein (ELP)1-6, CTU1-2, and ALKBH8-TRM112 illustrate the principal mechanism of codon-biased translation, with gene amplifications, increased RNA modifications, and enhanced tRNA stability promoting cancer proliferation. Furthermore, systems-level analyses of 34 tRNA writers and 493 tRNA genes highlight the theme of tRNA epitranscriptome dysregulation in many cancers and identify candidate tRNA writers, tRNA modifications, and tRNA molecules as drivers of pathological codon-biased translation.

Ribonucleic Acid Folding Prediction Based on Iterative Multiscale Simulation

RNA folding prediction is a challenge. Currently, many RNA folding models are coarse-grained (CG) with the potential derived from the known RNA structures. However, this potential is not suitable for modified and entirely new RNA. It is also not suitable for the folding simulation of RNA in the real cellular environment, including many kinds of molecular interactions. In contrast, our proposed model has the potential to address these issues, which is a multiscale simulation scheme based on all-atom (AA) force fields. We fit the CG force field using the trajectories generated by the AA force field and then iteratively perform molecular dynamics (MD) simulations of the two scales. The all-atom molecular dynamics (AAMD) simulation is mainly responsible for the correction of RNA structure, and the CGMD simulation is mainly responsible for efficient conformational sampling. On the basis of this scheme, we can successfully fold three RNAs belonging to a hairpin, a pseudoknot, and a four-way junction.

Rational Design and Optimization of m6A-RNA Demethylase FTO Inhibitors as Anticancer Agents

Aberrant regulation of N⁶-methyladenosine (m⁶A) RNA modification has been implicated in the progression of multiple diseases, including cancer. Previously, we identified a small molecule inhibitor of the m⁶A demethylase fat mass- and obesity-associated protein (FTO), which removes both m⁶A and N⁶,2′-O-dimethyladenosine (m⁶A_m) RNA modifications. In this work, we describe the rational design and optimization of a new class of FTO inhibitors derived from our previous lead FTO-04 with nanomolar potency and high selectivity against the homologous m⁶A RNA demethylase ALKBH5. The oxetanyl class of compounds comprise competitive inhibitors of FTO with potent antiproliferative effects in glioblastoma, acute myeloid leukemia, and gastric cancer models where lead FTO-43 demonstrated potency comparable to clinical chemotherapeutic 5-fluorouracil. Furthermore, FTO-43 increased m⁶A and m⁶A_m levels in a manner comparable to FTO knockdown in gastric cancer cells and regulated Wnt/PI3K-Akt signaling pathways. The oxetanyl class contains significantly improved anticancer agents with a variety of applications beyond glioblastoma.

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.

Epigenomic alterations in cancer: mechanisms and therapeutic potential

The human cell requires ways to specify its transcriptome without altering the essential sequence of DNA; this is achieved through mechanisms which govern the epigenetic state of DNA and epitranscriptomic state of RNA. These alterations can be found as modified histone proteins, cytosine DNA methylation, non-coding RNAs, and mRNA modifications, such as N6-methyladenosine (m6A). The different aspects of epigenomic and epitranscriptomic modifications require protein complexes to write, read, and erase these chemical alterations. Reflecting these important roles, many of these reader/writer/eraser proteins are either frequently mutated or differentially expressed in cancer. The disruption of epigenetic regulation in the cell can both contribute to cancer initiation and progression, and increase the likelihood of developing resistance to chemotherapies. Development of therapeutics to target proteins involved in epigenomic/epitranscriptomic modifications has been intensive, but further refinement is necessary to achieve ideal treatment outcomes without too many off-target effects for cancer patients. Therefore, further integration of clinical outcomes combined with large-scale genomic analyses is imperative for furthering understanding of epigenomic mechanisms in cancer.

Chemical biology and medicinal chemistry of RNA methyltransferases

RNA methyltransferases (MTases) are ubiquitous enzymes whose hitherto low profile in medicinal chemistry, contrasts with the surging interest in RNA methylation, the arguably most important aspect of the new field of epitranscriptomics. As MTases become validated as drug targets in all major fields of biomedicine, the development of small molecule compounds as tools and inhibitors is picking up considerable momentum, in academia as well as in biotech. Here we discuss the development of small molecules for two related aspects of chemical biology. Firstly, derivates of the ubiquitous cofactor S-adenosyl-l-methionine (SAM) are being developed as bioconjugation tools for targeted transfer of functional groups and labels to increasingly visible targets. Secondly, SAM-derived compounds are being investigated for their ability to act as inhibitors of RNA MTases. Drug development is moving from derivatives of cosubstrates towards higher generation compounds that may address allosteric sites in addition to the catalytic centre. Progress in assay development and screening techniques from medicinal chemistry have led to recent breakthroughs, e.g. in addressing human enzymes targeted for their role in cancer. Spurred by the current pandemic, new inhibitors against coronaviral MTases have emerged at a spectacular rate, including a repurposed drug which is now in clinical trial.

Bibliometric analysis of N6-methyladenosine in cancer: landscapes, highlights and trending topics

Background: N6-methyladenosine (m⁶A) is crucial in cancer and is being intensively studied. Aim: This bibliometric analysis seeks a broad picture of the role of m⁶A in cancer to guide and broaden future research. Methods: Publications were retrieved from Web of Science Core Collection and PubMed from 2000 to 2021, with keywords 'm⁶A' and 'cancer', and analyzed in biblioshiny and VOSviewer. Results: A total of 1013 documents were included, and China and the USA were the top countries with close collaboration. Mechanisms and predictive biomarkers of m⁶A regulator genes were highlighted. Cross-integration of m⁶A and other research hot spots, including 'immunotherapy', 'hypoxia' and 'polymorphism', were frontiers of m⁶A in cancer. Conclusion: This bibliometric study offered an updated perspective on m⁶A in cancer.

Targeted Profiling of Epitranscriptomic Reader, Writer, and Eraser Proteins Accompanied with Radioresistance in Breast Cancer Cells

Epitranscriptomic reader, writer, and eraser (RWE) proteins recognize, install, and remove modified nucleosides in RNA, which are known to play crucial roles in RNA processing, splicing, and stability. Here, we established a liquid chromatography-parallel-reaction monitoring (LC-PRM) method for high-throughput profiling of a total of 152 epitranscriptomic RWE proteins. We also applied the LC-PRM method, in conjunction with stable isotope labeling by amino acids in cell culture (SILAC), to quantify these proteins in two pairs of matched parental/radioresistant breast cancer cells (i.e., MDA-MB-231 and MCF-7 cells and their corresponding radioresistant C5 and C6 clones), with the goal of assessing the roles of these proteins in radioresistance. We found that eight epitranscriptomic RWE proteins were commonly altered by over 1.5-fold in the two pairs of breast cancer cells. Among them, TRMT1 (an m2,2G writer) may play a role in promoting breast cancer radioresistance due to its clinical relevance and its correlation with DNA repair gene sets. To our knowledge, this is the first report of a targeted proteomic method for comprehensive quantifications of epitranscriptomic RWE proteins. We envision that the LC-PRM method is applicable for studying the roles of these proteins in the metastatic transformation of cancer and therapeutic resistance of other types of cancer in the future.

MODOMICS: a database of RNA modification pathways. 2021 update

The MODOMICS database has been, since 2006, a manually curated and centralized resource, storing and distributing comprehensive information about modified ribonucleosides. Originally, it only contained data on the chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. Over the years, prompted by the accumulation of new knowledge and new types of data, it has been updated with new information and functionalities. In this new release, we have created a catalog of RNA modifications linked to human diseases, e.g., due to mutations in genes encoding modification enzymes. MODOMICS has been linked extensively to RCSB Protein Data Bank, and sequences of experimentally determined RNA structures with modified residues have been added. This expansion was accompanied by including nucleotide 5′-monophosphate residues. We redesigned the web interface and upgraded the database backend. In addition, a search engine for chemically similar modified residues has been included that can be queried by SMILES codes or by drawing chemical molecules. Finally, previously available datasets of modified residues, biosynthetic pathways, and RNA-modifying enzymes have been updated. Overall, we provide users with a new, enhanced, and restyled tool for research on RNA modification. MODOMICS is available at https://iimcb.genesilico.pl/modomics/.

Depletion of METTL3 alters cellular and extracellular levels of miRNAs containing m6A consensus sequences

Extracellular vesicles (EVs) are capable of transferring cargo from donor to recipient cells, but precisely how cargo content is regulated for export is mostly unknown. For miRNA cargo, we previously showed that when compared to isogenic colorectal cancer (CRC) cells expressing wild-type KRAS, a distinct subset of miRNAs are differentially enriched in EVs from KRAS mutant active CRC cells, with miR-100 being one of the most enriched. The mechanisms that could explain how miR-100 and other miRNAs are differentially exported into EVs have not been fully elucidated. Here, we tested the effect of N6-methyladenosine (m6A) modification on miRNA export into EVs by depletion of METTL3 and ALKBH5, a writer and eraser of m6A modification, respectively. While the effects of ALKBH5 knockdown were quite modest, decreased levels of METTL3 led to reduced cellular and extracellular levels of a subset of miRNAs that contain consensus sequences for m6A modification. Functional testing of EVs prepared from cells expressing shRNAs against METTL3 showed that they were less capable of conferring colony growth in 3D to wild-type KRAS cells and were also largely incapable of conferring the spread of cetuximab resistance. Our data support a role for METTL3 modification on cellular miRNA levels and export of specific miRNAs.