m(1)A and m(1)G disrupt A-RNA structure through the intrinsic instability of Hoogsteen base pairs

Clinician’s Guide to Epitranscriptomics: An Example of N1-Methyladenosine (m1A) RNA Modification and Cancer

Epitranscriptomics is the study of modifications of RNA molecules by small molecular residues, such as the methyl (-CH3) group. These modifications are inheritable and reversible. A specific group of enzymes called “writers” introduces the change to the RNA; “erasers” delete it, while “readers” stimulate a downstream effect. Epitranscriptomic changes are present in every type of organism from single-celled ones to plants and animals and are a key to normal development as well as pathologic processes. Oncology is a fast-paced field, where a better understanding of tumor biology and (epi)genetics is necessary to provide new therapeutic targets and better clinical outcomes. Recently, changes to the epitranscriptome have been shown to be drivers of tumorigenesis, biomarkers, and means of predicting outcomes, as well as potential therapeutic targets. In this review, we aimed to give a concise overview of epitranscriptomics in the context of neoplastic disease with a focus on N1-methyladenosine (m1A) modification, in layman’s terms, to bring closer this omics to clinicians and their future clinical practice.

In vitro selection of N 1-methyladenosine-sensitive RNA-cleaving deoxyribozymes with 105-fold selectivity over unmethylated RNA

RNA-cleaving DNAzymes (RCDs) are catalytically active DNA molecules that cleave a wide range of RNA targets with extremely high sequence-selectivity, but none is able to faithfully discriminate methylated from unmethylated RNA (typically <30-fold). We report the first efforts to isolate RCDs from a random-sequence DNA pool by in vitro selection that cleave RNA/DNA chimera containing N 1-methyladenosine (m1A), one of the most prevalent RNA modifications that plays important regulatory roles in gene expression and human cancers. A cis-acting deoxyribozyme, RCD1-S2m1A, exhibits an observed rate constant (k obs) of 5.3 × 10-2 min-1, resulting in up to 105-fold faster cleavage of the m1A-modified versus unmethylated RNA. Furthermore, a trans-acting fluorogenic deoxyribozyme was constructed by labeling a fluorophore and a quencher at the 5' and 3' ends of the chimeric substrate, respectively. It permits the synchronization of RNA-cleaving with real-time fluorescence signaling, thus allowing the selective monitoring of ALKBH3-mediated demethylation and inhibitor screening in living cells.

Dual genetic level modification engineering accelerate genome evolution of Corynebacterium glutamicum

High spontaneous mutation rate is crucial for obtaining ideal phenotype and exploring the relationship between genes and phenotype. How to break the genetic stability of organisms and increase the mutation frequency has become a research hotspot. Here, we present a practical and controllable evolutionary tool (oMut-Cgts) based on dual genetic level modification engineering for Corynebacterium glutamicum. Firstly, the modification engineering of transcription and replication levels based on RNA polymerase α subunit and DNA helicase Cgl0854 as the 'dock' of cytidine deaminase (pmCDA1) significantly increased the mutation rate, proving that the localization of pmCDA1 around transient ssDNA is necessary for genome mutation. Then, the combined modification and optimization of engineering at dual genetic level achieved 1.02 × 104-fold increased mutation rate. The genome sequencing revealed that the oMut-Cgts perform uniform and efficient C:G→T:A transitions on a genome-wide scale. Furthermore, oMut-Cgts-mediated rapid evolution of C. glutamicum with stress (acid, oxidative and ethanol) tolerance proved that the tool has powerful functions in multi-dimensional biological engineering (rapid phenotype evolution, gene function mining and protein evolution). The strategies for rapid genome evolution provided in this study are expected to be applicable to a variety of applications in all prokaryotic cells.

The role of RNA methylation in tumor immunity and its potential in immunotherapy

RNA methylation, a prevalent post-transcriptional modification, has garnered considerable attention in research circles. It exerts regulatory control over diverse biological functions by modulating RNA splicing, translation, transport, and stability. Notably, studies have illuminated the substantial impact of RNA methylation on tumor immunity. The primary types of RNA methylation encompass N6-methyladenosine (m6A), 5-methylcytosine (m5C), N1-methyladenosine (m1A), and N7-methylguanosine (m7G), and 3-methylcytidine (m3C). Compelling evidence underscores the involvement of RNA methylation in regulating the tumor microenvironment (TME). By affecting RNA translation and stability through the "writers", "erasers" and "readers", RNA methylation exerts influence over the dysregulation of immune cells and immune factors. Consequently, RNA methylation plays a pivotal role in modulating tumor immunity and mediating various biological behaviors, encompassing proliferation, invasion, metastasis, etc. In this review, we discussed the mechanisms and functions of several RNA methylations, providing a comprehensive overview of their biological roles and underlying mechanisms within the tumor microenvironment and among immunocytes. By exploring how these RNA modifications mediate tumor immune evasion, we also examine their potential applications in immunotherapy. This review aims to provide novel insights and strategies for identifying novel targets in RNA methylation and advancing cancer immunotherapy efficacy.

Comprehensive review on lipid metabolism and RNA methylation: Biological mechanisms, perspectives and challenges

Adipose tissue plays a crucial role in maintaining energy balance, regulating hormones, and promoting metabolic health. To address disorders related to obesity and develop effective therapies, it is essential to have a deep understanding of adipose tissue biology. In recent years, RNA methylation has emerged as a significant epigenetic modification involved in various cellular functions and metabolic pathways. Particularly in the realm of adipogenesis and lipid metabolism, extensive research is ongoing to uncover the mechanisms and functional importance of RNA methylation. Increasing evidence suggests that RNA methylation plays a regulatory role in adipocyte development, metabolism, and lipid utilization across different organs. This comprehensive review aims to provide an overview of common RNA methylation modifications, their occurrences, and regulatory mechanisms, focusing specifically on their intricate connections to fat metabolism. Additionally, we discuss the research methodologies used in studying RNA methylation and highlight relevant databases that can aid researchers in this rapidly advancing field.

RNA: The Unsuspected Conductor in the Orchestra of Macromolecular Crowding

This comprehensive Review delves into the chemical principles governing RNA-mediated crowding events, commonly referred to as granules or biological condensates. We explore the pivotal role played by RNA sequence, structure, and chemical modifications in these processes, uncovering their correlation with crowding phenomena under physiological conditions. Additionally, we investigate instances where crowding deviates from its intended function, leading to pathological consequences. By deepening our understanding of the delicate balance that governs molecular crowding driven by RNA and its implications for cellular homeostasis, we aim to shed light on this intriguing area of research. Our exploration extends to the methodologies employed to decipher the composition and structural intricacies of RNA granules, offering a comprehensive overview of the techniques used to characterize them, including relevant computational approaches. Through two detailed examples highlighting the significance of noncoding RNAs, NEAT1 and XIST, in the formation of phase-separated assemblies and their influence on the cellular landscape, we emphasize their crucial role in cellular organization and function. By elucidating the chemical underpinnings of RNA-mediated molecular crowding, investigating the role of modifications, structures, and composition of RNA granules, and exploring both physiological and aberrant phase separation phenomena, this Review provides a multifaceted understanding of the intriguing world of RNA-mediated biological condensates.

New directions for Ψ and m1A decoding in mRNA: deciphering the stoichiometry and function

Over the past decade, advancements in epitranscriptomics have significantly enhanced our understanding of mRNA metabolism and its role in human development and diseases. This period has witnessed breakthroughs in sequencing technologies and the identification of key proteins involved in RNA modification processes. Alongside the well-studied m6A, Ψ and m1A have emerged as key epitranscriptomic markers. Initially identified through transcriptome-wide profiling, these modifications are now recognized for their broad impact on RNA metabolism and gene expression. In this Perspective, we focus on the detections and functions of Ψ and m1A modifications in mRNA and discuss previous discrepancies and future challenges. We summarize recent advances and highlight the latest sequencing technologies for stoichiometric detection and their mechanistic investigations for functional unveiling in mRNA as the new research directions.

RNA modifications in cellular metabolism: implications for metabolism-targeted therapy and immunotherapy

Cellular metabolism is an intricate network satisfying bioenergetic and biosynthesis requirements of cells. Relevant studies have been constantly making inroads in our understanding of pathophysiology, and inspiring development of therapeutics. As a crucial component of epigenetics at post-transcription level, RNA modification significantly determines RNA fates, further affecting various biological processes and cellular phenotypes. To be noted, immunometabolism defines the metabolic alterations occur on immune cells in different stages and immunological contexts. In this review, we characterize the distribution features, modifying mechanisms and biological functions of 8 RNA modifications, including N6-methyladenosine (m6A), N6,2'-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), 5-methylcytosine (m5C), N4-acetylcytosine (ac4C), N7-methylguanosine (m7G), Pseudouridine (Ψ), adenosine-to-inosine (A-to-I) editing, which are relatively the most studied types. Then regulatory roles of these RNA modification on metabolism in diverse health and disease contexts are comprehensively described, categorized as glucose, lipid, amino acid, and mitochondrial metabolism. And we highlight the regulation of RNA modifications on immunometabolism, further influencing immune responses. Above all, we provide a thorough discussion about clinical implications of RNA modification in metabolism-targeted therapy and immunotherapy, progression of RNA modification-targeted agents, and its potential in RNA-targeted therapeutics. Eventually, we give legitimate perspectives for future researches in this field from methodological requirements, mechanistic insights, to therapeutic applications.

Emerging Roles and Mechanisms of RNA Modifications in Neurodegenerative Diseases and Glioma

Despite a long history of research, neurodegenerative diseases and malignant brain tumor gliomas are both considered incurable, facing challenges in the development of treatments. Recent evidence suggests that RNA modifications, previously considered as static components of intracellular RNAs, are in fact dynamically regulated across various RNA species in cells and play a critical role in major biological processes in the nervous system. Innovations in next-generation sequencing have enabled the accurate detection of modifications on bases and sugars within various RNA molecules. These RNA modifications influence the stability and transportation of RNA, and crucially affect its translation. This review delves into existing knowledge on RNA modifications to offer a comprehensive inventory of these modifications across different RNA species. The detailed regulatory functions and roles of RNA modifications within the nervous system are discussed with a focus on neurodegenerative diseases and gliomas. This article presents a comprehensive overview of the fundamental mechanisms and emerging roles of RNA modifications in these diseases, which can facilitate the creation of innovative diagnostics and therapeutics for these conditions.

Beyond reader proteins: RNA binding proteins and RNA modifications in conversation to regulate gene expression

Post-transcriptional mRNA modifications play diverse roles in gene expression and RNA function. In many cases, RNA modifications function by altering how cellular machinery such as RNA binding proteins (RBPs) interact with RNA substrates. For instance, N6-methyladenosine (m6A) is recognized by the well-characterized YTH domain-containing family of "reader" proteins. For other mRNA modifications, similar global readers of modification status have not been clearly defined. Rather, most interactions between RBPs and RNA modifications have a more complicated dependence on sequence context and binding modality. The current handful of studies that demonstrate modifications impacting protein binding likely represent only a fraction of the full landscape. In this review, we dissect the known instances of RNA modifications altering RBP binding, specifically m6A, N1-methyladenosine (m1A), 5-methylcytosine (m5C), pseudouridine (Ψ), and internal N7-methylguanosine. We then review the biochemical properties of these and other identified mRNA modifications including dihydrouridine (D), N4-acetylcytosine (ac4C), and 2'-O-Methylation (Nme). We focus on how these properties would be likely to impact RNA:RBP interactions, including by changes to hydrogen bond potential, base-stacking efficiency, and RNA conformational preferences. The effects of RNA modifications on secondary structure have been well-studied, and we briefly discuss how structural effects imparted by modifications can lead to protein binding changes. Finally, we discuss strategies for uncovering as-yet-to-be identified modification-sensitive RBP:RNA Interactions. Coordinating future efforts to intersect the epitranscriptome and the RNA-protein interactome will illuminate the rules governing RNA modification recognition and the mechanisms responsible for the biological consequences of mRNA modification. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Processing > RNA Editing and Modification.

Genetic Algorithm for Automated Parameterization of Network Hamiltonian Models of Amyloid Fibril Formation

The time scales of long-time atomistic molecular dynamics simulations are typically reported in microseconds, while the time scales for experiments studying the kinetics of amyloid fibril formation are typically reported in minutes or hours. This time scale deficit of roughly 9 orders of magnitude presents a major challenge in the design of computer simulation methods for studying protein aggregation events. Coarse-grained molecular simulations offer a computationally tractable path forward for exploring the molecular mechanism driving the formation of these structures, which are implicated in diseases such as Alzheimer's, Parkinson's, and type-II diabetes. Network Hamiltonian models of aggregation are centered around a Hamiltonian function that returns the total energy of a system of aggregating proteins, given the graph structure of the system as an input. In the graph, or network, representation of the system, each protein molecule is represented as a node, and noncovalent bonds between proteins are represented as edges. The parameter, i.e., a set of coefficients that determine the degree to which each topological degree of freedom is favored or disfavored, must be determined for each network Hamiltonian model, and is a well-known technical challenge. The methodology is first demonstrated by beginning with an initial set of randomly parametrized models of low fibril fraction (<5% fibrillar), and evolving to subsequent generations of models, ultimately leading to high fibril fraction models (>70% fibrillar). The methodology is also demonstrated by applying it to optimizing previously published network Hamiltonian models for the 5 key amyloid fibril topologies that have been reported in the Protein Data Bank (PDB). The models generated by the AI produced fibril fractions that surpass previously published fibril fractions in 3 of 5 cases, including the most naturally abundant amyloid fibril topology, the 1,2 2-ribbon, which features a steric zipper. The authors also aim to encourage more widespread use of the network Hamiltonian methodology for fitting a wide variety of self-assembling systems by releasing a free open-source implementation of the genetic algorithm introduced here.

Epitranscriptomics as a New Layer of Regulation of Gene Expression in Skeletal Muscle: Known Functions and Future Perspectives

Epitranscriptomics refers to post-transcriptional regulation of gene expression via RNA modifications and editing that affect RNA functions. Many kinds of modifications of mRNA have been described, among which are N6-methyladenosine (m6A), N1-methyladenosine (m1A), 7-methylguanosine (m7G), pseudouridine (Ψ), and 5-methylcytidine (m5C). They alter mRNA structure and consequently stability, localization and translation efficiency. Perturbation of the epitranscriptome is associated with human diseases, thus opening the opportunity for potential manipulations as a therapeutic approach. In this review, we aim to provide an overview of the functional roles of epitranscriptomic marks in the skeletal muscle system, in particular in embryonic myogenesis, muscle cell differentiation and muscle homeostasis processes. Further, we explored high-throughput epitranscriptome sequencing data to identify RNA chemical modifications in muscle-specific genes and we discuss the possible functional role and the potential therapeutic applications.

Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA

Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.

High-Order Quantum-Mechanical Analysis of Hydrogen Bonding in Hachimoji and Natural DNA Base Pairs

High-order quantum chemistry is applied to hydrogen-bonded natural DNA nucleobase pairs [adenine:thymine (A:T) and guanine:cytosine (G:C)] and non-natural Hachimoji nucleobase pairs [isoguanine:1-methylcytosine (B:S) and 2-aminoimidazo[1,2a][1,3,5]triazin-4(1H)-one:6-amino-5-nitropyridin-2-one (P:Z)] to see how the intermolecular interaction energies and their energetic components (electrostatics, exchange-repulsion, induction/polarization, and London dispersion interactions) vary among the base pairs. We examined the Hoogsteen (HG) geometries in addition to the traditional Watson-Crick (WC) geometries. Coupled-cluster theory through perturbative triples [CCSD(T)] extrapolated to the complete basis set (CBS) limit and high-order symmetry-adapted perturbation theory (SAPT) at the SAPT2+(3)(CCD)δMP2/aug-cc-pVTZ level are used to estimate highly accurate noncovalent interaction energies. Electrostatic interactions are the most attractive component of the interaction energies, but the sum of induction/polarization and London dispersion is nearly as large, for all base pairs and geometries considered. Interestingly, the non-natural Hachimoji base pairs interact more strongly than the corresponding natural base pairs, by -21.8 (B:S) and -0.3 (P:Z) kcal mol^-1 in the WC geometries, according to CCSD(T)/CBS. This is consistent with the H-bond distances being generally shorter in the non-natural base pairs. The natural base pairs are energetically more stabilized in their Hoogsteen geometries than in their WC geometries. The Hoogsteen geometry makes the A:T base pair slightly more stable, by -0.8 kcal mol^-1, and it greatly stabilizes the G:C⁺ base pair, by -15.3 kcal mol^-1. The G:C⁺ stabilization is mainly due to the fact that C has typically added a proton when found in Hoogsteen geometries. By contrast, Hoogsteen geometries are substantially less favorable than WC geometries for non-natural Hachimoji base pairs, by 17.3 (B:S) and 13.8 (P:Z) kcal mol^-1.

RNA methylation and cellular response to oxidative stress-promoting anticancer agents

Disruption of the complex network that regulates redox homeostasis often underlies resistant phenotypes, which hinder effective and long-lasting cancer eradication. In addition, the RNA methylome-dependent control of gene expression also critically affects traits of cellular resistance to anti-cancer agents. However, few investigations aimed at establishing whether the epitranscriptome-directed adaptations underlying acquired and/or innate resistance traits in cancer could be implemented through the involvement of redox-dependent or -responsive signaling pathways. This is unexpected mainly because: i) the effectiveness of many anti-cancer approaches relies on their capacity to promote oxidative stress (OS); ii) altered redox milieu and reprogramming of mitochondrial function have been acknowledged as critical mediators of the RNA methylome-mediated response to OS. Here we summarize the current state of understanding on this topic, as well as we offer new perspectives that might lead to original approaches and strategies to delay or prevent the problem of refractory cancer and tumor recurrence.

The Repertoire of RNA Modifications Orchestrates a Plethora of Cellular Responses

Although a plethora of DNA modifications have been extensively investigated in the last decade, recent breakthroughs in molecular biology, including high throughput sequencing techniques, have enabled the identification of post-transcriptional marks that decorate RNAs; hence, epitranscriptomics has arisen. This recent scientific field aims to decode the regulatory layer of the transcriptome and set the ground for the detection of modifications in ribose nucleotides. Until now, more than 170 RNA modifications have been reported in diverse types of RNA that contribute to various biological processes, such as RNA biogenesis, stability, and transcriptional and translational accuracy. However, dysfunctions in the RNA-modifying enzymes that regulate their dynamic level can lead to human diseases and cancer. The present review aims to highlight the epitranscriptomic landscape in human RNAs and match the catalytic proteins with the deposition or deletion of a specific mark. In the current review, the most abundant RNA modifications, such as N6-methyladenosine (m⁶A), N5-methylcytosine (m⁵C), pseudouridine (Ψ) and inosine (I), are thoroughly described, their functional and regulatory roles are discussed and their contributions to cellular homeostasis are stated. Ultimately, the involvement of the RNA modifications and their writers, erasers, and readers in human diseases and cancer is also discussed.

Force-Field-Dependent DNA Breathing Dynamics: A Case Study of Hoogsteen Base Pairing in A6-DNA

The Hoogsteen (HG) base pairing conformation, commonly observed in damaged and mutated DNA helices, facilitates DNA repair and DNA recognition. The free energy difference between HG and Watson-Crick (WC) base pairs has been computed in previous studies. However, the mechanism of the conformational transition is not well understood. A detailed understanding of the process of WC to HG base pair transition can provide a deeper understanding of DNA repair and recognition. In an earlier study, we explored the free energy landscape for this process using extensive computer simulation with the CHARMM36 force field. In this work, we study the impact of force field models in describing the WC to HG base pairing transition using meta-eABF enhanced sampling, quasi-harmonic entropy calculation, and nonbonded energy analysis. The secondary structures of both base pairing forms and the topology of the free energy landscapes were consistent over different force field models, although the relative free energy, entropy, and the interaction energies tend to vary. The relative stability of the WC and HG conformations is dictated by a delicate balance between the enthalpic stabilization and the reduced entropy of the structurally rigid HG structure. These findings highlight the impact that subtleties in force field models can have on accurately modeling DNA base pair dynamics and should stimulate further computational investigations into other dynamically important motions in DNA.

N1-methyladenosine modification in cancer biology: current status and future perspectives

Post-transcriptional modifications in RNAs regulate their biological behaviors and functions. N¹-methyladenosine (m¹A), which is dynamically regulated by writers, erasers and readers, has been found as a reversible modification in tRNA, mRNA, rRNA and long non-coding RNA (lncRNA). m¹A modification has impacts on the RNA processing, structure and functions of targets. Increasing studies reveal the critical roles of m¹A modification and its regulators in tumorigenesis. Due to the positive relevance between m¹A and cancer development, targeting m¹A modification and m¹A-related regulators has been of attention. In this review, we summarized the current understanding of m¹A in RNAs, covering the modulation of m¹A modification in cancer biology, as well as the possibility of targeting m¹A modification as a potential target for cancer diagnosis and therapy.

Alternative RNA Conformations: Companion or Combatant

RNA molecules, in one form or another, are involved in almost all aspects of cell physiology, as well as in disease development. The diversity of the functional roles of RNA comes from its intrinsic ability to adopt complex secondary and tertiary structures, rivaling the diversity of proteins. The RNA molecules form dynamic ensembles of many interconverting conformations at a timescale of seconds, which is a key for understanding how they execute their cellular functions. Given the crucial role of RNAs in various cellular processes, we need to understand the RNA molecules from a structural perspective. Central to this review are studies aimed at revealing the regulatory role of conformational equilibria in RNA in humans to understand genetic diseases such as cancer and neurodegenerative diseases, as well as in pathogens such as bacteria and viruses so as to understand the progression of infectious diseases. Furthermore, we also summarize the prior studies on the use of RNA structures as platforms for the rational design of small molecules for therapeutic applications.

RNA modifications: importance in immune cell biology and related diseases

RNA modifications have become hot topics recently. By influencing RNA processes, including generation, transportation, function, and metabolization, they act as critical regulators of cell biology. The immune cell abnormality in human diseases is also a research focus and progressing rapidly these years. Studies have demonstrated that RNA modifications participate in the multiple biological processes of immune cells, including development, differentiation, activation, migration, and polarization, thereby modulating the immune responses and are involved in some immune related diseases. In this review, we present existing knowledge of the biological functions and underlying mechanisms of RNA modifications, including N⁶-methyladenosine (m⁶A), 5-methylcytosine (m⁵C), N¹-methyladenosine (m¹A), N⁷-methylguanosine (m⁷G), N⁴-acetylcytosine (ac⁴C), pseudouridine (Ψ), uridylation, and adenosine-to-inosine (A-to-I) RNA editing, and summarize their critical roles in immune cell biology. Via regulating the biological processes of immune cells, RNA modifications can participate in the pathogenesis of immune related diseases, such as cancers, infection, inflammatory and autoimmune diseases. We further highlight the challenges and future directions based on the existing knowledge. All in all, this review will provide helpful knowledge as well as novel ideas for the researchers in this area.

Probing Watson-Crick and Hoogsteen base pairing in duplex DNA using dynamic nuclear polarization solid-state NMR spectroscopy

The majority of base pairs in double-stranded DNA exist in the canonical Watson-Crick geometry. However, they can also adopt alternate Hoogsteen conformations in various complexes of DNA with proteins and small molecules, which are key for biological function and mechanism. While detection of Hoogsteen base pairs in large DNA complexes and assemblies poses considerable challenges for traditional structural biology techniques, we show here that multidimensional dynamic nuclear polarization-enhanced solid-state NMR can serve as a unique spectroscopic tool for observing and distinguishing Watson-Crick and Hoogsteen base pairs in a broad range of DNA systems based on characteristic NMR chemical shifts and internuclear dipolar couplings. We illustrate this approach using a model 12-mer DNA duplex, free and in complex with the antibiotic echinomycin, which features two central adenine-thymine base pairs with Watson-Crick and Hoogsteen geometry, respectively, and subsequently extend it to the ∼200 kDa Widom 601 DNA nucleosome core particle.

RNA m1A methylation regulates glycolysis of cancer cells through modulating ATP5D

Studies on biological functions of RNA modifications such as N⁶-methyladenosine (m⁶A) in mRNA have sprung up in recent years, while the roles of N¹-methyladenosine (m¹A) in cancer progression remain largely unknown. We find m¹A demethylase ALKBH3 can regulate the glycolysis of cancer cells via a demethylation activity dependent manner. Specifically, sequencing and functional studies confirm that ATP5D, one of the most important subunit of adenosine 5'-triphosphate synthase, is involved in m¹A demethylase ALKBH3-regulated glycolysis of cancer cells. The m¹A modified A71 at the exon 1 of ATP5D negatively regulates its translation elongation via increasing the binding with YTHDF1/eRF1 complex, which facilitates the release of message RNA (mRNA) from ribosome complex. m¹A also regulates mRNA stability of E2F1, which directly binds with ATP5D promoter to initiate its transcription. Targeted specific demethylation of ATP5D m¹A by dm¹ACRISPR system can significantly increase the expression of ATP5D and glycolysis of cancer cells. In vivo data confirm the roles of m¹A/ATP5D in tumor growth and cancer progression. Our study reveals a crosstalk of mRNA m¹A modification and cell metabolism, which expands the understanding of such interplays that are essential for cancer therapeutic application.

Regulation and roles of RNA modifications in aging-related diseases

With the aging of the global population, accumulating interest is focused on manipulating the fundamental aging-related signaling pathways to delay the physiological aging process and eventually slow or prevent the appearance or severity of multiple aging-related diseases. Recently, emerging evidence has shown that RNA modifications, which were historically considered infrastructural features of cellular RNAs, are dynamically regulated across most of the RNA species in cells and thereby critically involved in major biological processes, including cellular senescence and aging. In this review, we summarize the current knowledge about RNA modifications and provide a catalog of RNA modifications on different RNA species, including mRNAs, miRNAs, lncRNA, tRNAs, and rRNAs. Most importantly, we focus on the regulation and roles of these RNA modifications in aging-related diseases, including neurodegenerative diseases, cardiovascular diseases, cataracts, osteoporosis, and fertility decline. This would be an important step toward a better understanding of fundamental aging mechanisms and thereby facilitating the development of novel diagnostics and therapeutics for aging-related diseases.

Measuring thermodynamic preferences to form non-native conformations in nucleic acids using ultraviolet melting

Thermodynamic preferences to form non-native conformations are crucial for understanding how nucleic acids fold and function. However, they are difficult to measure experimentally because this requires accurately determining the population of minor low-abundance (<10%) conformations in a sea of other conformations. Here, we show that melting experiments enable facile measurements of thermodynamic preferences to adopt nonnative conformations in DNA and RNA. The key to this "delta-melt" approach is to use chemical modifications to render specific minor non-native conformations the major state. The validity and robustness of delta-melt is established for four different non-native conformations under various physiological conditions and sequence contexts through independent measurements of thermodynamic preferences using NMR. Delta-melt is faster relative to NMR, simple, and cost-effective and enables thermodynamic preferences to be measured for exceptionally low-populated conformations. Using delta-melt, we obtained rare insights into conformational cooperativity, obtaining evidence for significant cooperativity (1.0 to 2.5 kcal/mol) when simultaneously forming two adjacent Hoogsteen base pairs. We also measured the thermodynamic preferences to form G-C⁺ and A-T Hoogsteen and A-T base open states for nearly all 16 trinucleotide sequence contexts and found distinct sequence-specific variations on the order of 2 to 3 kcal/mol. This rich landscape of sequence-specific non-native minor conformations in the DNA double helix may help shape the sequence specificity of DNA biochemistry. Thus, melting experiments can now be used to access thermodynamic information regarding regions of the free energy landscape of biomolecules beyond the native folded and unfolded conformations.

Development of Mild Chemical Catalysis Conditions for m1A-to-m6A Rearrangement on RNA

The conversion of N¹-methyladenosine (m¹A) to N⁶-methyladenosine (m⁶A) on RNA is an important step for both allowing efficient reverse transcription read-though for sequencing analysis and mapping modifications in the transcriptome. Enzymatic transformation is often used, but the efficiency of the removal can depend on local sequence context. Chemical conversion through the application of the Dimroth rearrangement, in which m¹A rearranges into m⁶A under heat and alkaline conditions, is an alternative, but the required alkaline conditions result in significant RNA degradation by hydrolysis of the phosphodiester backbone. Here, we report novel, mild pH conditions that catalyze m¹A-to-m⁶A arrangement using 4-nitrothiophenol as a catalyst. We demonstrate the efficient rearrangement in mononucleosides, synthetic RNA oligonucleotides, and RNAs isolated from human cell lines, thereby validating a new approach for converting m¹A-to-m⁶A in RNA samples for sequencing analyses.

m1A RNA Modification in Gene Expression Regulation

N1-methyladenosine (m1A) is a prevalent and reversible post-transcriptional RNA modification that decorates tRNA, rRNA and mRNA. Recent studies based on technical advances in analytical chemistry and high-throughput sequencing methods have revealed the crucial roles of m1A RNA modification in gene regulation and biological processes. In this review, we focus on progress in the study of m1A methyltransferases, m1A demethylases and m1A-dependent RNA-binding proteins and highlight the biological mechanisms and functions of m1A RNA modification, as well as its association with human disease. We also summarize the current understanding of detection approaches for m1A RNA modification.

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

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

TRMT6/61A-dependent base methylation of tRNA-derived fragments regulates gene-silencing activity and the unfolded protein response in bladder cancer

RNA modifications are important regulatory elements of RNA functions. However, most genome-wide mapping of RNA modifications has focused on messenger RNAs and transfer RNAs, but such datasets have been lacking for small RNAs. Here we mapped N1-methyladenosine (m1A) in the cellular small RNA space. Benchmarked with synthetic m1A RNAs, our workflow identified specific groups of m1A-containing small RNAs, which are otherwise disproportionally under-represented. In particular, 22-nucleotides long 3' tRNA-fragments are highly enriched for TRMT6/61A-dependent m1A located within the seed region. TRMT6/61A-dependent m1A negatively affects gene silencing by tRF-3s. In urothelial carcinoma of the bladder, where TRMT6/61A is over-expressed, higher m1A modification on tRFs is detected, correlated with a dysregulation of tRF targetome. Lastly, TRMT6/61A regulates tRF-3 targets involved in unfolded protein response. Together, our results reveal a mechanism of regulating gene expression via base modification of small RNA.

Calculations of pKa Values for a Series of Naturally Occurring Modified Nucleobases

Modified nucleobases are found in functionally important regions of RNA and are often responsible for essential structural roles. Many of these nucleobase modifications are dynamically regulated in nature, with each modification having a different biological role in RNA. Despite the high abundance of modifications, many of their characteristics are still poorly understood. One important property of a nucleobase is its pK_a value, which has been widely studied for unmodified nucleobases, but not for the modified versions. In this study, the pK_a values of modified nucleobases were determined by performing ab initio quantum mechanical calculations using a B3LYP density functional with the 6-31+G(d,p) basis set and a combination of implicit-explicit solvation systems. This method, which was previously employed to determine the pK_a values of unmodified nucleobases, is applicable to a variety of modified nucleobases. Comparisons of the pK_a values of modified nucleobases give insight into their structural and energetic impacts within nucleic acids.

Control of animal virus replication by RNA adenosine methylation

Methylation at the N⁶-position of either adenosine (m⁶A) or 2'-O-methyladenosine (m⁶Am) represents two of the most abundant internal modifications of coding and non-coding RNAs, influencing their maturation, stability and function. Additionally, although less abundant and less well-studied, monomethylation at the N¹-position (m¹A) can have profound effects on RNA folding. It has been known for several decades that RNAs produced by both DNA and RNA viruses can be m⁶A/m⁶Am modified and the list continues to broaden through advances in detection technologies and identification of the relevant methyltransferases. Recent studies have uncovered varied mechanisms used by viruses to manipulate the m⁶A pathway in particular, either to enhance virus replication or to antagonize host antiviral defenses. As such, RNA modifications represent an important frontier of exploration in the broader realm of virus-host interactions, and this new knowledge already suggests exciting opportunities for therapeutic intervention. In this review we summarize the principal mechanisms by which m⁶A/m⁶Am can promote or hinder viral replication, describe how the pathway is actively manipulated by biomedically important viruses, and highlight some remaining gaps in understanding how adenosine methylation of RNA controls viral replication and pathogenesis.

Structural Insights Into Tautomeric Dynamics in Nucleic Acids and in Antiviral Nucleoside Analogs

DNA (2'-deoxyribonucleic acid) and RNA (ribonucleic acid) play diverse functional roles in biology and disease. Despite being comprised primarily of only four cognate nucleobases, nucleic acids can adopt complex three-dimensional structures, and RNA in particular, can catalyze biochemical reactions to regulate a wide variety of biological processes. Such chemical versatility is due in part to the phenomenon of nucleobase tautomerism, whereby the bases can adopt multiple, yet distinct isomeric forms, known as tautomers. For nucleobases, tautomers refer to structural isomers that differ from one another by the position of protons. By altering the position of protons on nucleobases, many of which play critical roles for hydrogen bonding and base pairing interactions, tautomerism has profound effects on the biochemical processes involving nucleic acids. For example, the transient formation of minor tautomers during replication could generate spontaneous mutations. These mutations could arise from the stabilization of mismatches, in the active site of polymerases, in conformations involving minor tautomers that are indistinguishable from canonical base pairs. In this review, we discuss the evidence for tautomerism in DNA, and its consequences to the fidelity of DNA replication. Also reviewed are RNA systems, such as the riboswitches and self-cleaving ribozymes, in which tautomerism plays a functional role in ligand recognition and catalysis, respectively. We also discuss tautomeric nucleoside analogs that are efficacious as antiviral drug candidates such as molnupiravir for coronaviruses and KP1212 for HIV. The antiviral efficacy of these analogs is due, in part, to their ability to exist in multiple tautomeric forms and induce mutations in the replicating viral genomes. From a technical standpoint, minor tautomers of nucleobases are challenging to identify directly because they are rare and interconvert on a fast, millisecond to nanosecond, time scale. Nevertheless, many approaches including biochemical, structural, computational and spectroscopic methods have been developed to study tautomeric dynamics in RNA and DNA systems, and in antiviral nucleoside analogs. An overview of these methods and their applications is included here.

Revealing A-T and G-C Hoogsteen base pairs in stressed protein-bound duplex DNA

Watson–Crick base pairs (bps) are the fundamental unit of genetic information and the building blocks of the DNA double helix. However, A-T and G-C can also form alternative ‘Hoogsteen’ bps, expanding the functional complexity of DNA. We developed ‘Hoog-finder’, which uses structural fingerprints to rapidly screen Hoogsteen bps, which may have been mismodeled as Watson–Crick in crystal structures of protein–DNA complexes. We uncovered 17 Hoogsteen bps, 7 of which were in complex with 6 proteins never before shown to bind Hoogsteen bps. The Hoogsteen bps occur near mismatches, nicks and lesions and some appear to participate in recognition and damage repair. Our results suggest a potentially broad role for Hoogsteen bps in stressed regions of the genome and call for a community-wide effort to identify these bps in current and future crystal structures of DNA and its complexes.

A Small Molecule Exploits Hidden Structural Features within the RNA Repeat Expansion That Causes c9ALS/FTD and Rescues Pathological Hallmarks

The hexanucleotide repeat expansion GGGGCC [r(G₄C₂)^exp] within intron 1 of C9orf72 causes genetically defined amyotrophic lateral sclerosis and frontotemporal dementia, collectively named c9ALS/FTD. , the repeat expansion causes neurodegeneration via deleterious phenotypes stemming from r(G₄C₂)^exp RNA gain- and loss-of-function mechanisms. The r(G₄C₂)^exp RNA folds into both a hairpin structure with repeating 1 × 1 nucleotide GG internal loops and a G-quadruplex structure. Here, we report the identification of a small molecule (CB253) that selectively binds the hairpin form of r(G₄C₂)^exp. Interestingly, the small molecule binds to a previously unobserved conformation in which the RNA forms 2 × 2 nucleotide GG internal loops, as revealed by a series of binding and structural studies. NMR and molecular dynamics simulations suggest that the r(G₄C₂)^exp hairpin interconverts between 1 × 1 and 2 × 2 internal loops through the process of strand slippage. We provide experimental evidence that CB253 binding indeed shifts the equilibrium toward the 2 × 2 GG internal loop conformation, inhibiting mechanisms that drive c9ALS/FTD pathobiology, such as repeat-associated non-ATG translation formation of stress granules and defective nucleocytoplasmic transport in various cellular models of c9ALS/FTD.

Rapid assessment of Watson-Crick to Hoogsteen exchange in unlabeled DNA duplexes using high-power SELOPE imino 1H CEST

In duplex DNA, Watson-Crick A-T and G-C base pairs (bp's) exist in dynamic equilibrium with an alternative Hoogsteen conformation, which is low in abundance and short-lived. Measuring how the Hoogsteen dynamics varies across different DNA sequences, structural contexts and physiological conditions is key for identifying potential Hoogsteen hot spots and for understanding the potential roles of Hoogsteen base pairs in DNA recognition and repair. However, such studies are hampered by the need to prepare 13 C or 15 N isotopically enriched DNA samples for NMR relaxation dispersion (RD) experiments. Here, using SELective Optimized Proton Experiments (SELOPE) 1 H CEST experiments employing high-power radiofrequency fields (B 1  >  250 Hz) targeting imino protons, we demonstrate accurate and robust characterization of Watson-Crick to Hoogsteen exchange, without the need for isotopic enrichment of the DNA. For 13 residues in three DNA duplexes under different temperature and pH conditions, the exchange parameters deduced from high-power imino 1 H CEST were in very good agreement with counterparts measured using off-resonance 13 C /  15 N spin relaxation in the rotating frame (R 1 ρ ). It is shown that 1 H-1 H NOE effects which typically introduce artifacts in 1 H-based measurements of chemical exchange can be effectively suppressed by selective excitation, provided that the relaxation delay is short (≤  100 ms). The 1 H CEST experiment can be performed with ∼  10× higher throughput and ∼  100× lower cost relative to 13 C /  15 N R 1 ρ and enabled Hoogsteen chemical exchange measurements undetectable by R 1 ρ . The results reveal an increased propensity to form Hoogsteen bp's near terminal ends and a diminished propensity within A-tract motifs. The 1 H CEST experiment provides a basis for rapidly screening Hoogsteen breathing in duplex DNA, enabling identification of unusual motifs for more in-depth characterization.

Probing the Hydrogen-Bonding Environment of Individual Bases in DNA Duplexes with Isotope-Edited Infrared Spectroscopy

Measuring the strength of the hydrogen bonds between DNA base pairs is of vital importance for understanding how our genetic code is physically accessed and recognized in cells, particularly during replication and transcription. Therefore, it is important to develop probes for these key hydrogen bonds (H-bonds) that dictate events critical to cellular function, such as the localized melting of DNA. The vibrations of carbonyl bonds are well-known probes of their H-bonding environment, and their signals can be observed with infrared (IR) spectroscopy. Yet, pinpointing a single bond of interest in the complex IR spectrum of DNA is challenging due to the large number of carbonyl signals that overlap with each other. Here, we develop a method using isotope editing and infrared (IR) spectroscopy to isolate IR signals from the thymine (T) C2=O carbonyl. We use solvatochromatic studies to show that the TC2=O signal’s position in the IR spectrum is sensitive to the H-bonding capacity of the solvent. Our results indicate that C2=O of a single T base within DNA duplexes experiences weak H-bonding interactions. This finding is consistent with the existence of a third, noncanonical CH···O H-bond between adenine and thymine in both Watson–Crick and Hoogsteen base pairs in DNA.

Response to stress in biological disorders: Implications of stress granule assembly and function

It is indispensable for cells to adapt and respond to environmental stresses, in order for organisms to survive. Stress granules (SGs) are condensed membrane‐less organelles dynamically formed in the cytoplasm of eukaryotes cells to cope with diverse intracellular or extracellular stress factors, with features of liquid‐liquid phase separation. They are composed of multiple constituents, including translationally stalled mRNAs, translation initiation factors, RNA‐binding proteins and also non‐RNA‐binding proteins. SG formation is triggered by stress stimuli, viral infection and signal transduction, while aberrant assembly of SGs may contribute to tissue degenerative diseases. Recently, a growing body of evidence has emerged on SG response mechanisms for cells facing high temperatures, oxidative stress and osmotic stress. In this review, we aim to summarize factors affecting SGs assembly, present the impact of SGs on germ cell development and other biological processes. We particularly emphasize the significance of recently reported RNA modifications in SG stress responses. In parallel, we also review all current perspectives on the roles of SGs in male germ cells, with a particular focus on the dynamics of SG assembly.

Duplexes Formed by G4C2 Repeats Contain Alternate Slow- and Fast-Flipping G·G Base Pairs

Aberrant expansion of the hexanucleotide GGGGCC (or G4C2) repeat in the human C9ORF72 gene is the most common genetic factor found behind amyotrophic lateral sclerosis and frontotemporal dementia. The hypothesized pathways, through which the repeat expansions contribute to the pathology, involve one or more secondary structural forms of the DNA and/or RNA sequences, such as G-quadruplexes, duplexes, and hairpins. Here, we study the structures of DNA and RNA duplexes formed by G4C2 repeats, which contain G(syn)·G(anti) base pairs flanked by either G·C or C·G base pairs. We show that duplexes formed by G4C2 repeats contain alternately two types of G·G pair contexts exhibiting different syn-anti base flipping dynamics (∼100 ms vs ∼2 ms for DNA and ∼50 ms vs ∼20 ms for RNA at 10 °C, respectively) depending on the flanking bases, with the slow-flipping G·G pairs being flanked by a guanine at the 5'-end and the fast-flipping G·G pairs being flanked by a cytosine at the 5'-end. Our findings on the structures and dynamics of G·G base pairs in DNA and RNA duplexes formed by G4C2 repeats provide a foundation for further studies of the functions and targeting of such biologically relevant motifs.

Instrumental analysis of RNA modifications

Organisms from all domains of life invest a substantial amount of energy for the introduction of RNA modifications into nearly all transcripts studied to date. Instrumental analysis of RNA can focus on the modified residues and reveal the function of these epitranscriptomic marks. Here, we will review recent advances and breakthroughs achieved by NMR spectroscopy, sequencing, and mass spectrometry of the epitranscriptome.

The protective role of m1A during stress-induced granulation

Post-transcriptional methylation of N6-adenine and N1-adenine can affect transcriptome turnover and translation. Furthermore, the regulatory function of N6-methyladenine (m6A) during heat shock has been uncovered, including the enhancement of the phase separation potential of RNAs. In response to acute stress, e.g. heat shock, the orderly sequestration of mRNAs in stress granules (SGs) is considered important to protect transcripts from the irreversible aggregation. Until recently, the role of N1-methyladenine (m1A) on mRNAs during acute stress response remains largely unknown. Here we show that the methyltransferase complex TRMT6/61A, which generates the m1A tag, is involved in transcriptome protection during heat shock. Our bioinformatics analysis indicates that occurrence of the m1A motif is increased in mRNAs known to be enriched in SGs. Accordingly, the m1A-generating methyltransferase TRMT6/61A accumulated in SGs and mass spectrometry confirmed enrichment of m1A in the SG RNAs. The insertion of a single methylation motif in the untranslated region of a reporter RNA leads to more efficient recovery of protein synthesis from that transcript after the return to normal temperature. Our results demonstrate far-reaching functional consequences of a minimal RNA modification on N1-adenine during acute proteostasis stress.

Interrogation of highly structured RNA with multicomponent deoxyribozyme probes at ambient temperatures

Molecular analysis of RNA through hybridization with sequence-specific probes is challenging due to the intrinsic ability of RNA molecules to form stable secondary and tertiary structures. To overcome the energy barrier toward the probe-RNA complex formation, the probes are made of artificial nucleotides, which are more expensive than their natural counterparts and may still be inefficient. Here, we propose the use of a multicomponent probe based on an RNA-cleaving deoxyribozyme for the analysis of highly structured RNA targets. Efficient interrogation of two native RNA from Saccharomyces cerevisiae-a transfer RNA (tRNA) and 18S ribosomal RNA (rRNA)-was achieved at ambient temperature. We achieved detection limits of tRNA down to ∼0.3 nM, which is two orders of magnitude lower than that previously reported for molecular beacon probes. Importantly, no probe annealing to the target was required, with the hybridization assay performed at 37°C. Excess of nonspecific targets did not compromise the performance of the probe, and high interrogation efficiency was maintained by the probes even in complex matrices, such as cell lysate. A linear dynamic range of 0.3-150 nM tRNA was demonstrated. The probe can be adapted for differentiation of a single mismatch in the tRNA-probe complex. Therefore, this study opens a venue toward highly selective, sensitive, robust, and inexpensive assays for the interrogation of biological RNA.

Free Energy Landscape and Conformational Kinetics of Hoogsteen Base Pairing in DNA vs. RNA

Genetic information is encoded in the DNA double helix, which, in its physiological milieu, is characterized by the iconical Watson-Crick nucleo-base pairing. Recent NMR relaxation experiments revealed the transient presence of an alternative, Hoogsteen (HG) base pairing pattern in naked DNA duplexes, and estimated its relative stability and lifetime. In contrast with DNA, such structures were not observed in RNA duplexes. Understanding HG base pairing is important because the underlying "breathing" motion between the two conformations can significantly modulate protein binding. However, a detailed mechanistic insight into the transition pathways and kinetics is still missing. We performed enhanced sampling simulation (with combined metadynamics and adaptive force-bias method) and Markov state modeling to obtain accurate free energy, kinetics, and the intermediates in the transition pathway between Watson-Crick and HG base pairs for both naked B-DNA and A-RNA duplexes. The Markov state model constructed from our unbiased MD simulation data revealed previously unknown complex extrahelical intermediates in the seemingly simple process of base flipping in B-DNA. Extending our calculation to A-RNA, for which HG base pairing is not observed experimentally, resulted in relatively unstable, single-hydrogen-bonded, distorted Hoogsteen-like bases. Unlike B-DNA, the transition pathway primarily involved base paired and intrahelical intermediates with transition timescales much longer than that of B-DNA. The seemingly obvious flip-over reaction coordinate (i.e., the glycosidic torsion angle) is unable to resolve the intermediates. Instead, a multidimensional picture involving backbone dihedral angles and distance between hydrogen bond donor and acceptor atoms is required to gain insight into the molecular mechanism.

Hoogsteen base pairs increase the susceptibility of double-stranded DNA to cytotoxic damage

As the Watson-Crick faces of nucleobases are protected in dsDNA, it is commonly assumed that deleterious alkylation damage to the Watson-Crick faces of nucleobases predominantly occurs when DNA becomes single-stranded during replication and transcription. However, damage to the Watson-Crick faces of nucleobases has been reported in dsDNA in vitro through mechanisms that are not understood. In addition, the extent of protection from methylation damage conferred by dsDNA relative to ssDNA has not been quantified. Watson-Crick base pairs in dsDNA exist in dynamic equilibrium with Hoogsteen base pairs that expose the Watson-Crick faces of purine nucleobases to solvent. Whether this can influence the damage susceptibility of dsDNA remains unknown. Using dot-blot and primer extension assays, we measured the susceptibility of adenine-N1 to methylation by dimethyl sulfate (DMS) when in an A-T Watson-Crick versus Hoogsteen conformation. Relative to unpaired adenines in a bulge, Watson-Crick A-T base pairs in dsDNA only conferred ∼130-fold protection against adenine-N1 methylation, and this protection was reduced to ∼40-fold for A(syn)-T Hoogsteen base pairs embedded in a DNA-drug complex. Our results indicate that Watson-Crick faces of nucleobases are accessible to alkylating agents in canonical dsDNA and that Hoogsteen base pairs increase this accessibility. Given the higher abundance of dsDNA relative to ssDNA, these results suggest that dsDNA could be a substantial source of cytotoxic damage. The work establishes DMS probing as a method for characterizing A(syn)-T Hoogsteen base pairs in vitro and also lays the foundation for a sequencing approach to map A(syn)-T Hoogsteen and unpaired adenines genome-wide in vivo.

Impact of 5-formylcytosine on the melting kinetics of DNA by 1H NMR chemical exchange

5-Formylcytosine (5fC) is a chemically edited, naturally occurring nucleobase which appears in the context of modified DNA strands. The understanding of the impact of 5fC on dsDNA physical properties is to date limited. In this work, we applied temperature-dependent 1H Chemical Exchange Saturation Transfer (CEST) NMR experiments to non-invasively and site-specifically measure the thermodynamic and kinetic influence of formylated cytosine nucleobase on the melting process involving dsDNA. Incorporation of 5fC within symmetrically positioned CpG sites destabilizes the whole dsDNA structure-as witnessed from the ∼2°C decrease in the melting temperature and 5-10 kJ mol-1 decrease in ΔG°-and affects the kinetic rates of association and dissociation. We observed an up to ∼5-fold enhancement of the dsDNA dissociation and an up to ∼3-fold reduction in ssDNA association rate constants, over multiple temperatures and for several proton reporters. Eyring and van't Hoff analysis proved that the destabilization is not localized, instead all base-pairs are affected and the transition states resembles the single-stranded conformation. These results advance our knowledge about the role of 5fC as a semi-permanent epigenetic modification and assist in the understanding of its interactions with reader proteins.

Probing conformational transitions towards mutagenic Watson-Crick-like G·T mismatches using off-resonance sugar carbon R1ρ relaxation dispersion

NMR off-resonance R_1ρ relaxation dispersion measurements on base carbon and nitrogen nuclei have revealed that wobble G·T/U mismatches in DNA and RNA duplexes exist in dynamic equilibrium with short-lived, low-abundance, and mutagenic Watson-Crick-like conformations. As Watson-Crick-like G·T mismatches have base pairing geometries similar to Watson-Crick base pairs, we hypothesized that they would mimic Watson-Crick base pairs with respect to the sugar-backbone conformation as well. Using off-resonance R_1ρ measurements targeting the sugar C3' and C4' nuclei, a structure survey, and molecular dynamics simulations, we show that wobble G·T mismatches adopt sugar-backbone conformations that deviate from the canonical Watson-Crick conformation and that transitions toward tautomeric and anionic Watson-Crick-like G·T mismatches restore the canonical Watson-Crick sugar-backbone. These measurements also reveal kinetic isotope effects for tautomerization in D₂O versus H₂O, which provide experimental evidence in support of a transition state involving proton transfer. The results provide additional evidence in support of mutagenic Watson-Crick-like G·T mismatches, help rule out alternative inverted wobble conformations in the case of anionic G·T^-, and also establish sugar carbons as new non-exchangeable probes of this exchange process.

A molecular-level perspective on the frequency, distribution, and consequences of messenger RNA modifications

Cells use chemical modifications to alter the sterics, charge, and conformations of large biomolecules, modulating their biogenesis, function, and stability. Until recently post-transcriptional RNA modifications were thought to be largely limited to nonprotein coding RNA species. However, this dogma has rapidly transformed with the discovery of a host of modifications in protein coding messenger RNAs (mRNAs). Recent advancements in genome-wide sequencing technologies have enabled the identification of mRNA modifications as a potential new frontier in gene regulation-leading to the development of the epitranscriptome field. As a result, there has been a flurry of multiple groundbreaking discoveries, including new modifications, nucleoside modifying enzymes ("writers" and "erasers"), and RNA binding proteins that recognize chemical modifications ("readers"). These discoveries opened the door to understanding how post-transcriptional mRNA modifications can modulate the mRNA lifecycle, and established a link between the epitranscriptome and human health and disease. Despite a rapidly growing recognition of their importance, fundamental questions regarding the identity, prevalence, and functional consequences of mRNA modifications remain to be answered. Here, we highlight quantitative studies that characterize mRNA modification abundance, frequency, and interactions with cellular machinery. As the field progresses, we see a need for the further integration of quantitative and reductionist approaches to complement transcriptome wide studies in order to establish a molecular-level framework for understanding the consequences of mRNA chemical modifications on biological processes. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Processing > RNA Editing and Modification.

A Census and Categorization Method of Epitranscriptomic Marks

In the past few years, thorough investigation of chemical modifications operated in the cells on ribonucleic acid (RNA) molecules is gaining momentum. This new field of research has been dubbed “epitranscriptomics”, in analogy to best-known epigenomics, to stress the potential of ensembles of RNA modifications to constitute a post-transcriptional regulatory layer of gene expression orchestrated by writer, reader, and eraser RNA-binding proteins (RBPs). In fact, epitranscriptomics aims at identifying and characterizing all functionally relevant changes involving both non-substitutional chemical modifications and editing events made to the transcriptome. Indeed, several types of RNA modifications that impact gene expression have been reported so far in different species of cellular RNAs, including ribosomal RNAs, transfer RNAs, small nuclear RNAs, messenger RNAs, and long non-coding RNAs. Supporting functional relevance of this largely unknown regulatory mechanism, several human diseases have been associated directly to RNA modifications or to RBPs that may play as effectors of epitranscriptomic marks. However, an exhaustive epitranscriptome’s characterization, aimed to systematically classify all RNA modifications and clarify rules, actors, and outcomes of this promising regulatory code, is currently not available, mainly hampered by lack of suitable detecting technologies. This is an unfortunate limitation that, thanks to an unprecedented pace of technological advancements especially in the sequencing technology field, is likely to be overcome soon. Here, we review the current knowledge on epitranscriptomic marks and propose a categorization method based on the reference ribonucleotide and its rounds of modifications (“stages”) until reaching the given modified form. We believe that this classification scheme can be useful to coherently organize the expanding number of discovered RNA modifications.

Predicting sites of epitranscriptome modifications using unsupervised representation learning based on generative adversarial networks

Epitranscriptome is an exciting area that studies different types of modifications in transcripts and the prediction of such modification sites from the transcript sequence is of significant interest. However, the scarcity of positive sites for most modifications imposes critical challenges for training robust algorithms. To circumvent this problem, we propose MR-GAN, a generative adversarial network (GAN) based model, which is trained in an unsupervised fashion on the entire pre-mRNA sequences to learn a low dimensional embedding of transcriptomic sequences. MR-GAN was then applied to extract embeddings of the sequences in a training dataset we created for eight epitranscriptome modifications, including m⁶A, m¹A, m¹G, m²G, m⁵C, m⁵U, 2'-O-Me, Pseudouridine (Ψ) and Dihydrouridine (D), of which the positive samples are very limited. Prediction models were trained based on the embeddings extracted by MR-GAN. We compared the prediction performance with the one-hot encoding of the training sequences and SRAMP, a state-of-the-art m⁶A site prediction algorithm and demonstrated that the learned embeddings outperform one-hot encoding by a significant margin for up to 15% improvement. Using MR-GAN, we also investigated the sequence motifs for each modification type and uncovered known motifs as well as new motifs not possible with sequences directly. The results demonstrated that transcriptome features extracted using unsupervised learning could lead to high precision for predicting multiple types of epitranscriptome modifications, even when the data size is small and extremely imbalanced.

The roles of structural dynamics in the cellular functions of RNAs

RNAs fold into 3D structures that range from simple helical elements to complex tertiary structures and quaternary ribonucleoprotein assemblies. The functions of many regulatory RNAs depend on how their 3D structure changes in response to a diverse array of cellular conditions. In this Review, we examine how the structural characterization of RNA as dynamic ensembles of conformations, which form with different probabilities and at different timescales, is improving our understanding of RNA function in cells. We discuss the mechanisms of gene regulation by microRNAs, riboswitches, ribozymes, post-transcriptional RNA modifications and RNA-binding proteins, and how the cellular environment and processes such as liquid-liquid phase separation may affect RNA folding and activity. The emerging RNA-ensemble-function paradigm is changing our perspective and understanding of RNA regulation, from in vitro to in vivo and from descriptive to predictive.

Computational evolution of an RNA-binding protein towards enhanced oxidized-RNA binding

The oxidation of RNA has been implicated in the development of many diseases. Among the four ribonucleotides, guanosine is the most susceptible to oxidation, resulting in the formation of 8-oxo-7,8-dihydroguanosine (8-oxoG). Despite the limited knowledge about how cells regulate the detrimental effects of oxidized RNA, cellular factors involved in its regulation have begun to be identified. One of these factors is polynucleotide phosphorylase (PNPase), a multifunctional enzyme implicated in RNA turnover. In the present study, we have examined the interaction of PNPase with 8-oxoG in atomic detail to provide insights into the mechanism of 8-oxoG discrimination. We hypothesized that PNPase subunits cooperate to form a binding site using the dynamic SFF loop within the central channel of the PNPase homotrimer. We evolved this site using a novel approach that initially screened mutants from a library of beneficial mutations and assessed their interactions using multi-nanosecond Molecular Dynamics simulations. We found that evolving this single site resulted in a fold change increase in 8-oxoG affinity between 1.2 and 1.5 and/or selectivity between 1.5 and 1.9. In addition to the improvement in 8-oxoG binding, complementation of K12 Δpnp with plasmids expressing mutant PNPases caused increased cell tolerance to H₂O₂. This observation provides a clear link between molecular discrimination of RNA oxidation and cell survival. Moreover, this study provides a framework for the manipulation of modified-RNA protein readers, which has potential application in synthetic biology and epitranscriptomics.