Efficient Detection of Structure and Dynamics in Unlabeled RNAs: The SELOPE Approach

Quantitative and systematic NMR measurements of sequence-dependent A-T Hoogsteen dynamics uncovers unique conformational specificity in the DNA double helix

The propensities to form lowly-populated short-lived conformations of DNA could vary with sequence, providing an important source of sequence-specificity in biochemical reactions. However, comprehensively measuring how these dynamics vary with sequence is challenging. Using 1H CEST and 13C R 1 ρ NMR, we measured Watson-Crick to Hoogsteen dynamics for an A-T base pair in thirteen trinucleotide sequence contexts. The Hoogsteen population and exchange rate varied 4-fold and 16-fold, respectively, and were dependent on both the 3'- and 5'-neighbors but only weakly dependent on monovalent ion concentration (25 versus 100 mM NaCl) and pH (6.8 versus 8.0). Flexible TA and CA dinucleotide steps exhibited the highest Hoogsteen populations, and their kinetics rates strongly depended on the 3'-neighbor. In contrast, the stiffer AA and GA steps had the lowest Hoogsteen population, and their kinetics were weakly dependent on the 3'-neighbor. The Hoogsteen lifetime was especially short when G-C neighbors flanked the A-T base pair. The Hoogsteen dynamics had a distinct sequence-dependence compared to duplex stability and minor groove width. Thus, our results uncover a unique source of sequence-specificity hidden within the DNA double helix in the form of A-T Hoogsteen dynamics and establish the utility of 1H CEST to quantitively measure sequence-dependent DNA dynamics.

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.

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.

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.

1H NMR Chemical Exchange Techniques Reveal Local and Global Effects of Oxidized Cytosine Derivatives

5-Carboxycytosine (5caC) is a rare epigenetic modification found in nucleic acids of all domains of life. Despite its sparse genomic abundance, 5caC is presumed to play essential regulatory roles in transcription, maintenance and base-excision processes in DNA. In this work, we utilize nuclear magnetic resonance (NMR) spectroscopy to address the effects of 5caC incorporation into canonical DNA strands at multiple pH and temperature conditions. Our results demonstrate that 5caC has a pH-dependent global destabilizing and a base-pair mobility enhancing local impact on dsDNA, albeit without any detectable influence on the ground-state B-DNA structure. Measurement of hybridization thermodynamics and kinetics of 5caC-bearing DNA duplexes highlighted how acidic environment (pH 5.8 and 4.7) destabilizes the double-stranded structure by ∼10–20 kJ mol^–1 at 37 °C when compared to the same sample at neutral pH. Protonation of 5caC results in a lower activation energy for the dissociation process and a higher barrier for annealing. Studies on conformational exchange on the microsecond time scale regime revealed a sharply localized base-pair motion involving exclusively the modified site and its immediate surroundings. By direct comparison with canonical and 5-formylcytosine (5fC)-edited strands, we were able to address the impact of the two most oxidized naturally occurring cytosine derivatives in the genome. These insights on 5caC’s subtle sensitivity to acidic pH contribute to the long-standing questions of its capacity as a substrate in base excision repair processes and its purpose as an independent, stable epigenetic mark.

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.

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.

Production of Structured RNA Fragments by In Vitro Transcription and HPLC Purification

The understanding of the functional importance of RNA has increased enormously in the last decades. This has required research on the RNA molecules themselves, with the concomitant need for obtaining purified RNA samples, such as for structural studies by NMR or other methods. The main method to create labeled and unlabeled RNA, T7 in vitro transcription, suffers from sequence-dependent yield and often low homogeneity for short constructs (<100 nt) and requires laborious purification. Additionally, the design of structured RNA fragments mimicking the structure of a larger biological RNA is often not straightforward. Secondary structure simulations can be used to make reliable predictions about the folding of a particular RNA fragment. In this article, we describe how to design an RNA construct of interest from a larger sequence, and we combine several previously published improvements of the in vitro transcription method, such as the use of 2'-methoxy modifications and dimethyl sulfoxide or the use of tandem repeats, to increase yield and purity of in vitro-transcribed RNA. Together with a high-performance liquid chromatography (HPLC) purification procedure using both reversed-phase ion-pairing and ion-exchange HPLC, we provide a robust protocol to obtain highly pure RNA of short to intermediate length in large quantities. The protocol optimizes yield, especially for RNA starting with nucleotides other than G. At the same time, it is simplified, and the required time is reduced. The protocols described here constitute a versatile pipeline for the production of purified RNA samples and are suitable for users with little experience in liquid chromatography. © 2021 The Authors. Basic Protocol 1: RNA construct design Basic Protocol 2: DNA template production and in vitro transcription Alternate Protocol: Tandem transcription and RNase H cleavage Basic Protocol 3: Reversed-phase ion-pairing HPLC purification Basic Protocol 4: Ion-exchange HPLC purification.

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.

Chemical shift prediction of RNA imino groups: application toward characterizing RNA excited states

NH groups in proteins or nucleic acids are the most challenging target for chemical shift prediction. Here we show that the RNA base pair triplet motif dictates imino chemical shifts in its central base pair. A lookup table is established that links each type of base pair triplet to experimental chemical shifts of the central base pair, and can be used to predict imino chemical shifts of RNAs to remarkable accuracy. Strikingly, the semiempirical method can well interpret the variations of chemical shifts for different base pair triplets, and is even applicable to non-canonical motifs. This finding opens an avenue for predicting chemical shifts of more complicated RNA motifs. Furthermore, we combine the imino chemical shift prediction with NMR relaxation dispersion experiments targeting both ¹⁵N and ¹H^N of the imino group, and verify a previously characterized excited state of P5abc subdomain including an earlier speculated non-native G•G mismatch.

Prediction of chemical shifts is critical for extracting structural and dynamic information from biomolecular NMR data. Here the authors report an RNA imino group chemical shift predictor, showing that the imino chemical shifts of a residue are dictated by the surrounding base pair triplet.

Structural dynamics of double-stranded DNA with epigenome modification

Modification of cytosine plays an important role in epigenetic regulation of gene expression and genome stability. Cytosine is converted to 5-methylcytosine (5mC) by DNA methyltransferase; in turn, 5mC may be oxidized to 5-hydroxymethylcytosine (5hmC) by ten-eleven translocation enzyme. The structural flexibility of DNA is known to affect the binding of proteins to methylated DNA. Here, we have carried out a semi-quantitative analysis of the dynamics of double-stranded DNA (dsDNA) containing various epigenetic modifications by combining data from imino 1H exchange and imino 1H R1ρ relaxation dispersion NMR experiments in a complementary way. Using this approach, we characterized the base-opening (kopen) and base-closing (kclose) rates, facilitating a comparison of the base-opening and -closing process of dsDNA containing cytosine in different states of epigenetic modification. A particularly striking result is the increase in the kopen rate of hemi-methylated dsDNA 5mC/C relative to unmodified or fully methylated dsDNA, indicating that the Watson-Crick base pairs undergo selective destabilization in 5mC/C. Collectively, our findings imply that the epigenetic modulation of cytosine dynamics in dsDNA mediates destabilization of the GC Watson-Crick base pair to allow base-flipping in living cells.

Impact of 5-formylcytosine on the melting kinetics of DNA by 1H NMR chemical exchange

5-Formylcytosine (5fC) is a chemically edited, naturally occurring nucleobase which appears in the context of modified DNA strands. The understanding of the impact of 5fC on dsDNA physical properties is to date limited. In this work, we applied temperature-dependent 1H Chemical Exchange Saturation Transfer (CEST) NMR experiments to non-invasively and site-specifically measure the thermodynamic and kinetic influence of formylated cytosine nucleobase on the melting process involving dsDNA. Incorporation of 5fC within symmetrically positioned CpG sites destabilizes the whole dsDNA structure-as witnessed from the ∼2°C decrease in the melting temperature and 5-10 kJ mol-1 decrease in ΔG°-and affects the kinetic rates of association and dissociation. We observed an up to ∼5-fold enhancement of the dsDNA dissociation and an up to ∼3-fold reduction in ssDNA association rate constants, over multiple temperatures and for several proton reporters. Eyring and van't Hoff analysis proved that the destabilization is not localized, instead all base-pairs are affected and the transition states resembles the single-stranded conformation. These results advance our knowledge about the role of 5fC as a semi-permanent epigenetic modification and assist in the understanding of its interactions with reader proteins.

A robust and versatile method for production and purification of large-scale RNA samples for structural biology

Recent findings in genome-wide transcriptomics revealed that RNAs are involved in almost every biological process, across all domains of life. The characterization of native RNAs of unknown function and structure is particularly challenging due to their typical low abundance in the cell and the inherent sensitivity toward ubiquitous RNA degrading enzymes. Therefore, robust in vitro synthesis and extensive work-up methods are often needed to obtain samples amenable for biochemical, biophysical, and structural studies. Here, we present a protocol that combines the most recent advances in T7 in vitro transcription methodology with reverse phase ion pairing and ion exchange HPLC purification of RNAs for the production of yield-optimized large-scale samples. The method is easy to follow, robust and suitable for users with little or no experience within the field of biochemistry or chromatography. The complete execution of this method, for example, for production of isotopically labeled NMR samples, can be performed in less than a week.

Base-pair conformational switch modulates miR-34a targeting of Sirt1 mRNA

MicroRNAs (miRNAs) regulate the levels of translation of messenger RNAs (mRNAs). At present, the major parameter that can explain the selection of the target mRNA and the efficiency of translation repression is the base pairing between the 'seed' region of the miRNA and its counterpart mRNA¹. Here we use R_1ρ relaxation-dispersion nuclear magnetic resonance² and molecular simulations³ to reveal a dynamic switch-based on the rearrangement of a single base pair in the miRNA-mRNA duplex-that elongates a weak five-base-pair seed to a complete seven-base-pair seed. This switch also causes coaxial stacking of the seed and supplementary helix fitting into human Argonaute 2 protein (Ago2), reminiscent of an active state in prokaryotic Ago^4,5. Stabilizing this transient state leads to enhanced repression of the target mRNA in cells, revealing the importance of this miRNA-mRNA structure. Our observations tie together previous findings regarding the stepwise miRNA targeting process from an initial 'screening' state to an 'active' state, and unveil the role of the RNA duplex beyond the seed in Ago2.

Zinc-finger protein CNBP alters the 3-D structure of lncRNA Braveheart in solution

Long non-coding RNAs (lncRNAs) constitute a significant fraction of the transcriptome, playing important roles in development and disease. However, our understanding of structure-function relationships for this emerging class of RNAs has been limited to secondary structures. Here, we report the 3-D atomistic structural study of epigenetic lncRNA, Braveheart (Bvht), and its complex with CNBP (Cellular Nucleic acid Binding Protein). Using small angle X-ray scattering (SAXS), we elucidate the ensemble of Bvht RNA conformations in solution, revealing that Bvht lncRNA has a well-defined, albeit flexible 3-D structure that is remodeled upon CNBP binding. Our study suggests that CNBP binding requires multiple domains of Bvht and the RHT/AGIL RNA motif. We show that RHT/AGIL, previously shown to interact with CNBP, contains a highly flexible loop surrounded by more ordered helices. As one of the largest RNA-only 3-D studies, the work lays the foundation for future structural studies of lncRNA-protein complexes.

RNA Dynamics by NMR Spectroscopy

An ever-increasing number of functional RNAs require a mechanistic understanding. RNA function relies on changes in its structure, so-called dynamics. To reveal dynamic processes and higher energy structures, new NMR methods have been developed to elucidate these dynamics in RNA with atomic resolution. In this Review, we provide an introduction to dynamics novices and an overview of methods that access most dynamic timescales, from picoseconds to hours. Examples are provided as well as insight into theory, data acquisition and analysis for these different methods. Using this broad spectrum of methodology, unprecedented detail and invisible structures have been obtained and are reviewed here. RNA, though often more complicated and therefore neglected, also provides a great system to study structural changes, as these RNA structural changes are more easily defined-Lego like-than in proteins, hence the numerous revelations of RNA excited states.

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

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

Efficient Detection of Structure and Dynamics in Unlabeled RNAs: The SELOPE Approach

The knowledge of structure and dynamics is crucial to explain the function of RNAs. While nuclear magnetic resonance (NMR) is well suited to probe these for complex biomolecules, it requires expensive, isotopically labeled samples, and long measurement times. Here we present SELOPE, a new robust, proton-only NMR method that allows us to obtain site-specific overview of structure and dynamics in an entire RNA molecule using an unlabeled sample. SELOPE simplifies assignment and allows for cost-effective screening of the response of nucleic acids to physiological changes (e.g. ion concentration) or screening of drugs in a high throughput fashion. This single technique allows us to probe an unprecedented range of exchange time scales (the whole μs to ms motion range) with increased sensitivity, surpassing all current experiments to detect chemical exchange. For the first time we could describe an RNA excited state using an unlabeled RNA.