Atomistic and Thermodynamic Analysis of N6-Methyladenosine (m6A) Recognition by the Reader Domain of YTHDC1

Thermodynamic insights into N6-methyladenosine-modified ribonucleic acids and their interactions with the RNA recognition motif of heterogeneous nuclear ribonucleoprotein C

N6-methyladenosine (m6A) is a prevalent RNA modification that regulates key functions such as splicing, transport, translation, and stability across various RNA types, including mRNA, tRNA, rRNA, and lncRNA. Transcriptome-wide studies reveal that approximately one-third of mammalian mRNAs carry 3-5 m6A modifications, enriched in the consensus motif RRA*CH. While some studies suggest m6A induces structural changes in RNA to facilitate protein binding through an "m6A switch" mechanism, others propose it primarily primes RNA for enhanced protein interactions, emphasizing the need for further exploration of m6A's role. Here, we investigated how m6A influences the binding of heterogeneous nuclear ribonucleoprotein C (hnRNPC), which recognizes poly (U) tracts via its RNA recognition motif (RRM). Using naturally occurring m6A-modified RNAs, including lncRNA MALAT1, we examined the effects of m6A on RNA folding and protein binding. Biophysical experiments (UV melting, circular dichroism, and molecular dynamics simulations) revealed that m6A subtly alters RNA stability and folding. Binding studies using EMSA, Microscale Thermophoresis (MST), and Isothermal Titration Calorimetry (ITC) showed m6A primes RNA for hnRNPC recognition rather than inducing structural switches. These findings refine our understanding of m6A's role in RNA-protein interactions, highlighting its regulatory importance in RNA metabolism.

Systematic molecular profiling of non-native N6-substitution effects on m6A binding to the YTH domains of human RNA m6A readers in diabetes

The RNA N6-adenosine methylation, resulting in N6-methyl adenosine (m6A), is one of the most important post-transcriptional modification events in the eukaryotic transcriptome, which is dynamically regulated by methyltransferases (writers), recognition proteins (readers) and demethylases (erasers). Human has five m6A readers namely YTHDC1, YTHDC2, YTHDF1, YTHDF2 and YTHDF3 that specifically recognize and bind to the methylated m6A residue of RNA through their YT521-B homology (YTH) domains, which have been involved in the pathogenesis of diabetes mellitus and its diverse complications such as diabetic nephropathy. Instead of the native N6-methylation, we herein attempted to explore the molecular effect of various non-native N6-substitutions on adenosine (A) binding behavior to YTH domains. A systematic interaction profile of 40 reported N6-substituted adenosine (x6A) mononucleotides with 5 human reader YTH domains was created computationally. Heuristic clustering of the profile divided these YTH domains and these x6A mononucleotides into two subfamilies and three classes, respectively; they represent distinct intrinsic interaction modes between the domains and mononucleotides. Statistical survey unraveled that the volume (Vg) and hydrophobicity (Hg) of N6-substituted chemical groups exhibit linear and nonlinear correlations with the binding energy (ΔGttl) of x6A mononucleotides to YTH domains, respectively; N6-substitutions with moderate size and weak polarity are favorable for the x6A binding. From the profile the N6-bromomethyl adenosine (brm6A) was identified as a potent binder of YTHDF2 YTH domain; its affinity was improved significantly by 77.2-fold from A and considerably by 19.5-fold from m6A. Structural modeling observed that the N6-bromomethyl group of brm6A is tightly packed against an aromatic cage defined by the Trp432-Trp486-Trp491 triad of YTHDF2 YTH domain. Electron-correlation analysis revealed that the bromine atom can form geometrically and energetically satisfactory halogen-π interactions with the aromatic cage, thus conferring considerable affinity and specificity to the domain-brm6A interaction.

On the specificity of the recognition of m6A-RNA by YTH reader domains

Most processes of life are the result of polyvalent interactions between macromolecules, often of heterogeneous types and sizes. Frequently, the times associated with these interactions are prohibitively long for interrogation using atomistic simulations. Here, we study the recognition of N6-methylated adenine (m6A) in RNA by the reader domain YTHDC1, a prototypical, cognate pair that challenges simulations through its composition and required timescales. Simulations of RNA pentanucleotides in water reveal that the unbound state can impact (un)binding kinetics in a manner that is both model- and sequence-dependent. This is important because there are two contributions to the specificity of the recognition of the Gm6AC motif: from the sequence adjacent to the central adenine and from its methylation. Next, we establish a reductionist model consisting of an RNA trinucleotide binding to the isolated reader domain in high salt. An adaptive sampling protocol allows us to quantitatively study the dissociation of this complex. Through joint analysis of a data set including both the cognate and control sequences (GAC, Am6AA, and AAA), we derive that both contributions to specificity, sequence, and methylation, are significant and in good agreement with experimental numbers. Analysis of the kinetics suggests that flexibility in both the RNA and the YTHDC1 recognition loop leads to many low-populated unbinding pathways. This multiple-pathway mechanism might be dominant for the binding of unstructured polymers, including RNA and peptides, to proteins when their association is driven by polyvalent, electrostatic interactions.

Molecular Simulations to Investigate the Impact of N6-Methylation in RNA Recognition: Improving Accuracy and Precision of Binding Free Energy Prediction

N6-Methyladenosine (m6A) is a prevalent RNA post-transcriptional modification that plays crucial roles in RNA stability, structural dynamics, and interactions with proteins. The YT521-B (YTH) family of proteins, which are notable m6A readers, functions through its highly conserved YTH domain. Recent structural investigations and molecular dynamics (MD) simulations have shed light on the mechanism of recognition of m6A by the YTHDC1 protein. Despite advancements, using MD to predict the stabilization induced by m6A on the free energy of binding between RNA and YTH proteins remains challenging due to inaccuracy of the employed force field and limited sampling. For instance, simulations often fail to sufficiently capture the hydration dynamics of the binding pocket. This study addresses these challenges through an innovative methodology that integrates metadynamics, alchemical simulations, and force-field refinement. Importantly, our research identifies hydration of the binding pocket as giving only a minor contribution to the binding free energy and emphasizes the critical importance of precisely tuning force-field parameters to experimental data. By employing a fitting strategy built on alchemical calculations, we refine the m6A partial charge parameters, thereby enabling the simultaneous reproduction of N6 methylation on both the protein binding free energy and the thermodynamic stability of nine RNA duplexes. Our findings underscore the sensitivity of binding free energies to partial charges, highlighting the necessity for thorough parametrization and validation against experimental observations across a range of structural contexts.

Structure-Based Design of a Potent and Selective YTHDC1 Ligand

N6-Adenosine methylation (m6A) is a prevalent post-transcriptional modification of mRNA, with YTHDC1 being the reader protein responsible for recognizing this modification in the cell nucleus. Here, we present a protein structure-based medicinal chemistry campaign that resulted in the YTHDC1 inhibitor 40, which shows an equilibrium dissociation constant (Kd) of 49 nM. The crystal structure of the complex (1.6 Å resolution) validated the design. Compound 40 is selective against the cytoplasmic m6A-RNA readers YTHDF1-3 and YTHDC2 and shows antiproliferative activity against the acute myeloid leukemia (AML) cell lines THP-1, MOLM-13, and NOMO-1. For the series of compounds that culminated into ligand 40, the good correlation between the affinity in the biochemical assay and antiproliferative activity in the THP-1 cell line provides evidence of YTHDC1 target engagement in the cell. The binding to YTHDC1 in the cell is further supported by the cellular thermal shift assay. Thus, ligand 40 is a tool compound for studying the role of YTHDC1 in AML.

Adenine Methylation Enhances the Conformational Flexibility of an RNA Hairpin Tetraloop

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

Methylation across the central dogma in health and diseases: new therapeutic strategies

The proper transfer of genetic information from DNA to RNA to protein is essential for cell-fate control, development, and health. Methylation of DNA, RNAs, histones, and non-histone proteins is a reversible post-synthesis modification that finetunes gene expression and function in diverse physiological processes. Aberrant methylation caused by genetic mutations or environmental stimuli promotes various diseases and accelerates aging, necessitating the development of therapies to correct the disease-driver methylation imbalance. In this Review, we summarize the operating system of methylation across the central dogma, which includes writers, erasers, readers, and reader-independent outputs. We then discuss how dysregulation of the system contributes to neurological disorders, cancer, and aging. Current small-molecule compounds that target the modifiers show modest success in certain cancers. The methylome-wide action and lack of specificity lead to undesirable biological effects and cytotoxicity, limiting their therapeutic application, especially for diseases with a monogenic cause or different directions of methylation changes. Emerging tools capable of site-specific methylation manipulation hold great promise to solve this dilemma. With the refinement of delivery vehicles, these new tools are well positioned to advance the basic research and clinical translation of the methylation field.

Structure-Based Design of Inhibitors of the m6A-RNA Writer Enzyme METTL3

Methyltransferase-like 3 (METTL3) and METTL14 form a heterodimeric complex that catalyzes the most abundant internal mRNA modification, N6-methyladenosine (m6A). METTL3 is the catalytic subunit that binds the co-substrate S-adenosyl methionine (SAM), while METTL14 is involved in mRNA binding. The m6A modification provides post-transcriptional level control over gene expression as it affects almost all stages of the mRNA life cycle, including splicing, nuclear export, translation, and decay. There is increasing evidence for an oncogenic role of METTL3 in acute myeloid leukemia. Here, we use structural and dynamic details of the catalytic subunit METTL3 for developing small-molecule inhibitors that compete with SAM. Starting from a hit identified by high-throughput docking, protein crystallography and molecular dynamics simulations were employed to guide the optimization of inhibitory activity. The potency was successfully improved by 8000-fold as measured by a homogeneous time-resolved fluorescence assay. The optimized compound is selective against the off-targets RNA methyltransferases METTL1 and METTL16.

Characterization of epitranscriptome reader proteins experimentally and in silico: Current knowledge and future perspectives beyond the YTH domain

To date, over 150 chemical modifications to the four canonical RNA bases have been discovered, known collectively as the epitranscriptome. Many of these modifications have been implicated in a variety of cellular processes and disease states. Additional work has been done to identify proteins known as "readers" that selectively interact with RNAs that contain specific chemical modifications. Protein interactomes with N6-methyladenosine (m⁶A), N1-methyladenosine (m¹A), N5-methylcytosine (m⁵C), and 8-oxo-7,8-dihydroguanosine (8-oxoG) have been determined, mainly through experimental advances in proteomics techniques. However, relatively few proteins have been confirmed to bind directly to RNA containing these modifications. Furthermore, for many of these protein readers, the exact binding mechanisms as well as the exclusivity for recognition of modified RNA species remain elusive, leading to questions regarding their roles within different cellular processes. In the case of the YT-521B homology (YTH) family of proteins, both experimental and in silico techniques have been leveraged to provide valuable biophysical insights into the mechanisms of m⁶A recognition at atomic resolution. To date, the YTH family is one of the best characterized classes of readers. Here, we review current knowledge about epitranscriptome recognition of the YTH domain proteins from previously published experimental and computational studies. We additionally outline knowledge gaps for proteins beyond the well-studied human YTH domains and the current in silico techniques and resources that can enable investigation of protein interactions with modified RNA outside of the YTH-m⁶A context.

Recent Advances of RNA m6A Modifications in Cancer Immunoediting and Immunotherapy

Cancer immunotherapy, which modulates immune responses against tumors using immune-checkpoint inhibitors or adoptive cell transfer, has emerged as a novel and promising therapy for tumors. However, only a minority of patients demonstrate durable responses, while the majority of patients are resistant to immunotherapy. The immune system can paradoxically constrain and promote tumor development and progression. This process is referred to as cancer immunoediting. The mechanisms of resistance to immunotherapy seem to be that cancer cells undergo immunoediting to evade recognition and elimination by the immune system. RNA modifications, specifically N6-methyladenosine (m6A) methylation, have emerged as a key regulator of various post-transcriptional gene regulatory processes, such as RNA export, splicing, stability, and degradation, which play unappreciated roles in various physiological and pathological processes, including immune system development and cancer pathogenesis. Therefore, a deeper understanding of the mechanisms by which RNA modifications impact the cancer immunoediting process can provide insight into the mechanisms of resistance to immunotherapies and the strategies that can be used to overcome such resistance. In this chapter, we briefly introduce the background of cancer immunoediting and immunotherapy. We also review and discuss the roles and mechanisms of RNA m6A modifications in fine-tuning the innate and adaptive immune responses, as well as in regulating tumor escape from immunosurveillance. Finally, we summarize the current strategies targeting m6A regulators for cancer immunotherapy.

Interaction preferences between protein side chains and key epigenetic modifications 5-methylcytosine, 5-hydroxymethycytosine and N6-methyladenine

Covalent modifications of standard DNA/RNA nucleobases affect epigenetic regulation of gene expression by modulating interactions between nucleic acids and protein readers. We derive here the absolute binding free energies and analyze the binding modalities between key modified nucleobases 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) and N⁶-methyladenine (m⁶A) and all non-prolyl/non-glycyl protein side chains using molecular dynamics simulations and umbrella sampling in both water and methanol, the latter mimicking the low dielectric environment at the dehydrated nucleic-acid/protein interfaces. We verify the derived affinities by comparing against a comprehensive set of high-resolution structures of nucleic-protein complexes involving 5mC. Our analysis identifies protein side chains that are highly tuned for detecting cytosine methylation as a function of the environment and can thus serve as microscopic readers of epigenetic marks. Conversely, we show that the relative ordering of sidechain affinities for 5hmC and m⁶A does not differ significantly from those for their precursor bases, cytosine and adenine, respectively, especially in the low dielectric environment. For those two modified bases, the effect is more nuanced and manifests itself primarily at the level of absolute changes in the binding free energy. Our results contribute towards establishing a quantitative foundation for understanding, predicting and modulating the interactions between modified nucleic acids and proteins at the atomistic level.

Molecular Simulations Matching Denaturation Experiments for N6-Methyladenosine

Post-transcriptional modifications are crucial for RNA function and can affect its structure and dynamics. Force-field-based classical molecular dynamics simulations are a fundamental tool to characterize biomolecular dynamics, and their application to RNA is flourishing. Here, we show that the set of force-field parameters for N⁶-methyladenosine (m⁶A) developed for the commonly used AMBER force field does not reproduce duplex denaturation experiments and, specifically, cannot be used to describe both paired and unpaired states. Then, we use reweighting techniques to derive new parameters matching available experimental data. The resulting force field can be used to properly describe paired and unpaired m⁶A in both syn and anti conformation, which thus opens the way to the use of molecular simulations to investigate the effects of N6 methylations on RNA structural dynamics.

Cancer epitranscriptomics in a nutshell

Remarkable technological progress has led to breakthrough discoveries in epitranscriptomics, reshaping our understanding of modifications decorating RNA. The past decade has seen a tremendous endeavor to describe the nature, functions, and biological roles of messenger RNA (mRNA) modifications, positioning epitranscriptomics as a crucial pillar in tumor biology. Like DNA and histone modifications, mRNA marks have been increasingly linked to cancer pathogenesis. Here, we summarize the latest research in cancer epitranscriptomics with emphasis on N6-methyladenosine, untangling its contribution to five prime oncogenic features: tumor growth, activating invasion and metastasis, stemness, metabolic reprogramming, and tumor microenvironment. We discuss mRNA-modifying enzymes, their impact on biological processes, and contribution to cancer hallmarks. We spotlight epitranscriptomics as a promising bonanza for forthcoming targeting approaches in cancer therapy.

RNA N6-Methyladenosine Modifications and Its Roles in Alzheimer's Disease

The importance of epitranscriptomics in regulating gene expression has received widespread attention. Recently, RNA methylation modifications, particularly N6-methyladenosine (m6A), have received marked attention. m6A, the most common and abundant type of eukaryotic methylation modification in RNAs, is a dynamic reversible modification that regulates nuclear splicing, stability, translation, and subcellular localization of RNAs. These processes are involved in the occurrence and development of many diseases. An increasing number of studies have focused on the role of m6A modification in Alzheimer's disease, which is the most common neurodegenerative disease. This review focuses on the general features, mechanisms, and functions of m6A methylation modification and its role in Alzheimer's disease.

Mutant NPM1-Regulated FTO-Mediated m6A Demethylation Promotes Leukemic Cell Survival via PDGFRB/ERK Signaling Axis

Acute myeloid leukemia (AML) with nucleophosmin 1 (NPM1) mutations exhibits distinct biological and clinical features, accounting for approximately one-third of AML. Recently, the N⁶-methyladenosine (m⁶A) RNA modification has emerged as a new epigenetic modification to contribute to tumorigenesis and development. However, there is limited knowledge on the role of m⁶A modifications in NPM1-mutated AML. In this study, the decreased m⁶A level was first detected and high expression of fat mass and obesity-associated protein (FTO) was responsible for the m⁶A suppression in NPM1-mutated AML. FTO upregulation was partially induced by NPM1 mutation type A (NPM1-mA) through impeding the proteasome pathway. Importantly, FTO promoted leukemic cell survival by facilitating cell cycle and inhibiting cell apoptosis. Mechanistic investigations demonstrated that FTO depended on its m⁶A RNA demethylase activity to activate PDGFRB/ERK signaling axis. Our findings indicate that FTO-mediated m⁶A demethylation plays an oncogenic role in NPM1-mutated AML and provide a new layer of epigenetic insight for future treatments of this distinctly leukemic entity.

Structural effects of m6A modification of the Xist A-repeat AUCG tetraloop and its recognition by YTHDC1

The A-repeat region of the lncRNA Xist is critical for X inactivation and harbors several N⁶-methyladenosine (m⁶A) modifications. How the m⁶A modification affects the conformation of the conserved AUCG tetraloop hairpin of the A-repeats and how it can be recognized by the YTHDC1 reader protein is unknown. Here, we report the NMR solution structure of the (m⁶A)UCG hairpin, which reveals that the m⁶A base extends 5′ stacking of the A-form helical stem, resembling the unmethylated AUCG tetraloop. A crystal structure of YTHDC1 bound to the (m⁶A)UCG tetraloop shows that the (m⁶A)UC nucleotides are recognized by the YTH domain of YTHDC1 in a single-stranded conformation. The m⁶A base inserts into the aromatic cage and the U and C bases interact with a flanking charged surface region, resembling the recognition of single-stranded m⁶A RNA ligands. Notably, NMR and fluorescence quenching experiments show that the binding requires local unfolding of the upper stem region of the (m⁶A)UCG hairpin. Our data show that m⁶A can be readily accommodated in hairpin loop regions, but recognition by YTH readers requires local unfolding of flanking stem regions. This suggests how m⁶A modifications may regulate lncRNA function by modulating RNA structure.

Recognition of N6-Methyladenosine by the YTHDC1 YTH Domain Studied by Molecular Dynamics and NMR Spectroscopy: The Role of Hydration

The YTH domain of YTHDC1 belongs to a class of protein "readers", recognizing the N6-methyladenosine (m⁶A) chemical modification in mRNA. Static ensemble-averaged structures revealed details of N6-methyl recognition via a conserved aromatic cage. Here, we performed molecular dynamics (MD) simulations along with nuclear magnetic resonance (NMR) and isothermal titration calorimetry (ITC) to examine how dynamics and solvent interactions contribute to the m⁶A recognition and negative selectivity toward an unmethylated substrate. The structured water molecules surrounding the bound RNA and the methylated substrate's ability to exclude bulk water molecules contribute to the YTH domain's preference for m⁶A. Intrusions of bulk water deep into the binding pocket disrupt binding of unmethylated adenosine. The YTHDC1's preference for the 5'-Gm⁶A-3' motif is partially facilitated by a network of water-mediated interactions between the 2-amino group of the guanosine and residues in the m⁶A binding pocket. The 5'-Im⁶A-3' (where I is inosine) motif can be recognized too, but disruption of the water network lowers affinity. The D479A mutant also disrupts the water network and destabilizes m⁶A binding. Our interdisciplinary study of the YTHDC1 protein-RNA complex reveals an unusual physical mechanism by which solvent interactions contribute toward m⁶A recognition.