m6A minimally impacts the structure, dynamics, and Rev ARM binding properties of HIV-1 RRE stem IIB

Authors: Chia-Chieh Chu ^aff001; Bei Liu ^aff001; Raphael Plangger ^aff002; Christoph Kreutz ^aff002; Hashim M. Al-Hashimi ^aff001
Authors place of work: Department of Biochemistry, Duke University School of Medicine, Durham, NC, United States of America ^aff001; Institute of Organic Chemistry and Center for Molecular Biosciences CMBI, Universität Innsbruck, Innsbruck, Austria ^aff002; Department of Chemistry, Duke University, Durham, NC, United States of America ^aff003
Published in the journal: PLoS ONE 14(12)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0224850

Summary

N⁶-methyladenosine (m⁶A) is a ubiquitous RNA post-transcriptional modification found in coding as well as non-coding RNAs. m⁶A has also been found in viral RNAs where it is proposed to modulate host-pathogen interactions. Two m⁶A sites have been reported in the HIV-1 Rev response element (RRE) stem IIB, one of which was shown to enhance binding to the viral protein Rev and viral RNA export. However, because these m⁶A sites have not been observed in other studies mapping m⁶A in HIV-1 RNA, their significance remains to be firmly established. Here, using optical melting experiments, NMR spectroscopy, and in vitro binding assays, we show that m⁶A minimally impacts the stability, structure, and dynamics of RRE stem IIB as well as its binding affinity to the Rev arginine-rich-motif (ARM) in vitro. Our results indicate that if present in stem IIB, m⁶A is unlikely to substantially alter the conformational properties of the RNA. Our results add to a growing view that the impact of m⁶A on RNA depends on sequence context and Mg²⁺.

Keywords:

Methylation – HIV-1 – Viral replication – NMR spectroscopy – Sodium phosphate – RNA structure – Melting – Fluorescence polarization

Introduction

N⁶-methyladenosine (m⁶A) is an abundant reversible epitranscriptomic modification found in coding and non-coding RNAs [1–4]. It plays important roles in RNA metabolism [5–8] and is implicated in a growing number of cellular processes [9–15]. m⁶A has also been found in viral RNAs where it is proposed to modulate host-pathogen interactions [16–19].

A number of studies have reported m⁶A in the HIV-1 RNA genome [18, 20, 21]. One study examining HIV-1 infected human T cells [20] reported two m⁶A sites (A68 and A62, Fig 1) in the Rev response element (RRE) stem IIB. RRE is a ~350 nt cis-acting RNA element that is recognized by viral Rev protein to promote the export of unspliced or partially spliced viral RNA to express the structural proteins required for viral replication [22–24]. The two m⁶A sites were found in stem IIB (RREIIB), which is the primary binding site for Rev [23, 25–27]. Knocking down the methyltransferase complex (METTL3/METTL14) was shown to suppress Rev-RRE mediated RNA export and viral replication, and point substitution mutation of one of the two highly conserved adenines (A68) strongly suppressed viral replication (> 90%) [20]. It was proposed [20] that the methyl group of m⁶A68 may interact with Rev protein to stabilize Rev-RRE binding, and/or that m⁶A68 may alter the conformational properties of stem IIB to facilitate Rev recognition. Another different study employing distinct cell lines and mapping methods did not observe these m⁶A sites on RREIIB suggesting that m⁶A can enhance HIV-1 replication and mRNA expression through recruitment of the m⁶A reader proteins YTH-domain containing family (YTHDF) [21]. A third study found m⁶A in RRE but did not identify the specific site [18]. The study proposed that the YTHDF proteins inhibit HIV-1 replication and infection by blocking viral reverse transcription.

m<sup>6</sup>A has limited impact on the thermal stability of RREIIB stem loop. — **Fig. 1. m⁶A has limited impact on the thermal stability of RREIIB stem loop.**
(A) The secondary structure with the proposed N⁶-methyl modified A68 bulge (highlighted in red) of RREIIB native GS and two non-native ES1 and ES2. (B) Impact of m⁶A68 on the thermal stability of the RREIIB hairpins with and without 3 mM Mg²⁺. Shown are the free energy of annealing of the methylated and unmodified RREIIB with and without Mg²⁺. (C) RNA thermodynamic parameters from UV melting. Data are presented as the mean ± standard deviation of at least three independent measurements. Uncertainty in the calculated thermodynamic parameters were determined by error propagation as previously described [48].

Here, we asked whether methylation of RREIIB with m⁶A leads to changes in its conformational and Rev arginine rich motif (ARM) binding properties. We were driven to test this hypothesis because our recent studies showed that the internal loop region of RREIIB near the A68 bulge is highly flexible, and can adopt conformations with alternative secondary structures (Fig 1A) that have different export activities [28, 29]. By redistributing this dynamic ensemble of RREIIB, m⁶A could potentially impact the Rev-RRE interaction and RNA export. Prior studies have shown that m⁶A can reshape RNA-protein and RNA-RNA interactions by modulating RNA structure [30]. For example, a single m⁶A was shown to destabilize RNA duplexes by 0.5–1.7 kcal/mol [31, 32] thus enhancing the binding affinity of proteins to their single-stranded RNA targets [33]. The modification destabilizes A-U base pairs because hydrogen bonding requires that the N⁶-methyl group adopt the unfavorable anti conformation [34, 35]. m⁶A has also been shown to disrupt the non-canonical sheared G-A base pairs to block the assembly of the box C/D snoRNP complexes [36].

We find that modification of the A68 bulge has little effect on the stability, structure, dynamics, as well as the Rev-ARM binding properties of RREIIB. The results indicate that if RREIIB is m⁶A modified, it is unlikely to substantially change the conformational properties of RRE although we cannot rule out that small changes in the conformational properties could modulate Rev-RRE binding in vivo. The results also add to a growing view that the impact of m⁶A on RNA structure depends on sequence context and Mg²⁺ [31, 32, 36, 37]. For example, while m⁶A has been shown to destabilize canonical duplexes [31, 32], it can stabilize junctional A-U base pairs in a Mg²⁺ and secondary structure dependent manner [37] as well as in contexts in which the m⁶A is in a dangling end [32].

Materials and methods

Preparation of RNA samples for NMR studies

RREIIB and m⁶A modified RREIIB (RREIIB^m6A68 and RREIIB^m6A26m6A68) were chemically synthesized using an in-house oligo synthesizer (MerMade 6, BioAutomation) with solid-phase RNA synthesis using N-acetyl protected 2’-tBDSilyl-phosphoramidites (ChemGenes Corporation) and 1 μmol standard columns (1000 Å, BioAutomation) with 4,4’-dimethoxytrityl (DMT)-off synthesis followed by base and 2’-O deprotection (100 μmol DMSO, 125 μL TEA•3HF, heat at 65ºC for 2.5hrs), and ethanol precipitation[28]. A similar approach was used to synthesize site-labeled ¹⁵N3-U72-RREII^m6A68 three-way junction using ¹⁵N3-uridine phosphoramidites [38] and using DMT-off 2’-O deprotection to obtain cleaner NMR spectra.

All RNA samples were purified using 20% (w/v) denaturing polyacrylamide (29:1) gel within 8M urea, 20 mM Tris Borate and 1 mM ethylene-diaminetetraacetate (EDTA) TBE buffer followed by Elutrap electro-elution system (Whatmann, GE healthcare) with 40 mM Tris Acetate and 1 mM EDTA (TAE) buffer then ethanol precipitation. The RNA pellets were dissolved in water and annealed by heating at 95°C for 10 mins then rapidly cooling on ice. After measuring the concentration, the RNA samples were buffer-exchanged into NMR buffer (15 mM sodium phosphate, 25 mM NaCl, 0.1 mM EDTA, with or without 3 mM MgCl₂ at pH = 6.4) three times using 3kDa Amicon Ultra centrifugal filters (EMD Millipore). The final RNA concentrations were ~1.3 mM and ~0.4 mM for RREIIB and RREII, respectively.

NMR experiments

NMR experiments were performed on Bruker Avance III 600 MHz or 700 MHz NMR spectrometers equipped with 5 mm triple-resonance cryogenic probes at 25°C. 2D [¹H,¹³C] (C2/C6/C8-H2/H6/H8) and 2D [¹H,¹⁵N] (N1/N3-H1/H3) heteronuclear single quantum coherence (HSQC) spectra were recorded in the absence or presence of 3 mM Mg²⁺ at 25 ^oC in NMR buffer with 10% D₂O. All NMR data were analyzed using NMRPipe [39] and SPARKY (T.D. Goddard and D.G. Kneller, SPARKY 3, University of California, San Francisco). Unless stated otherwise, all NMR spectra were collected in the presence of Mg²⁺. Resonance intensities were measured using 2D [¹H–¹³C] HSQC experiments. The intensity for each type of C–H spin was normalized to a value of 0.1 using A52-C8H8, A52-C2H2 and U56-C6H6, respectively.

UV melting experiments

Thermal melting experiments were performed on m⁶A modified and unmodified RREIIB using a PerkinElmer Lambda 25 UV/VIS spectrometer equipped with an RTP 6 Peltier Temperature Programmer and a PCB 1500 Water Peltier System. All RNA samples were buffer exchanged into NMR buffer. RNA samples (concentration ~ 1 mM) were then diluted (with NMR buffer) to 3 μM prior to UV melting measurement, which were performed in triplicate (or more) using a sample volume of 400 μL in a Teflon-stoppered 1 cm path length quartz cell. The absorbance at 260 nm was monitored while the temperature was varied between 15 and 95°C at a rate of 1C/min. Thermodynamic parameters were obtained by fitting the UV melting curves using nonlinear model fitting in Mathematica 10.0 (Wolfram Research) as previously described [37].

Fluorescence polarization binding assays

Fluorescence polarization (A) was measured using a CLARIOstar plate reader (BMG LABTECH) using 480 nm excitation and a 540 nm emission filter [28, 40]. Fluorescence polarization binding assays were carried out using 3’-end fluorescein labeled Rev-ARM peptide (Rev-Fl, TRQARRNRRRRWRERQRAAAACK-FITC, LifeTein LLC). The serially diluted RNA in the reaction buffer (30 mM HEPES, pH = 7.0, 100 mM KCl, 10mM sodium phosphate, 10 mM ammonium acetate, 10 mM guanidinium chloride, 2 mM MgCl₂, 20 mM NaCl, 0.5 mM EDTA, and 0.001% (v/v) Triton-X100) was incrementally added into a 384-well plate containing 10 nM Rev-Fl [40]. The binding curves were fitted to single-site binding equation using least-squares methods implemented in Mathematica 10.0 (Wolfram Research).

where A is the measured fluorescence polarization; A_free is the polarization without Rev-Fl binding; A_bound is the polarization with saturated Rev-Fl binding; R_T is the total RNA concentration; L_T is the total Rev-Fl concentration; K_d is the dissociation constant. The uncertainty in (A) was deduced based on the standard deviation over triplicate measurements.

Results

m⁶A68 does not alter the thermal stability of RREIIB

We first used optical melting experiments to examine whether methylation of A68 impacts the thermal stability of RREIIB. All experiments were performed in the presence of 3 mM Mg²⁺ unless stated otherwise. This is important given the impact of Mg²⁺ on RNA folding and dynamics [41] and also given recent studies showing that the impact of m⁶A on RNA structural dynamics and stability can depend on Mg²⁺ [37]. For these experiments, we used a stem IIB construct (Fig 1B) containing the wild-type sequence that was recently shown to recapitulate the conformation of stem IIB in the larger three-way junction context [28]. Prior X-ray [42, 43], NMR [44, 45], SAXS [46] as well as chemical probing data on larger fragments of the HIV genome (~350 nt) [47] indicate that RREII stem IIB adopts the predominant secondary structure shown in Fig 1A. In the dominant RREIIB ground state (GS) conformation, A68 adopts an unpaired bulged conformation. When located in bulged nucleotides, m⁶A has previously been shown to slightly destabilize RNA hairpins by 0.4–0.7 kcal/mol [37] most likely due to the disruption of stacking interactions in the flipped-in conformation. Based on the UV melting data, m⁶A68 had a negligible effect on RREIIB stability (Fig 1B and 1C). Interestingly, a small degree of destabilization (by ~0.4 kcal/mol) was observed in the absence of Mg²⁺, most likely because A68 adopts a partially flipped in conformation in the absence of Mg²⁺, and such a conformation could be more susceptible to destabilization by m⁶A [28]. Thus, the m⁶A effect is overridden by Mg²⁺. Indeed, Mg²⁺ stabilizes RREIIB by ~2 kcal/mol while m⁶A only destabilizes it by ~0.4 kcal/mol (Fig 1C).

m⁶A68 minimally affects structural and dynamic properties of the RREIIB ground state

We used NMR to examine whether methylation of A68 affects the conformation of RREIIB. A single resonance was observed for the N⁶-methyl group in 1D ¹H NMR spectra of RREIIB^m6A68, consistent with a single dominant conformation for m⁶A68 (Fig 2A). The 1D ¹H imino spectrum of RREIIB^m6A68 is virtually identical to its unmodified counterpart, indicating that the methylation does not alter the RREIIB secondary structure (Fig 2B). Very good agreement was also observed between the 2D [¹³C,¹H] HSQC spectra of modified and unmodified RREIIB with only few residues in the internal loop region (A68, U66, G50, G46) showing minor chemical shift perturbations (Fig 2C). These results indicate that m⁶A68 modification does not substantially affect the structure of RREIIB GS.

The m<sup>6</sup>A modification has minor effects on RREIIB structure and dynamics. — **Fig. 2. The m⁶A modification has minor effects on RREIIB structure and dynamics.**
(A) Secondary structure of RRE2Bm^6A68, with the residues showing line-broadening with m⁶A, highlighted in red (with Mg²⁺) and blue (no Mg²⁺). Comparison of the 1D ¹H NMR spectrum of RREIIB^m6A68 with or without Mg²⁺ with the methyl peak indicated by arrows. The comparison of 1D imino spectra (B) and 2D [¹H,¹³C]-HSQC spectra (C) of RREIIB^m6A68 and RREIIB in the presence (red) and absence (blue) of 3 mM Mg²⁺. Resonances exhibiting shifting are indicated using arrows, and those with ambiguous assignments denoted using an asterisk. (D) Normalized resonance intensities in 2D [¹H,¹³C]-HSQC spectra of RREIIB^m6A68 and RREIIB in the presence (red) and absence (blue) of 3 mM Mg²⁺. A52-C8H8, A52-C2H2 and U56-C6H6 were used as a reference and normalized to 0.1. The sample conditions were 1.2–1.5 mM RREIIB^m6A68 or RREIIB in 15 mM sodium phosphate, 25 mM NaCl, 0.1 mM EDTA, pH 6.4 with or without 3 mM MgCl₂.

Prior studies showed that m⁶A can impact RNA conformation in a Mg²⁺ dependent manner [37]. Interestingly, the modification had a larger impact on NMR spectra of RREIIB recorded in the absence of Mg²⁺. Multiple resonances are observed for the N⁶-methyl group, indicating the co-existence of multiple conformations (Fig 2A). The methylation also induced larger perturbations at A52-U66 base pair near the m⁶A68 bulge (Fig 2B). In 2D [¹H, ¹³C] HSQC spectra, the modification induced severe line broadening [49] for resonances belonging to residues G47, G48, G50, A68, C69, G71 and U72 in the internal loop, consistent with enhanced dynamics at the micro-to-millisecond timescales, where only G50 resonance shows broadening in the presence of Mg²⁺ (Fig 2C and 2D). It is likely that the m⁶A induced line broadening in the absence of Mg²⁺ is due to a changes in ES kinetics or populations. The UV melting data showing that m⁶A68 destabilizes RREIIB GS in the absence of Mg²⁺ suggests that the enhanced ES dynamics could arise from destabilization of the GS in the absence of Mg²⁺. Similar NMR results were obtained for a double m⁶A modified RREIIB (RREIIB^m6A62m6A68) (S1 Fig).

m⁶A68 minimally affects the structural and dynamic properties of the excited states in the more native RREII three-way junction

We recently showed that RREII transiently adopts two low-abundance alternative secondary structures (‘excited states’, ES) referred to as ES1 and ES2 (Fig 3A). A68 remains as a bulge in ES1 while forms a G50-A68 mismatch in ES2. The combined m⁶A effects to A68 bulge in ES1 and to the specific G(anti)-A(anti) mismatch in ES2 remain unknown. We therefore examined whether methylation of A68 impacts the dynamics between the GS and ESs in more native RREII three-way junction. This conformational exchange can in principle be measured quantitatively with the use of NMR R_1ρ relaxation dispersion (RD) data [28, 50–52]. However, in practice, severe line-broadening and resonance overlap in NMR spectra of the more native three-way junction (RREII) in the presence of Mg²⁺ complicates measurements [28]. We therefore turned to an alternative approach which we recently developed which uses site-specific stable isotope labeling to directly observe the low-populated ES that form under slow exchange kinetics [28]. With this approach, we were previously able to directly observe imino resonances belonging to ES1 and ES2 in RREII, which slows down the exchange kinetics relative to RREIIB (Fig 3A) [28]. In particular, by site-specifically labeling ¹⁵N3-U72, we observed the G48-U72 mismatch, which uniquely forms in both ES1 and ES2, based on the characteristic imino chemical shift of G-U mismatches (Fig 3A) [28].

Selective site-labeling strategy to probe conformational exchange in larger RREII<sup>m6A68</sup> three-way junction. — **Fig. 3. Selective site-labeling strategy to probe conformational exchange in larger RREII^m6A68 three-way junction.**
(A) Proposed secondary structure for ES1 and ES2 in RREII based on ESs observed in RREIIB. Nucleotides that experience exchange due to ES1 and ES2 are colored green and blue, respectively. (B) The comparison of 1D imino spectrum of RREII^m6A68 and RREII. (C) 2D [¹H,¹⁵N]-HSQC spectra of site-specifically labeled ¹⁵N3-U72-RREII^m6A68 and ¹⁵N3-U72-RREII at 25°C showing a single imino resonance at the characteristic chemical shift region (~11.6 ppm) expected for a G-U in ESs. The sample conditions were 0.3–0.5 mM RNA in 15 mM sodium phosphate, 25 mM NaCl, 0.1 mM EDTA at pH 6.4 with 3 mM MgCl₂.

We used the above strategy to examine how methylation of A68 impacts the GS-ES exchange in RREII. Samples of m⁶A68 modified and unmodified RREII were chemically synthesized with site-labeled ¹⁵N3-U72 (¹⁵N3-U72-RREII^m6A68 and ¹⁵N3-U72-RREII, respectively). The 1D ¹H imino spectra of modified and unmodified RREII were very similar (Fig 3B) indicating that the modification minimally impacts the GS secondary structure even in the more native three-way junction context. The 2D [¹H,¹⁵N] HSQC spectrum of RREII^m6A68 includes a resonance characteristic of the ES G48-U72 mismatch, and it shows excellent overlap with the corresponding resonance observed in unmodified RREII (Fig 3C). No other resonances were observed indicating that m⁶A does not lead to stabilization of alternative conformations in which U72 is base paired. Similar results were obtained in the absence of Mg²⁺ (S2 Fig). These results indicate that m⁶A68 does not significantly impact structural properties of the GS or the populations of the ESs in the native three-way junction both in the presence and absence of Mg²⁺. Note that we cannot entirely rule out that m⁶A destabilizes ES2 or ES1 and that the observed G48-U72 mismatch reflects either ES1 or ES2, respectively.

m⁶A68 has a negligible effect on Rev-RRM binding to RREIIB

We used a fluorescence polarization binding assay to examine whether methylation of A68 impacts binding of fluorescein labeled Rev-ARM peptide (Rev-Fl) to RREIIB (Fig 4) [28, 40]. Unmodified RREIIB binds to Rev-ARM peptide with apparent K_d = 30.6 ± 5.8 nM, in agreement with prior studies [44, 53]. The binding affinity decreased two-fold (K_d = 62.2 ± 23.8 nM) for the modified RREIIB^m6A68, indicating that the binding affinity for m⁶A modified was only slightly weakened relative to unmodified RREIIB especially when considering the uncertainty.

Measurement of the binding affinity of Rev-ARM peptide to m<sup>6</sup>A modified or unmodified RREIIB using fluorescence polarization. — **Fig. 4. Measurement of the binding affinity of Rev-ARM peptide to m⁶A modified or unmodified RREIIB using fluorescence polarization.**
Normalized anisotropy values measured for RREIIB^m6A68 and RREIIB fitted with one-site binding model (see ‘*Materials and Methods*‘ section). The anisotropy value observed in the absence of RNA was normalized to 0.1. Uncertainty reflects the standard deviation from three independent measurements. Buffer conditions: 30 mM HEPES, pH = 7.0, 100 mM KCl, 10mM sodium phosphate, 10 mM ammonium acetate, 10 mM guanidinium chloride, 2 mM MgCl₂, 20 mM NaCl, 0.5 mM EDTA and 0.001% (v/v) Triton-X100. Concentration of Rev-Fl peptide was 10 nM.

Discussions

Methylation of RREIIB can promote Rev-RRE interaction by changing the RREIIB conformation so as to favorably bind to Rev or by promoting Rev-RRE binding through direct interaction involving the methyl group. Our results argue against a significant impact on RREIIB structure, as might be expected given placement of m⁶A68 in the bulge. In addition, the modification slightly weakened binding of the Rev-ARM, and this is consistent with prior NMR [54] and X-ray crystallography [55] studies showing that A68 bulge interacts with Rev protein through non-specific Van Der Waals contacts, and that deleting A68 has minor effects on Rev or Rev-ARM binding affinity to RRE stem IIB [26, 56]. Our data cannot rule out that the exposed and accessible m⁶A enhances binding in the context of full-length Rev and RRE or promotes RNA export through the recruitment of other host export factors, and that this in turn gives rise to the reported ~2–3 fold decrease in HIV viral RNA pull down using Rev when knocking down METTL3/METTL14 responsible for producing m⁶A [20]. It is also possible that the structural and dynamic properties of m⁶A modified RREIIB differ in vivo relative to in vitro. A recent study [29] showed that the key dynamic properties of unmodified RREIIB are similar in vitro and in cells. Further studies are needed to address these possibilities.

The modification had a greater effect on the stability and conformation of RREIIB in the absence of Mg²⁺. We previously showed that Mg²⁺ redistributes the RREIIB ensemble by stabilizing the GS and ES1 relative to ES2 [28]. Our results indicate that Mg²⁺ has a larger effect on the relative stability of these different conformations as compared to the methylation. This underscores the importance of studying the impact of m⁶A modifications in the presence of Mg²⁺ as in a more physiological conditions.

Supporting information

S1 Fig [a]
The double mA modifications have minor effects on RREIIB structure and dynamics in the presence of Mg.

S2 Fig [a]
Selective site-labeling strategy to probe conformational exchange in larger 3-way junction RREII in the absence of Mg.

Zdroje

1. Meyer KD, Saletore Y, Zumbo P, Elemento O, Mason CE, Jaffrey SR. Comprehensive analysis of mRNA methylation reveals enrichment in 3' UTRs and near stop codons. Cell. 2012;149(7):1635–46. doi: 10.1016/j.cell.2012.05.003 22608085; PubMed Central PMCID: PMC3383396.

2. Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, et al. Topology of the human and mouse m(6)A RNA methylomes revealed by m(6)A-seq. Nature. 2012;485(7397):201–U84. doi: 10.1038/nature11112 WOS:000303799800033. 22575960

3. Desrosiers R, Friderici K, Rottman F. Identification of Methylated Nucleosides in Messenger-Rna from Novikoff Hepatoma-Cells. Proceedings of the National Academy of Sciences of the United States of America. 1974;71(10):3971–5. doi: 10.1073/pnas.71.10.3971 WOS:A1974U614000040. 4372599

4. Li X, Xiong X, Yi C. Epitranscriptome sequencing technologies: decoding RNA modifications. Nat Methods. 2016;14(1):23–31. doi: 10.1038/nmeth.4110 28032622.

5. Wang X, Zhao BS, Roundtree IA, Lu ZK, Han DL, Ma HH, et al. N-6-methyladenosine Modulates Messenger RNA Translation Efficiency. Cell. 2015;161(6):1388–99. doi: 10.1016/j.cell.2015.05.014 WOS:000355935000017. 26046440

6. Wang X, Lu ZK, Gomez A, Hon GC, Yue YN, Han DL, et al. N-6-methyladenosine-dependent regulation of messenger RNA stability. Nature. 2014;505(7481):117–+. doi: 10.1038/nature12730 WOS:000329163300035. 24284625

7. Zhao X, Yang Y, Sun BF, Shi Y, Yang X, Xiao W, et al. FTO-dependent demethylation of N6-methyladenosine regulates mRNA splicing and is required for adipogenesis. Cell Res. 2014;24(12):1403–19. doi: 10.1038/cr.2014.151 WOS:000345894800006. 25412662

8. Xiao W, Adhikari S, Dahal U, Chen YS, Hao YJ, Sun BF, et al. Nuclear m(6)A Reader YTHDC1 Regulates mRNA Splicing. Molecular cell. 2016;61(4):507–19. doi: 10.1016/j.molcel.2016.01.012 WOS:000372326300005. 26876937

9. Zhao BXS, Roundtree IA, He C. Post-transcriptional gene regulation by mRNA modifications. Nat Rev Mol Cell Bio. 2017;18(1):31–42. doi: 10.1038/nrm.2016.132 WOS:000393267500008. 27808276

10. Roundtree IA, Evans ME, Pan T, He C. Dynamic RNA Modifications in Gene Expression Regulation. Cell. 2017;169(7):1187–200. doi: 10.1016/j.cell.2017.05.045 WOS:000403332400009. 28622506

11. Meyer KD, Jaffrey SR. Rethinking m(6)A Readers, Writers, and Erasers. Annu Rev Cell Dev Bi. 2017;33:319–42. doi: 10.1146/annurev-cellbio-100616-060758 WOS:000413584700015. 28759256

12. Deng XL, Su R, Feng XS, Wei MJ, Chen JJ. Role of N-6-methyladenosine modification in cancer. Curr Opin Genet Dev. 2018;48:1–7. doi: 10.1016/j.gde.2017.10.005 WOS:000428232100002. 29040886

13. Chen T, Hao Y-J, Zhang Y, Li M-M, Wang M, Han W, et al. m6A RNA Methylation Is Regulated by MicroRNAs and Promotes Reprogramming to Pluripotency. Cell stem cell. 2015;16(3):289–301. doi: 10.1016/j.stem.2015.01.016 25683224

14. Weng Y-L, Wang X, An R, Cassin J, Vissers C, Liu Y, et al. Epitranscriptomic m6A Regulation of Axon Regeneration in the Adult Mammalian Nervous System. Neuron. 2018;97(2):313–25.e6. doi: 10.1016/j.neuron.2017.12.036 29346752

15. Gokhale NS, McIntyre ABR, McFadden MJ, Roder AE, Kennedy EM, Gandara JA, et al. N6-Methyladenosine in Flaviviridae Viral RNA Genomes Regulates Infection. Cell Host Microbe. 2016;20(5):654–65. Epub 2016/10/25. doi: 10.1016/j.chom.2016.09.015 27773535; PubMed Central PMCID: PMC5123813.

16. Gokhale NS, McIntyre ABR, McFadden MJ, Roder AE, Kennedy EM, Gandara JA, et al. N6-Methyladenosine in Flaviviridae Viral RNA Genomes Regulates Infection. Cell Host Microbe. 2016;20(5):654–65. doi: 10.1016/j.chom.2016.09.015 WOS:000389472000014. 27773535

17. Lichinchi G, Zhao BS, Wu YA, Lu ZK, Qin Y, He C, et al. Dynamics of Human and Viral RNA Methylation during Zika Virus Infection. Cell Host Microbe. 2016;20(5):666–73. doi: 10.1016/j.chom.2016.10.002 WOS:000389472000015. 27773536

18. Tirumuru N, Zhao BS, Lu WX, Lu ZK, He C, Wu L. N-6-methyladenosine of HIV-1 RNA regulates viral infection and HIV-1 Gag protein expression. eLife. 2016;5. ARTN e15528 10.7554/eLife.15528. WOS:000380858000001.

19. Hesser CR, Karijolich J, Dominissini D, He C, Glaunsinger BA. N-6-methyladenosine modification and the YTHDF2 reader protein play cell type specific roles in lytic viral gene expression during Kaposi's sarcoma-associated herpesvirus infection. Plos Pathog. 2018;14(4). ARTN e1006995 10.1371/journal.ppat.1006995. WOS:000431135400036.

20. Lichinchi G, Gao S, Saletore Y, Gonzalez GM, Bansal V, Wang Y, et al. Dynamics of the human and viral m6A RNA methylomes during HIV-1 infection of T cells. Nature Microbiology. 2016;1. doi: 10.1038/nmicrobiol.2016.11 27572442

21. Kennedy EM, Bogerd HP, Kornepati AVR, Kang D, Ghoshal D, Marshall JB, et al. Posttranscriptional m(6)A Editing of HIV-1 mRNAs Enhances Viral Gene Expression. Cell Host Microbe. 2016;19(5):675–85. doi: 10.1016/j.chom.2016.04.002 WOS:000375595500017. 27117054

22. Malim MH, Hauber J, Le SY, Maizel JV, Cullen BR. The HIV-1 rev trans-activator acts through a structured target sequence to activate nuclear export of unspliced viral mRNA. Nature. 1989;338(6212):254–7. doi: 10.1038/338254a0 2784194.

23. Malim MH, Tiley LS, McCarn DF, Rusche JR, Hauber J, Cullen BR. HIV-1 structural gene expression requires binding of the Rev trans-activator to its RNA target sequence. Cell. 1990;60(4):675–83. doi: 10.1016/0092-8674(90)90670-a 2406030.

24. Malim MH, Cullen BR. HIV-1 structural gene expression requires the binding of multiple Rev monomers to the viral RRE: implications for HIV-1 latency. Cell. 1991;65(2):241–8. doi: 10.1016/0092-8674(91)90158-u 2015625.

25. Bartel DP, Zapp ML, Green MR, Szostak JW. HIV-1 Rev regulation involves recognition of non-Watson-Crick base pairs in viral RNA. Cell. 1991;67(3):529–36. doi: 10.1016/0092-8674(91)90527-6 1934059.

26. Heaphy S, Finch JT, Gait MJ, Karn J, Singh M. Human immunodeficiency virus type 1 regulator of virion expression, rev, forms nucleoprotein filaments after binding to a purine-rich "bubble" located within the rev-responsive region of viral mRNAs. Proceedings of the National Academy of Sciences of the United States of America. 1991;88(16):7366–70. doi: 10.1073/pnas.88.16.7366 1871141; PubMed Central PMCID: PMC52296.

27. Huang XJ, Hope TJ, Bond BL, McDonald D, Grahl K, Parslow TG. Minimal Rev-response element for type 1 human immunodeficiency virus. Journal of virology. 1991;65(4):2131–4. 2002556; PubMed Central PMCID: PMC240087.

28. Chu CC, Plangger R, Kreutz C, Al-Hashimi HM. Dynamic ensemble of HIV-1 RRE stem IIB reveals non-native conformations that disrupt the Rev-binding site. Nucleic acids research. 2019;47(13):7105–17. doi: 10.1093/nar/gkz498 31199872; PubMed Central PMCID: PMC6649712.

29. Ganser LR, Chu C-C, Bogerd HP, Kelly ML, Cullen BR, Al-Hashimi HM. Probing RNA conformational equilibria within the functional cellular context. bioRxiv. 2019:634576. doi: 10.1101/634576

30. Spitale RC, Flynn RA, Zhang QC, Crisalli P, Lee B, Jung JW, et al. Structural imprints in vivo decode RNA regulatory mechanisms. Nature. 2015;519(7544):486–90. doi: 10.1038/nature14263 25799993; PubMed Central PMCID: PMC4376618.

31. Kierzek E, Kierzek R. The thermodynamic stability of RNA duplexes and hairpins containing N6-alkyladenosines and 2-methylthio-N6-alkyladenosines. Nucleic acids research. 2003;31(15):4472–80. doi: 10.1093/nar/gkg633 12888507; PubMed Central PMCID: PMC169893.

32. Roost C, Lynch SR, Batista PJ, Qu K, Chang HY, Kool ET. Structure and thermodynamics of N6-methyladenosine in RNA: a spring-loaded base modification. Journal of the American Chemical Society. 2015;137(5):2107–15. doi: 10.1021/ja513080v 25611135; PubMed Central PMCID: PMC4405242.

33. Liu N, Dai Q, Zheng G, He C, Parisien M, Pan T. N6-methyladenosine-dependent RNA structural switches regulate RNA–protein interactions. Nature. 2015;518:560. doi: 10.1038/nature14234 25719671

34. Roost C, Lynch SR, Batista PJ, Qu K, Chang HY, Kool ET. Structure and Thermodynamics of N6-Methyladenosine in RNA: A Spring-Loaded Base Modification. Journal of the American Chemical Society. 2015;137(5):2107–15. doi: 10.1021/ja513080v 25611135

35. Kierzek E, Kierzek R. The thermodynamic stability of RNA duplexes and hairpins containing N-6-alkyladenosines and 2-methylthio-N-6-alkyladenosines. Nucleic acids research. 2003;31(15):4472–80. doi: 10.1093/nar/gkg633 WOS:000184532900031. 12888507

36. Huang L, Ashraf S, Wang J, Lilley DM. Control of box C/D snoRNP assembly by N(6)-methylation of adenine. EMBO Rep. 2017;18(9):1631–45. doi: 10.15252/embr.201743967 28623187; PubMed Central PMCID: PMC5579392.

37. Liu B, Merriman DK, Choi SH, Schumacher MA, Plangger R, Kreutz C, et al. A potentially abundant junctional RNA motif stabilized by m(6)A and Mg(2). Nat Commun. 2018;9(1):2761. doi: 10.1038/s41467-018-05243-z 30018356; PubMed Central PMCID: PMC6050335.

38. Neuner S, Santner T, Kreutz C, Micura R. The "Speedy" Synthesis of Atom-Specific (15)N Imino/Amido-Labeled RNA. Chemistry. 2015;21(33):11634–43. doi: 10.1002/chem.201501275 26237536; PubMed Central PMCID: PMC4946632.

39. Delaglio F, Grzesiek S, Vuister GW, Zhu G, Pfeifer J, Bax A. NMRPipe: a multidimensional spectral processing system based on UNIX pipes. Journal of biomolecular NMR. 1995;6(3):277–93. doi: 10.1007/bf00197809 8520220.

40. Luedtke NW, Tor Y. Fluorescence-based methods for evaluating the RNA affinity and specificity of HIV-1 Rev-RRE inhibitors. Biopolymers. 2003;70(1):103–19. doi: 10.1002/bip.10428 12925996.

41. Draper DE, Grilley D, Soto AM. Ions and RNA folding. Annu Rev Bioph Biom. 2005;34:221–43. doi: 10.1146/annurev.biophys.34.040204.144511 WOS:000230099600010. 15869389

42. Ippolito JA, Steitz TA. The structure of the HIV-1 RRE high affinity rev binding site at 1.6 A resolution. Journal of molecular biology. 2000;295(4):711–7. doi: 10.1006/jmbi.1999.3405 10656783.

43. Hung LW, Holbrook EL, Holbrook SR. The crystal structure of the Rev binding element of HIV-1 reveals novel base pairing and conformational variability. Proceedings of the National Academy of Sciences of the United States of America. 2000;97(10):5107–12. doi: 10.1073/pnas.090588197 10792052; PubMed Central PMCID: PMC25789.

44. Battiste JL, Tan R, Frankel AD, Williamson JR. Binding of an HIV Rev peptide to Rev responsive element RNA induces formation of purine-purine base pairs. Biochemistry. 1994;33(10):2741–7. doi: 10.1021/bi00176a001 8130185.

45. Peterson RD, Bartel DP, Szostak JW, Horvath SJ, Feigon J. 1H NMR studies of the high-affinity Rev binding site of the Rev responsive element of HIV-1 mRNA: base pairing in the core binding element. Biochemistry. 1994;33(18):5357–66. doi: 10.1021/bi00184a001 8180157.

46. Fang X, Wang J, O'Carroll IP, Mitchell M, Zuo X, Wang Y, et al. An unusual topological structure of the HIV-1 Rev response element. Cell. 2013;155(3):594–605. doi: 10.1016/j.cell.2013.10.008 24243017; PubMed Central PMCID: PMC3918456.

47. Bai Y, Tambe A, Zhou K, Doudna JA. RNA-guided assembly of Rev-RRE nuclear export complexes. eLife. 2014;3:e03656. doi: 10.7554/eLife.03656 25163983; PubMed Central PMCID: PMC4142337.

48. Siegfried NA, Metzger SL, Bevilacqua PC. Folding cooperativity in RNA and DNA is dependent on position in the helix. Biochemistry. 2007;46(1):172–81. doi: 10.1021/bi061375l WOS:000243157300018. 17198387

49. Getz MM, Andrews AJ, Fierke CA, Al-Hashimi HM. Structural plasticity and Mg2+ binding properties of RNase P P4 from combined analysis of NMR residual dipolar couplings and motionally decoupled spin relaxation. Rna. 2007;13(2):251–66. Epub 2006/12/30. doi: 10.1261/rna.264207 17194721; PubMed Central PMCID: PMC1781369.

50. Palmer AG 3rd, Massi F. Characterization of the dynamics of biomacromolecules using rotating-frame spin relaxation NMR spectroscopy. Chem Rev. 2006;106(5):1700–19. doi: 10.1021/cr0404287 16683750.

51. Korzhnev DM, Orekhov VY, Kay LE. Off-resonance R(1rho) NMR studies of exchange dynamics in proteins with low spin-lock fields: an application to a Fyn SH3 domain. Journal of the American Chemical Society. 2005;127(2):713–21. doi: 10.1021/ja0446855 15643897.

52. Rangadurai A, Szymaski ES, Kimsey IJ, Shi H, Al-Hashimi HM. Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1rho relaxation dispersion. Prog Nucl Magn Reson Spectrosc. 2019;112–113:55–102. Epub 2019/09/05. doi: 10.1016/j.pnmrs.2019.05.002 31481159; PubMed Central PMCID: PMC6727989.

53. Daugherty MD, D'Orso I, Frankel AD. A solution to limited genomic capacity: using adaptable binding surfaces to assemble the functional HIV Rev oligomer on RNA. Molecular cell. 2008;31(6):824–34. doi: 10.1016/j.molcel.2008.07.016 18922466; PubMed Central PMCID: PMC2651398.

54. Battiste JL, Mao H, Rao NS, Tan R, Muhandiram DR, Kay LE, et al. Alpha helix-RNA major groove recognition in an HIV-1 rev peptide-RRE RNA complex. Science. 1996;273(5281):1547–51. doi: 10.1126/science.273.5281.1547 8703216.

55. Jayaraman B, Crosby DC, Homer C, Ribeiro I, Mavor D, Frankel AD. RNA-directed remodeling of the HIV-1 protein Rev orchestrates assembly of the Rev-Rev response element complex. eLife. 2014;3:e04120. doi: 10.7554/eLife.04120 25486594; PubMed Central PMCID: PMC4360532.

56. Kjems J, Calnan BJ, Frankel AD, Sharp PA. Specific binding of a basic peptide from HIV-1 Rev. The EMBO journal. 1992;11(3):1119–29. 1547776; PubMed Central PMCID: PMC556554.