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.

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.

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.

Hydrogen bonding in duplex DNA probed by DNP enhanced solid-state NMR N-H bond length measurements

Numerous biological processes and mechanisms depend on details of base pairing and hydrogen bonding in DNA. Hydrogen bonds are challenging to quantify by X-ray crystallography and cryo-EM due to difficulty of visualizing hydrogen atom locations but can be probed with site specificity by NMR spectroscopy in solution and the solid state with the latter particularly suited to large, slowly tumbling DNA complexes. Recently, we showed that low-temperature dynamic nuclear polarization (DNP) enhanced solid-state NMR is a valuable tool for distinguishing Hoogsteen base pairs (bps) from canonical Watson-Crick bps in various DNA systems under native-like conditions. Here, using a model 12-mer DNA duplex containing two central adenine-thymine (A-T) bps in either Watson-Crick or Hoogsteen confirmation, we demonstrate DNP solid-state NMR measurements of thymine N3-H3 bond lengths, which are sensitive to details of N-H···N hydrogen bonding and permit hydrogen bonds for the two bp conformers to be systematically compared within the same DNA sequence context. For this DNA duplex, effectively identical TN3-H3 bond lengths of 1.055 ± 0.011 Å and 1.060 ± 0.011 Å were found for Watson-Crick A-T and Hoogsteen A (syn)-T base pairs, respectively, relative to a reference amide bond length of 1.015 ± 0.010 Å determined for N-acetyl-valine under comparable experimental conditions. Considering that prior quantum chemical calculations which account for zero-point motions predict a somewhat longer effective peptide N-H bond length of 1.041 Å, in agreement with solution and solid-state NMR studies of peptides and proteins at ambient temperature, to facilitate direct comparisons with these earlier studies TN3-H3 bond lengths for the DNA samples can be readily scaled appropriately to yield 1.083 Å and 1.087 Å for Watson-Crick A-T and Hoogsteen A (syn)-T bps, respectively, relative to the 1.041 Å reference peptide N-H bond length. Remarkably, in the context of the model DNA duplex, these results indicate that there are no significant differences in N-H···N A-T hydrogen bonds between Watson-Crick and Hoogsteen bp conformers. More generally, high precision measurements of N-H bond lengths by low-temperature DNP solid-state NMR based methods are expected to facilitate detailed comparative analysis of hydrogen bonding for a range of DNA complexes and base pairing environments.

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.

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.

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.

Turing, von Neumann, and the computational architecture of biological machines

In the mid-1930s, the English mathematician and logician Alan Turing invented an imaginary machine which could emulate the process of manipulating finite symbolic configurations by human computers. His machine launched the field of computer science and provided a foundation for the modern-day programmable computer. A decade later, building on Turing's machine, the American-Hungarian mathematician John von Neumann invented an imaginary self-reproducing machine capable of open-ended evolution. Through his machine, von Neumann answered one of the deepest questions in Biology: Why is it that all living organisms carry a self-description in the form of DNA? The story behind how two pioneers of computer science stumbled on the secret of life many years before the discovery of the DNA double helix is not well known, not even to biologists, and you will not find it in biology textbooks. Yet, the story is just as relevant today as it was eighty years ago: Turing and von Neumann left a blueprint for studying biological systems as if they were computing machines. This approach may hold the key to answering many remaining questions in Biology and could even lead to advances in computer science.

RNA conformational propensities determine cellular activity

Cellular processes are the product of interactions between biomolecules, which associate to form biologically active complexes1. These interactions are mediated by intermolecular contacts, which if disrupted, lead to alterations in cell physiology. Nevertheless, the formation of intermolecular contacts nearly universally requires changes in the conformations of the interacting biomolecules. As a result, binding affinity and cellular activity crucially depend both on the strength of the contacts and on the inherent propensities to form binding-competent conformational states2,3. Thus, conformational penalties are ubiquitous in biology and must be known in order to quantitatively model binding energetics for protein and nucleic acid interactions4,5. However, conceptual and technological limitations have hindered our ability to dissect and quantitatively measure how conformational propensities affect cellular activity. Here we systematically altered and determined the propensities for forming the protein-bound conformation of HIV-1 TAR RNA. These propensities quantitatively predicted the binding affinities of TAR to the RNA-binding region of the Tat protein and predicted the extent of HIV-1 Tat-dependent transactivation in cells. Our results establish the role of ensemble-based conformational propensities in cellular activity and reveal an example of a cellular process driven by an exceptionally rare and short-lived RNA conformational state.

mRNA interactions with disordered regions control protein activity

The cytoplasm is compartmentalized into different translation environments. mRNAs use their 3'UTRs to localize to distinct cytoplasmic compartments, including TIS granules (TGs). Many transcription factors, including MYC, are translated in TGs. It was shown that translation of proteins in TGs enables the formation of protein complexes that cannot be established when these proteins are translated in the cytosol, but the mechanism is poorly understood. Here we show that MYC protein complexes that involve binding to the intrinsically disordered region (IDR) of MYC are only formed when MYC is translated in TGs. TG-dependent protein complexes require TG-enriched mRNAs for assembly. These mRNAs bind to a new and widespread RNA-binding domain in neutral or negatively charged IDRs in several transcription factors, including MYC. RNA-IDR interaction changes the conformational ensemble of the IDR, enabling the formation of MYC protein complexes that act in the nucleus and control functions that cannot be accomplished by cytosolically-translated MYC. We propose that certain mRNAs have IDR chaperone activity as they control IDR conformations. In addition to post-translational modifications, we found a novel mode of protein activity regulation. Since RNA-IDR interactions are prevalent, we suggest that mRNA-dependent control of protein functional states is widespread.

RBM45 is an m6A-binding protein that affects neuronal differentiation and the splicing of a subset of mRNAs

N⁶-methyladenosine (m⁶A) is deposited co-transcriptionally on thousands of cellular mRNAs and plays important roles in mRNA processing and cellular function. m⁶A is particularly abundant within the brain and is critical for neurodevelopment. However, the mechanisms through which m⁶A contributes to brain development are incompletely understood. RBM45 acts as an m⁶A-binding protein that is highly expressed during neurodevelopment. We find that RBM45 binds to thousands of cellular RNAs, predominantly within intronic regions. Rbm45 depletion disrupts the constitutive splicing of a subset of target pre-mRNAs, leading to altered mRNA and protein levels through both m⁶A-dependent and m⁶A-independent mechanisms. Finally, we find that RBM45 is necessary for neuroblastoma cell differentiation and that its depletion impacts the expression of genes involved in several neurodevelopmental signaling pathways. Altogether, our findings show a role for RBM45 in controlling mRNA processing and neuronal differentiation, mediated in part by the recognition of methylated RNA.

Choi et al. identify RBM45 as an m⁶A-binding protein enriched in the developing brain. RBM45 binds to thousands of cellular RNAs, primarily within introns, and regulates constitutive splicing of target transcripts. Loss of RBM45 causes altered expression of neurodevelopmental genes and defects in the proliferation and differentiation of neuroblastoma cells.

Probing Watson-Crick and Hoogsteen base pairing in duplex DNA using dynamic nuclear polarization solid-state NMR spectroscopy

The majority of base pairs in double-stranded DNA exist in the canonical Watson-Crick geometry. However, they can also adopt alternate Hoogsteen conformations in various complexes of DNA with proteins and small molecules, which are key for biological function and mechanism. While detection of Hoogsteen base pairs in large DNA complexes and assemblies poses considerable challenges for traditional structural biology techniques, we show here that multidimensional dynamic nuclear polarization-enhanced solid-state NMR can serve as a unique spectroscopic tool for observing and distinguishing Watson-Crick and Hoogsteen base pairs in a broad range of DNA systems based on characteristic NMR chemical shifts and internuclear dipolar couplings. We illustrate this approach using a model 12-mer DNA duplex, free and in complex with the antibiotic echinomycin, which features two central adenine-thymine base pairs with Watson-Crick and Hoogsteen geometry, respectively, and subsequently extend it to the ∼200 kDa Widom 601 DNA nucleosome core particle.

Revealing A-T and G-C Hoogsteen base pairs in stressed protein-bound duplex DNA

Watson–Crick base pairs (bps) are the fundamental unit of genetic information and the building blocks of the DNA double helix. However, A-T and G-C can also form alternative ‘Hoogsteen’ bps, expanding the functional complexity of DNA. We developed ‘Hoog-finder’, which uses structural fingerprints to rapidly screen Hoogsteen bps, which may have been mismodeled as Watson–Crick in crystal structures of protein–DNA complexes. We uncovered 17 Hoogsteen bps, 7 of which were in complex with 6 proteins never before shown to bind Hoogsteen bps. The Hoogsteen bps occur near mismatches, nicks and lesions and some appear to participate in recognition and damage repair. Our results suggest a potentially broad role for Hoogsteen bps in stressed regions of the genome and call for a community-wide effort to identify these bps in current and future crystal structures of DNA and its complexes.

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.

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.

Understanding the characteristics of nonspecific binding of drug-like compounds to canonical stem-loop RNAs and their implications for functional cellular assays

Identifying small molecules that selectively bind an RNA target while discriminating against all other cellular RNAs is an important challenge in RNA-targeted drug discovery. Much effort has been directed toward identifying drug-like small molecules that minimize electrostatic and stacking interactions that lead to nonspecific binding of aminoglycosides and intercalators to many stem-loop RNAs. Many such compounds have been reported to bind RNAs and inhibit their cellular activities. However, target engagement and cellular selectivity assays are not routinely performed, and it is often unclear whether functional activity directly results from specific binding to the target RNA. Here, we examined the propensities of three drug-like compounds, previously shown to bind and inhibit the cellular activities of distinct stem-loop RNAs, to bind and inhibit the cellular activities of two unrelated HIV-1 stem-loop RNAs: the transactivation response element (TAR) and the rev response element stem IIB (RREIIB). All compounds bound TAR and RREIIB in vitro, and two inhibited TAR-dependent transactivation and RRE-dependent viral export in cell-based assays while also exhibiting off-target interactions consistent with nonspecific activity. A survey of X-ray and NMR structures of RNA-small molecule complexes revealed that aminoglycosides and drug-like molecules form hydrogen bonds with functional groups commonly accessible in canonical stem-loop RNA motifs, in contrast to ligands that specifically bind riboswitches. Our results demonstrate that drug-like molecules can nonspecifically bind stem-loop RNAs most likely through hydrogen bonding and electrostatic interactions and reinforce the importance of assaying for off-target interactions and RNA selectivity in vitro and in cells when assessing novel RNA-binders.

DNA mismatches reveal conformational penalties in protein-DNA recognition

Transcription factors recognize specific genomic sequences to regulate complex gene-expression programs. Although it is well-established that transcription factors bind to specific DNA sequences using a combination of base readout and shape recognition, some fundamental aspects of protein-DNA binding remain poorly understood1,2. Many DNA-binding proteins induce changes in the structure of the DNA outside the intrinsic B-DNA envelope. However, how the energetic cost that is associated with distorting the DNA contributes to recognition has proven difficult to study, because the distorted DNA exists in low abundance in the unbound ensemble3-9. Here we use a high-throughput assay that we term SaMBA (saturation mismatch-binding assay) to investigate the role of DNA conformational penalties in transcription factor-DNA recognition. In SaMBA, mismatched base pairs are introduced to pre-induce structural distortions in the DNA that are much larger than those induced by changes in the Watson-Crick sequence. Notably, approximately 10% of mismatches increased transcription factor binding, and for each of the 22 transcription factors that were examined, at least one mismatch was found that increased the binding affinity. Mismatches also converted non-specific sites into high-affinity sites, and high-affinity sites into 'super sites' that exhibit stronger affinity than any known canonical binding site. Determination of high-resolution X-ray structures, combined with nuclear magnetic resonance measurements and structural analyses, showed that many of the DNA mismatches that increase binding induce distortions that are similar to those induced by protein binding-thus prepaying some of the energetic cost incurred from deforming the DNA. Our work indicates that conformational penalties are a major determinant of protein-DNA recognition, and reveals mechanisms by which mismatches can recruit transcription factors and thus modulate replication and repair activities in the cell10,11.

Hoogsteen base pairs increase the susceptibility of double-stranded DNA to cytotoxic damage

As the Watson-Crick faces of nucleobases are protected in dsDNA, it is commonly assumed that deleterious alkylation damage to the Watson-Crick faces of nucleobases predominantly occurs when DNA becomes single-stranded during replication and transcription. However, damage to the Watson-Crick faces of nucleobases has been reported in dsDNA in vitro through mechanisms that are not understood. In addition, the extent of protection from methylation damage conferred by dsDNA relative to ssDNA has not been quantified. Watson-Crick base pairs in dsDNA exist in dynamic equilibrium with Hoogsteen base pairs that expose the Watson-Crick faces of purine nucleobases to solvent. Whether this can influence the damage susceptibility of dsDNA remains unknown. Using dot-blot and primer extension assays, we measured the susceptibility of adenine-N1 to methylation by dimethyl sulfate (DMS) when in an A-T Watson-Crick versus Hoogsteen conformation. Relative to unpaired adenines in a bulge, Watson-Crick A-T base pairs in dsDNA only conferred ∼130-fold protection against adenine-N1 methylation, and this protection was reduced to ∼40-fold for A(syn)-T Hoogsteen base pairs embedded in a DNA-drug complex. Our results indicate that Watson-Crick faces of nucleobases are accessible to alkylating agents in canonical dsDNA and that Hoogsteen base pairs increase this accessibility. Given the higher abundance of dsDNA relative to ssDNA, these results suggest that dsDNA could be a substantial source of cytotoxic damage. The work establishes DMS probing as a method for characterizing A(syn)-T Hoogsteen base pairs in vitro and also lays the foundation for a sequencing approach to map A(syn)-T Hoogsteen and unpaired adenines genome-wide in vivo.

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.

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.

Environmental Effects on Guanine-Thymine Mispair Tautomerization Explored with Quantum Mechanical/Molecular Mechanical Free Energy Simulations

DNA bases can adopt energetically unfavorable tautomeric forms that enable the formation of Watson-Crick-like (WC-like) mispairs, which have been proposed to give rise to spontaneous mutations in DNA and misincorporation errors in DNA replication and translation. Previous NMR and computational studies have indicated that the population of WC-like guanine-thymine (G-T) mispairs depends on the environment, such as the local nucleic acid sequence and solvation. To investigate these environmental effects, herein G-T mispair tautomerization processes are studied computationally in aqueous solution, in A-form and B-form DNA duplexes, and within the active site of a DNA polymerase λ variant. The wobble G-T (wG-T), WC-like G-T*, and WC-like G*-T forms are considered, where * indicates the enol tautomer of the base. The minimum free energy paths for the tautomerization from the wG-T to the WC-like G-T* and from the WC-like G-T* to the WC-like G*-T are computed with mixed quantum mechanical/molecular mechanical (QM/MM) free energy simulations. The reaction free energies and free energy barriers are found to be significantly influenced by the environment. The wG-T→G-T* tautomerization is predicted to be endoergic in aqueous solution and the DNA duplexes but slightly exoergic in the polymerase, with Arg517 and Asn513 providing electrostatic stabilization of G-T*. The G-T*→G*-T tautomerization is also predicted to be slightly more thermodynamically favorable in the polymerase relative to these DNA duplexes. These simulations are consistent with an experimentally driven kinetic misincorporation model suggesting that G-T mispair tautomerization occurs in the ajar polymerase conformation or concertedly with the transition from the ajar to the closed polymerase conformation. Furthermore, the order of the associated two proton transfer reactions is predicted to be different in the polymerase than in aqueous solution and the DNA duplexes. These studies highlight the impact of the environment on the thermodynamics, kinetics, and fundamental mechanisms of G-T mispair tautomerization, which plays a role in a wide range of biochemically important processes.

The roles of structural dynamics in the cellular functions of RNAs

RNAs fold into 3D structures that range from simple helical elements to complex tertiary structures and quaternary ribonucleoprotein assemblies. The functions of many regulatory RNAs depend on how their 3D structure changes in response to a diverse array of cellular conditions. In this Review, we examine how the structural characterization of RNA as dynamic ensembles of conformations, which form with different probabilities and at different timescales, is improving our understanding of RNA function in cells. We discuss the mechanisms of gene regulation by microRNAs, riboswitches, ribozymes, post-transcriptional RNA modifications and RNA-binding proteins, and how the cellular environment and processes such as liquid-liquid phase separation may affect RNA folding and activity. The emerging RNA-ensemble-function paradigm is changing our perspective and understanding of RNA regulation, from in vitro to in vivo and from descriptive to predictive.

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.

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.

Hidden Structural Modules in a Cooperative RNA Folding Transition

Large-scale, cooperative rearrangements underlie the functions of RNA in RNA-protein machines and gene regulation. To understand how such rearrangements are orchestrated, we used high-throughput chemical footprinting to dissect a seemingly concerted rearrangement in P5abc RNA, a paradigm of RNA folding studies. With mutations that systematically disrupt or restore putative structural elements, we found that this transition reflects local folding of structural modules, with modest and incremental cooperativity that results in concerted behavior. First, two distant secondary structure changes are coupled through a bridging three-way junction and Mg²⁺-dependent tertiary structure. Second, long-range contacts are formed between modules, resulting in additional cooperativity. Given the sparseness of RNA tertiary contacts after secondary structure formation, we expect that modular folding and incremental cooperativity are generally important for specifying functional structures while also providing productive kinetic paths to these structures. Additionally, we expect our approach to be useful for uncovering modularity in other complex RNAs.

The RRM of the kRNA-editing protein TbRGG2 uses multiple surfaces to bind and remodel RNA

Kinetoplastid RNA (kRNA) editing takes place in the mitochondria of kinetoplastid protists and creates translatable mRNAs by uridine insertion/deletion. Extensively edited (pan-edited) transcripts contain quadruplex forming guanine stretches, which must be remodeled to promote uridine insertion/deletion. Here we show that the RRM domain of the essential kRNA-editing factor TbRGG2 binds poly(G) and poly(U) RNA and can unfold both. A region C-terminal to the RRM mediates TbRGG2 dimerization, enhancing RNA binding. A RRM-U4 RNA structure reveals a unique RNA-binding mechanism in which the two RRMs of the dimer employ aromatic residues outside the canonical RRM RNA-binding motifs to encase and wrench open the RNA, while backbone atoms specify the uridine bases. Notably, poly(G) RNA is bound via a different binding surface. Thus, these data indicate that TbRGG2 RRM can bind and remodel several RNA substrates suggesting how it might play multiple roles in the kRNA editing process.

Why are Hoogsteen base pairs energetically disfavored in A-RNA compared to B-DNA?

A(syn)-U/T and G(syn)-C⁺ Hoogsteen (HG) base pairs (bps) are energetically more disfavored relative to Watson–Crick (WC) bps in A-RNA as compared to B-DNA by >1 kcal/mol for reasons that are not fully understood. Here, we used NMR spectroscopy, optical melting experiments, molecular dynamics simulations and modified nucleotides to identify factors that contribute to this destabilization of HG bps in A-RNA. Removing the 2′-hydroxyl at single purine nucleotides in A-RNA duplexes did not stabilize HG bps relative to WC. In contrast, loosening the A-form geometry using a bulge in A-RNA reduced the energy cost of forming HG bps at the flanking sites to B-DNA levels. A structural and thermodynamic analysis of purine-purine HG mismatches reveals that compared to B-DNA, the A-form geometry disfavors syn purines by 1.5–4 kcal/mol due to sugar-backbone rearrangements needed to sterically accommodate the syn base. Based on MD simulations, an additional penalty of 3–4 kcal/mol applies for purine-pyrimidine HG bps due to the higher energetic cost associated with moving the bases to form hydrogen bonds in A-RNA versus B-DNA. These results provide insights into a fundamental difference between A-RNA and B-DNA duplexes with important implications for how they respond to damage and post-transcriptional modifications.

Dynamic ensemble of HIV-1 RRE stem IIB reveals non-native conformations that disrupt the Rev-binding site

The HIV-1 Rev response element (RRE) RNA element mediates the nuclear export of intron containing viral RNAs by forming an oligomeric complex with the viral protein Rev. Stem IIB and nearby stem II three-way junction nucleate oligomerization through cooperative binding of two Rev molecules. Conformational flexibility at this RRE region has been shown to be important for Rev binding. However, the nature of the flexibility has remained elusive. Here, using NMR relaxation dispersion, including a new strategy for directly observing transient conformational states in large RNAs, we find that stem IIB alone or when part of the larger RREII three-way junction robustly exists in dynamic equilibrium with non-native excited state (ES) conformations that have a combined population of ∼20%. The ESs disrupt the Rev-binding site by changing local secondary structure, and their stabilization via point substitution mutations decreases the binding affinity to the Rev arginine-rich motif (ARM) by 15- to 80-fold. The ensemble clarifies the conformational flexibility observed in stem IIB, reveals long-range conformational coupling between stem IIB and the three-way junction that may play roles in cooperative Rev binding, and also identifies non-native RRE conformational states as new targets for the development of anti-HIV therapeutics.

Infrared Spectroscopic Observation of a G-C+ Hoogsteen Base Pair in the DNA:TATA-Box Binding Protein Complex Under Solution Conditions

Hoogsteen DNA base pairs (bps) are an alternative base pairing to canonical Watson-Crick bps and are thought to play important biochemical roles. Hoogsteen bps have been reported in a handful of X-ray structures of protein-DNA complexes. However, there are several examples of Hoogsteen bps in crystal structures that form Watson-Crick bps when examined under solution conditions. Furthermore, Hoogsteen bps can sometimes be difficult to resolve in DNA:protein complexes by X-ray crystallography due to ambiguous electron density and by solution-state NMR spectroscopy due to size limitations. Here, using infrared spectroscopy, we report the first direct solution-state observation of a Hoogsteen (G-C⁺ ) bp in a DNA:protein complex under solution conditions with specific application to DNA-bound TATA-box binding protein. These results support a previous assignment of a G-C⁺ Hoogsteen bp in the complex, and indicate that Hoogsteen bps do indeed exist under solution conditions in DNA:protein complexes.

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.

Structural and Functional Characterization of Sulfonium Carbon-Oxygen Hydrogen Bonding in the Deoxyamino Sugar Methyltransferase TylM1

The N-methyltransferase TylM1 from Streptomyces fradiae catalyzes the final step in the biosynthesis of the deoxyamino sugar mycaminose, a substituent of the antibiotic tylosin. The high-resolution crystal structure of TylM1 bound to the methyl donor S-adenosylmethionine (AdoMet) illustrates a network of carbon-oxygen (CH···O) hydrogen bonds between the substrate's sulfonium cation and residues within the active site. These interactions include hydrogen bonds between the methyl and methylene groups of the AdoMet sulfonium cation and the hydroxyl groups of Tyr14 and Ser120 in the enzyme. To examine the functions of these interactions, we generated Tyr14 to phenylalanine (Y14F) and Ser120 to alanine (S120A) mutations to selectively ablate the CH···O hydrogen bonding to AdoMet. The TylM1 S120A mutant exhibited a modest decrease in its catalytic efficiency relative to that of the wild type (WT) enzyme, whereas the Y14F mutation resulted in an approximately 30-fold decrease in catalytic efficiency. In contrast, site-specific substitution of Tyr14 by the noncanonical amino acid p-aminophenylalanine partially restored activity comparable to that of the WT enzyme. Correlatively, quantum mechanical calculations of the activation barrier energies of WT TylM1 and the Tyr14 mutants suggest that substitutions that abrogate hydrogen bonding with the AdoMet methyl group impair methyl transfer. Together, these results offer insights into roles of CH···O hydrogen bonding in modulating the catalytic efficiency of TylM1.

Characterizing Watson-Crick versus Hoogsteen Base Pairing in a DNA-Protein Complex Using Nuclear Magnetic Resonance and Site-Specifically 13C- and 15N-Labeled DNA

A( syn)-T and G( syn)-C⁺ Hoogsteen base pairs in protein-bound DNA duplexes can be difficult to resolve by X-ray crystallography due to ambiguous electron density and by nuclear magnetic resonance (NMR) spectroscopy due to poor chemical shift dispersion and size limitations with solution-state NMR spectroscopy. Here we describe an NMR strategy for characterizing Hoogsteen base pairs in protein-DNA complexes, which relies on site-specifically incorporating ¹³C- and ¹⁵N-labeled nucleotides into DNA duplexes for unambiguous resonance assignment and to improve spectral resolution. The approach was used to resolve the conformation of an A-T base pair in a crystal structure of an ∼43 kDa complex between a 34 bp duplex DNA and the integration host factor (IHF) protein. In the crystal structure (Protein Data Bank entry 1IHF ), this base pair adopts an unusual Hoogsteen conformation with a distorted sugar backbone that is accommodated by a nearby nick used to aid in crystallization. The NMR chemical shifts and interproton nuclear Overhauser effects indicate that this base pair predominantly adopts a Watson-Crick conformation in the intact DNA-IHF complex under solution conditions. Consistent with these NMR findings, substitution of 7-deazaadenine at this base pair resulted in only a small (∼2-fold) decrease in the IHF-DNA binding affinity. The NMR strategy provides a new approach for resolving crystallographic ambiguity and more generally for studying the structure and dynamics of protein-DNA complexes in solution.

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

A potentially abundant junctional RNA motif stabilized by m6A and Mg2

N6-Methyladenosine (m6A) is an abundant post-transcriptional RNA modification that influences multiple aspects of gene expression. In addition to recruiting proteins, m6A can modulate RNA function by destabilizing base pairing. Here, we show that when neighbored by a 5' bulge, m6A stabilizes m6A-U base pairs, and global RNA structure by ~1 kcal mol-1. The bulge most likely provides the flexibility needed to allow optimal stacking between the methyl group and 3' neighbor through a conformation that is stabilized by Mg2+. A bias toward this motif can help explain the global impact of methylation on RNA structure in transcriptome-wide studies. While m6A embedded in duplex RNA is poorly recognized by the YTH domain reader protein and m6A antibodies, both readily recognize m6A in this newly identified motif. The results uncover potentially abundant and functional m6A motifs that can modulate the epitranscriptomic structure landscape with important implications for the interpretation of transcriptome-wide data.

Atomic structures of excited state A-T Hoogsteen base pairs in duplex DNA by combining NMR relaxation dispersion, mutagenesis, and chemical shift calculations

NMR relaxation dispersion studies indicate that in canonical duplex DNA, Watson-Crick base pairs (bps) exist in dynamic equilibrium with short-lived low abundance excited state Hoogsteen bps. N1-methylated adenine (m¹A) and guanine (m¹G) are naturally occurring forms of damage that stabilize Hoogsteen bps in duplex DNA. NMR dynamic ensembles of DNA duplexes with m¹A-T Hoogsteen bps reveal significant changes in sugar pucker and backbone angles in and around the Hoogsteen bp, as well as kinking of the duplex towards the major groove. Whether these structural changes also occur upon forming excited state Hoogsteen bps in unmodified duplexes remains to be established because prior relaxation dispersion probes provided limited information regarding the sugar-backbone conformation. Here, we demonstrate measurements of C3' and C4' spin relaxation in the rotating frame (R1ρ) in uniformly ¹³C/¹⁵N labeled DNA as sensitive probes of the sugar-backbone conformation in DNA excited states. The chemical shifts, combined with structure-based predictions using an automated fragmentation quantum mechanics/molecular mechanics method, show that the dynamic ensemble of DNA duplexes containing m¹A-T Hoogsteen bps accurately model the excited state Hoogsteen conformation in two different sequence contexts. Formation of excited state A-T Hoogsteen bps is accompanied by changes in sugar-backbone conformation that allow the flipped syn adenine to form hydrogen-bonds with its partner thymine and this in turn results in overall kinking of the DNA toward the major groove. Results support the assignment of Hoogsteen bps as the excited state observed in canonical duplex DNA, provide an atomic view of DNA dynamics linked to formation of Hoogsteen bps, and lay the groundwork for a potentially general strategy for solving structures of nucleic acid excited states.

Dynamic basis for dG•dT misincorporation via tautomerization and ionization

Tautomeric and anionic Watson-Crick-like mismatches have important roles in replication and translation errors through mechanisms that are not fully understood. Here, using NMR relaxation dispersion, we resolve a sequence-dependent kinetic network connecting G•T/U wobbles with three distinct Watson-Crick mismatches: two rapidly exchanging tautomeric species (Genol•T/UG•Tenol/Uenol; population less than 0.4%) and one anionic species (G•T-/U-; population around 0.001% at neutral pH). The sequence-dependent tautomerization or ionization step was inserted into a minimal kinetic mechanism for correct incorporation during replication after the initial binding of the nucleotide, leading to accurate predictions of the probability of dG•dT misincorporation across different polymerases and pH conditions and for a chemically modified nucleotide, and providing mechanisms for sequence-dependent misincorporation. Our results indicate that the energetic penalty for tautomerization and/or ionization accounts for an approximately 10-2 to 10-3-fold discrimination against misincorporation, which proceeds primarily via tautomeric dGenol•dT and dG•dTenol, with contributions from anionic dG•dT- dominant at pH 8.4 and above or for some mutagenic nucleotides.

The Mechanism of HdeA Unfolding and Chaperone Activation

HdeA is a periplasmic chaperone that is rapidly activated upon shifting the pH to acidic conditions. This activation is thought to involve monomerization of HdeA. There is evidence that monomerization and partial unfolding allow the chaperone to bind to proteins denatured by low pH, thereby protecting them from aggregation. We analyzed the acid-induced unfolding of HdeA using NMR spectroscopy and fluorescence measurements, and obtained experimental evidence suggesting a complex mechanism in HdeA's acid-induced unfolding pathway, as previously postulated from molecular dynamics simulations. Counterintuitively, dissociation constant measurements show a stabilization of the HdeA dimer upon exposure to mildly acidic conditions. We provide experimental evidence that protonation of Glu37, a glutamate residue embedded in a hydrophobic pocket of HdeA, is important in controlling HdeA stabilization and thus the acid activation of this chaperone. Our data also reveal a sharp transition from folded dimer to unfolded monomer between pH3 and pH 2, and suggest the existence of a low-populated, partially folded intermediate that could assist in chaperone activation or function. Overall, this study provides a detailed experimental investigation into the mechanism by which HdeA unfolds and activates.

Insights into Watson-Crick/Hoogsteen breathing dynamics and damage repair from the solution structure and dynamic ensemble of DNA duplexes containing m1A

In the canonical DNA double helix, Watson–Crick (WC) base pairs (bps) exist in dynamic equilibrium with sparsely populated (∼0.02–0.4%) and short-lived (lifetimes ∼0.2–2.5 ms) Hoogsteen (HG) bps. To gain insights into transient HG bps, we used solution-state nuclear magnetic resonance spectroscopy, including measurements of residual dipolar couplings and molecular dynamics simulations, to examine how a single HG bp trapped using the N1-methylated adenine (m¹A) lesion affects the structural and dynamic properties of two duplexes. The solution structure and dynamic ensembles of the duplexes reveals that in both cases, m¹A forms a m¹A•T HG bp, which is accompanied by local and global structural and dynamic perturbations in the double helix. These include a bias toward the BI backbone conformation; sugar repuckering, major-groove directed kinking (∼9°); and local melting of neighboring WC bps. These results provide atomic insights into WC/HG breathing dynamics in unmodified DNA duplexes as well as identify structural and dynamic signatures that could play roles in m¹A recognition and repair.

Resolving sugar puckers in RNA excited states exposes slow modes of repuckering dynamics

Recent studies have shown that RNAs exist in dynamic equilibrium with short-lived low-abundance 'excited states' that form by reshuffling base pairs in and around non-canonical motifs. These conformational states are proposed to be rich in non-canonical motifs and to play roles in the folding and regulatory functions of non-coding RNAs but their structure proves difficult to characterize given their transient nature. Here, we describe an approach for determining sugar pucker conformation in RNA excited states through nuclear magnetic resonance measurements of C1΄ and C4΄ rotating frame spin relaxation (R1ρ) in uniformly 13C/15N labeled RNA samples. Application to HIV-1 TAR exposed slow modes of sugar repuckering dynamics at the μs and ms timescale accompanying transitions between non-helical (C2΄-endo) to helical (C3΄-endo) conformations during formation of two distinct excited states. In contrast, we did not obtain any evidence for slow sugar repuckering dynamics for nucleotides in a variety of structural contexts that do not undergo non-helical to helical transitions. Our results outline a route for significantly improving the conformational characterization of RNA excited states and suggest that slow modes of repuckering dynamics gated by transient changes in secondary structure are quite common in RNA.

Robust IR-based detection of stable and fractionally populated G-C+ and A-T Hoogsteen base pairs in duplex DNA

Noncanonical G-C+ and A-T Hoogsteen base pairs can form in duplex DNA and play roles in recognition, damage repair, and replication. Identifying Hoogsteen base pairs in DNA duplexes remains challenging due to difficulties in resolving syn versus antipurine bases with X-ray crystallography; and size limitations and line broadening can make them difficult to characterize by NMR spectroscopy. Here, we show how infrared (IR) spectroscopy can identify G-C+ and A-T Hoogsteen base pairs in duplex DNA across a range of different structural contexts. The utility of IR-based detection of Hoogsteen base pairs is demonstrated by characterizing the first example of adjacent A-T and G-C+ Hoogsteen base pairs in a DNA duplex where severe broadening complicates detection with NMR.

Pseudouridine and N6-methyladenosine modifications weaken PUF protein/RNA interactions

RNA modifications are ubiquitous in biology, with over 100 distinct modifications. While the vast majority were identified and characterized on abundant noncoding RNA such as tRNA and rRNA, the advent of sensitive sequencing-based approaches has led to the discovery of extensive and regulated modification of eukaryotic messenger RNAs as well. The two most abundant mRNA modifications-pseudouridine (Ψ) and N⁶-methyladenosine (m⁶A)-affect diverse cellular processes including mRNA splicing, localization, translation, and decay and modulate RNA structure. Here, we test the hypothesis that RNA modifications directly affect interactions between RNA-binding proteins and target RNA. We show that Ψ and m⁶A weaken the binding of the human single-stranded RNA binding protein Pumilio 2 (hPUM2) to its consensus motif, with individual modifications having effects up to approximately threefold and multiple modifications giving larger effects. While there are likely to be some cases where RNA modifications essentially fully ablate protein binding, here we see modest responses that may be more common. Such modest effects could nevertheless profoundly alter the complex landscape of RNA:protein interactions, and the quantitative rather than qualitative nature of these effects underscores the need for quantitative, systems-level accounting of RNA:protein interactions to understand post-transcriptional regulation.

Amiloride as a new RNA-binding scaffold with activity against HIV-1 TAR

Diversification of RNA-targeted scaffolds offers great promise in the search for selective ligands of therapeutically relevant RNA such as HIV-1 TAR. We herein report the establishment of amiloride as a novel RNA-binding scaffold along with synthetic routes for combinatorial C(5)- and C(6)-diversification. Iterative modifications at the C(5)- and C(6)- positions yielded derivative 24, which demonstrated a 100-fold increase in activity over the parent dimethylamiloride in peptide displacement assays. NMR chemical shift mapping was performed using the 2D SOFAST- [1H-13C] HMQC NMR method, which allowed for facile and rapid evaluation of binding modes for all library members. Cheminformatic analysis revealed distinct differences between selective and non-selective ligands. In this study, we evolved dimethylamiloride from a weak TAR ligand to one of the tightest binding selective TAR ligands reported to date through a novel combination of synthetic methods and analytical techniques. We expect these methods to allow for rapid library expansion and tuning of the amiloride scaffold for a range of RNA targets and for SOFAST NMR to allow unprecedented evaluation of small molecule:RNA interactions.

Direct NMR Evidence that Transient Tautomeric and Anionic States in dG·dT Form Watson-Crick-like Base Pairs

The replicative and translational machinery utilizes the unique geometry of canonical G·C and A·T/U Watson-Crick base pairs to discriminate against DNA and RNA mismatches in order to ensure high fidelity replication, transcription, and translation. There is growing evidence that spontaneous errors occur when mismatches adopt a Watson-Crick-like geometry through tautomerization and/or ionization of the bases. Studies employing NMR relaxation dispersion recently showed that wobble dG·dT and rG·rU mismatches in DNA and RNA duplexes transiently form tautomeric and anionic species with probabilities (≈0.01-0.40%) that are in concordance with replicative and translational errors. Although computational studies indicate that these exceptionally short-lived and low-abundance species form Watson-Crick-like base pairs, their conformation could not be directly deduced from the experimental data, and alternative pairing geometries could not be ruled out. Here, we report direct NMR evidence that the transient tautomeric and anionic species form hydrogen-bonded Watson-Crick-like base pairs. A guanine-to-inosine substitution, which selectively knocks out a Watson-Crick-type (G)N2H2···O2(T) hydrogen bond, significantly destabilized the transient tautomeric and anionic species, as assessed by lack of any detectable chemical exchange by imino nitrogen rotating frame spin relaxation (R1ρ) experiments. An 15N R1ρ NMR experiment targeting the amino nitrogen of guanine (dG-N2) provides direct evidence for Watson-Crick (G)N2H2···O2(T) hydrogen bonding in the transient tautomeric state. The strategy presented in this work can be generally applied to examine hydrogen-bonding patterns in nucleic acid transient states including in other tautomeric and anionic species that are postulated to play roles in replication and translational errors.