A quantitative model predicts how m6A reshapes the kinetic landscape of nucleic acid hybridization and conformational transitions.

N6-methyladenosine (m6A) is a post-transcriptional modification that controls gene expression by recruiting proteins to RNA sites. The modification also slows biochemical processes through mechanisms that are not understood. Using temperature-dependent (20°C-65°C) NMR relaxation dispersion, we show that m6A pairs with uridine with the methylamino group in the anti conformation to form a Watson-Crick base pair that transiently exchanges on the millisecond timescale with a singly hydrogen-bonded low-populated (1%) mismatch-like conformation in which the methylamino group is syn. This ability to rapidly interchange between Watson-Crick or mismatch-like forms, combined with different syn:anti isomer preferences when paired (~1:100) versus unpaired (~10:1), explains how m6A robustly slows duplex annealing without affecting melting at elevated temperatures via two pathways in which isomerization occurs before or after duplex annealing. Our model quantitatively predicts how m6A reshapes the kinetic landscape of nucleic acid hybridization and conformational transitions, and provides an explanation for why the modification robustly slows diverse cellular processes.

NMR Chemical Exchange Measurements Reveal That N6-Methyladenosine Slows RNA Annealing.

N6-Methyladenosine (m6A) is an abundant epitranscriptomic modification that plays important roles in many aspects of RNA metabolism. While m6A is thought to mainly function by recruiting reader proteins to specific RNA sites, the modification can also reshape RNA-protein and RNA-RNA interactions by altering RNA structure mainly by destabilizing base pairing. Little is known about how m6A and other epitranscriptomic modifications might affect the kinetic rates of RNA folding and other conformational transitions that are also important for cellular activity. Here, we used NMR R1ρ relaxation dispersion and chemical exchange saturation transfer to noninvasively and site-specifically measure nucleic acid hybridization kinetics. The methodology was validated on two DNA duplexes and then applied to examine how a single m6A alters the hybridization kinetics in two RNA duplexes. The results show that m6A minimally impacts the rate constant for duplex dissociation, changing koff by ∼1-fold but significantly slows the rate of duplex annealing, decreasing kon by ∼7-fold. A reduction in the annealing rate was observed robustly for two different sequence contexts at different temperatures, both in the presence and absence of Mg2+. We propose that rotation of the N6-methyl group from the preferred syn conformation in the unpaired nucleotide to the energetically disfavored anti conformation required for Watson-Crick pairing is responsible for the reduced annealing rate. The results help explain why in mRNA m6A slows down tRNA selection and more generally suggest that m6A may exert cellular functions by reshaping the kinetics of RNA conformational transitions.

Branch site bulge conformations in domain 6 determine functional sugar puckers in group II intron splicing.

Although group II intron ribozymes are intensively studied the question how structural dynamics affects splicing catalysis has remained elusive. We report for the first time that the group II intron domain 6 exists in a secondary structure equilibrium between a single- and a two-nucleotide bulge conformation, which is directly linked to a switch between sugar puckers of the branch site adenosine. Our study determined a functional sugar pucker equilibrium between the transesterification active C2'-endo conformation of the branch site adenosine in the 1nt bulge and an inactive C3'-endo state in the 2nt bulge fold, allowing the group II intron to switch its activity from the branching to the exon ligation step. Our detailed NMR spectroscopic investigation identified magnesium (II) ions and the branching reaction as regulators of the equilibrium populations. The tuneable secondary structure/sugar pucker equilibrium supports a conformational selection mechanism to up- and downregulate catalytically active and inactive states of the branch site adenosine to orchestrate the multi-step splicing process. The conformational dynamics of group II intron domain 6 is also proposed to be a key aspect for the directionality selection in reversible splicing.

Eukaryotic Translation Elongation is Modulated by Single Natural Nucleotide Derivatives in the Coding Sequences of mRNAs.

RNA modifications are crucial factors for efficient protein synthesis. All classes of RNAs that are involved in translation are modified to different extents. Recently, mRNA modifications and their impact on gene regulation became a focus of interest because they can exert a variety of effects on the fate of mRNAs. mRNA modifications within coding sequences can either directly or indirectly interfere with protein synthesis. In order to investigate the roles of various natural occurring modified nucleotides, we site-specifically introduced them into the coding sequence of reporter mRNAs and subsequently translated them in HEK293T cells. The analysis of the respective protein products revealed a strong position-dependent impact of RNA modifications on translation efficiency and accuracy. Whereas a single 5-methylcytosine (m5C) or pseudouridine () did not reduce product yields, N¹-methyladenosine (m¹A) generally impeded the translation of the respective modified mRNA. An inhibitory effect of 2'O-methlyated nucleotides (Nm) and N6-methyladenosine (m6A) was strongly dependent on their position within the codon. Finally, we could not attribute any miscoding potential to the set of mRNA modifications tested in HEK293T cells.

Studying sparsely populated conformational states in RNA combining chemical synthesis and solution NMR spectroscopy.

Using chemical synthesis and solution NMR spectroscopy, RNA structural ensembles including a major ground state and minor populated excited states can be studied at atomic resolution. In this work, atom-specific 13C labeled RNA building blocks - a 5-13C-uridine and a 2,8-13C2-adenosine building block - are used to introduce isolated 13C-1H-spin topologies into a target RNA to probe such structural ensembles via NMR spectroscopy. First, the 5-13C-uridine 2'-O-TBDMS-phosphoramidite building block was introduced into a 21 nucleotide (nt) tP5c stem construct of the tP5abc subdomain of the Tetrahymena group I ribozyme. Then, the 2,8-13C2-adenosine 2'-O-TBDMS-phosphoramidite building block was incorporated into a 9 kDa and a 15 kD construct derived from the epsilon (ε) RNA element of the duck Hepatitis B virus. The 2,8-13C2-adenosine resonances of the 9 kDa 28 nt sequence could be mapped to the full-length 53 nt construct. The isolated NMR active nuclei pairs were used to probe for low populated excited states (<10%) via 13C-Carr-Purcell-Meiboom-Gill (CPMG)-relaxation dispersion NMR spectroscopy. The 13C-CPMG relaxation dispersion experiment recapitulated a secondary structure switching event in the P5c hairpin of the group I intron construct previously revealed by 15N relaxation dispersion experiments. In the ε-HBV RNA an unfolding event occurring on the millisecond time scale was found in the upper stem in-line with earlier observations. This unpaired conformational state is presumed to be important for the binding of the epsilon reverse transcriptase (RT) enzyme. Thus, a full description of an RNA's folding landscape helps to obtain a deeper understanding of its function, as these high energy conformational states often represent functionally important intermediates involved in (un)folding or ribozyme catalysis.

Synthesis and incorporation of 13C-labeled DNA building blocks to probe structural dynamics of DNA by NMR.

We report the synthesis of atom-specifically 13C-modified building blocks that can be incorporated into DNA via solid phase synthesis to facilitate investigations on structural and dynamic features via NMR spectroscopy. In detail, 6-13C-modified pyrimidine and 8-13C purine DNA phosphoramidites were synthesized and incorporated into a polypurine tract DNA/RNA hybrid duplex to showcase the facile resonance assignment using site-specific labeling. We also addressed micro- to millisecond dynamics in the mini-cTAR DNA. This DNA is involved in the HIV replication cycle and our data points toward an exchange process in the lower stem of the hairpin that is up-regulated in the presence of the HIV-1 nucleocapsid protein 7. As another example, we picked a G-quadruplex that was earlier shown to exist in two folds. Using site-specific 8-13C-2'deoxyguanosine labeling we were able to verify the slow exchange between the two forms on the chemical shift time scale. In a real-time NMR experiment the re-equilibration of the fold distribution after a T-jump could be monitored yielding a rate of 0.012 min-1. Finally, we used 13C-ZZ-exchange spectroscopy to characterize the kinetics between two stacked X-conformers of a Holliday junction mimic. At 25°C, the refolding process was found to occur at a forward rate constant of 3.1 s-1 and with a backward rate constant of 10.6 s-1.