Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

Direct Measurement of 8OG Syn-Anti Flips in Mutagenic 8OG·A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter their propensities to form low-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn·Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly populated (pop. ∼ 5%) and short-lived (lifetime ∼50 ms) nonmutagenic 8OGanti·Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on the labeled partner adenine. The Watson-Crick-like 8OGsyn·Aanti mismatch also rescued the kinetics of Hoogsteen motions at distant A-T base pairs, which the G·A mismatch had slowed down. The results lend further support for 8OGanti·Aanti as a minor conformational state of 8OG·A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G·T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R1ρ experiments, we show that dG·rU and dT·rG mismatches in two RNA:DNA hybrids transiently form tautomeric (Genol·T/U $ \mathbin{\lower.3ex\hbox{$\buildrel\textstyle\rightarrow\over {\smash{\leftarrow}\vphantom{_{\vbox to.5ex{\vss}}}}$}}$ G·Tenol/Uenol) and anionic (G·T-/U-) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG·rU- formed with a ten-fold higher propensity relative to dT-·rG and dG·dT- and this could be attributed to the lower pKa (ΔpKa ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G·T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G·U versus G·T, and indicate that anionic Watson-Crick-like G·U- could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

Beyond slow two-state protein conformational exchange using CEST: applications to three-state protein interconversion on the millisecond timescale

Although NMR spectroscopy is routinely used to study the conformational dynamics of biomolecules, robust analyses of the data are challenged in cases where exchange is more complex than two-state, such as when a 'visible' major conformer exchanges with two 'invisible' minor states on the millisecond timescale. It is becoming increasingly clear that chemical exchange saturation transfer (CEST) NMR experiments that were initially developed to study systems undergoing slow interconversion are also sensitive to intermediate-fast timescale biomolecular conformational exchange. Here we investigate the utility of the amide 15N CEST experiment to characterise protein three-state exchange occurring on the millisecond timescale by studying the interconversion between the folded (F) state of the FF domain from human HYPA/FBP11 (WT FF) and two of its folding intermediates I1 and I2. Although 15N CPMG experiments are consistent with the F state interconverting with a single minor state on the millisecond timescale, 15N CEST data clearly establish an exchange process between F and a pair of minor states. A unique three-state exchange model cannot be obtained by analysis of 15N CEST data recorded at a single temperature. However, including the relative sign of the difference in the chemical shifts of the two minor states based on a simple two-state analysis of CEST data recorded at multiple temperatures, results in a robust three-state model in which the F, I1 and I2 states interconvert with each other on the millisecond timescale ( k e x , F I 1 ~ 550 s-1, k e x , F I 2 ~ 1200 s-1, k e x , I 1 I 2 ~ 5000 s-1), with I1 and I2 sparsely populated at ~ 0.15% and ~ 0.35%, respectively, at 15 °C. A computationally demanding grid-search of exchange parameter space is not required to extract the best-fit exchange parameters from the CEST data. The utility of the CEST experiment, thus, extends well beyond studies of conformers in slow exchange on the NMR chemical shift timescale, to include systems with interconversion rates on the order of thousands/second.

Direct Measurement of 8OG syn-anti Flips in Mutagenic 8OG•A and Long-Range Damage-Dependent Hoogsteen Breathing Dynamics Using 1H CEST NMR

Elucidating how damage impacts DNA dynamics is essential for understanding the mechanisms of damage recognition and repair. Many DNA lesions alter the propensities to form lowly-populated and short-lived conformational states. However, NMR methods to measure these dynamics require isotopic enrichment, which is difficult for damaged nucleotides. Here, we demonstrate the utility of the 1H chemical exchange saturation transfer (CEST) NMR experiment in measuring the dynamics of oxidatively damaged 8-oxoguanine (8OG) in the mutagenic 8OGsyn•Aanti mismatch. Using 8OG-H7 as an NMR probe of the damaged base, we directly measured 8OG syn-anti flips to form a lowly-populated (pop. ~ 5%) and short-lived (lifetime ~ 50 ms) non-mutagenic 8OGanti•Aanti. These exchange parameters were in quantitative agreement with values from 13C off-resonance R1ρ and CEST on a labeled partner adenine. The Watson-Crick-like 8OGsyn•Aanti mismatch also rescued the kinetics of Hoogsteen motions at distance A-T base pairs, which the G•A mismatch had slowed down. The results lend further support for 8OGanti•Aanti as a minor conformational state of 8OG•A, reveal that 8OG damage can impact Hoogsteen dynamics at a distance, and demonstrate the utility of 1H CEST for measuring damage-dependent dynamics in unlabeled DNA.

Z-Form Adoption of Nucleic Acid is a Multi-Step Process Which Proceeds through a Melted Intermediate

The left-handed Z-conformation of nucleic acids can be adopted by both DNA and RNA when bound by Zα domains found within a variety of innate immune response proteins. Zα domains stabilize this higher-energy conformation by making specific interactions with the unique geometry of Z-DNA/Z-RNA. However, the mechanism by which a right-handed helix contorts to become left-handed in the presence of proteins, including the intermediate steps involved, is poorly understood. Through a combination of nuclear magnetic resonance (NMR) and other biophysical measurements, we have determined that in the absence of Zα, under low salt conditions at room temperature, d(CpG) and r(CpG) constructs show no observable evidence of transient Z-conformations greater than 0.5% on either the intermediate or slow NMR time scales. At higher temperatures, we observed a transient unfolded intermediate. The ease of melting a nucleic acid duplex correlates with Z-form adoption rates in the presence of Zα. The largest contributing factor to the activation energies of Z-form adoption as calculated by Arrhenius plots is the ease of flipping the sugar pucker, as required for Z-DNA and Z-RNA. Together, these data validate the previously proposed "zipper model" for Z-form adoption in the presence of Zα. Overall, Z-conformations are more likely to be adopted by double-stranded DNA and RNA regions flanked by less stable regions and by RNAs experiencing torsional/mechanical stress.

An RNA excited conformational state at atomic resolution

Sparse and short-lived excited RNA conformational states are essential players in cell physiology, disease, and therapeutic development, yet determining their 3D structures remains challenging. Combining mutagenesis, NMR spectroscopy, and computational modeling, we determined the 3D structural ensemble formed by a short-lived (lifetime ~2.1 ms) lowly-populated (~0.4%) conformational state in HIV-1 TAR RNA. Through a strand register shift, the excited conformational state completely remodels the 3D structure of the ground state (RMSD from the ground state = 7.2 ± 0.9 Å), forming a surprisingly more ordered conformational ensemble rich in non-canonical mismatches. The structure impedes the formation of the motifs recognized by Tat and the super elongation complex, explaining why this alternative TAR conformation cannot activate HIV-1 transcription. The ability to determine the 3D structures of fleeting RNA states using the presented methodology holds great promise for our understanding of RNA biology, disease mechanisms, and the development of RNA-targeting therapeutics.

Kinetic Resolution of the Atomic 3D Structures Formed by Ground and Excited Conformational States in an RNA Dynamic Ensemble

Knowing the 3D structures formed by the various conformations populating the RNA free-energy landscape, their relative abundance, and kinetic interconversion rates is required to obtain a quantitative and predictive understanding of how RNAs fold and function at the atomic level. While methods integrating ensemble-averaged experimental data with computational modeling are helping define the most abundant conformations in RNA ensembles, elucidating their kinetic rates of interconversion and determining the 3D structures of sparsely populated short-lived RNA excited conformational states (ESs) remains challenging. Here, we developed an approach integrating Rosetta-FARFAR RNA structure prediction with NMR residual dipolar couplings and relaxation dispersion that simultaneously determines the 3D structures formed by the ground-state (GS) and ES subensembles, their relative abundance, and kinetic rates of interconversion. The approach is demonstrated on HIV-1 TAR, whose six-nucleotide apical loop was previously shown to form a sparsely populated (∼13%) short-lived (lifetime ∼ 45 μs) ES. In the GS, the apical loop forms a broad distribution of open conformations interconverting on the pico-to-nanosecond time scale. Most residues are unpaired and preorganized to bind the Tat-superelongation protein complex. The apical loop zips up in the ES, forming a narrow distribution of closed conformations, which sequester critical residues required for protein recognition. Our work introduces an approach for determining the 3D ensemble models formed by sparsely populated RNA conformational states, provides a rare atomic view of an RNA ES, and kinetically resolves the atomic 3D structures of RNA conformational substates, interchanging on time scales spanning 6 orders of magnitude, from picoseconds to microseconds.

Probing excited state 1Hα chemical shifts in intrinsically disordered proteins with a triple resonance-based CEST experiment: Application to a disorder-to-order switch

Over 40% of eukaryotic proteomes and 15% of bacterial proteomes are predicted to be intrinsically disordered based on their amino acid sequence. Intrinsically disordered proteins (IDPs) exist as heterogeneous ensembles of interconverting conformations and pose a challenge to the structure-function paradigm by apparently functioning without possessing stable structural elements. IDPs play a prominent role in biological processes involving extensive intermolecular interaction networks and their inherently dynamic nature facilitates their promiscuous interaction with multiple structurally diverse partner molecules. NMR spectroscopy has made pivotal contributions to our understanding of IDPs because of its unique ability to characterize heterogeneity at atomic resolution. NMR methods such as Chemical Exchange Saturation Transfer (CEST) and relaxation dispersion have enabled the detection of 'invisible' excited states in biomolecules which are transiently and sparsely populated, yet central for function. Here, we develop a 1Hα CEST pulse sequence which overcomes the resonance overlap problem in the 1Hα-13Cα plane of IDPs by taking advantage of the superior resolution in the 1H-15N correlation spectrum. In this sequence, magnetization is transferred after 1H CEST using a triple resonance coherence transfer pathway from 1Hα (i) to 1HN(i + 1) during which the 15N(t1) and 1HN(t2) are frequency labelled. This approach is integrated with spin state-selective CEST for eliminating spurious dips in CEST profiles resulting from dipolar cross-relaxation. We apply this sequence to determine the excited state 1Hα chemical shifts of the intrinsically disordered DNA binding domain (CytRN) of the bacterial cytidine repressor (CytR), which transiently acquires a functional globally folded conformation. The structure of the excited state, calculated using 1Hα chemical shifts in conjunction with other excited state NMR restraints, is a three-helix bundle incorporating a helix-turn-helix motif that is vital for binding DNA.

NMR measurements of transient low-populated tautomeric and anionic Watson-Crick-like G·T/U in RNA:DNA hybrids: Implications for the fidelity of transcription and CRISPR/Cas9 gene editing

Many biochemical processes use the Watson-Crick geometry to distinguish correct from incorrect base pairing. However, on rare occasions, mismatches such as G•T/U can transiently adopt Watson-Crick-like conformations through tautomerization or ionization of the bases, giving rise to replicative and translational errors. The propensities to form Watson-Crick-like mismatches in RNA:DNA hybrids remain unknown, making it unclear whether they can also contribute to errors during processes such as transcription and CRISPR/Cas editing. Here, using NMR R 1ρ experiments, we show that dG•rU and dT•rG mismatches in two RNA:DNA hybrids transiently form tautomeric (G enol •T/U ⇄G•T enol /U enol ) and anionic (G•T - /U - ) Watson-Crick-like conformations. The tautomerization dynamics were like those measured in A-RNA and B-DNA duplexes. However, anionic dG•rU - formed with a ten-fold higher propensity relative to dT - •rG and dG•dT - and this could be attributed to the lower pK a (Δ pK a ∼0.4-0.9) of U versus T. Our findings suggest plausible roles for Watson-Crick-like G•T/U mismatches in transcriptional errors and CRISPR/Cas9 off-target gene editing, uncover a crucial difference between the chemical dynamics of G•U versus G•T, and indicate that anionic Watson-Crick-like G•U - could play a significant role evading Watson-Crick fidelity checkpoints in RNA:DNA hybrids and RNA duplexes.

Zero- to low-field relaxometry of chemical and biological fluids

Nuclear magnetic resonance (NMR) relaxometry is an analytical method that provides information about molecular environments, even for NMR "silent" molecules (spin-0), by analyzing the properties of NMR signals versus the magnitude of the longitudinal field. Conventionally, this technique is performed at fields much higher than Earth's magnetic field, but our work focuses on NMR relaxometry at zero and ultra-low magnetic fields (ZULFs). Operating under such conditions allows us to investigate slow (bio)chemical processes occurring on a timescale from milliseconds to seconds, which coincide with spin evolution. ZULFs also minimize T2 line broadening in heterogeneous samples resulting from magnetic susceptibility. Here, we use ZULF NMR relaxometry to analyze (bio)chemical compounds containing 1H-13C, 1H-15N, and 1H-31P spin pairs. We also detected high-quality ULF NMR spectra of human whole-blood at 0.8 μT, despite a shortening of spin relaxation by blood proteomes (e.g., hemoglobin). Information on proton relaxation times of blood, a potential early biomarker of inflammation, can be acquired in under a minute using inexpensive, portable/small-size NMR spectrometers based on atomic magnetometers.

Studying micro to millisecond protein dynamics using simple amide 15N CEST experiments supplemented with major-state R2 and visible peak-position constraints

Over the last decade amide 15N CEST experiments have emerged as a popular tool to study protein dynamics that involves exchange between a 'visible' major state and sparsely populated 'invisible' minor states. Although initially introduced to study exchange between states that are in slow exchange with each other (typical exchange rates of, 10 to 400 s-1), they are now used to study interconversion between states on the intermediate to fast exchange timescale while still using low to moderate (5 to 350 Hz) 'saturating' B1 fields. The 15N CEST experiment is very sensitive to exchange as the exchange delay TEX can be quite long (~0.5 s) allowing for a large number of exchange events to occur making it a very powerful tool to detect minor sates populated ([Formula: see text]) to as low as 1%. When systems are in fast exchange and the 15N CEST data has to be described using a model that contains exchange, the exchange parameters are often poorly defined because the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus exchange rate ([Formula: see text]) plots can be quite flat with shallow or no minima and the analysis of such 15N CEST data can lead to wrong estimates of the exchange parameters due to the presence of 'spurious' minima. Here we show that the inclusion of experimentally derived constraints on the intrinsic transverse relaxation rates and the inclusion of visible state peak-positions during the analysis of amide 15N CEST data acquired with moderate B1 values (~50 to ~350 Hz) results in convincing minima in the [Formula: see text] versus [Formula: see text] and the [Formula: see text] versus [Formula: see text] plots even when exchange occurs on the 100 μs timescale. The utility of this strategy is demonstrated on the fast-folding Bacillus stearothermophilus peripheral subunit binding domain that folds with a rate constant ~104 s-1. Here the analysis of 15N CEST data alone results in [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots that contain shallow minima, but the inclusion of visible-state peak positions and restraints on the intrinsic transverse relaxation rates of both states during the analysis of the 15N CEST data results in pronounced minima in the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots and precise exchange parameters even in the fast exchange regime ([Formula: see text]~5). Using this strategy we find that the folding rate constant of PSBD is invariant (~10,500 s-1) from 33.2 to 42.9 °C while the unfolding rates (~70 to ~500 s-1) and unfolded state populations (~0.7 to ~4.3%) increase with temperature. The results presented here show that protein dynamics occurring on the 10 to 104 s-1 timescale can be studied using amide 15N CEST experiments.

Dynamic basis for dA•dGTP and dA•d8OGTP misincorporation via Hoogsteen base pairs

Replicative errors contribute to the genetic diversity needed for evolution but in high frequency can lead to genomic instability. Here, we show that DNA dynamics determine the frequency of misincorporating the A•G mismatch, and altered dynamics explain the high frequency of 8-oxoguanine (8OG) A•8OG misincorporation. NMR measurements revealed that Aanti•Ganti (population (pop.) of >91%) transiently forms sparsely populated and short-lived Aanti+•Gsyn (pop. of ~2% and kex = kforward + kreverse of ~137 s-1) and Asyn•Ganti (pop. of ~6% and kex of ~2,200 s-1) Hoogsteen conformations. 8OG redistributed the ensemble, rendering Aanti•8OGsyn the dominant state. A kinetic model in which Aanti+•Gsyn is misincorporated quantitatively predicted the dA•dGTP misincorporation kinetics by human polymerase β, the pH dependence of misincorporation and the impact of the 8OG lesion. Thus, 8OG increases replicative errors relative to G because oxidation of guanine redistributes the ensemble in favor of the mutagenic Aanti•8OGsyn Hoogsteen state, which exists transiently and in low abundance in the A•G mismatch.

19F NMR Untersuchung des Konformationsaustauschs mehrerer Zustände im synthetischen Neomycin-bindenden Riboschalter

Der synthetische Neomycin‐bindende Riboschalter interagiert mit seinem Liganden Neomycin sowie mit den verwandten Antibiotika Ribostamycin und Paromomycin. Die Bindung dieser Aminoglykoside induziert sehr ähnliche Grundzustandsstrukturen in der RNA, allerdings kann nur Neomycin die Initiierung der Translation effizient unterdrücken. Der molekulare Ursprung dieser Unterschiede wurde auf Unterschiede in der Dynamik der Ligand‐Riboschalter‐Komplexe zurückgeführt. In diesem Artikel kombinieren wir fünf komplementäre fluorbasierte NMR‐Methoden, um die Dynamik der drei Riboschalter‐Komplexe im Sekunden‐ bis Mikrosekundenbereich genau zu quantifizieren. Unsere Daten offenbaren komplexe Austauschprozesse mit bis zu vier strukturell unterschiedlichen Zuständen. Wir interpretieren unsere Ergebnisse in einem Modell, das ein Zusammenspiel zwischen verschiedenen chemischen Gruppen in den Antibiotika und spezifischen Basen im Riboschalter zeigt. Allgemeiner unterstreichen unsere Daten das Potenzial von 19F NMR‐Methoden, komplexe Austauschprozesse mit mehreren angeregten Zuständen zu charakterisieren.

Multi-Site Conformational Exchange in the Synthetic Neomycin-Sensing Riboswitch Studied by 19 F NMR

The synthetic neomycin-sensing riboswitch interacts with its cognate ligand neomycin as well as with the related antibiotics ribostamycin and paromomycin. Binding of these aminoglycosides induces a very similar ground state structure in the RNA, however, only neomycin can efficiently repress translation initiation. The molecular origin of these differences has been traced back to differences in the dynamics of the ligand:riboswitch complexes. Here, we combine five complementary fluorine based NMR methods to accurately quantify seconds to microseconds dynamics in the three riboswitch complexes. Our data reveal complex exchange processes with up to four structurally different states. We interpret our findings in a model that shows an interplay between different chemical groups in the antibiotics and specific bases in the riboswitch. More generally, our data underscore the potential of 19 F NMR methods to characterize complex exchange processes with multiple excited states.

SPICY: a method for single scan rotating frame relaxometry

T _1ρ is an NMR relaxation mode that is sensitive to low frequency molecular motions, making it an especially valuable tool in biomolecular research. Here, we introduce a new method, SPICY, for measuring T_1ρ relaxation times. In contrast to conventional T_1ρ experiments, in which the sequence is repeated many times to determine the T_1ρ time, the SPICY sequence allows determination of T_1ρ within a single scan, shortening the experiment time remarkably. We demonstrate the method using ¹H T_1ρ relaxation dispersion experiments. Additionally, we combine the sequence with spatial encoding to produce 1D images in a single scan. We show that T_1ρ relaxation times obtained using the single scan approach are in good agreement with those obtained using the traditional experiments.

Measuring the Enthalpy of an Individual Hydrogen Bond in a DNA Duplex with Nucleobase Isotope Editing and Variable-Temperature Infrared Spectroscopy

The level of interest in probing the strength of noncovalent interactions in DNA duplexes is high, as these weak forces dictate the range of suprastructures the double helix adopts under different conditions, in turn directly impacting the biological functions and industrial applications of duplexes that require making and breaking them to access the genetic code. However, few experimental tools can measure these weak forces embedded within large biological suprastructures in the native solution environment. Here, we develop experimental methods for detecting the presence of a single noncovalent interaction [a hydrogen bond (H-bond)] within a large DNA duplex in solution and measure its formation enthalpy (ΔH_f). We report that introduction of a H-bond into the TC2═O group from the noncanonical nucleobase 2-aminopurine produces an expected decrease ∼10 ± 0.76 cm^-1 (from ∼1720 cm^-1 in Watson-Crick to ∼1710 cm^-1 in 2-aminopurine), which correlates with an enthalpy of ∼0.93 ± 0.066 kcal/mol for this interaction.

Site-based description of R1ρ relaxation in local reference frames

Nuclear magnetic spin relaxation in the presence of an applied radiofrequency field depends critically on chemical exchange processes that transfer nuclear spins between chemical or conformational environments with distinct resonance frequencies. Characterization of chemical exchange processes in R1ρ relaxation dispersion, CEST, and DEST experiments provides powerful insights into chemical and conformational kinetics of biological macromolecules. The present work reformulates expressions for magnetization evolution and the R1ρ relaxation rate constant by focussing on the orientations of the tilted rotating frames of reference for magnetization components in different sites, by treating the spin-locking field strength as a perturbation to free-precession evolution, and by applying the Homotopy Analysis and Laplace transform methods to approximate solutions to the Bloch-McConnell equations. The results provide an expression for R1ρ that is invariant to the topology of the kinetic scheme, an approximate equation for evolution of spin-locked z-magnetization, and an approach for effective simplification of chemical exchange topologies for 2- and N-site chemical exchange processes. The theoretical approach also provides an accurate approximation for relaxation during a constant-amplitude radiofrequency field in the absence of exchange.

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

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

Isotope Labels Combined with Solution NMR Spectroscopy Make Visible the Invisible Conformations of Small-to-Large RNAs

RNA is central to the proper function of cellular processes important for life on earth and implicated in various medical dysfunctions. Yet, RNA structural biology lags significantly behind that of proteins, limiting mechanistic understanding of RNA chemical biology. Fortunately, solution NMR spectroscopy can probe the structural dynamics of RNA in solution at atomic resolution, opening the door to their functional understanding. However, NMR analysis of RNA, with only four unique ribonucleotide building blocks, suffers from spectral crowding and broad linewidths, especially as RNAs grow in size. One effective strategy to overcome these challenges is to introduce NMR-active stable isotopes into RNA. However, traditional uniform labeling methods introduce scalar and dipolar couplings that complicate the implementation and analysis of NMR measurements. This challenge can be circumvented with selective isotope labeling. In this review, we outline the development of labeling technologies and their application to study biologically relevant RNAs and their complexes ranging in size from 5 to 300 kDa by NMR spectroscopy.

Synthesis of [7-15N]-GTPs for RNA structure and dynamics by NMR spectroscopy

Several isotope-labeling strategies have been developed for the study of RNA by nuclear magnetic resonance (NMR) spectroscopy. Here, we report a combined chemical and enzymatic synthesis of [7-¹⁵N]-guanosine-5′-triphosphates for incorporation into RNA via T7 RNA polymerase-based in vitro transcription. We showcase the utility of these labels to probe both structure and dynamics in two biologically important RNAs.

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

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

The A39G FF domain folds on a volcano-shaped free energy surface via separate pathways

Conformational dynamics play critical roles in protein folding, misfolding, function, misfunction, and aggregation. While detecting and studying the different conformational states populated by protein molecules on their free energy surfaces (FESs) remain a challenge, NMR spectroscopy has emerged as an invaluable experimental tool to explore the FES of a protein, as conformational dynamics can be probed at atomic resolution over a wide range of timescales. Here, we use chemical exchange saturation transfer (CEST) to detect "invisible" minor states on the energy landscape of the A39G mutant FF domain that exhibited "two-state" folding kinetics in traditional experiments. Although CEST has mostly been limited to studies of processes with rates between ∼5 to 300 s^-1 involving sparse states with populations as low as ∼1%, we show that the line broadening that is often associated with minor state dips in CEST profiles can be exploited to inform on additional conformers, with lifetimes an order of magnitude shorter and populations close to 10-fold smaller than what typically is characterized. Our analysis of CEST profiles that exploits the minor state linewidths of the 71-residue A39G FF domain establishes a folding mechanism that can be described in terms of a four-state exchange process between interconverting states spanning over two orders of magnitude in timescale from ∼100 to ∼15,000 μs. A similar folding scheme is established for the wild-type domain as well. The study shows that the folding of this small domain proceeds through a pair of sparse, partially structured intermediates via two discrete pathways on a volcano-shaped FES.

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

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

Mutate-and-chemical-shift-fingerprint (MCSF) to characterize excited states in RNA using NMR spectroscopy

It is important to understand the dynamics and higher energy structures of RNA, called excited states, to achieve better understanding of RNA function. R_1ρ relaxation dispersion NMR spectroscopy (RD) determines chemical shift differences between the most stable, ground state and the short-lived, low-populated excited states. We describe a procedure for deducing the excited state structure from these chemical shift differences using the mutate-and-chemical-shift-fingerprint (MCSF) method, which requires ~2-6 weeks and moderate understanding of NMR and RNA structure. We recently applied the MCSF methodology to elucidate the excited state of microRNA 34a targeting the SIRT1 mRNA and use this example to demonstrate the analysis. The protocol comprises the following steps: (i) determination of the secondary structure of the excited state from RD chemical shift data, (ii) design of trapped excited state RNA, (iii) validation of the excited state structure by NMR, and (iv) MCSF analysis comparing the chemical shifts of the trapped excited state with the RD-derived chemical shift differences. MCSF enables observation of the short-lived RNA structures, which can be functionally and structurally characterized by entrapment.

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

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

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.

Subsecond Time-Resolved Mass Spectrometry in Dynamic Structural Biology

Life at the molecular level is a dynamic world, where the key players-proteins, oligonucleotides, lipids, and carbohydrates-are in a perpetual state of structural flux, shifting rapidly between local minima on their conformational free energy landscapes. The techniques of classical structural biology, X-ray crystallography, structural NMR, and cryo-electron microscopy (cryo-EM), while capable of extraordinary structural resolution, are innately ill-suited to characterize biomolecules in their dynamically active states. Subsecond time-resolved mass spectrometry (MS) provides a unique window into the dynamic world of biological macromolecules, offering the capacity to directly monitor biochemical processes and conformational shifts with a structural dimension provided by the electrospray charge-state distribution, ion mobility, covalent labeling, or hydrogen-deuterium exchange. Over the past two decades, this suite of techniques has provided important insights into the inherently dynamic processes that drive function and pathogenesis in biological macromolecules, including (mis)folding, complexation, aggregation, ligand binding, and enzyme catalysis, among others. This Review provides a comprehensive account of subsecond time-resolved MS and the advances it has enabled in dynamic structural biology, with an emphasis on insights into the dynamic drivers of protein function.

Developments in solution-state NMR yield broader and deeper views of the dynamic ensembles of nucleic acids

Nucleic acids do not fold into a single conformation, and dynamic ensembles are needed to describe their propensities to cycle between different conformations when performing cellular functions. We review recent advances in solution-state nuclear magnetic resonance (NMR) methods and their integration with computational techniques that are improving the ability to probe the dynamic ensembles of DNA and RNA. These include computational approaches for predicting chemical shifts from structure and generating conformational libraries from sequence, measurements of exact nuclear Overhauser effects, development of new probes to study chemical exchange using relaxation dispersion, faster and more sensitive real-time NMR techniques, and new NMR approaches to tackle large nucleic acid assemblies. We discuss how these advances are leading to new mechanistic insights into gene expression and regulation.

Probing RNA Conformational Equilibria within the Functional Cellular Context

Low-abundance short-lived non-native conformations referred to as excited states (ESs) are increasingly observed in vitro and implicated in the folding and biological activities of regulatory RNAs. We developed an approach for assessing the relative abundance of RNA ESs within the functional cellular context. Nuclear magnetic resonance (NMR) spectroscopy was used to estimate the degree to which substitution mutations bias conformational equilibria toward the inactive ES in vitro. The cellular activity of the ES-stabilizing mutants was used as an indirect measure of the conformational equilibria within the functional cellular context. Compensatory mutations that restore the ground-state conformation were used to control for changes in sequence. Using this approach, we show that the ESs of two regulatory RNAs from HIV-1, the transactivation response element (TAR) and the Rev response element (RRE), likely form in cells with abundances comparable to those measured in vitro, and their targeted stabilization may provide an avenue for developing anti-HIV therapeutics.

2'-O-Methylation can increase the abundance and lifetime of alternative RNA conformational states

2'-O-Methyl (Nm) is a highly abundant post-transcriptional RNA modification that plays important biological roles through mechanisms that are not entirely understood. There is evidence that Nm can alter the biological activities of RNAs by biasing the ribose sugar pucker equilibrium toward the C3'-endo conformation formed in canonical duplexes. However, little is known about how Nm might more broadly alter the dynamic ensembles of flexible RNAs containing bulges and internal loops. Here, using NMR and the HIV-1 transactivation response (TAR) element as a model system, we show that Nm preferentially stabilizes alternative secondary structures in which the Nm-modified nucleotides are paired, increasing both the abundance and lifetime of low-populated short-lived excited states by up to 10-fold. The extent of stabilization increased with number of Nm modifications and was also dependent on Mg2+. Through phi-value analysis, the Nm modification also provided rare insights into the structure of the transition state for conformational exchange. Our results suggest that Nm could alter the biological activities of Nm-modified RNAs by modulating their secondary structural ensembles as well as establish the utility of Nm as a tool for the discovery and characterization of RNA excited state conformations.

Slow motions in A·T rich DNA sequence

In free B-DNA, slow (microsecond-to-millisecond) motions that involve equilibrium between Watson-Crick (WC) and Hoogsteen (HG) base-pairing expand the DNA dynamic repertoire that could mediate DNA-protein assemblies. R_1ρ relaxation dispersion NMR methods are powerful tools to capture such slow conformational exchanges in solution using ¹³C/¹⁵ N labelled DNA. Here, these approaches were applied to a dodecamer containing a TTAAA element that was assumed to facilitate nucleosome formation. NMR data and inferred exchange parameters assign HG base pairs as the minor, transient conformers specifically observed in three successive A·T base pairs forming the TAA·TTA segment. The abundance of these HG A·T base pairs can be up to 1.2% which is high compared to what has previously been observed. Data analyses support a scenario in which the three adenines undergo non-simultaneous motions despite their spatial proximity, thus optimising the probability of having one HG base pair in the TAA·TTA segment. Finally, revisiting previous NMR data on H2 resonance linewidths on the basis of our results promotes the idea of there being a special propensity of A·T base pairs in TAA·TTA tracts to adopt HG pairing. In summary, this study provides an example of a DNA functional element submitted to slow conformational exchange. More generally, it strengthens the importance of the role of the DNA sequence in modulating its dynamics, over a nano- to milli-second time scale.

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

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

Reduction in Dynamics of Base pair Opening upon Ligand Binding by the Cocaine-Binding Aptamer

We have used magnetization transfer NMR experiments to measure the exchange rate constant (k_ex) of the imino protons in the unbound, cocaine-bound, and quinine-bound forms of the cocaine-binding DNA aptamer. Both long-stem 1 (MN4) and short-stem 1 (MN19) variants were analyzed, corresponding to structures with a prefolded secondary structure and ligand-induced-folding versions of this aptamer, respectively. The k_ex values were measured as a function of temperature from 5 to 45°C to determine the thermodynamics of the base pair opening for MN4. We find that the base pairs close to the ligand-binding site become stronger upon ligand binding, whereas those located away from the binding site do not strengthen. With the buffer conditions used in this study, we observe imino ¹H signals in MN19 not previously seen, which leads us to conclude that in the free form, both stem 2 and parts of stem 3 are formed and that the base pairs in stem 1 become structured or more rigid upon binding. This is consistent with the k_ex values for MN19 decreasing in both stem 1 and at the ligand-binding site. Based on the temperature dependence of the k_ex values, we find that MN19 is more dynamic than MN4 in the free and both ligand-bound forms. For MN4, ligand-binding results in the reduction of dynamics that are localized to the binding site. These results demonstrate that an aptamer in which the base pairs are preformed also experiences a reduction in dynamics with ligand binding.

Recent advances in solid-state relaxation dispersion techniques

This review describes two rotating-frame (R_1ρ) relaxation dispersion methods, namely the Bloch-McConnell Relaxation Dispersion and the Near-rotary Resonance Relaxation Dispersion, which enable the study of microsecond time-scale conformational fluctuations in the solid state using magic-angle-spinning nuclear magnetic resonance spectroscopy. The goal is to provide the reader with key ideas, experimental descriptions, and practical considerations associated with R_1ρ measurements that are needed for analyzing relaxation dispersion and quantifying conformational exchange. While the focus is on protein motion, many presented concepts can be equally well adapted to study the microsecond time-scale dynamics of other bio- (e.g. lipids, polysaccharides, nucleic acids), organic (e.g. pharmaceutical compounds), or inorganic molecules (e.g., metal organic frameworks). This article summarizes the essential contributions made by recent theoretical and experimental solid-state NMR studies to our understanding of protein motion. Here we discuss recent advances in fast MAS applications that enable the observation and atomic level characterization of sparsely populated conformational states which are otherwise inaccessible for other experimental methods. Such high-energy states are often associated with protein functions such as molecular recognition, ligand binding, or enzymatic catalysis, as well as with disease-related properties such as misfolding and amyloid formation.

Extending the Sensitivity of CEST NMR Spectroscopy to Micro-to-Millisecond Dynamics in Nucleic Acids Using High-Power Radio-Frequency Fields

Biomolecules undergo motions on the micro-to-millisecond timescale to adopt low-populated transient states that play important roles in folding, recognition, and catalysis. NMR techniques, such as Carr-Purcell-Meiboom-Gill (CPMG), chemical exchange saturation transfer (CEST), and R1ρ are the most commonly used methods for characterizing such transitions at atomic resolution under solution conditions. CPMG and CEST are most effective at characterizing motions on the millisecond timescale. While some implementations of the R1ρ experiment are more broadly sensitive to motions on the micro-to-millisecond timescale, they entail the use of selective irradiation schemes and inefficient 1D data acquisition methods. Herein, we show that high-power radio-frequency fields can be used in CEST experiments to extend the sensitivity to faster motions on the micro-to-millisecond timescale. Given the ease of implementing high-power fields in CEST, this should make it easier to characterize micro-to-millisecond dynamics in biomolecules.

Exposing Hidden High-Affinity RNA Conformational States

RNA recognition frequently results in conformational changes that optimize intermolecular binding. As a consequence, the overall binding affinity of RNA to its binding partners depends not only on the intermolecular interactions formed in the bound state but also on the energy cost associated with changing the RNA conformational distribution. Measuring these "conformational penalties" is, however, challenging because bound RNA conformations tend to have equilibrium populations in the absence of the binding partner that fall outside detection by conventional biophysical methods. In this study we employ as a model system HIV-1 TAR RNA and its interaction with the ligand argininamide (ARG), a mimic of TAR's cognate protein binding partner, the transactivator Tat. We use NMR chemical shift perturbations and relaxation dispersion in combination with Bayesian inference to develop a detailed thermodynamic model of coupled conformational change and ligand binding. Starting from a comprehensive 12-state model of the equilibrium, we estimate the energies of six distinct detectable thermodynamic states that are not accessible by currently available methods. Our approach identifies a minimum of four RNA intermediates that differ in terms of the TAR conformation and ARG occupancy. The dominant bound TAR conformation features two bound ARG ligands and has an equilibrium population in the absence of ARG that is below detection limit. Consequently, even though ARG binds to TAR with an apparent overall weak affinity (Kdapp ≈ 0.2 mM), it binds the prefolded conformation with a Kd in the nM range. Our results show that conformational penalties can be major determinants of RNA-ligand binding affinity as well as a source of binding cooperativity, with important implications for a predictive understanding of how RNA is recognized and for RNA-targeted drug discovery.

Probing excited conformational states of nucleic acids by nitrogen CEST NMR spectroscopy

Characterizing low-populated and short-lived excited conformational states has become increasingly important for understanding mechanisms of RNA function. Interconversion between RNA ground and excited conformational states often involves base pairing rearrangements that lead to changes in the hydrogen-bond network. Here, we present two 15N chemical exchange saturation transfer (CEST) NMR experiments that utilize protonated and non-protonated nitrogens, which are key hydrogen-bond donors and acceptors, for characterizing excited conformational states in RNA. We demonstrated these approaches on the B. Cereus fluoride riboswitch, where 15N CEST profiles complement 13C CEST profiles in depicting a potential pathway for ligand-dependent allosteric regulation of the excited conformational state of the fluoride riboswitch.

Demonstration that Small Molecules can Bind and Stabilize Low-abundance Short-lived RNA Excited Conformational States

Many promising RNA drug targets have functions that require the formation of RNA-protein complexes, but inhibiting RNA-protein interactions can prove difficult using small molecules. Regulatory RNAs have been shown to transiently form excited conformational states (ESs) that remodel local aspects of secondary structure. In some cases, the ES conformation has been shown to be inactive and to be poorly recognized by protein binding partners. In these cases, specifically targeting and stabilizing the RNA ES using a small molecule provides a rational structure-based strategy for inhibiting RNA activity. However, this requires that a small molecule discriminates between two conformations of the same RNA to preferentially bind and stabilize the short-lived low-abundance ES relative to the long-lived more abundant ground state (GS). Here, we tested the feasibility of this approach by designing a mutant that inverts the conformational equilibrium of the HIV-1 transactivation response element (TAR) RNA, such that the native GS conformation becomes a low-abundance ES. Using this mutant and NMR chemical shift mapping experiments, we show that argininamide, a ligand mimic of TAR's cognate protein binding partner Tat, is able to restore a native-like conformation by preferentially binding and stabilizing the transient and low-populated ES. A synthetic small molecule optimized to bind the TAR GS also partially stabilized the ES, whereas an aminoglycoside molecule that binds RNAs nonspecifically did not preferentially stabilize the ES to a similar extent. These results support the feasibility of inhibiting RNA activity using small molecules that preferentially bind and stabilize the ES.

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

N⁶-Methyladenosine (m⁶A) is an abundant epitranscriptomic modification that plays important roles in many aspects of RNA metabolism. While m⁶A 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 m⁶A 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 R_1ρ 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 m⁶A alters the hybridization kinetics in two RNA duplexes. The results show that m⁶A minimally impacts the rate constant for duplex dissociation, changing k_off by ∼1-fold but significantly slows the rate of duplex annealing, decreasing k_on 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 Mg²⁺. We propose that rotation of the N⁶-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 m⁶A slows down tRNA selection and more generally suggest that m⁶A may exert cellular functions by reshaping the kinetics of RNA conformational transitions.

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

N6-methyladenosine (m6A) is a ubiquitous RNA post-transcriptional modification found in coding as well as non-coding RNAs. m6A has also been found in viral RNAs where it is proposed to modulate host-pathogen interactions. Two m6A 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 m6A sites have not been observed in other studies mapping m6A 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 m6A 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, m6A is unlikely to substantially alter the conformational properties of the RNA. Our results add to a growing view that the impact of m6A on RNA depends on sequence context and Mg2+.

Direct evidence for (G)O6···H2-N4(C)+ hydrogen bonding in transient G(syn)-C+ and G(syn)-m5C+ Hoogsteen base pairs in duplex DNA from cytosine amino nitrogen off-resonance R1ρ relaxation dispersion measurements

NMR relaxation dispersion studies have shown that Watson-Crick G-C and A-T base pairs in duplex DNA exist in dynamic equilibrium with their Hoogsteen counterparts. Hoogsteen base pairs form through concurrent rotation of the purine base about the glycosidic bond from an anti to a syn conformation and constriction of the C1'-C1' distance across the base pair by ∼2 Å to allow Hoogsteen type hydrogen bonding. Owing to their unique structure, Hoogsteen base pairs can play important roles in DNA recognition, the accommodation, recognition, and repair of DNA damage, and in DNA replication. NMR relaxation dispersion experiments targeting imino nitrogen and protonated base and sugar carbons have provided insights into many structural features of transient Hoogsteen base pairs, including one of two predicted hydrogen bonds involving (G)N7···H-N3(C)⁺ and (A)N7···H-N3(T). Here, through measurement of cytosine amino (N4) R_1ρ relaxation dispersion, we provide direct evidence for the second (G)O6···H₂-N4(C)⁺ hydrogen bond in G(syn)-C⁺ transient Hoogsteen base pairs. The utility of cytosine N4 R_1ρ relaxation dispersion as a new sensitive probe of transient Hoogsteen base pairs, and cytosine dynamics in general, is further demonstrated by measuring G(syn)-C⁺ Hoogsteen exchange near neutral pH and in the context of the naturally occurring DNA modification 5-methyl cytosine (m⁵C), in DNA samples prepared using chemical synthesis and a ¹⁵N labeled m⁵C phosphoramidite.