RNA binding to human METTL3-METTL14 restricts N6-deoxyadenosine methylation of DNA in vitro
Abstract
Methyltransferase like-3 (METTL3) and METTL14 complex transfers a methyl group from S-adenosyl-L-methionine to N6 amino group of adenosine bases in RNA (m6A) and DNA (m6dA). Emerging evidence highlights a role of METTL3-METTL14 in the chromatin context, especially in processes where DNA and RNA are held in close proximity. However, a mechanistic framework about specificity for substrate RNA/DNA and their interrelationship remain unclear. By systematically studying methylation activity and binding affinity to a number of DNA and RNA oligos with different propensities to form inter- or intra-molecular duplexes or single-stranded molecules in vitro, we uncover an inverse relationship for substrate binding and methylation and show that METTL3-METTL14 preferentially catalyzes the formation of m6dA in single-stranded DNA (ssDNA), despite weaker binding affinity to DNA. In contrast, it binds structured RNAs with high affinity, but methylates the target adenosine in RNA (m6A) much less efficiently than it does in ssDNA. We also show that METTL3-METTL14-mediated methylation of DNA is largely restricted by structured RNA elements prevalent in long noncoding and other cellular RNAs.
Editor's evaluation
This manuscript will be of interest to researchers in the fields of nucleic acid chemical biology in general and diseases related to nucleic acid methylation in particular. The data presented support the conclusions of the paper within the current context, provide new evidence and plausible explanations to previously inexplicable mechanisms.
https://doi.org/10.7554/eLife.67150.sa0Introduction
N6-methyladenosine (m6A) is considered a major covalent modification of the adenosine (A) base in coding and noncoding (nc) RNAs (Desrosiers et al., 1974; Bokar et al., 1994; Rottman et al., 1994; Dominissini et al., 2012; Meyer et al., 2012). It is linked to diverse physiological processes, including – but not limited to – RNA turnover (Ke et al., 2017; Knuckles et al., 2017; Liu et al., 2020), stem cell differentiation (Batista et al., 2014; Geula et al., 2015,) oncogenic translation (Lin et al., 2016; Barbieri et al., 2017; Vu et al., 2017; Choe et al., 2018), and DNA damage repair (Xiang et al., 2017; Zhang et al., 2020). In human cells, most m6A in cellular RNAs is installed by methyltransferase like-3 (METTL3) and METTL14 methyltransferases (Bokar et al., 1994; Dominissini et al., 2012; Liu et al., 2014), both of which belong to the β-class of S-adenosyl methionine (SAM)-dependent methyltransferases (N6-MTases) (Bujnicki et al., 2002; Woodcock et al., 2020) that also include N6-deoxyadenosine DNA methyltransferases (m6dA MTases), especially those from Type III restriction-modification (R-M) systems in bacteria (e.g., ModA/B dimer of EcoP15I) (Gupta et al., 2015). The first structures of METTL3-METTL14 complexes suggested common features in EcoP15I and METTL3-METTL14 (Wang et al., 2016b; Sledz and Jinek, 2016). We analyzed these MTases and found high structural similarity of the cores of human METTL3, METTL14, and EcoP15I, suggesting their evolutionary origin from a common ancestor (Figure 1a–e and Figure 1—figure supplement 1).
While m6dA in bacterial DNA is deposited in a strictly sequence-dependent manner in double-stranded DNA (dsDNA) (Gupta et al., 2015; Malone et al., 1995), m6A in RNA by human METTL3-METTL14 shows less stringent sequence dependency for sequences flanking the target A within the recognition motif DRACH (D = A/G/U, R = A/G, H = A/C/U) (Bokar et al., 1994; Liu et al., 2014; Linder et al., 2015). Such relaxed specificity maybe required to expand the repertoire of target transcripts of METTL3-METTL14. However, this does not correlate with low levels of m6A in cellular RNA, which remains at ~0.1% of total pool of ribonucleotides in cultured mammalian cells (Dubin and Taylor, 1975). Of note, recent evidence suggests that all m6A is deposited co-transcriptionally, with most occurring on chromatin-associated RNAs, for example, nascent pre-mRNA (CA-RNA), promoter-associated RNA, enhancer RNA, and repeat RNA elements (Ke et al., 2017; Knuckles et al., 2017; Liu et al., 2020; Slobodin et al., 2017). Consistently, the METTL3-METTL14 m6A writers can be recruited to chromatin by modified histone tails (e.g., H3K36me3) (Huang et al., 2019b) or the transcription machinery (Barbieri et al., 2017; Slobodin et al., 2017). Certain long noncoding (lnc) RNAs that are essential for regulated transcription, such as chromatin-associated XIST (Patil et al., 2016) and NEAT2 (MALAT1) (Liu et al., 2015), also harbor m6A marks. A role of m6A in DNA damage response and repair has also been suggested (Zhang et al., 2020; Xiang et al., 2017). In those contexts, METTL3 localizes at sites of dsDNA breaks and installs the m6A on DNA-damage-associated RNAs (Zhang et al., 2020). METTL3 is also emerged as a pro-viral host factor in SARS-CoV-2 pathogenesis. Consistently, genetic depletion or inhibition of METTL3 by small molecules triggers innate immune response and limits viral growth (Burgess et al., 2021; Li et al., 2021).
The increasingly appreciated roles of m6A in nucleic acid metabolism thus seem to converge to a common point: the occurrence of m6A at chromatin, a predominantly DNA- and RNA-rich environment inside the nucleus. A recent study showed that METTL3-METTL14 can convert dA to m6dA within equivalent DRACH (D = A/G/U, R = A/G, H = A/C/U) motif in a single-stranded DNA (ssDNA) with a 2.5-fold preference over the equivalent RNA in vitro (Woodcock et al., 2019). Since m6dA levels in mammalian DNA are extremely low (0.006–0.01%) and thought to originate from RNA catabolism rather than from SAM-dependent MTase reactions (Musheev et al., 2020), fundamental questions about substrate(s) specificity and mechanisms of action of METTL3-METTL14 have emerged. By using classical biochemical and mass spectrometry assays, we show that METTL3-METTL14 preferentially methylates N6dA on ssDNAs to m6dA, and this activity is largely regulated by structured RNA elements prevalent in lnc RNAs, that are also found in other cellular RNAs. Thus, our work provides a framework to explore a new regulatory axis of RNA-mediated restriction of N6-deoxyadenosine (m6dA) methylation in mammalian genomes.
Results
We compared the MTase cores of (ModB) of EcoP15I (aa 90–132, 169–261, 385–511, PDB: 4ZCF Gupta et al., 2015) and METTL14 (aa 116–402, PDB: 5IL0 Wang et al., 2016b). A secondary structure-based superposition (Krissinel and Henrick, 2004) of these two structures revealed common features and similar arrangement of canonical MTase motifs (motifs I and IV-X), including the motif IV (D/EPPY/W/L) that surrounds the Watson-Crick edge of the flipped target adenine base (Figure 1a–e and Figure 1—figure supplement 1). Previous studies identified three loops encompassing the sequence intervening MTase motifs IV and V (loop 1), VIII and VIII’ (loop 2), and motifs IX and X (loop 3) contribute to substrate binding (Wang et al., 2016b; Sledz and Jinek, 2016; Wang et al., 2016a). Due to lack of experimental structures of METTL3 or ModA (in complex with a flipped adenine base), we chose ModB and METTL14 for this comparison. While the MTase domain of METTL14 is enzymatically inactive, it harbors features that contribute to substrate binding (Wang et al., 2016b; Sledz and Jinek, 2016; Wang et al., 2016a). These observations prompted us to address the prevailing question about substrate specificity of the human full-length METTL3-METTL14 enzyme complex. We co-purified the full-length human METTL3-METTL14 complex and its truncated form devoid of the RGG motif of METTL14 (METL3-METL14-RGG) from insect cells (see details in Materials and methods section) (Figure 1f). The RGG motifs are clustered sequences of arginines and glycines. These motifs are commonly present in diverse set of RNA-binding proteins that play important roles in physiological processes, such as RNA synthesis and processing (Thandapani et al., 2013). In human METTL14, a total of seven RGG triplets and two RG motifs are present at its C-terminus tail (aa 408–457). The region also harbors a number of aromatic amino acids, which can further stabilize the bound nucleic acids via hydrophobic or stacking interactions. Consistently, this region contributes to substrate binding and m6A activity of the METTL3-METTL14 enzyme complex (Schöller et al., 2018). The sequence encompassing the RGG motifs in METTL14 is well conserved in higher vertebrates (Liu et al., 2014).
Next, we designed RNA substrates of varied lengths (5–30 nucleotides) with at least one canonical m6A signature motif (RRACH), including an RNA oligo (rPal-Top+ bat) in which two RRACH sites were arranged in a palindromic fashion (Figure 1g). We included a 14-mer DNA oligo (d6T) according to Woodcock et al., 2019, that reported that METTL3-METTL14 exhibited highest methyltransferase (MTase) activity on d6T. We also included two structured RNAs with propensity to form a perfect stem-loop (23-mer rTCE23 or also known as SRE; Oberstrass et al., 2006) without a DRACH motif, and a 30-mer rNEAT2 (or MALAT1) encompassing one DRACH motif and predicted to form a stem-loop with a bulged stem (Figure 1h–i).
The m6A mark in mammalian RNAs is widespread and expands beyond internal adenines in mRNAs to lncRNAs (XIST Patil et al., 2016, NEAT2 or MALAT1 Liu et al., 2015), pri-miRNAs (Alarcón et al., 2015), and snoRNAs (Linder et al., 2015) – all of which are crucial to gene regulation and maintenance of cellular homeostasis. Of note, NEAT2 lncRNA can relocate a large subset of transcription units across cellular compartments from polycomb bodies to interchromatin granules (Yang et al., 2011). According to this model, a 49-nt long fragment of NEAT2 (nt 4917–4966) facilitates the recruitment of a mega transcription complex (Yang et al., 2011). We observed that this NEAT2 fragment harbors one RRACH motif (AAACA) in its loop region and fulfils the criteria of m6A peak signature (RRACH Bokar et al., 1994; Liu et al., 2014 or DRACH Linder et al., 2015 or ACA Garcia-Campos et al., 2019; Figure 1g). Thus, we decided to test the activity of METTL3-METTL14 on this regulatory module of the lncRNA. For simplicity of synthesis and purification, we focused on a slightly shorter version (30-nt; rNEAT2-30) of NEAT2-49 in which both the m6A signature motif and bulged stem regions were well preserved (Figure 1g–h). The rTCE23 could serve as a control – it forms a perfect stem-loop structure in solution – confirmed by solution NMR previously (Oberstrass et al., 2006) – but lacks an m6A signature motif (Figure 1g and i).
We observed higher enzymatic activity of METTL3-METTL14 on the d6T DNA oligo compared with its RNA counterpart (r6T) (Figure 2a and c). By closely examining the sequence of the 14-mer d6T/r6T, we found a propensity of these oligos to form duplex DNA (for d6T) or RNA (for r6T): 6 of the 14 nucleotides at 5’-end will pair with a complementary sequence to create a duplex with an 8-nt overhang at the 3’-end of each strand (Figures 1g and 2d). As a result of this self-annealing, two adenine bases are available for methylation (m6dA) at the end of the short palindromic sequence (5’-CCG•CGG-3’). We thus hypothesized that d6T (or 6T, as referred to by Woodcock et al., 2019) may not truly represent an ssDNA substrate. Thus, we designed a new DNA oligo (d6T*) in which two nucleotides, a deoxycytidine and a deoxyguanosine at –3 and –2 positions (relative to the signature motif RRACH), respectively, were replaced by two thymidines to disrupt base pairing and palindrome formation without affecting the core RRACH motif required for m6dA installation (Figure 1g). We reasoned that if METTL3-METTL14 indeed prefers a ssDNA for methylation, it should show higher activity for d6T*, a true ssDNA substrate, over d6T. Indeed, we found that the methyltransferase activity of METTL3-METTL14 toward d6T* increased by 2-fold compared with d6T (and >12-fold compared with r6T) (Figure 2a). Consistent with previous observations (Woodcock et al., 2019), METTL3-METTL14 showed no activity in the presence of a perfect duplex version of d6T DNA (d6T-ds). These results confirm a robust activity of METTL3-METTL14 on ssDNA substrates.
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/67150/elife-67150-fig2-data1-v3.xlsx
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/67150/elife-67150-fig2-data2-v3.xlsx
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/67150/elife-67150-fig2-data3-v3.xlsx
-
Figure 2—source data 4
- https://cdn.elifesciences.org/articles/67150/elife-67150-fig2-data4-v3.xlsx
-
Figure 2—source data 5
- https://cdn.elifesciences.org/articles/67150/elife-67150-fig2-data5-v3.xlsx
We observed no methyltransferase activity on rTCE23 RNA, as expected; but surprisingly, some low (but consistently detectable) methyltransferase activity on rNEAT2 (Figure 2a). We then asked what attributes (sequence, length, and shape) in an RNA substrate are critical for METTL3-METTL14 to yield a high methyltransferase activity. To answer this question, we designed eight RNA oligos with varied length (5- to 22-nt) and sequences (Figure 1g). We also varied the position of the RRACH motif and the flanking sequences in these oligos. All oligos comprised one RRACH motif, except the RNA duplex Pal-Top+ bot, where two RRACH sites were arranged in a palindromic fashion within each 10-mer strand. We designed another variant of probe rC2-flip (22-mer, RRACH in the center), rC-flip wherein the RRACH motif resided at the 3’-end. The other three RNA oligos, FoxM1-p1 (13-nt, RRACH at 3’-end), c-Myc-p5 (8-nt), and c-Myc-p6 (11-nt, RRACH at 5’-end), were derived from the FOXM1 and MYC genes because of their involvement in m6A-mediated processes in glioblastoma and acute myeloid leukemia (Vu et al., 2017; Zhang et al., 2017). We covalently attached a fluorescein moiety at the 5’-end of all these oligos to facilitate the quantitative measurement of their binding affinities (equilibrium dissociation constant or Kd) to METTL3-METTL14 (see Materials and methods section for details).
We used a radiometric assay and fluorescence polarization (FP) to determine methylation activity and RNA binding, respectively. As shown in Figure 2a and c, RNA oligos with high propensity toward forming secondary structures (rNEAT2, rTCE23, Pal-Top+ bot) showed low to negligible levels of methylation by METTL3-METTL14. The other RNA oligos showed some levels of methylation, though they could not be methylated with as much efficiency as the ssDNA substrate, d6T*. The results of binding experiments are shown in Figure 2b and Table 1. The full-length METTL3-METTL14 binds to rNEAT2 with the highest affinity (Kd = 2 nM), followed by the hairpin rTCE23 (Kd = 25 nM), the 3’-overhang r6T (Kd = 114 nM), the palindromic Pal-Top+ bot (Kd = 127 nM), and the ssRNA FoxM1-p1 (Kd = 294 nM). Surprisingly, the DNA substrates (d6T and d6T*) that showed highest methyltransferase activity were the poorest binders, with Kd values of 370 and 509 nM, respectively.
If METTL3-METTL14 were to methylate ssDNA in vivo, and the evidence to date does not support its direct enzymatic origin (Musheev et al., 2020), then what could be the significance of high affinity RNA binding to structured RNA oligos (rNEAT2 and rTCE23)? One explanation is that structural motifs in ncRNA or mRNAs can serve as a recruitment platform for METTL3-METTL14 in a specific context. But how would the RNA-binding activity then affect DNA methylation activity? To answer this question, we measured the methyltransferase activity of METTL3-METTL14 in the presence of roughly equimolar mixtures of d6T* DNA with each RNA oligonucleotide. The methyltransferase activity of METTL3-METTL14 was significantly reduced in the presence of RNA (Figure 2c). Of note, RNAs with propensity to form secondary structures and highest binding affinity to METTL3-METTL14 (rNEAT2, rTCE23) attenuated methyltransferase activity to the most negligible levels (Figure 2c). To further explore the relationship between methyltransferase activity and substrate-binding affinity, we synthesized a variant of the stem-loop oligo (rNEAT2*) with the DRACH motif (5’-AAACA-3’) changed to an ideal m6A target sequence, that is, 5’-GGACU-3’ (Figure 1g). Although higher methylation levels were observed for rNEAT2* compared to rNEAT2, most likely attributed to a perfect GGACU motif in rNEAT2*, yet the methylation was significantly lower than that of d6T* DNA (last panel of Figure 2c). Moreover, rNEAT2* showed ~6 fold reduction in affinity compared to rNEAT2 (Figure 2b). Importantly, rNEAT2* also attenuated methylation of the single-stranded d6T* DNA by METTL3-METL14 but to a lesser extent than the rNEAT2 (Figure 2c). These results suggests that: (a) the sequence and shape of the target RNAs dictate the m6A activity of METTL3-METTL14, and (b) that higher affinity substrates tend to show lower methylation levels, possibly due to slow off rates of m6A RNA.
To further investigate RNA-mediated inhibition, we showed that the two structured RNAs (rTCE23 and rNEAT2) inhibited the DNA methylation activity of METTL3-METTL14 in a dose-dependent manner, with rNEAT2 (IC50 = 216 nM) showing 4-fold stronger inhibition than rTCE23 (IC50 = 887 nM) (Figure 2e). Next, we employed an LC/MS-based oligonucleotide intact mass analysis to further validate the results of the radiometric assay. As shown in Figure 2f, an equimolar (with respect to d6T* DNA) addition of rNEAT2 or rTCE23 significantly reduced the methyltransferase activity from 80.5–14.5% to 7.3%, respectively. We then performed nucleoside analysis of the METTL3-METTL14 reactions to unambiguously determine the identity of methylated base and efficiency of methylation (Figure 2g). We confirmed that deoxyadenosine on d6T* and adenosine on rNEAT2 were the only modified bases. In the absence of rNEAT2 or rTCE23, 86.3% of dAs on d6T* were modified. When rNEAT2 or rTCE23 were present, only 14.4% and 7.4%, respectively, of dAs on d6T* were modified. Additionally, 5.8% of the rAs on rNEAT2 were modified. As expected, no methylation on rTCE23 was detected (no RRACH motif). These results are consistent with those obtained by oligonucleotide intact mass analysis. Altogether, the methyltransferase and binding assays confirm that the binding to structured RNAs, especially those lacking the GGACU sequence, almost completely abolishes the methyltransferase activity METTL3-METTL14.
The C-terminal RGG repeat motif of METTL14 contributes to RNA binding and activity of METTL3-METTL14 (Schöller et al., 2018; Figure 3a). Thus, we hypothesized that the absence of RGG motif in METTL3-METTL14 should diminish its ability to bind RNA and consequently attenuate the effect of RNA on DNA methylation. In fact, we observed a 2- to 10-fold decrease in binding affinity by METTL3-METTL14-RGG (Figure 3b and Table 1). This suggests that the RGG motifs play a major role in RNA and DNA substrates, while other parts of the enzyme (e.g., CCCH-type zinc finger domains in METTL3 Huang et al., 2019a and MTase core of METTL14 Wang et al., 2016a) may also contribute to the overall binding (Figure 3b). As expected, the deletion of RGG caused 60% reduction in DNA methylation activity. In the presence of rNEAT2, the activity of METTL3-METTL14-RGG on d6T* DNA was reduced to about 75% of its activity in absence of the RNA; comparatively, a >90% activity reduction was observed for the full-length METTL3-METTL14 in presence of the same RNA (Figure 3c). These results suggest that the RGG motif, as a general RNA-binding domain, dictates the shape-specific RNA recognition and promotes an RNA-mediated restriction of the DNA methylation activity of METTL3-METTL14 (Table 1).
Discussion
Our results are noteworthy given the relevance of METTL3-METTL14 in regions where RNA and DNA are held in closed proximity, as supported by recent studies suggesting: (a) localization of METTL3 to sites of dsDNA breaks, and its role in DNA/RNA hybrid accumulation (Zhang et al., 2020), (b) recruitment of METTL3 to chromatin by promoter-bound transcription factors (Barbieri et al., 2017), (c) UV-induced DNA damage signaling (Xiang et al., 2017), and (d) metabolism of R-loops (Abakir et al., 2020). Earlier studies suggest the existence of a large complex of multiple subunits of ~1 MDa size, possessing RNA N6-adenosine methylation activity distributed over two sub-complexes of approximate molecular mass of ~200 and ~800 kDa (Bokar et al., 1994Bokar et al., 1997). One such subunit is WTAP, which alters both the pattern and location of the m6A transcriptome (Schwartz et al., 2014). Future studies could examine the role of WTAP and other partners of METTL3-METTL14 in substrate recognition and specificity, especially within the context of the regulatory role of RNA in DNA methylation.
In summary, we unveiled a new N6-deoxyadenosine methylation activity on ssDNA, which is regulated by binding of METTL3-METTL14 to shape-dependent structured regions of ncRNAs and chromatin-associated nascent pre-mRNAs (CA-RNA) (Figure 3d and e). The structured elements in RNA could facilitate the recruitment of METTL3-METTL14 to chromatin and/or keep the writer complex engaged during co-transcriptional occurrence of m6A. Such a topological arrangement of METTL3-METTL14/RNA interaction could also restrict m6dA deposition, especially to the single-stranded regions of DNA exposed during transcription (Slobodin et al., 2017), DNA recombination (Zhang et al., 2020), damage and repair (Xiang et al., 2017), and R-loop metabolism (Abakir et al., 2020).
A major difference in prokaryotic and eukaryotic genomes is the contrasting occurrence of the base methylation. The m6dA is a predominant modification mark in bacterial DNA, whereas the m5dC is more prevalent in eukaryotic genomes. How is the eukaryotic genome then protected from m6dA deposition despite having an efficient enzyme machinery (e.g., METTL3-METTL14) to accomplish such task? A comparative sequence/structural analysis of the bacterial EcoP15I and human METTL3 and METTL14 MTases provides some hints (Figures 1a–c–3a). Of note, there are two major differences that exist in the sequences of METTL3 and METTL14: (a) loss of the target recognition domain (TRD) and (b) acquisition of the RGG motif in METTL14. The TRD motif (in EcoP15I MTase) is a major contributor to binding and specificity for a dsDNA (Gupta et al., 2015), whereas the RGG motif is known for RNA binding (Thandapani et al., 2013). It is likely that gene fusion/genome arrangement event(s) could have occurred during evolution (Bujnicki et al., 2002; Woodcock et al., 2020) to ensure very low/no occurrence (Douvlataniotis et al., 2020) of m6dA in mammalian genome. The high affinity to structured RNAs as demonstrated here could also ensure efficient recruitment of METTL3-METTL14 to specific nuclear compartment(s) in human cells.
Interestingly, a recent study showed that a comprehensive network of RNA-RNA and RNA-DNA interactions mediated by ncRNAs (including MALAT1) creates distinct nuclear compartments and facilitates recruitment of proteins for RNA processing, heterochromatin assembly, and gene regulation (Quinodoz et al., 2020). Thus, the secondary structure and 3D shape of the RNAs appears to be crucial for these important biological processes. In line with these observations, our results suggest that structured elements in RNAs may play an important role in regulating the ssDNA and ssRNA methylation activity of METTL3-METTL14.
Materials and methods
Protein expression, purification, and oligonucleotide preparation
Request a detailed protocolThe coding sequence of full-length human METTL3 (NCBI reference sequence GI: 33301371) was used in this study. The gene was cloned into a plasmid suitable for expression of proteins in insect cells with an N-terminal poly histidine tag followed by a tobacco etch virus (TEV) protease site. This plasmid (5TEY (METTL3)) was a kind gift from Dr Cheryl Arrowsmith (Addgene plasmid # 101892; http://n2t.net/addgene:101892; RRID:Addgene_101892). The coding sequence of full-length human METTL14 (NCBI reference number GI: 172045930) was cloned into the pFastBac1 vector between BamHI and NotI restriction enzyme sites. For expression of recombinant proteins, we used ExpiSf Expression System (Thermo Fisher). First, the viral DNA bacmids for each gene were prepared from individual plasmids transformed into MAX Efficiency DH10Bac competent cells (Thermo Fisher). We confirmed the identity of the inserted genes by PCR amplification and DNA sequencing. The recombinant P0 virus for each gene were generated in ExpiSf9 insect cells following the manufacturer’s recommendations (Thermo Fisher). The amount of each virus needed for infection has been optimized for maximal production of the two proteins. Infected cells were grown in ExpiSf CD medium (Thermo Fisher) at 27°C on an orbital shaker (speed 125 rpm) in a controlled (non-humidified, air-regulated) environment. Cells were harvested at 72 hr post-infection by spinning at 300× g for 5 min, washed with phosphate buffered saline (PBS) twice, and resuspended in a buffer A (0.025 M Tris-HCl pH 8.0, 0.5 M NaCl, 0.005 M imidazole, 0.1 mM TCEP, 10% [(v/v)] glycerol) supplemented with 0.5% (v/v) Igepal, two mini-protease inhibitor cocktail tablets (Roche), DNase I, and lysozyme. Cells were lysed using a microfluidizer (Analytik, UK) and the soluble fraction was separated by centrifuging the lysate at 40,000 rpm for 40 min. The supernatant was passed through a 0.22 µm filter.
The clarified soluble fraction was loaded on to a Nuvia IMAC column (Bio-Rad) pre-equilibrated in buffer A. The proteins were eluted by increasing the concentration of imidazole. The 6XHis tag was then proteolytically removed using TEV protease. Any uncleaved fractions and remaining protein impurities were removed by second passage through an IMAC column. The METTL3-METTL14 complex was further purified by successive passage through HiTrap heparin and HiLoad Superdex 200 columns (GE Healthcare). The protein complex was eluted as a single homogenous species in a final buffer containing 0.02 M Tris-HCl pH 8.0, 0.15 M NaCl, 0.1 mM TCEP. The purified protein complex was concentrated to 3–5 mg/mL and used immediately, or flash-frozen in liquid nitrogen and stored at –80°C. The METTL3-METT14-RGG complex lacks the C-terminus RGG repeats motif (amino acids 400–457) of METTL14. It was expressed and purified using same method as the full-length METTL3-METTL14 complex.
All RNA and DNA oligonucleotides used in this study were synthesized, purified by HPLC, and received in deprotected (for RNAs) and desalted form from Dharmacon and Integrated DNA Technologies, respectively. Oligos were dissolved in 1× buffer containing 0.01 M Tris-HCl pH 8.0, 0.05 M NaCl. To prepare double-stranded oligos, the two strands were mixed in roughly equimolar amounts. Annealing was carried out by heating the mixture to 95°C for 2 min followed by gradual cooling to 25°C over 45 min in a thermocycler (Eppendorf).
Methyltransferase assays
Radiometric assay
Request a detailed protocolEach reaction was carried out in a 5 µL mixture containing 50 mM HEPES pH 7.5, 5 mM NaCl, 1 mM dithiothreitol (DTT), 5 µM [methyl-3H] SAM (PerkinElmer), 10 µM substrate RNA/DNA probe, and 1 µM co-purified full-length METTL3-METTL14 enzyme complex. The reactions were incubated at 37°C for 1 hr and 3 µL of each reaction was quenched by blotting on Hybond-N+ membrane (Amersham). RNA/DNA probes were crosslinked to membrane by exposure to ultraviolet light (254 nm) for 2 min. The membranes were successively washed three times by 1× PBS followed by three ethanol washes for 5 min each. The membranes were air-dried in the hood for 15 min and the counts per minute (c.p.m.) of the RNA/DNA probes were measured using a scintillation counter (Beckman LS6500). We also performed this assay in a reaction buffer containing higher salt (50.0 mM NaCl) and observed consistent results (dose-dependent inhibition of m6dA activity of METTL3-METTL14 by small, structured RNA oligos) as presented in Figure 2e. Results presented here are an average of three independent experiments (n = 3) with one standard deviation (s.d.) for each oligonucleotide shown as error bars. Source data are provided as a Source Data file. The results of two groups in Figure 3c were derived from three independent experiments and were analyzed using two-tailed Student’s t-test (unpaired).
Quantitation of methylation and validation of identity of methylated base by mass spectrometry
Request a detailed protocolMethylation reactions were carried out in a 10 μL mixture containing 50 mM Tris-HCl, pH 7.5, 5 mM NaCl, 1 mM DTT, 200 μM SAM, 10 μM substrate RNA/DNA probe, and 1 μM METTL3-METTL14 enzyme complex. The reactions were incubated at 37°C for 1 hr, followed by incubation with 0.8 units of proteinase K (New England Biolabs) at 37°C for 10 min. Two microliters of each reaction were subjected to oligonucleotide intact mass analysis by liquid chromatography-mass spectrometry (LC-MS) as described previously (Viswanathan et al., 2020). The raw data was deconvoluted using ProMass HR (Novatia, LCC). The deconvoluted-mass peak area ratio between reactants and expected products was used to estimate the percentage of methylation. Results from three independent experiments (n = 3) are shown. To verify the identity of the modified nucleotides, METTL3-METTL14 reactions were also subjected to nucleoside analysis. Briefly, 8 μL of each of reaction were purified using Oligo Clean-Up and Concentration Kit (Norgen Biotek) according to manufacturer’s instructions. The nucleic acids were eluted in 20 μL nuclease-free water. To seventeen microliters of the eluates were added 2 μL of 10× Nucleoside Digestion Mix reaction buffer and 1 μL of Nucleoside Digestion Mix (New England Biolabs). Reactions were incubated at 37°C for 1 hr. The entirety of the nucleoside digestion reaction was subjected to ultra-high-performance liquid chromatography (UHPLC) or UHPLC-MS analysis without purification. UHPLC analysis was performed on an Agilent 1290 Infinity II system equipped with a G7117A diode array detector (DAD). UHPLC-MS analysis was performed an Agilent 1290 Infinity II system equipped with a G7117A DAD and a 6135 XT mass detector. LC separations were carried out using a Waters XSelect HSS T3 XP column (2.1 × 100 mm, 2.5 µm) with the gradient mobile phase consisting of methanol and 10 mM ammonium acetate buffer (pH 4.5). The identity of each nucleoside was confirmed by analysis of the LC retention times relative to authentic standards and by mass spectrometric analysis. The relative abundance of each nucleoside was determined by UV absorbance.
Determining the nucleic acid-binding affinities of METTL3-METTL14 and METTL3-METTL14-RGG
Request a detailed protocolFP-based assay was used to measure equilibrium dissociation-binding constants (Kd). The reactions were carried out in a buffer 0.01 M HEPES pH 7.5 and 0.05 M KCl. All oligonucleotides contain a fluorescein moiety covalently attached to the 5’-end. A constant concentration of the fluorescein labeled oligo (1 nM) was used with increasing concentrations of full-length METTL3-METTL14 protein or its truncated version (METTL3-METTL14-RGG (0–2000 nM)) in a 384-well plate. Significant changes in FP upon increasing concentrations of protein were indicative of direct binding. The FP (emission wavelength = 530 nm, excitation wavelength = 485 nm) value for each dilution was measured using PHERAstar FS (BMG Labtech). The buffer corrected values were used to calculate the equilibrium dissociation constant (Kd) for protein-DNA/RNA binding using a simple 1:1 specific binding model (Y = Bmax*X/(Kd+ X), where X denotes the concentration of protein ligand, Y = specific binding, Bmax = maximum binding, Kd = equilibrium dissociation constant in same units as X). Data were fitted by a single-site binding model in GraphPad Prism (GraphPad Software, San Diego, CA). Results presented here are an average of three independent experiments (n = 3) with one standard deviation (s.d.) for each oligonucleotide shown as error bars. Source data are provided as a Source Data file.
Data availability
The information about coding sequences of human METTL3 (NCBI reference sequence GI: 33301371) and METTL14 (NCBI reference sequence GI: 172045930) used in this study is available at NCBI. Source data are provided as a separate Source Data file. Correspondence and requests for material should be addressed to Y.K.G. (guptay@uthscsa.edu).
References
-
Characterization and partial purification of mRNA N6-adenosine methyltransferase from HeLa cell nucleiThe Journal of Biological Chemistry 269:17697–17704.
-
Purification and cDNA cloning of the AdoMet-binding subunit of the human mRNA (N6-adenosine)-methyltransferaseRNA 3:1233–1247.
-
Targeting the m6A RNA modification pathway blocks SARS-CoV-2 and HCoV-OC43 replicationGenes & Development 35:1005–1019.https://doi.org/10.1101/gad.348320.121
-
No evidence for DNA N (6)-methyladenine in mammalsScience Advances 6:eaay3335.https://doi.org/10.1126/sciadv.aay3335
-
The methylation state of poly A-containing messenger RNA from cultured hamster cellsNucleic Acids Research 2:1653–1668.https://doi.org/10.1093/nar/2.10.1653
-
Stem cells m6A mRNA methylation facilitates resolution of naïve pluripotency toward differentiationScience (New York, N.Y.) 347:1002–1006.https://doi.org/10.1126/science.1261417
-
RNA fate determination through cotranscriptional adenosine methylation and microprocessor bindingNature Structural & Molecular Biology 24:561–569.https://doi.org/10.1038/nsmb.3419
-
Secondary-structure matching (SSM), a new tool for fast protein structure alignment in three dimensionsActa Crystallographica. Section D, Biological Crystallography 60:2256–2268.https://doi.org/10.1107/S0907444904026460
-
A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylationNature Chemical Biology 10:93–95.https://doi.org/10.1038/nchembio.1432
-
The origin of genomic N6-methyl-deoxyadenosine in mammalian cellsNature Chemical Biology 16:630–634.https://doi.org/10.1038/s41589-020-0504-2
-
Shape-specific recognition in the structure of the Vts1p SAM domain with RNANature Structural & Molecular Biology 13:160–167.https://doi.org/10.1038/nsmb1038
-
Structural basis of RNA cap modification by SARS-CoV-2Nature Communications 11:3718.https://doi.org/10.1038/s41467-020-17496-8
-
Beta class amino methyltransferases from bacteria to humans: evolution and structural consequencesNucleic Acids Research 48:10034–10044.https://doi.org/10.1093/nar/gkaa446
Article and author information
Author details
Funding
Max and Minnie Tomerlin Voelcker Fund
- Yogesh K Gupta
IIMS/CTSA pilot award (UL1 TR002645)
- Yogesh K Gupta
Greehey Children's Cancer Research Institute
- Shan Qi
- Yogesh K Gupta
National Institute of Allergy and Infectious Diseases (1R01AI161363-01)
- Yogesh K Gupta
Cancer Prevention & Research Institute of Texas (RP190534)
- Yogesh K Gupta
CPRIT (RP200110)
- Yogesh K Gupta
- Manjeet K Rao
NIH (1R01AI161363)
- Yogesh K Gupta
New England Biolabs
- Siu-Hong Chan
- Nan Dai
- Ivan R Corrêa
University of Texas System (Rising STARs award)
- Yogesh K Gupta
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was partly supported by funding from the Max and Minnie Tomerlin Voelcker Foundation, San Antonio Partnership for Precision Therapeutics, IIMS/CTSA pilot award, and laboratory startup funds from the Greehey Children’s Cancer Research Institute (GCCRI) of UT Health San Antonio to YKG. SQ is supported by the GCCRI. YKG is a recipient of a high impact/high risk award (RP190534) from the Cancer Prevention and Research Institute of Texas (CPRIT), and a Rising STARs award from the University of Texas System. YKG and MKR are supported by CPRIT grant (RP200110). YKG is grateful for the support from NIH grant 1R01AI161363-01. S-HC, ND, and IRC are grateful to Donald Comb, Jim Ellard, and Rich Roberts for their support of research at New England Biolabs.
Copyright
© 2022, Qi et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,708
- views
-
- 391
- downloads
-
- 16
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Cell Biology
Activation of the Wnt/β-catenin pathway crucially depends on the polymerization of dishevelled 2 (DVL2) into biomolecular condensates. However, given the low affinity of known DVL2 self-interaction sites and its low cellular concentration, it is unclear how polymers can form. Here, we detect oligomeric DVL2 complexes at endogenous protein levels in human cell lines, using a biochemical ultracentrifugation assay. We identify a low-complexity region (LCR4) in the C-terminus whose deletion and fusion decreased and increased the complexes, respectively. Notably, LCR4-induced complexes correlated with the formation of microscopically visible multimeric condensates. Adjacent to LCR4, we mapped a conserved domain (CD2) promoting condensates only. Molecularly, LCR4 and CD2 mediated DVL2 self-interaction via aggregating residues and phenylalanine stickers, respectively. Point mutations inactivating these interaction sites impaired Wnt pathway activation by DVL2. Our study discovers DVL2 complexes with functional importance for Wnt/β-catenin signaling. Moreover, we provide evidence that DVL2 condensates form in two steps by pre-oligomerization via high-affinity interaction sites, such as LCR4, and subsequent condensation via low-affinity interaction sites, such as CD2.
-
- Biochemistry and Chemical Biology
The development of proteolysis targeting chimeras (PROTACs), which induce the degradation of target proteins by bringing them into proximity with cellular E3 ubiquitin ligases, has revolutionized drug development. While the human genome encodes more than 600 different E3 ligases, current PROTACs use only a handful of them, drastically limiting their full potential. Furthermore, many PROTAC development campaigns fail because the selected E3 ligase candidates are unable to induce degradation of the particular target of interest. As more and more ligands for novel E3 ligases are discovered, the chemical effort to identify the best E3 ligase for a given target is exploding. Therefore, a genetic system to identify degradation-causing E3 ligases and suitable target/E3 ligase pairs is urgently needed. Here, we used the well-established dimerization of the FKBP12 protein and FRB domain by rapamycin to bring the target protein WDR5 into proximity with candidate E3 ligases. Strikingly, this rapamycin-induced proximity assay (RiPA) revealed that VHL, but not Cereblon, is able to induce WDR5 degradation - a finding previously made by PROTACs, demonstrating its predictive power. By optimizing the steric arrangement of all components and fusing the target protein with a minimal luciferase, RiPA can identify the ideal E3 for any target protein of interest in living cells, significantly reducing and focusing the chemical effort in the early stages of PROTAC development.