RNA Editing Research Articles

WTAP-mediated m6A modification of lncRNA Snhg1 improves myocardial ischemia-reperfusion injury via miR-361-5p/OPA1-dependent mitochondrial fusion

BackgroundMyocardial ischemia-reperfusion injury (MIRI) is caused by reperfusion after ischemic heart disease. LncRNA Snhg1 regulates the progression of various diseases. N6-methyladenosine (m6A) is the frequent RNA modification and plays a critical role in MIRI. However, it is unclear whether lncRNA Snhg1 regulates MIRI progression and whether the lncRNA Snhg1 was modified by m6A methylation.MethodsMouse cardiomyocytes HL-1 cells were utilized to construct the hypoxia/reoxygenation (H/R) injury model. HL-1 cell viability was evaluated utilizing CCK-8 method. Cell apoptosis, mitochondrial reactive oxygen species (ROS), and mitochondrial membrane potential (MMP) were quantitated utilizing flow cytometry. RNA immunoprecipitation and dual-luciferase reporter assays were applied to measure the m6A methylation and the interactions between lncRNA Snhg1 and targeted miRNA or target miRNAs and its target gene. The I/R mouse model was constructed with adenovirus expressing lncRNA Snhg1. HE and TUNEL staining were used to evaluate myocardial tissue damage and apoptosis.ResultsLncRNA Snhg1 was down-regulated after H/R injury, and overexpressed lncRNA Snhg1 suppressed H/R-stimulated cell apoptosis, mitochondrial ROS level and polarization. Besides, lncRNA Snhg1 could target miR-361-5p, and miR-361-5p targeted OPA1. Overexpressed lncRNA Snhg1 suppressed H/R-stimulated cell apoptosis, mitochondrial ROS level and polarization though the miR-361-5p/OPA1 axis. Furthermore, WTAP induced lncRNA Snhg1 m6A modification in H/R-stimulated HL-1 cells. Moreover, enforced lncRNA Snhg1 repressed I/R-stimulated myocardial tissue damage and apoptosis and regulated the miR-361-5p and OPA1 levels.ConclusionWTAP-mediated m6A modification of lncRNA Snhg1 regulated MIRI progression through modulating myocardial apoptosis, mitochondrial ROS production, and mitochondrial polarization via miR-361-5p/OPA1 axis, providing the evidence for lncRNA as the prospective target for alleviating MIRI progression.

ADR-2 regulates fertility and oocyte fate in C. elegans.

RNA binding proteins play essential roles in coordinating germline gene expression and development in all organisms. Here, we report that loss of ADR-2, a member of the Adenosine DeAminase acting on RNA (ADAR) family of RNA binding proteins and the sole adenosine-to-inosine RNA editing enzyme in C. elegans, can improve fertility in multiple genetic backgrounds. First, we show that loss of RNA editing by ADR-2 restores normal embryo production to subfertile animals that transgenically express a vitellogenin (yolk protein) fusion to green fluorescent protein. Using this phenotype, a high-throughput screen was designed to identify RNA binding proteins that when depleted yield synthetic phenotypes with loss of adr-2. The screen uncovered a genetic interaction between ADR-2 and SQD-1, a member of the heterogenous nuclear ribonucleoprotein (hnRNP) family of RNA binding proteins. Microscopy, reproductive assays, and high-throughput sequencing reveal that sqd-1 is essential for the onset of oogenesis and oogenic gene expression in young adult animals, and that loss of adr-2 can counteract the effects of loss of sqd-1 on gene expression and rescue the switch from spermatogenesis to oogenesis. Together, these data demonstrate that ADR-2 can contribute to the suppression of fertility and suggest novel roles for both RNA editing-dependent and independent mechanisms in regulating embryogenesis.

Comprehensive analysis of the complete mitochondrial genome of Lilium tsingtauense reveals a novel multichromosome structure.

Lilium tsingtauense mitogenome comprises 27 independent chromosome molecules, it undergoes frequent genomic recombination, and the rate of recombination and mutation between different repetitive sequences affects the formation of multichromosomal structures. Given the extremely large genome of Lily, which likely harbors additional genetic resources, it serves as an ideal material for studying the phylogenetic evolution of organisms. Although the Lilium chloroplast genome has been documented, the sequence of its mitochondrial genome (mitogenome) remains uncharted. Using BGI short reads and Nanopore long reads, we sequenced, assembled, and annotated the mitogenome of Lilium tsingtauense. This effort culminated in the characterization of Lilium's first complete mitogenome. Comparative analysis with other angiosperms revealed the unique multichromosomal structure of the L. tsingtauense mitogenome, spanning 1,125,108 bp and comprising 27 independent circular chromosomes. It contains 36 protein-coding genes, 12 tRNA genes, and 3 rRNA genes, with a GC content of 44.90%. Notably, three chromosomes in the L. tsingtauense mitogenome lack identifiable genes, hinting at the potential existence of novel genes and noncoding elements. The high degree of observed genome fragmentation implies frequent reorganization, with recombination and mutation rates among diverse repetitive sequences likely driving the formation of multichromosomal structures. Our comprehensive analysis, covering genome size, coding genes, structure, RNA editing, repetitive sequences, and sequence migration, sheds light on the evolutionary and molecular biology of multichromosomal mitochondria in Lilium. This high-quality mitogenome of L. tsingtauense not only enriches our understanding of multichromosomal mitogenomes but also establishes a solid foundation for future genome breeding and germplasm innovation in Lilium.

Dynamic changes in RNA m6A and 5 hmC influence gene expression programs during macrophage differentiation and polarisation

RNA modifications are essential for the establishment of cellular identity. Although increasing evidence indicates that RNA modifications regulate the innate immune response, their role in monocyte-to-macrophage differentiation and polarisation is unclear. While m6A has been widely studied, other RNA modifications, including 5 hmC, remain poorly characterised. We profiled m6A and 5 hmC epitranscriptomes, transcriptomes, translatomes and proteomes of monocytes and macrophages at rest and pro- and anti-inflammatory states. Transcriptome-wide mapping of m6A and 5 hmC reveals enrichment of m6A and/or 5 hmC on specific categories of transcripts essential for macrophage differentiation. Our analyses indicate that m6A and 5 hmC modifications are present in transcripts with critical functions in pro- and anti-inflammatory macrophages. Notably, we also discover the co-occurrence of m6A and 5 hmC on alternatively-spliced isoforms and/or opposing ends of the untranslated regions (UTR) of mRNAs with key roles in macrophage biology. In specific examples, RNA 5 hmC controls the decay of transcripts independently of m6A. This study provides (i) a comprehensive dataset to interrogate the role of RNA modifications in a plastic system (ii) a resource for exploring different layers of gene expression regulation in the context of human monocyte-to-macrophage differentiation and polarisation, (iii) new insights into RNA modifications as central regulators of effector cells in innate immunity.

Antiretroviral Therapy Suppresses RNA N6-Methyladenosine Modification in Peripheral Blood Mononuclear Cells from HIV-1-Infected Individuals.

RNA N6-methyladenosine (m6A) modification is important for regulating gene expression and innate immune responses to viral infection. HIV-1 in vitro infection induces a significant increase in m6A modification of cellular RNA; however, whether m6A levels of cellular RNA are affected by HIV-1 replication or by antiretroviral therapy (ART) in infected individuals remains unknown. Using dot blot or enzyme-linked immunosorbent assay, we measured RNA m6A levels of peripheral blood mononuclear cells (PBMCs) from healthy donors or HIV-1-infected individuals with or without ART. Using a reverse transcription-quantitative polymerase chain reaction array, we quantified expression levels of 84 type-I interferon (IFN-I)-responsive genes in PBMCs from some individuals of these three groups. RNA m6A levels in PBMCs from HIV-1 viremic patients (n = 10) were significantly higher (p ≤ .0001) compared with ART-treated individuals (n = 22) or 1.5-fold higher compared with healthy donors (n = 14). However, the increase in RNA m6A levels did not correlate with changes in the expression of 10 m6A-regulatory genes. We found significant upregulation and downregulation in the expression of several IFN-I-responsive genes from HIV-1 viremic patients (n = 4) and ART-treated patients (n = 6) compared with healthy donors (n = 5), respectively. Our results suggest that post-transcriptional m6A modification may contribute to the regulation of IFN-I-responsive gene expression during HIV-1 infection and ART.

The adenosine deaminase family acting on RNA 1 can be a useful diagnostic biomarker in chorioamnionitis

IntroductionChorioamnionitis (CAM) involves infection and inflammation of the chorion and amniotic membrane, but there are still no effective diagnostic biomarkers for CAM. MethodsWe investigated the correlation between RNA editing enzyme Adenosine deaminase family acting on RNA 1 (ADAR1) and CAM in chorion and amniotic membrane specimens derived from premature rupture of the membrane (PROM), CAM (pathologically diagnosed), and clinical CAM (clinically diagnosed) patients using reverse transcription polymerase chain reaction (RT-PCR). ResultsADAR1 was upregulated in the chorion and amniotic membrane specimens of CAM and clinical CAM patients (p < 0.001 and p = 0.005). ADAR1 had a significantly higher area under the curve (AUC) (0.735 and 0.828) than markers of inflammation characteristics in diagnosing CAM and clinical CAM patients. ADAR1 also had significantly higher AUC (0.701 and 0.837) than clinical characteristics for CAM and clinical CAM patients. DiscussionADAR1 can be a useful diagnostic biomarker in CAM patients.

IGF2BP2 inhibits invasion and migration of clear cell renal cell carcinoma via targeting Netrin-4 in an m6A-dependent manner.

Clear cell renal cell carcinoma (ccRCC), the most common subtype of renal cell carcinoma, often leads to a poor prognosis due to metastasis. The investigation of N6-methyladenosine (m6A) methylation, a crucial RNA modification, and its role in ccRCC, particularly through the m6A reader insulin-like growth factor 2 mRNA-binding protein 2(IGF2BP2), revealed significant insights. We found that IGF2BP2 was notably downregulated in ccRCC, which correlated with tumor aggressiveness and poor prognosis. Thus, IGFBP2 has emerged as an independent prognostic factor of ccRCC. Moreover, a strong positive correlation was observed between the expression of IGF2BP2 and Netrin-4. Netrin-4 was also downregulated in ccRCC, and its lower levels were associated with increased malignancy and poor prognosis. Overexpression of IGF2BP2 and Netrin-4 suppressed the invasion and migration of ccRCC cells, while Netrin-4 knockdown reversed these effects in ccRCC cell lines. RNA immunoprecipitation (RIP)-quantitative polymerase chain reaction validated the robust enrichment of Netrin-4 mRNA in anti-IGF2BP2 antibody immunoprecipitates. MeRlP showed significantly increased Netrin4 m6A levels after lGF2BP2 overexpression. Moreover, we found that IGF2BP2 recognized and bound to the m6A site within the coding sequence of Netrin-4, enhancing its mRNA stability. Collectively, these results showed that IGF2BP2 plays a suppressive role in the invasion and migration of ccRCC cells by targeting Netrin-4 in an m6A-dependent manner. These findings underscore the potential of IGF2BP2/Netrin-4 as a promising prognostic biomarker and therapeutic target in patients with ccRCC metastasis.

Morphine-Driven m6A Epitranscriptomic Neuroadaptations in Primary Cortical Cultures.

Opioid overdose is the leading cause of accidental death in the United States and remains a major public health concern, despite significant resources aimed at combating opioid misuse. Neurobiological research to elucidate molecular and cellular consequences of opioid exposure is required to define avenues to explore for reversal of opioid-induced neuroadaptations. Opioids impart well-documented regulation of the transcriptome and epigenetic modifications in the brain, but opioid-induced epitranscriptomic posttranscriptional regulation of RNA is vastly understudied. N6-methyladenosine (m6A) RNA methylation is significantly enriched in the brain and involved in learning, memory, and reward. m6A modifications have not been studied in opioid use disorder, despite being the most common RNA modification. We detected significant regulation of m6A-modifying enzymes in rat primary cortical cultures following morphine treatment, including AlkB Homolog 5 (Alkbh5). The m6a demethylaseALKBH5 functions as an m6A eraser, removing m6A modifications from mRNA. We hypothesized that chronic opioid treatment regulates m6A modifications through modulation of Alkbh5 and profiled m6A modifications in primary cortical cultures following chronic morphine treatment and Alkbh5 knock-down. We observed differential regulation of m6A modifications for a common set of transcripts following morphine or Alkbh5 knock-down, and the two treatments elicited concordant m6A epitranscriptomic profiles, suggesting that a subset of morphine-driven m6A modifications may be mediated through downregulation of Alkbh5 in cortical cultures. Gene Ontology terms of commonly regulated transcripts included serotonin secretion, synapse disassembly, neuron remodeling, and immune response. Thus, we conclude that morphine can drive epitranscriptomic changes, a subset of which may occur in an Alkbh5-dependent manner.

Hypoxia-induced downregulation of RNA m6A protein machinery in the naked mole-rat heart

Naked mole-rats, Heterocephalus glaber, are champion hypoxia-tolerant rodents that live under low oxygen conditions in their subterranean burrows. Detrimental effects of low oxygen can be mitigated through metabolic rate depression (MRD), metabolic reorganization, and global downregulation of nonessential cellular processes. Recent research has progressively implicated epigenetic modifications – rapid, reversible changes to gene expression that do not alter the DNA sequence itself – as major players in implementing and maintaining MRD. N6-adenosine (m6A) methylation is the most prevalent mammalian RNA modification and is responsible for pre-mRNA processing and mRNA export from the nucleus. Hence, m6A -mediated conformational changes alter the cellular fate of transcripts. The present study investigated the role of m6A RNA methylation responses to 24 h of hypoxia exposure in H. glaber cardiac tissue. Total protein levels of m6A writers/readers/erasers, m6A demethylase activity, and total m6A quantification were measured under normoxic vs. hypoxic conditions in H. glaber heart. While there was no change in either demethylase activity or total m6A content, many proteins of the m6A pathway were downregulated during hypoxia. Overall, m6A may not be a signature hypoxia-responsive characteristic in H. glaber heart, but downregulation of the protein machinery involved in m6A cycling points to an alternate biological involvement. Further research will explore other forms of RNA modifications and other epigenetic mechanisms to determine the controls on hypoxia endurance in this subterranean mammal.

IGF2BP2-modified circular RNA circCHD7 promotes endometrial cancer progression via stabilizing PDGFRB and activating JAK/STAT signaling pathway.

Circular RNAs (circRNAs) represent a class of covalently closed, single-stranded RNAs and have been linked to cancer progression. N6-methyladenosine (m6A) methylation is a ubiquitous RNA modification in cancer cells. Increasing evidence suggests that m6A can mediate the effects of circRNAs in cancer biology. In contrast, the post-transcriptional systems of m6A and circRNA in the progression of endometrial cancer (EC) remain obscure. The current study identified a novel circRNA with m6A modification, hsa_circ_0084582 (circCHD7), which was upregulated in EC tissues. Functionally, circCHD7 was found to promote the proliferation of EC cells. Mechanistically, circCHD7 interacted with insulin-like growth factor 2 mRNA-binding protein (IGF2BP2) to amplify its enrichment. Moreover, circCHD7 increased the mRNA stability of platelet-derived growth factor receptor beta (PDGFRB) in an m6A-dependent manner, thereby enhancing its expression. In addition, the circCHD7/IGF2BP2/PDGFRB axis activated the JAK/STAT signaling pathway and promoted EC cell proliferation. In conclusion, these findings provide new insights into the regulation of circRNA-mediated m6A modification, and the new "circCHD7-PDGFRB" model of regulation offers new perspectives on circCHD7 as a potential target for EC therapy.

Identification and validation of a novel risk model based on cuproptosis‑associated m6A for head and neck squamous cell carcinoma

BackgroundHead and neck squamous cell carcinoma (HNSCC) is a prevalent cancer with a poor survival rate due to anatomical limitations of the head and a lack of reliable biomarkers. Cuproptosis represents a novel cellular regulated death pathway, and N6-methyladenosine (m6A) is the most common internal RNA modification in mRNA. They are intricately connected to tumor formation, progression, and prognosis. This study aimed to construct a risk model for HNSCC using a set of mRNAs associated with m6A regulators and cuproptosis genes (mcrmRNA).MethodsRNA-seq and clinical data of HNSCC patients from The Cancer Genome Atlas (TCGA) database were analyzed to develop a risk model through the least absolute shrinkage and selection operator (LASSO) analysis. Survival analysis and receiver operating characteristic (ROC) analysis were performed for the high- and low-risk groups. Additionally, the model was validated using the GSE41613 dataset from the Gene Expression Omnibus (GEO) database. GSEA and CIBERSORT were applied to investigate the immune microenvironment of HNSCC.ResultsA risk model consisting of 32 mcrmRNA was developed using the LASSO analysis. The risk score of patients was confirmed to be an independent prognostic indicator by multivariate Cox analysis. The high-risk group exhibited a higher tumor mutation burden. Additionally, CIBERSORT analysis indicated varying levels of immune cell infiltration between the two groups. Significant disparities in drug sensitivity to common medications were also observed. Enrichment analysis further unveiled significant differences in metabolic pathways and RNA processing between the two groups.ConclusionsOur risk model can predict outcomes for HNSCC patients and offers valuable insights for personalized therapeutic approaches.

ADAR1 promotes cisplatin resistance in intrahepatic cholangiocarcinoma by regulating BRCA2 expression through A-to-I editing manner.

Aberrant A-to-I RNA editing, mediated by ADAR1 has been found to be associated with increased tumourigenesis and the development of chemotherapy resistance in various types of cancer. Intrahepatic cholangiocarcinoma (iCCA) is a highly aggressive malignancy with a poor prognosis, and overcoming chemotherapy resistance poses a significant clinical challenge. This study aimed to clarify the roles of ADAR1 in tumour resistance to cisplatin in iCCA. We discovered that ADAR1 expression is elevated in iCCA patients, particularly in those resistant to cisplatin, and associated with poor clinical outcomes. Downregulation of ADAR1 can increase the sensitivity of iCCA cells to cisplatin treatment, whereas its overexpression has the inverse effect. By integrating RNA sequencing and Sanger sequencing, we identified BRCA2, a critical DNA damage repair gene, as a downstream target of ADAR1 in iCCA. ADAR1 mediates the A-to-I editing in BRCA2 3'UTR, inhibiting miR-3157-5p binding, consequently increasing BRCA2 mRNA and protein levels. Furthermore, ADAR1 enhances cellular DNA damage repair ability and facilitates cisplatin resistance in iCCA cells. Combining ADAR1 targeting with cisplatin treatment markedly enhances the anticancer efficacy of cisplatin. In conclusion, ADAR1 promotes tumour progression and cisplatin resistance of iCCA. ADAR1 targeting could inform the development of innovative combination therapies for iCCA.

Discovery of Photoexcited 2-Chloro-3,5-Dinitrobenzoic Acid as a Chemical Deprenylase of i6A RNA.

RNA modifications play crucial roles in regulating gene expression and cellular homeostasis. Modulating RNA modifications, particularly by targeting the enzymes responsible for their catalysis, has emerged as a promising therapeutic strategy. However, limitations, such as the lack of identified modifying enzymes and compensatory mechanisms, hinder targeted interventions. Chemical approaches independent of enzymatic activity offer an alternative strategy for RNA modification modulation. Here, we present the identification of 2-chloro-3,5-dinitrobenzoic acid as a highly effective photochemical deprenylase of i6A RNA. This method demonstrates exceptional selectivity towards i6A, converting its substituent into a "N-doped" ozonide, which upon hydrolysis releases natural adenine. We believe that this chemical approach will pave the way for a better understanding of RNA modification biology and the development of novel therapeutic modalities.

Clonal Transgenerational Effects of Parental Grazing Environment on Offspring Shade Avoidance

Grassland plants that endure livestock grazing exhibit a dwarf phenotype, which can be transmitted to clonal offspring. Yet to date, it remains poorly understood whether such transgenerational dwarf effects alter the plants’ response to shade. Here, we conducted a common garden experiment under sunlight and shade conditions with clonal Leymus chinensis offspring, the parents of which had endured livestock overgrazing (OG) and non-grazing (NG) in the field, respectively. Plant morphological, physiological, and transcriptomic analyses were carried out. The results indicated that NG offspring showed greater shade avoidance than OG offspring. That is, NG offspring exhibited greater plasticity of vegetative height and leaf width, which may be contributed to their greater photosynthetic capacity and gibberellin (GA3) content compared with OG offspring when treated with shade. In addition, RNA-Seq profiling showed that differentially expressed genes in NG offspring were mainly enriched in RNA modification and metabolic processes, which facilitated rapid response to shade. Phytochrome interacting factors (PIFs) promoted downstream shade marker genes in NG offspring by significantly downregulating the expression of PHYC, SPY, and DELLA. Our findings suggest that light conditions should be taken into account to better understand transgenerational dwarf effects induced by livestock grazing on grassland ecosystems. These results provide new insights into the inducible factors of phenotypic variations in grassland plants that experience grazing.

Metabolic Recoding of NSUN2-Mediated m5C Modification Promotes the Progression of Colorectal Cancer via the NSUN2/YBX1/m5C-ENO1 Positive Feedback Loop.

The RNA modification, 5-methylcytosine (m5C), has recently gained prominence as a pivotal post-transcriptional regulator of gene expression, intricately intertwined with various tumorigenic processes. However,the exact mechanisms governing m5C modifications during the onset and progressionof colorectal cancer (CRC) remain unclear. Here, it is determined that the m5C methyltransferase NSUN2 exhibits significantly elevated expression and exerts an oncogenic function in CRC. Mechanistically, NSUN2 and YBX1 are identified as the "writer" and "reader" of ENO1, culminating in the reprogramming of the glucose metabolism and increased production of lactic acid in an m5C-dependent manner. The accumulation of lactic acid derived from CRC cells, in turn, activates the transcription of NSUN2 through histone H3K18 lactylation (H3K18la), and induces the lactylation of NSUN2 at the Lys356 residue (K356), which is crucial for capturing target RNAs. Together,these findings reveal an intriguing positive feedback loop involving the NSUN2/YBX1/m5C-ENO1 signaling axis, thereby bridging the connection between metabolic reprogramming and epigenetic remodeling, which may shed light on the therapeutic potential of combining an NSUN2 inhibitor with immunotherapy for CRC.

Mitoribosome structure with cofactors and modifications reveals mechanism of ligand binding and interactions with L1 stalk

The mitoribosome translates mitochondrial mRNAs and regulates energy conversion that is a signature of aerobic life forms. We present a 2.2 Å resolution structure of human mitoribosome together with validated mitoribosomal RNA (rRNA) modifications, including aminoacylated CP-tRNAVal. The structure shows how mitoribosomal proteins stabilise binding of mRNA and tRNA helping to align it in the decoding center, whereas the GDP-bound mS29 stabilizes intersubunit communication. Comparison between different states, with respect to tRNA position, allowed us to characterize a non-canonical L1 stalk, and molecular dynamics simulations revealed how it facilitates tRNA transitions in a way that does not require interactions with rRNA. We also report functionally important polyamines that are depleted when cells are subjected to an antibiotic treatment. The structural, biochemical, and computational data illuminate the principal functional components of the translation mechanism in mitochondria and provide a description of the structure and function of the human mitoribosome.

Unveiling the Role of RNA Recognition Motif Proteins in Orchestrating Nucleotide-Binding Site and Leucine-Rich Repeat Protein Gene Pairs and Chloroplast Immunity Pathways: Insights into Plant Defense Mechanisms.

In plants, nucleotide-binding site and leucine-rich repeat proteins (NLRs) play pivotal roles in effector-triggered immunity (ETI). However, the precise mechanisms underlying NLR-mediated disease resistance remain elusive. Previous studies have demonstrated that the NLR gene pair Pik-H4 confers resistance to rice blast disease by interacting with the transcription factor OsBIHD1, consequently leading to the upregulation of hormone pathways. In the present study, we identified an RNA recognition motif (RRM) protein, OsRRM2, which interacted with Pik1-H4 and Pik2-H4 in vesicles and chloroplasts. OsRRM2 exhibited a modest influence on Pik-H4-mediated rice blast resistance by upregulating resistance genes and genes associated with chloroplast immunity. Moreover, the RNA-binding sequence of OsRRM2 was elucidated using systematic evolution of ligands by exponential enrichment. Transcriptome analysis further indicated that OsRRM2 promoted RNA editing of the chloroplastic gene ndhB. Collectively, our findings uncovered a chloroplastic RRM protein that facilitated the translocation of the NLR gene pair and modulated chloroplast immunity, thereby bridging the gap between ETI and chloroplast immunity.

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design.

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

Queuosine biosynthetic enzyme, QueE moonlights as a cell division regulator.

In many organisms, stress responses to adverse environments can trigger secondary functions of certain proteins by altering protein levels, localization, activity, or interaction partners. Escherichia coli cells respond to the presence of specific cationic antimicrobial peptides by strongly activating the PhoQ/PhoP two-component signaling system, which regulates genes important for growth under this stress. As part of this pathway, a biosynthetic enzyme called QueE, which catalyzes a step in the formation of queuosine (Q) tRNA modification is upregulated. When cellular QueE levels are high, it co-localizes with the central cell division protein FtsZ at the septal site, blocking division and resulting in filamentous growth. Here we show that QueE affects cell size in a dose-dependent manner. Using alanine scanning mutagenesis of amino acids in the catalytic active site, we pinpoint residues in QueE that contribute distinctly to each of its functions-Q biosynthesis or regulation of cell division, establishing QueE as a moonlighting protein. We further show that QueE orthologs from enterobacteria like Salmonella typhimurium and Klebsiella pneumoniae also cause filamentation in these organisms, but the more distant counterparts from Pseudomonas aeruginosa and Bacillus subtilis lack this ability. By comparative analysis of E. coli QueE with distant orthologs, we elucidate a unique region in this protein that is responsible for QueE's secondary function as a cell division regulator. A dual-function protein like QueE is an exception to the conventional model of "one gene, one enzyme, one function", which has divergent roles across a range of fundamental cellular processes including RNA modification and translation to cell division and stress response.

Alteration of RNA m6A methylation mediates aberrant RNA binding protein expression and alternative splicing in condyloma acuminatum.

Condyloma acuminatum (CA) is caused by low-risk human papillomavirus, and is characterized by high recurrence after treatment. The RNA modification N6-methyladenosine (m6A) plays an important role during diverse viral infections, including high-risk HPV infection in cervical cancer. However, it is unclear whether low-risk HPV infection changes the RNA m6A methylation in CA. High-throughputm6A-sequencing was performed to profile the transcriptome-wide mRNA modifications of CA tissues infected by LR-HPVs and the paired normal tissues from CA patients. We further investigated the regulation of alternative splicing by RNA binding proteins (RBPs) with altered m6A modification and constructed a regulatory network among these RBPs, regulated alternative splicing events (RASEs) and regulated alternative splicing genes (RASGs) in CA. The results show that the m6A level in CA tissues differed from that in the paired controls. Furthermore, cell cycle- and cell adhesion- associated genes with m6A modification were differentially expressed in CA tissues compared to the paired controls. In particular, seven RNA binding protein genes with specific m6A methylated sites, showed a higher or lower expression at the mRNA level in CA tissues than in the paired normal tissues. In addition, these differentially expressed RNA binding protein genes would regulate the alternative splicing pattern of apoptotic process genes in CA tissue. Our study reveals a sophisticated m6A modification profile in CA tissue that affects the response of host cells to HPV infection, and provides cues for the further exploration of the roles of m6A and the development of a novel treatment strategy for CA.