Gadolinium-Based NMR Spin Relaxation Measurements of Near-Surface Electrostatic Potentials of Biomolecules.

NMR spectroscopy is an important tool for the measurement of the electrostatic properties of biomolecules. To this point, paramagnetic relaxation enhancements (PREs) of 1H nuclei arising from nitroxide cosolutes in biomolecular solutions have been used to measure effective near-surface electrostatic potentials (ϕENS) of proteins and nucleic acids. Here, we present a gadolinium (Gd)-based NMR method, exploiting Gd chelates with different net charges, for measuring ϕENS values and demonstrate its utility through applications to a number of biomolecular systems. The use of Gd-based cosolutes offers several advantages over nitroxides for ϕENS measurements. First, unlike nitroxide compounds, Gd chelates enable electrostatic potential measurements on oxidation-sensitive proteins that require reducing agents. Second, the large electron spin quantum number of Gd (7/2) results in notably larger PREs for Gd chelates when used at the same concentrations as nitroxide radicals. Thus, it is possible to measure ϕENS values exclusively from + and - charged compounds even for highly charged biomolecules, avoiding the use of neutral cosolutes that, as we further establish here, limits the accuracy of the measured electrostatic potentials. In addition, the smaller concentrations of cosolutes required minimize potential binding to sites on macromolecules. Fourth, the closer proximity of the paramagnetic center and charged groups within Gd chelates, in comparison to the corresponding nitroxide compounds, enables more accurate predictions of ϕENS potentials for cross-validation of the experimental results. The Gd-based method described here, thus, broadens the applicability of studies of biomolecular electrostatics using solution NMR spectroscopy.

Phase Separation Modulates the Thermodynamics and Kinetics of RNA Hybridization.

Biomolecular condensates can influence cellular function in a number of ways, including by changing the structural dynamics and conformational equilibria of the molecules partitioned within them. Here we use methyl transverse relaxation optimized spectroscopy (methyl-TROSY) NMR in conjunction with 2'-O-methyl labeling of RNA to characterize the thermodynamics and kinetics of RNA-RNA base pairing in condensates formed by the C-terminal intrinsically disordered region of CAPRIN1, an RNA-binding protein involved in RNA transport, translation, and stability. CAPRIN1 condensates destabilize RNA-RNA base pairing, resulting from a ∼270-fold decrease and a concomitant ∼15-fold increase in the on- and off-rates for duplex formation, respectively. The ∼30-fold slower diffusion of RNA single strands within the condensed phase partially accounts for the reduced on-rate, but the further ∼9-fold reduction likely reflects shedding of CAPRIN1 chains that are interacting with the RNA prior to hybridization. Our study emphasizes the important role of protein solvation in modulating nucleic acid recognition processes inside condensates.

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 R1ρ 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.

A delayed decoupling methyl-TROSY pulse sequence for atomic resolution studies of folded proteins and RNAs in condensates.

Solution NMR spectroscopy has tremendous potential for providing atomic resolution insights into the interactions between proteins and nucleic acids partitioned into condensed phases of phase-separated systems. However, the highly viscous nature of the condensed phase challenges applications, and in particular, the extraction of quantitative, site-specific information. Here, we present a delayed decoupling-based HMQC pulse sequence for methyl-TROSY studies of 'client' proteins and nucleic acids partitioned into 'scaffold' proteinaceous phase-separated solvents. High sensitivity and excellent quality spectra are recorded of a nascent form of superoxide dismutase and of a small RNA fragment partitioned into CAPRIN1 condensates.

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.

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.

Lifetime of ground conformational state determines the activity of structured RNA.

Biomolecules continually sample alternative conformations. Consequently, even the most energetically favored ground conformational state has a finite lifetime. Here, we show that, in addition to the 3D structure, the lifetime of a ground conformational state determines its biological activity. Using hydrogen-deuterium exchange nuclear magnetic resonance spectroscopy, we found that Zika virus exoribonuclease-resistant RNA (xrRNA) encodes a ground conformational state with a lifetime that is ~105-107 longer than that of canonical base pairs. Mutations that shorten the apparent lifetime of the ground state without affecting its 3D structure decreased exoribonuclease resistance in vitro and impaired virus replication in cells. Additionally, we observed this exceptionally long-lived ground state in xrRNAs from diverse infectious mosquito-borne flaviviruses. These results demonstrate the biological significance of the lifetime of a preorganized ground state and further suggest that elucidating the lifetimes of dominant 3D structures of biomolecules may be crucial for understanding their behaviors and functions.

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.

Sometimes pulses just have to be perfect - An example based on the measurement of amide proton transverse relaxation rates in proteins.

The measurement of spin relaxation rates provides a unique avenue for quantifying dynamic processes in biomolecules. In order to simplify analysis of the measurements so that a few key intuitive parameters can be extracted, it is often the case that experiments are designed to eliminate interference effects between different classes of spin relaxation. One example emerges in the measurement of amide proton (1HN) transverse relaxation rates in 15N labeled proteins, where 15N inversion pulses are applied during a relaxation element to eliminate cross-correlated spin relaxation between 1HN-15N dipole-1HN CSA interactions. We show that unless these pulses are essentially perfect, significant oscillations in magnetization decay profiles can be obtained, due to the excitation of multiple-quantum coherences, leading potentially to errors in measured R2 rates. With the recent development of experiments for quantifying electrostatic potentials via amide proton relaxation rates, the need for highly accurate measurement schemes becomes critical. Straightforward modifications to existing pulse sequences are suggested to achieve this goal.

Practical considerations for the measurement of near-surface electrostatics based on solvent paramagnetic relaxation enhancements.

Electrostatic interactions can play important roles in regulating various biological processes. Quantifying surface electrostatics of biomolecules is, therefore, of significant interest. Recent advances in solution NMR spectroscopy have enabled site-specific measurements of de novo near-surface electrostatic potentials (ϕENS) based on a comparison of solvent paramagnetic relaxation enhancements generated from differently charged paramagnetic co-solutes with similar structures. Although the NMR-derived near-surface electrostatic potentials have been shown to agree with theoretical calculations in the context of folded proteins and nucleic acids, such benchmark comparisons may not always be possible, particularly in cases where high-resolution structural models are lacking, such as in the study of intrinsically disordered proteins. Cross-validation of ϕENS potentials can be achieved by comparing values obtained using three pairs of paramagnetic co-solutes, each with a different net charge. Notably we have found cases where agreement of ϕENS potentials between the three pairs is poor and herein we investigate the source of this discrepancy in some detail. We show that for the systems considered here ϕENS potentials obtained from cationic and anionic co-solutes are accurate and that the use of paramagnetic co-solutes with different structures can be a viable option for validation, although the optimal choice of paramagnetic compounds depends on the system of interest.

Surface electrostatics dictate RNA-binding protein CAPRIN1 condensate concentration and hydrodynamic properties.

Biomolecular condensates concentrate proteins, nucleic acids, and small molecules and play an essential role in many biological processes. Their formation is tuned by a balance between energetically favorable and unfavorable contacts, with charge-charge interactions playing a central role in some systems. The positively charged intrinsically disordered carboxy-terminal region of the RNA-binding protein CAPRIN1 is one such example, phase separating upon addition of negatively charged ATP or high concentrations of sodium chloride (NaCl). Using solution NMR spectroscopy, we measured residue-specific near-surface electrostatic potentials (ϕENS) of CAPRIN1 along its NaCl-induced phase separation trajectory to compare with those obtained using ATP. In both cases, electrostatic shielding decreases ϕENS values, yet surface potentials of CAPRIN1 in the two condensates can be different, depending on the amount of NaCl or ATP added. Our results establish that even small differences in ϕENS can significantly affect the level of protein enrichment and the mechanical properties of the condensed phase, leading, potentially, to the regulation of biological processes.

Measurement of 1Hα transverse relaxation rates in proteins: application to solvent PREs.

It has recently been demonstrated that accurate near surface electrostatic potentials can be calculated for proteins from solvent paramagnetic relaxation enhancements (PREs) of amide protons measured using spin labels of similar structures but different charges (Yu et al. in Proc Natl Acad Sci 118(25):e2104020118, 2021). Here we develop methodology for extending such measurements to intrinsically disordered proteins at neutral pH where amide spectra are of very poor quality. Under these conditions it is shown that accurate PRE values can be measured using the haCONHA experiment that has been modified for recording 1Hα transverse relaxation rates. The optimal pulse scheme includes a spin-lock relaxation element for suppression of homonuclear scalar coupled evolution for all 1Hα protons, except those derived from Ser and Thr residues, and minimizes the radiation damping field from water magnetization that would otherwise increase measured relaxation rates. The robustness of the experiment is verified by developing a second approach using a band selective adiabatic decoupling scheme for suppression of scalar coupling modulations during 1Hα relaxation and showing that the measured PRE values from the two methods are in excellent agreement. The near surface electrostatic potential of a 103-residue construct comprising the C-terminal intrinsically disordered region of the RNA-binding protein CAPRIN1 is obtained at pH 5.5 using both 1HN and 1Hα-based relaxation rates, and at pH 7.4 where only 1Hα rates can be quantified, with very good agreement between potentials obtained under all experimental conditions.

Mapping the per-residue surface electrostatic potential of CAPRIN1 along its phase-separation trajectory.

Electrostatic interactions and charge balance are important for the formation of biomolecular condensates involving proteins and nucleic acids. However, a detailed, atomistic picture of the charge distribution around proteins during the phase-separation process is lacking. Here, we use solution NMR spectroscopy to measure residue-specific near-surface electrostatic potentials (ϕENS) of the positively charged carboxyl-terminal intrinsically disordered 103 residues of CAPRIN1, an RNA-binding protein localized to membraneless organelles playing an important role in messenger RNA (mRNA) storage and translation. Measured ϕENS values have been mapped along the adenosine triphosphate (ATP)-induced phase-separation trajectory. In the absence of ATP, ϕENS values for the mixed state of CAPRIN1 are positive and large and progressively decrease as ATP is added. This is coupled to increasing interchain interactions, particularly between aromatic-rich and arginine-rich regions of the protein. Upon phase separation, CAPRIN1 molecules in the condensed phase are neutral (ϕENS [Formula: see text] 0 mV), with ∼five molecules of ATP associated with each CAPRIN1 chain. Increasing the ATP concentration further inverts the CAPRIN1 electrostatic potential, so that molecules become negatively charged, especially in aromatic-rich regions, leading to re-entrance into a mixed phase. Our results collectively show that a subtle balance between electrostatic repulsion and interchain attractive interactions regulates CAPRIN1 phase separation and provides insight into how nucleotides, such as ATP, can induce formation of and subsequently dissolve protein condensates.

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

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

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.

Probing the Hydrogen-Bonding Environment of Individual Bases in DNA Duplexes with Isotope-Edited Infrared Spectroscopy.

Measuring the strength of the hydrogen bonds between DNA base pairs is of vital importance for understanding how our genetic code is physically accessed and recognized in cells, particularly during replication and transcription. Therefore, it is important to develop probes for these key hydrogen bonds (H-bonds) that dictate events critical to cellular function, such as the localized melting of DNA. The vibrations of carbonyl bonds are well-known probes of their H-bonding environment, and their signals can be observed with infrared (IR) spectroscopy. Yet, pinpointing a single bond of interest in the complex IR spectrum of DNA is challenging due to the large number of carbonyl signals that overlap with each other. Here, we develop a method using isotope editing and infrared (IR) spectroscopy to isolate IR signals from the thymine (T) C2═O carbonyl. We use solvatochromatic studies to show that the TC2═O signal's position in the IR spectrum is sensitive to the H-bonding capacity of the solvent. Our results indicate that C2═O of a single T base within DNA duplexes experiences weak H-bonding interactions. This finding is consistent with the existence of a third, noncanonical CH···O H-bond between adenine and thymine in both Watson-Crick and Hoogsteen base pairs in DNA.

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 13C or 15N isotopically enriched DNA samples for NMR relaxation dispersion (RD) experiments. Here, using SELective Optimized Proton Experiments (SELOPE) 1H CEST experiments employing high-power radiofrequency fields (B1 > 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 1H CEST were in very good agreement with counterparts measured using off-resonance 13C / 15N spin relaxation in the rotating frame (R1ρ). It is shown that 1H-1H NOE effects which typically introduce artifacts in 1H-based measurements of chemical exchange can be effectively suppressed by selective excitation, provided that the relaxation delay is short (≤ 100 ms). The 1H CEST experiment can be performed with ∼ 10× higher throughput and ∼ 100× lower cost relative to 13C / 15N R1ρ and enabled Hoogsteen chemical exchange measurements undetectable by R1ρ. The results reveal an increased propensity to form Hoogsteen bp's near terminal ends and a diminished propensity within A-tract motifs. The 1H CEST experiment provides a basis for rapidly screening Hoogsteen breathing in duplex DNA, enabling identification of unusual motifs for more in-depth characterization.

Rapid and accurate determination of atomistic RNA dynamic ensemble models using NMR and structure prediction.

Biomolecules form dynamic ensembles of many inter-converting conformations which are key for understanding how they fold and function. However, determining ensembles is challenging because the information required to specify atomic structures for thousands of conformations far exceeds that of experimental measurements. We addressed this data gap and dramatically simplified and accelerated RNA ensemble determination by using structure prediction tools that leverage the growing database of RNA structures to generate a conformation library. Refinement of this library with NMR residual dipolar couplings provided an atomistic ensemble model for HIV-1 TAR, and the model accuracy was independently supported by comparisons to quantum-mechanical calculations of NMR chemical shifts, comparison to a crystal structure of a substate, and through designed ensemble redistribution via atomic mutagenesis. Applications to TAR bulge variants and more complex tertiary RNAs support the generality of this approach and the potential to make the determination of atomic-resolution RNA ensembles routine.

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.

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

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