Base flexibility in HIV-2 TAR RNA mapped by solution (15)N, (13)C NMR relaxation

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

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

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

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

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s00706-022-02892-1.

Quantifying the effects of long-range 13C-13C dipolar coupling on measured relaxation rates in RNA

Selective stable isotope labeling has transformed structural and dynamics analysis of RNA by NMR spectroscopy. These methods can remove ¹³C-¹³C dipolar couplings that complicate ¹³C relaxation analyses. While these phenomena are well documented for sites with adjacent ¹³C nuclei (e.g. ribose C1'), less is known about so-called isolated sites (e.g. adenosine C2). To investigate and quantify the effects of long-range (> 2 Å) ¹³C-¹³C dipolar interactions on RNA dynamics, we simulated adenosine C2 relaxation rates in uniformly [U-¹³C/¹⁵N]-ATP or selectively [2-¹³C]-ATP labeled RNAs. Our simulations predict non-negligible ¹³C-¹³C dipolar contributions from adenosine C4, C5, and C6 to C2 longitudinal (R₁) relaxation rates in [U-¹³C/¹⁵N]-ATP labeled RNAs. Moreover, these contributions increase at higher magnetic fields and molecular weights to introduce discrepancies that exceed 50%. This will become increasingly important at GHz fields. Experimental R₁ measurements in the 61 nucleotide human hepatitis B virus encapsidation signal ε RNA labeled with [U-¹³C/¹⁵N]-ATP or [2-¹³C]-ATP corroborate these simulations. Thus, in the absence of selectively labeled samples, long-range ¹³C-¹³C dipolar contributions must be explicitly taken into account when interpreting adenosine C2 R₁ rates in terms of motional models for large RNAs.

Targeting structural features of viral genomes with a nano-sized supramolecular drug

RNA targeting is an exciting frontier for drug design. Intriguing targets include functional RNA structures in structurally-conserved untranslated regions (UTRs) of many lethal viruses. However, computational docking screens, valuable in protein structure targeting, fail for inherently flexible RNA. Herein we harness MD simulations with Markov state modeling to enable nanosize metallo-supramolecular cylinders to explore the dynamic RNA conformational landscape of HIV-1 TAR untranslated region RNA (representative for many viruses) replicating experimental observations. These cylinders are exciting as they have unprecedented nucleic acid binding and are the first supramolecular helicates shown to have anti-viral activity in cellulo: the approach developed in this study provides additional new insight about how such viral UTR structures might be targeted with the cylinder binding into the heart of an RNA-bulge cavity, how that reduces the conformational flexibility of the RNA and molecular details of the insertion mechanism. The approach and understanding developed represents a new roadmap for design of supramolecular drugs to target RNA structural motifs across biology and nucleic acid nanoscience.

In vivo, large-scale preparation of uniformly (15)N- and site-specifically (13)C-labeled homogeneous, recombinant RNA for NMR studies

Knowledge of how ribonucleic acid (RNA) structures fold to form intricate, three-dimensional structures has provided fundamental insights into understanding the biological functions of RNA. Nuclear magnetic resonance (NMR) spectroscopy is a particularly useful high-resolution technique to investigate the dynamic structure of RNA. Effective study of RNA by NMR requires enrichment with isotopes of (13)C or (15)N or both. Here, we present a method to produce milligram quantities of uniformly (15)N- and site-specifically (13)C-labeled RNAs using wild-type K12 and mutant tktA Escherichia coli in combination with a tRNA-scaffold approach. The method includes a double selection protocol to obtain an E. coli clone with consistently high expression of the recombinant tRNA-scaffold. We also present protocols for the purification of the tRNA-scaffold from a total cellular RNA extract and the excision of the RNA of interest from the tRNA-scaffold using DNAzymes. Finally, we showcase NMR applications to demonstrate the benefit of using in vivo site-specifically (13)C-labeled RNA.

Shortening the HIV-1 TAR RNA Bulge by a Single Nucleotide Preserves Motional Modes over a Broad Range of Time Scales

Helix-junction-helix (HJH) motifs are flexible building blocks of RNA architecture that help define the orientation and dynamics of helical domains. They are also frequently involved in adaptive recognition of proteins and small molecules and in the formation of tertiary contacts. Here, we use a battery of nuclear magnetic resonance techniques to examine how deleting a single bulge residue (C24) from the human immunodeficiency virus type 1 (HIV-1) transactivation response element (TAR) trinucleotide bulge (U23-C24-U25) affects dynamics over a broad range of time scales. Shortening the bulge has an effect on picosecond-to-nanosecond interhelical and local bulge dynamics similar to that casued by increasing the Mg(2+) and Na(+) concentration, whereby a preexisting two-state equilibrium in TAR is shifted away from a bent flexible conformation toward a coaxial conformation, in which all three bulge residues are flipped out and flexible. Surprisingly, the point deletion minimally affects microsecond-to-millisecond conformational exchange directed toward two low-populated and short-lived excited conformational states that form through reshuffling of bases pairs throughout TAR. The mutant does, however, adopt a slightly different excited conformational state on the millisecond time scale, in which U23 is intrahelical, mimicking the expected conformation of residue C24 in the excited conformational state of wild-type TAR. Thus, minor changes in HJH topology preserve motional modes in RNA occurring over the picosecond-to-millisecond time scales but alter the relative populations of the sampled states or cause subtle changes in their conformational features.

Elucidating molecular motion through structural and dynamic filters of energy-minimized conformer ensembles

Complex RNA structures are constructed from helical segments connected by flexible loops that move spontaneously and in response to binding of small molecule ligands and proteins. Understanding the conformational variability of RNA requires the characterization of the coupled time evolution of interconnected flexible domains. To elucidate the collective molecular motions and explore the conformational landscape of the HIV-1 TAR RNA, we describe a new methodology that utilizes energy-minimized structures generated by the program “Fragment Assembly of RNA with Full-Atom Refinement (FARFAR)”. We apply structural filters in the form of experimental residual dipolar couplings (RDCs) to select a subset of discrete energy-minimized conformers and carry out principal component analyses (PCA) to corroborate the choice of the filtered subset. We use this subset of structures to calculate solution T₁ and T_1ρ relaxation times for ¹³C spins in multiple residues in different domains of the molecule using two simulation protocols that we previously published. We match the experimental T₁ times to within 2% and the T_1ρ times to within less than 10% for helical residues. These results introduce a protocol to construct viable dynamic trajectories for RNA molecules that accord well with experimental NMR data and support the notion that the motions of the helical portions of this small RNA can be described by a relatively small number of discrete conformations exchanging over time scales longer than 1 μs.

Deriving quantitative dynamics information for proteins and RNAs using ROTDIF with a graphical user interface

To facilitate rigorous analysis of molecular motions in proteins, DNA, and RNA, we present a new version of ROTDIF, a program for determining the overall rotational diffusion tensor from single- or multiple-field nuclear magnetic resonance relaxation data. We introduce four major features that expand the program's versatility and usability. The first feature is the ability to analyze, separately or together, (13)C and/or (15)N relaxation data collected at a single or multiple fields. A significant improvement in the accuracy compared to direct analysis of R2/R1 ratios, especially critical for analysis of (13)C relaxation data, is achieved by subtracting high-frequency contributions to relaxation rates. The second new feature is an improved method for computing the rotational diffusion tensor in the presence of biased errors, such as large conformational exchange contributions, that significantly enhances the accuracy of the computation. The third new feature is the integration of the domain alignment and docking module for relaxation-based structure determination of multi-domain systems. Finally, to improve accessibility to all the program features, we introduced a graphical user interface that simplifies and speeds up the analysis of the data. Written in Java, the new ROTDIF can run on virtually any computer platform. In addition, the new ROTDIF achieves an order of magnitude speedup over the previous version by implementing a more efficient deterministic minimization algorithm. We not only demonstrate the improvement in accuracy and speed of the new algorithm for synthetic and experimental (13)C and (15)N relaxation data for several proteins and nucleic acids, but also show that careful analysis required especially for characterizing RNA dynamics allowed us to uncover subtle conformational changes in RNA as a function of temperature that were opaque to previous analysis.

Independent alignment of RNA for dynamic studies using residual dipolar couplings

Molecular motion and dynamics play an essential role in the biological function of many RNAs. An important source of information on biomolecular motion can be found in residual dipolar couplings which contain dynamics information over the entire ms-ps timescale. However, these methods are not fully applicable to RNA because nucleic acid molecules tend to align in a highly collinear manner in different alignment media. As a consequence, information on dynamics that can be obtained with this method is limited. In order to overcome this limitation, we have generated a chimeric RNA containing both the wild type TAR RNA, the target of our investigation of dynamics, as well as the binding site for U1A protein. When U1A protein was bound to the portion of the chimeric RNA containing its binding site, we obtained independent alignment of TAR by exploiting the physical chemical characteristics of this protein. This technique can allow the extraction of new information on RNA dynamics, which is particularly important for time scales not covered by relaxation methods where important RNA motions occur.

Examining the relationship between RNA function and motion using nuclear magnetic resonance

The biological function of proteins and nucleic acids relies on their complex structures, yet dynamics provides an additional layer of functional adaptability. Numerous studies have demonstrated that RNA is only able to perform the multitude of functions for which it is responsible by readily changing its conformation in response to binding of proteins or small molecules. Examination of RNA dynamics is therefore essential to understanding its biological function. Nuclear magnetic resonance (NMR) has emerged as a leading technique for the examination of RNA motion and conformational transitions. It can examine domain motions as well as motion with atomic level resolution over a wide range of time scales. This review examines how NMR spectroscopy can be applied to examine the relationship between function and dynamics in RNA.

Selective 13C labeling of nucleotides for large RNA NMR spectroscopy using an E. coli strain disabled in the TCA cycle

Escherichia coli (E. coli) is an ideal organism to tailor-make labeled nucleotides for biophysical studies of RNA. Recently, we showed that adding labeled formate enhanced the isotopic enrichment at protonated carbon sites in nucleotides. In this paper, we show that growth of a mutant E. coli strain DL323 (lacking succinate and malate dehydrogenases) on (13)C-2-glycerol and (13)C-1,3-glycerol enables selective labeling at many useful sites for RNA NMR spectroscopy. For DL323 E. coli grown in (13)C-2-glycerol without labeled formate, all the ribose carbon atoms are labeled except the C3' and C5' carbon positions. Consequently the C1', C2' and C4' positions remain singlet. In addition, only the pyrimidine base C6 atoms are substantially labeled to ~96% whereas the C2 and C8 atoms of purine are labeled to ~5%. Supplementing the growth media with (13)C-formate increases the labeling at C8 to ~88%, but not C2. Not unexpectedly, addition of exogenous formate is unnecessary for attaining the high enrichment levels of ~88% for the C2 and C8 purine positions in a (13)C-1,3-glycerol based growth. Furthermore, the ribose ring is labeled in all but the C4' carbon position, such that the C2' and C3' positions suffer from multiplet splitting but the C5' position remains singlet and the C1' position shows a small amount of residual C1'-C2' coupling. As expected, all the protonated base atoms, except C6, are labeled to ~90%. In addition, labeling with (13)C-1,3-glycerol affords an isolated methylene ribose with high enrichment at the C5' position (~90%) that makes it particularly attractive for NMR applications involving CH(2)-TROSY modules without the need for decoupling the C4' carbon. To simulate the tumbling of large RNA molecules, perdeuterated glycerol was added to a mixture of the four nucleotides, and the methylene TROSY experiment recorded at various temperatures. Even under conditions of slow tumbling, all the expected carbon correlations were observed, which indicates this approach of using nucleotides obtained from DL323 E. coli will be applicable to high molecular weight RNA systems.

Probing adenine rings and backbone linkages using base specific isotope-edited Raman spectroscopy: application to group II intron ribozyme domain V

Raman difference spectroscopy is used to probe the properties of a 36-nt RNA molecule, "D5", which lies at the heart of the catalytic apparatus in group II introns. For D5 that has all of its adenine residues labeled with (13)C and (15)N and utilizing Raman difference spectroscopy, we identify the conformationally sensitive -C-O-P-O-C- stretching modes of the unlabeled bonds adjacent to adenine bases, as well as the adenine ring modes themselves. The phosphodiester modes can be assigned to individual adenine residues based on earlier NMR data. The effect of Mg(2+) binding was explored by analyzing the Raman difference spectra for [D5 + Mg(2+)] minus [D5 no Mg(2+)], for D5 unlabeled, or D5 labeled with (13)C/(15)N-enriched adenine. In both sets of data we assign differential features to G ring modes perturbed by Mg(2+) binding at the N7 position. In the A-labeled spectra we attribute a Raman differential near 1450 cm(-1) and changes of intensity at 1296 cm(-1) to Mg binding at the N7 position of adenine bases. The A and G bases involved in Mg(2+) binding again can be identified using earlier NMR results. For the unlabeled D5, a change in the C-O-P-O-C stretch profile at 811 cm(-1) upon magnesium binding is due to a "tightening up" (in the sense of a more rigid molecule with less dynamic interchange among competing ribose conformers) of the D5 structure. For adenine-labeled D5, small changes in the adenine backbone bond signatures in the 810-830 cm(-1) region suggest that small conformational changes occur in the tetraloop and bulge regions upon binding of Mg(2+). The PO(2)(-) stretching vibration, near 1100 cm(-1), from the nonbridging phosphate groups, probes the effect of Mg(2+)-hydrate inner-sphere interactions that cause an upshift. In turn, the upshift is modulated by the presence of monovalent cations since in the presence of Na(+) and Li(+) the upshift is 23 +/- 2 cm(-1) while in the presence of K(+) and Cs(+) it is 13 +/- 3 cm(-1), a finding that correlates with the differences in hydration radii. These subtle differences in electrostatic interactions may be related to observed variations in catalytic activity. For a reconstructed ribozyme comprising domains 1-3 (D123) connected in cis plus domain 5 (D5) supplied in trans, cleavage of spliced exon substrates in the presence of magnesium and K(+) or Cs(+) is more efficient than that in the presence of magnesium with Na(+) or Li(+).

Structures of HIV TAR RNA-ligand complexes reveal higher binding stoichiometries

Target TAR by NMR: Tripeptides containing arginines as terminal residues and non-natural amino acids as central residues are good leads for drug design to target the HIV trans-activation response element (TAR). The structural characterization of the RNA-ligand complex by NMR spectroscopy reveals two specific binding sites that are located at bulge residue U23 and around the pyrimidine-stretch U40-C41-U42 directly adjacent to the bulge.

Hydration dependent dynamics in RNA

The essential role played by local and collective motions in RNA function has led to a growing interest in the characterization of RNA dynamics. Recent investigations have revealed that even relatively simple RNAs experience complex motions over multiple time scales covering the entire ms-ps motional range. In this work, we use deuterium solid-state NMR to systematically investigate motions in HIV-1 TAR RNA as a function of hydration. We probe dynamics at three uridine residues in different structural environments ranging from helical to completely unrestrained. We observe distinct and substantial changes in (2)H solid-state relaxation times and lineshapes at each site as hydration levels increase. By comparing solid-state and solution state (13)C relaxation measurements, we establish that ns-micros motions that may be indicative of collective dynamics suddenly arise in the RNA as hydration reaches a critical point coincident with the onset of bulk hydration. Beyond that point, we observe smaller changes in relaxation rates and lineshapes in these highly hydrated solid samples, compared to the dramatic activation of motion occurring at moderate hydration.

Constructing RNA dynamical ensembles by combining MD and motionally decoupled NMR RDCs: new insights into RNA dynamics and adaptive ligand recognition

We describe a strategy for constructing atomic resolution dynamical ensembles of RNA molecules, spanning up to millisecond timescales, that combines molecular dynamics (MD) simulations with NMR residual dipolar couplings (RDC) measured in elongated RNA. The ensembles are generated via a Monte Carlo procedure by selecting snap-shot from an MD trajectory that reproduce experimentally measured RDCs. Using this approach, we construct ensembles for two variants of the transactivation response element (TAR) containing three (HIV-1) and two (HIV-2) nucleotide bulges. The HIV-1 TAR ensemble reveals significant mobility in bulge residues C24 and U25 and to a lesser extent U23 and neighboring helical residue A22 that give rise to large amplitude spatially correlated twisting and bending helical motions. Omission of bulge residue C24 in HIV-2 TAR leads to a significant reduction in both the local mobility in and around the bulge and amplitude of inter-helical bending motions. In contrast, twisting motions of the helices remain comparable in amplitude to HIV-1 TAR and spatial correlations between them increase significantly. Comparison of the HIV-1 TAR dynamical ensemble and ligand bound TAR conformations reveals that several features of the binding pocket and global conformation are dynamically preformed, providing support for adaptive recognition via a ‘conformational selection’ type mechanism.

How binding of small molecule and peptide ligands to HIV-1 TAR alters the RNA motional landscape

The HIV-1 TAR RNA represents a well-known paradigm to study the role of dynamics and conformational change in RNA function. This regulatory RNA changes conformation in response to binding of Tat protein and of a variety of peptidic and small molecule ligands, indicating that its conformational flexibility and intrinsic dynamics play important roles in molecular recognition. We have used ¹³C NMR relaxation experiments to examine changes in the motional landscape of HIV-1 TAR in the presence of three ligands of different affinity and specificity. The ligands are argininamide, a linear peptide mimic of the Tat basic domain and a cyclic peptide that potently inhibits Tat-dependent activation of transcription. All three molecules induce the same motional characteristics within the three nucleotides bulge that represents the Tat-binding site. However, the cyclic peptide has a unique motional signature in the apical loop, which represents a binding site for the essential host co-factor cyclin T1. These results suggest that all peptidic mimics of Tat induce the same dynamics in TAR within this protein binding site. However, the new cyclic peptide mimic of Tat represents a new class of ligands with a unique effect on the dynamics and the structure of the apical loop.

Extensive backbone dynamics in the GCAA RNA tetraloop analyzed using 13C NMR spin relaxation and specific isotope labeling

Conformational dynamics play a key role in the properties and functions of proteins and nucleic acids. Heteronuclear NMR spin relaxation is a uniquely powerful site-specific probe of dynamics in proteins and has found increasing applications to nucleotide base side chains and anomeric sites in RNA. Applications to the nucleic acid ribose backbone, however, have been hampered by strong magnetic coupling among ring carbons in uniformly 13C-labeled samples. In this work, we apply a recently developed, metabolically directed isotope labeling scheme that places 13C with high efficiency and specificity at the nucleotide ribose C2' and C4' sites. We take advantage of this scheme to explore backbone dynamics in the well-studied GCAA RNA tetraloop. Using a combination of CPMG (Carr-Purcell-Meiboom-Gill) and R(1rho) relaxation dispersion spectroscopy to explore exchange processes on the microsecond to millisecond time scale, we find an extensive pattern of dynamic transitions connecting a set of relatively well-defined conformations. In many cases, the observed transitions appear to be linked to C3'-endo/C2'-endo sugar pucker transitions of the corresponding nucleotides, and may also be correlated across multiple nucleotides within the tetraloop. These results demonstrate the power of NMR spin relaxation based on alternate-site isotope labeling to open a new window into the dynamic properties of ribose backbone groups in RNA.

RNA in motion

Although RNA duplex regions are highly structured and inflexible, other elements of an RNA molecule are capable of dynamic motions. These flexible regions are the sites of interactions with small molecules, proteins, and other RNAs, yet there are few descriptions of these regions that include the timescale and amplitude of their motions. No one technique is sufficient to accurately describe these motions, but the combination of in vitro methods, particularly NMR relaxation methods, and more robust in silico methods, is beginning to yield the type of data that can be used to understand RNA function. Very few RNAs have been described by both techniques, and here one such RNA and one RNA:protein complex are reviewed.

Strong correlation between SHAPE chemistry and the generalized NMR order parameter (S2) in RNA

The functions of most RNA molecules are critically dependent on the distinct local dynamics that characterize secondary structure and tertiary interactions and on structural changes that occur upon binding by proteins and small molecule ligands. Measurements of RNA dynamics at nucleotide resolution set the foundation for understanding the roles of individual residues in folding, catalysis, and ligand recognition. In favorable cases, local order in small RNAs can be quantitatively analyzed by NMR in terms of a generalized order parameter, S2. Alternatively, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) chemistry measures local nucleotide flexibility in RNAs of any size using structure-sensitive reagents that acylate the 2'-hydroxyl position. In this work, we compare per-residue RNA dynamics, analyzed by both S2 and SHAPE, for three RNAs: the HIV-1 TAR element, the U1A protein binding site, and the Tetrahymena telomerase stem loop 4. We find a very strong correlation between the two measurements: nucleotides with high SHAPE reactivities consistently have low S2 values. We conclude that SHAPE chemistry quantitatively reports local nucleotide dynamics and can be used with confidence to analyze dynamics in large RNAs, RNA-protein complexes, and RNAs in vivo.

Characterizing complex dynamics in the transactivation response element apical loop and motional correlations with the bulge by NMR, molecular dynamics, and mutagenesis

The HIV-1 transactivation response element (TAR) RNA binds a variety of proteins and is a target for developing anti-HIV therapies. TAR has two primary binding sites: a UCU bulge and a CUGGGA apical loop. We used NMR residual dipolar couplings, carbon spin relaxation (R(1) and R(2)), and relaxation dispersion (R(1rho)) in conjunction with molecular dynamics and mutagenesis to characterize the dynamics of the TAR apical loop and investigate previously proposed long-range interactions with the distant bulge. Replacement of the wild-type apical loop with a UUCG loop did not significantly affect the structural dynamics at the bulge, indicating that the apical loop and the bulge act largely as independent dynamical recognition centers. The apical loop undergoes complex dynamics at multiple timescales that are likely important for adaptive recognition: U31 and G33 undergo limited motions, G32 is highly flexible at picosecond-nanosecond timescales, and G34 and C30 form a dynamic Watson-Crick basepair in which G34 and A35 undergo a slow (approximately 30 mus) likely concerted looping in and out motion, with A35 also undergoing large amplitude motions at picosecond-nanosecond timescales. Our study highlights the power of combining NMR, molecular dynamics, and mutagenesis in characterizing RNA dynamics.

Solid-state deuterium NMR studies reveal micros-ns motions in the HIV-1 transactivation response RNA recognition site

Solution and solid-state NMR measurements were used together to examine motion in three sites in the HIV-1 TAR RNA. We wished to investigate the dynamics facilitating the conformational rearrangements the TAR RNA must undergo for tat binding, and in particular to characterize the full range of motional timescales accessible to this RNA. Our results demonstrate that the dynamics in TAR involving residues essential to tat binding include not only the faster motions detected by solution relaxation measurements, but also a significant component in the μs-ns timescale.

Tripeptides from synthetic amino acids block the Tat-TAR association and slow down HIV spread in cell cultures

Non-natural amino acids with aromatic or heteroaromatic side chains were incorporated into tripeptides of the general structure Arg-X-Arg and tested as ligands of the HIV RNA element TAR. Some of these compounds could compete efficiently with the association of TAR and Tat and downregulated a TAR-controlled reporter gene in HeLa cells. Peptide 7, which contains a 2-pyrimidinyl-alkyl chain, also inhibited the spread of HIV-1 in cell cultures. NMR studies of 7 bound to HIV-2-TAR gave evidence for contacts in the bulge region.

Conserved nucleotides in an RNA essential for hepatitis B virus replication show distinct mobility patterns

The number of regulatory RNAs with identified non-canonical structures is increasing, and structural transitions often play a role in their biological function. This stimulates interest in internal motions of RNA, which can underlie structural transitions. Heteronuclear NMR relaxation measurements, which are commonly used to study internal motion, only report on local motions of few sites within the molecule. Here we have studied a 27-nt segment of the human hepatitis B virus (HBV) pregenomic RNA, which is essential for viral replication. We combined heteronuclear relaxation with the new off-resonance ROESY technique, which reports on internal motions of H,H contacts. Using off-resonance ROESY, we could for the first time detect motion of through-space H,H contacts, such as in intra-residue base-ribose contacts or inter-nucleotide contacts, both essential for NMR structure determination. Motions in non-canonical structure elements were found primarily on the sub-nanosecond timescale. Different patterns of mobility were observed among several mobile nucleotides. The most mobile nucleotides are highly conserved among different HBV strains, suggesting that their mobility patterns may be necessary for the RNA’s biological function.

SHARP-TACSY: triple-band tailored correlated spectroscopy for base-to-sugar transfer in nucleic acid residues with intermediate time scale motions

Established experiments to identify the sugar-to-base connectivity in isotopically labeled RNA require long transfer periods and are inefficient for residues undergoing intermediate time scale motions (microsecond to millisecond). Here, an alternative transfer experiment is introduced, whereby the C1'-N1/9-C6/8 spin system is selectively brought to the so-called Hartmann-Hahn condition using selectiveheteronuclear planar triple-band tailored correlated spectroscopy (SHARP-TACSY). Results are shown for the fully labeled 30-mer oligonucleotide TAR RNA with particular attention placed on residues from and close to the bulge and the loop. For these residues, the faster relaxation can be attributed to exchange contributions stemming from transient stacking and unstacking of the bases and/or from the isomerization of the ribose sugar pucker. The new experiment shows improved signal-to-noise for residues exhibiting large microsecond-millisecond time scale motions with respect to established experiments, thus providing a valid alternative for resonance assignment in mobile RNA regions.

NMR studies of dynamics in RNA and DNA by 13C relaxation

RNA and DNA molecules experience motions on a wide range of time scales, ranging from rapid localized motions to much slower collective motions of entire helical domains. The many functions of RNA in biology very often require this molecule to change its conformation in response to biological signals in the form of small molecules, proteins or other nucleic acids, whereas local motions in DNA may facilitate protein recognition and allow enzymes acting on DNA to access functional groups on the bases that would otherwise be buried in Watson-Crick base pairs. Although these statements make a compelling case to study the sequence dependent dynamics in nucleic acids, there are few residue-specific studies of nucleic acid dynamics. Fortunately, NMR studies of dynamics of nucleic acids and nucleic acids-protein complexes are gaining increased attention. The aim of this review is to provide an update of the recent progress in studies of nucleic acid dynamics by NMR based on the application of solution relaxation techniques.

Single-Molecule Study of the Inhibition of HIV-1 Transactivation Response Region DNA/DNA Annealing by Argininamide

Single-molecule spectroscopy was used to examine how a model inhibitor of HIV-1, argininamide, modulates the nucleic acid chaperone activity of the nucleocapsid protein (NC) in the minus-strand transfer step of HIV-1 reverse transcription, in vitro. In minus-strand transfer, the transactivation response region (TAR) RNA of the genome is annealed to the complementary "TAR DNA" generated during minus-strand strong-stop DNA synthesis. Argininamide and its analogs are known to bind to the hairpin bulge region of TAR RNA as well as to various DNA loop structures, but its ability to inhibit the strand transfer process has only been implied. Here, we explore how argininamide modulates the annealing kinetics and secondary structure of TAR DNA. The studies reveal that the argininamide inhibitory mechanism involves a shift of the secondary structure of TAR, away from the NC-induced "Y" form, an intermediate in reverse transcription, and toward the free closed or "C" form. In addition, more potent inhibition of the loop-mediated annealing pathway than stem-mediated annealing is observed. Taken together, these data suggest a molecular mechanism wherein argininamide inhibits NC-facilitated TAR RNA/DNA annealing in vitro by interfering with the formation of key annealing intermediates.

Probing Na(+)-induced changes in the HIV-1 TAR conformational dynamics using NMR residual dipolar couplings: new insights into the role of counterions and electrostatic interactions in adaptive recognition

Many regulatory RNAs undergo large changes in structure upon recognition of proteins and ligands, but the mechanism by which this occurs remains poorly understood. Using NMR residual dipolar coupling (RDCs), we characterized Na+-induced changes in the structure and dynamics of the bulge-containing HIV-1 transactivation response element (TAR) RNA that mirrors changes induced by small molecules bearing a different number of cationic groups. Increasing the Na+ concentration from 25 to 320 mM led to a continuous reduction in the average inter-helical bend angle (from 46 degrees to 22 degrees ), inter-helical twist angle (from 66 degrees to -18 degrees ), and inter-helix flexibility (as measured by an increase in the internal generalized degree of order from 0.56 to 0.74). Similar conformational changes were observed with Mg2+, indicating that nonspecific electrostatic interactions drive the conformational transition, although results also suggest that Na+ and Mg2+ may associate with TAR in distinct modes. The transition can be rationalized on the basis of a population-weighted average of two ensembles comprising an electrostatically relaxed bent and flexible TAR conformation that is weakly associated with counterions and a globally rigid coaxial conformation that has stronger electrostatic potential and association with counterions. The TAR inter-helical orientations that are stabilized by small molecules fall around the metal-induced conformational pathway, indicating that counterions may help predispose the TAR conformation for target recognition. Our results underscore the intricate sensitivity of RNA conformational dynamics to environmental conditions and demonstrate the ability to detect subtle conformational changes using NMR RDCs.

Strong coupling effects during X-pulse CPMG experiments recorded on heteronuclear ABX spin systems: artifacts and a simple solution

Simulation and experiment have been used to establish that significant artifacts can be generated in X-pulse CPMG relaxation dispersion experiments recorded on heteronuclear ABX spin-systems, such as (13)C( i )-(13)C( j )-(1)H, where (13)C( i ) and (13)C( j ) are strongly coupled. A qualitative explanation of the origin of these artifacts is presented along with a simple method to significantly reduce them. An application to the measurement of (1)H CPMG relaxation dispersion profiles in an HIV-2 TAR RNA molecule where all ribose sugars are protonated at the 2' position, deuterated at all other sugar positions and (13)C labeled at all sugar carbons is presented to illustrate the problems that strong (13)C-(13)C coupling introduces and a simple solution is proposed.

Binding of U1A protein changes RNA dynamics as observed by 13C NMR relaxation studies

Recognition of RNA by proteins and small molecules often involves large changes in RNA structure and dynamics, yet very few studies have so far characterized these motional changes. Here we extend to the protein-bound RNA recent 13C relaxation studies of motions in the RNA recognized by human U1A protein, a well-known model for protein-RNA recognition. Changes in relaxation observed upon complex formation demonstrate that the protein-binding site becomes rigid in the complex, but the upper stem-loop that defines the secondary structure of this RNA experiences unexpected motional freedom. By using a helix elongation strategy, we observe that the upper stem-loop moves independently of the remainder of the structure also in the absence of U1A. Surprisingly, RNA residues making important intermolecular contacts in the structure of the complex exhibit increased flexibility in the presence of the protein. Both of these results support the hypothesis that RNA-binding proteins select a structure that optimizes intermolecular contacts in the manifold of conformations sampled by the free RNA and that protein binding quenches these motions. Together with previous studies of the RNA-bound protein, they also demonstrate that protein-RNA interfaces experience complex motions that modulate the strength of individual interactions.

Internal Bulge and Tetraloop of the Catalytic Domain 5 of a Group II Intron Ribozyme Are Flexible: Implications for Catalysis

RNA molecules have an inherent flexibility that enables recognition of other interacting partners through potential disorder-order transitions, yet studies to quantify such motional dynamics remain few. With an increasing database of three-dimensional structures of biologically important RNA molecules, quantifying such motions becomes important to link structural deformations with function. One such system studied intensely is domain 5 (D5) from the self-splicing group II introns, which is at the heart of its catalytic machinery. We report the dynamics of a 36 nucleotide D5 from the Pylaiella littoralis group II intron in the presence and absence of magnesium ions, and at a range of temperatures (298K-318 K). Using high-resolution NMR experiments of heteronuclear nuclear Overhauser enhancement (NOE), spin-lattice (R(1)), and spin-spin (R(2)) (13)C relaxation rates, we determined the rotational diffusion tensor of D5 using the ROTDIF program modified for RNA dynamic analysis (ROTDIF_RNA). The D5 rotational diffusion tensor has an axial symmetric ratio (D(||)/D(perpendicular)) of 1.7+/-0.3, consistent with an estimated overall rotational correlation time of tau(m)=(2D(||)+4D(perpendicular))(-1) of 6.1(+/-0.3) ns at 298 K and 4.1(+/-0.2) ns at 318 K. The measured relaxation data were analyzed with the reduced spectral density mapping formalism using assumed values of the chemical shift anisotropy of the (13)C spins. Both the relaxation data and the values of the spectral density function reveal that the functional groups in D5 implicated in magnesium ion binding and catalysis (catalytic triad, internal bulge, and tetraloop regions) exhibit thermally induced motion on a wide variety of timescales. Because these motions parallel those observed in the intramolecular stem-loop of the U6 element within the spliceosome, we hypothesize that such extensive dynamic disorder likely facilitates D5 engaging both binding and catalytic regions of the ribozyme, and these may be a conserved feature of the catalytic machinery essential for catalysis.

Alternate-site isotopic labeling of ribonucleotides for NMR studies of ribose conformational dynamics in RNA

Heteronuclear NMR spin relaxation studies of conformational dynamics are coming into increasing use to help understand the functions of ribozymes and other RNAs. Due to strong 13C-13C magnetic interactions within the ribose ring, however, these studies have thus far largely been limited to (13)C and (15)N resonances on the nucleotide base side chains. We report here the application of the alternate-site (13)C isotopic labeling scheme, pioneered by LeMaster for relaxation studies of amino acid side chains, to nucleic acid systems. We have used different strains of E. coli to prepare mononucleotides containing (13)C label in one of two patterns: Either C1' or C2' in addition to C4', termed (1'/2',4') labeling, or nearly complete labeling at the C2' and C4' sites only, termed (2',4') labeling. These patterns provide isolated 13C-1H spin systems on the labeled carbon atoms and thus allow spin relaxation studies without interference from 13C-13C scalar or dipolar coupling. Using relaxation studies of AMP dissolved in glycerol at varying temperature to produce systems with correlation times characteristic of different size RNAs, we demonstrate the removal of errors due to 13C-13C interaction in T (1) measurements of larger nucleic acids and in T (1rho) measurements in RNA molecules. By extending the applicability of spin relaxation measurements to backbone ribose groups, this technology should greatly improve the flexibility and completeness of NMR analyses of conformational dynamics in RNA.

Structural study of elements of Tetrahymena telomerase RNA stem-loop IV domain important for function

Tetrahymena telomerase RNA (TER) contains several regions in addition to the template that are important for function. Central among these is the stem-loop IV domain, which is involved in both catalysis and RNP assembly, and includes binding sites for both the holoenzyme assembly protein p65 and telomerase reverse transcriptase (TERT). Stem-loop IV contains two regions with high evolutionary sequence conservation: a central GA bulge between helices, and a terminal loop. We solved the solution structure of loop IV and modeled the structure of the helical region containing the GA bulge, using NMR and residual dipolar couplings. The central GA bulge with flanking C-G base pairs induces a approximately 50 degrees semi-rigid bend in the helix. Loop IV is highly structured, and contains a conserved C-U base pair at the top of the helical stem. Analysis of new and previous biochemical data in light of the structure provides a rationale for some of the sequence conservation in this region of TER. The results suggest that during holoenzyme assembly the protein p65 recognizes a bend in stem IV, and this binding to central stem IV helps to position the structured loop IV for interaction with TERT and other region(s) of TER.

Decoding RNA motional codes

When proteins and small molecules bind to RNA, they often alter its conformation. These structural changes are an essential aspect of the ability of RNA to sense signaling molecules and modulate gene expression. Thus far, few studies have been dedicated to understanding how RNA moves at a residue level and how these motions change upon complex formation. A recent report highlights how intrinsic motions in RNA correlate with its ability to bind to cognate ligands.

New RNA Labeling Methods Offer Dramatic Sensitivity Enhancements in 2H NMR Relaxation Spectra

A new labeling strategy is presented that greatly facilitates the measurement of 2H spin relaxation rates in RNA molecules as a probe of pico- to nanosecond time scale dynamics. In this labeling scheme the sugar positions are uniformly 13C-labeled, with position 2' protonated and all other sites on the sugar deuterated. Pulse sequences are presented for measurement of 2H R1 and R2 relaxation rates at positions 1', 3', and 4' with sensitivity gains that are on the order of 5-fold relative to previous methods that employed random fractional deuteration. The improved sensitivity is transformative and facilitates the study of motion in moderately sized RNA molecules with good sensitivity. The utility of the approach is demonstrated with an application to HIV-2 TAR, where the site-specific measures of molecular dynamics at sugar positions obtained here complement previous studies of dynamics at aromatic sites in the molecule.

Measurement of long-range 1H-19F scalar coupling constants and their glycosidic torsion dependence in 5-fluoropyrimidine-substituted RNA

Long-range scalar 5J(H1',F) couplings were observed in 5-fluoropyrimidine-substituted RNA. We developed a novel S3E-19F-alpha,beta-edited NOESY experiment for quantitation of these long-range scalar 5J(H1',F) couplings, where the J-couplings can be extracted from inspection of intraresidual (H1',H6) NOE cross-peaks. Quantum chemical calculations were exploited to investigate the relation between scalar couplings and conformations around the glycosidic bond in oligonucleotides. The theoretical dependence of the observed 5J(H1',F) couplings on the torsion angle chi can be described by a generalized Karplus relationship. The corresponding density functional theory (DFT) analysis is outlined. Additional NMR experiments facilitating the resonance assignments of 5-fluoropyrimidine-substituted RNAs are described, and chemical shift changes due to altered shielding in the presence of fluorine-19 (19F) are presented.

The structure of a mixed LNA/DNA:RNA duplex is driven by conformational coupling between LNA and deoxyribose residues as determined from 13C relaxation measurements

A study of the internal dynamics of an LNA/DNA:RNA duplex has been performed to further characterize the conformational changes associated with the incorporation of locked nucleic acid (LNA) nucleotides in a DNA:RNA duplex. In general, it was demonstrated that the LNA/DNA:RNA duplex has a very high degree of order compared to dsDNA and dsRNA duplexes. The order parameters of the aromatic carbon atoms in the LNA/DNA strand are uniformly high, whereas a sharp drop in the degree of order was seen in the RNA strand in the beginning of the AUAU stretch in the middle of the strand. This can be related to a return to normal dsRNA dynamics for the central A:U base pair. The high order of the heteroduplex is consistent with preorganization of the chimera strand for an A-form duplex conformation. These results partly explain the dramatic increase in T(m) of the chimeric heteroduplex over dsDNA and DNA:RNA hybrids of the same sequence.

Dynamics-based amplification of RNA function and its characterization by using NMR spectroscopy

The ever-increasing cellular roles ascribed to RNA raise fundamental questions regarding how a biopolymer composed of only four chemically similar building-block nucleotides achieves such functional diversity. Here, I discuss how RNA achieves added mechanistic and chemical complexity by undergoing highly controlled conformational changes in response to a variety of cellular signals. I examine pathways for achieving selectivity in these conformational changes that rely to different extents on the structure and dynamics of RNA. Finally, I review solution-state NMR techniques that can be used to characterize RNA structural dynamics and its relationship to function.

NMR methods for studying the structure and dynamics of RNA

Proper functioning of RNAs requires the formation of complex three-dimensional structures combined with the ability to rapidly interconvert between multiple functional states. This review covers recent advances in isotope-labeling strategies and NMR experimental approaches that have promise for facilitating solution structure determinations and dynamics studies of biologically active RNAs. Improved methods for the production of isotopically labeled RNAs combined with new multidimensional heteronuclear NMR experiments make it possible to dramatically reduce spectral crowding and simplify resonance assignments for RNAs. Several novel applications of experiments that directly detect hydrogen-bonding interactions are discussed. These studies demonstrate how NMR spectroscopy can be used to distinguish between possible secondary structures and identify mechanisms of ligand binding in RNAs. A variety of recently developed methods for measuring base and sugar residual dipolar couplings are described. NMR residual dipolar coupling techniques provide valuable data for determining the long-range structure and orientation of helical regions in RNAs. A number of studies are also presented where residual dipolar coupling constraints are used to determine the global structure and dynamics of RNAs. NMR relaxation data can be used to probe the dynamics of macromolecules in solution. The power dependence of transverse rotating-frame relaxation rates was used here to study dynamics in the minimal hammerhead ribozyme. Improved methods for isotopically labeling RNAs combined with new types of structural data obtained from a growing repertoire of NMR experiments are facilitating structural and dynamic studies of larger RNAs.

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.

Dynamics and metal ion binding in the U6 RNA intramolecular stem-loop as analyzed by NMR

The U6 RNA intramolecular stem-loop (ISL) is a conserved component of the spliceosome, and contains an essential metal ion binding site centered between a protonated adenine, A79, and U80. Correlated with protonation of A79, U80 undergoes a base-flipping conformational change accompanied by significant helical movement. We have investigated the dynamics of the U6 ISL by analyzing the power dependence of 13C NMR relaxation rates in the rotating frame. The data provide evidence that the conformational transition is centered around an exchange lifetime of 84 micros. The U80 nucleotide displays low internal mobility on the picosecond time-scale at pH 7.0 but high internal mobility at pH 6.0, in agreement with the global transition resulting in the base of U80 adopting a looped-out conformation with increased dynamic disorder. A kinetic analysis suggests that the conformational change, rather than adenine protonation, is the rate-limiting step in the pathway of the conformational transition. Two nucleotides, U70 and U80, were found from chemical shift perturbation mapping to interact with the magnesium ion, with apparent K(d) values in the micromolar to millimolar range. These nucleotides also displayed metal ion-induced elevation of R1 rates, which can be explained by a model that assumes dynamic metal ion coordination concomitant with an induced higher shielding anisotropy for the base 13C nuclei. Addition of Mg2+ shifts the conformational equilibrium toward the high-pH (base-stacked) structure, accompanied by a significant drop in the apparent pK(a) of A79.

Site-directed spin labeling studies reveal solution conformational changes in a GAAA tetraloop receptor upon Mg(2+)-dependent docking of a GAAA tetraloop

The Mg(2+)-dependent GAAA tetraloop interaction with its 11 nucleotide receptor is one of the most frequently occurring long-range tertiary interactions in RNAs. To explore conformational changes in the receptor during tetraloop docking, nitroxide spin labels were attached at each of four uridine bases, one at a time, within an RNA molecule containing the receptor sequence. In the presence of Mg2+ and the tetraloop, the electron paramagnetic resonance (EPR) spectrum of one of the labeled bases reflected a large increase in mobility, indicating unstacking of the base upon tetraloop docking. This provides direct evidence that base unstacking is an intrinsic feature of the solution tetraloop-receptor complex formed in the presence of Mg2+. Additional evidence suggests that in solution the bound receptor conformation is similar to that observed in the crystal structure of a group I intron ribozyme domain. In Mg2+ alone, a receptor conformation with an unstacked base was not detectable, suggesting that this conformation is of higher standard state free energy than that of the free receptor. This leads to the conclusion that the extensive RNA-RNA interactions observed in the crystal structure of the tetraloop-receptor complex provide larger interaction energy than the measured apparent affinity between the tetraloop and the free receptor. This is compatible with a high specificity of the tetraloop-receptor interaction.

Identification of amino acids that promote specific and rigid TAR RNA-tat protein complex formation

The Tat protein and the transactivation responsive (TAR) RNA form an essential complex in the HIV lifecycle, and mutations in the basic region of the Tat protein alter this RNA-protein molecular recognition. Here, EPR spectroscopy was used to identify amino acids, flanking an essential arginine of the Tat protein, which contribute to specific and rigid TAR-Tat complex formation by monitoring changes in the mobility of nitroxide spin-labeled TAR RNA nucleotides upon binding. Arginine to lysine N-terminal mutations did not affect TAR RNA interfacial dynamics. In contrast, C-terminal point mutations, R56 in particular, affected the mobility of nucleotides U23 and U38, which are involved in a base-triple interaction in the complex. This report highlights the role of dynamics in specific molecular complex formation and demonstrates the ability of EPR spectroscopy to study interfacial dynamics of macromolecular complexes.

Iron responsive element RNA flexibility described by NMR and isotropic reorientational eigenmode dynamics

The first example of the application of reorientational eigenmode dynamics (RED) to RNA is shown here for the small and floppy Iron Responsive Element (IRE) RNA hairpin. Order parameters calculated for bases and riboses from a 12 ns molecular dynamics trajectory are compared to experimentally determined order parameters from 13C-1H NMR relaxation experiments, and shown to be in qualitative agreement. Given the small size of the IRE hairpin and its very flexible loop, isotropic RED (iRED) was also used to analyze the trajectory in order to describe its dynamic motions. iRED analysis shows that the global and internal dynamics of the IRE are not rigorously separable, which will result in inaccurate experimental order parameters. In addition, the iRED analysis described the many correlated motions that comprise the dynamics of the IRE RNA. The combined use of NMR relaxation, RED, and iRED provide a uniquely detailed description of IRE RNA dynamics.

A suite of 2H NMR spin relaxation experiments for the measurement of RNA dynamics

A suite of (2)H-based spin relaxation NMR experiments is presented for the measurement of molecular dynamics in a site-specific manner in uniformly (13)C, randomly fractionally deuterated ( approximately 50%) RNA molecules. The experiments quantify (2)H R(1) and R(2) relaxation rates that can subsequently be analyzed to obtain information about dynamics on a pico- to nanosecond time scale. Sensitivity permitting, the consistency of the data can be evaluated by measuring all five rates that are accessible for a spin 1 particle and establishing that the rates obey relations that are predicted from theory. The utility of the methodology is demonstrated with studies of the dynamics of a 14-mer RNA containing the UUCG tetraloop at temperatures of 25 and 5 degrees C. The high quality of the data, even at 5 degrees C, suggests that the experiments will be of use for the study of RNA molecules that are as large as 30 nucleotides.

Monitoring tat peptide binding to TAR RNA by solid-state 31P-19F REDOR NMR

Complexes of the HIV transactivation response element (TAR) RNA with the viral regulatory protein tat are of special interest due in particular to the plasticity of the RNA at this binding site and to the potential for therapeutic targeting of the interaction. We performed REDOR solid-state NMR experiments on lyophilized samples of a 29 nt HIV-1 TAR construct to measure conformational changes in the tat-binding site concomitant with binding of a short peptide comprising the residues of the tat basic binding domain. Peptide binding was observed to produce a nearly 4 A decrease in the separation between phosphorothioate and 2'F labels incorporated at A27 in the upper helix and U23 in the bulge, respectively, consistent with distance changes observed in previous solution NMR studies, and with models showing significant rearrangement in position of bulge residue U23 in the bound-form RNA. In addition to providing long-range constraints on free TAR and the TAR-tat complex, these results suggest that in RNAs known to undergo large deformations upon ligand binding, 31P-19F REDOR measurements can also serve as an assay for complex formation in solid-state samples. To our knowledge, these experiments provide the first example of a solid-state NMR distance measurement in an RNA-peptide complex.

A single intermolecular contact mediates intramolecular stabilization of both RNA and protein

An arginine-rich peptide from the Jembrana disease virus (JDV) Tat protein is a structural "chameleon" that binds bovine immunodeficiency virus (BIV) or HIV TAR RNAs in two different binding modes, with an affinity for BIV TAR even higher than the cognate BIV peptide. We determined the NMR structure of the JDV Tat-BIV TAR high-affinity complex and found that the C-terminal tyrosine in JDV Tat forms a network of inter- and intramolecular hydrogen bonding and stacking interactions that simultaneously stabilize the beta-hairpin conformation of the peptide and a base triple in the RNA. A neighboring histidine also appears to help stabilize the peptide conformation. Induced fit binding is recurrent in protein-protein and protein-nucleic acid interactions, and the JDV Tat complex demonstrates how high affinity can be achieved not only by optimization of the binding interface but also by inducing new intramolecular contacts that stabilize each binding partner. Comparison to the cognate BIV Tat peptide-TAR complex shows how such a costabilization mechanism can evolve with only small changes to the peptide sequence. In addition, the bound structure of BIV TAR in the chameleon peptide complex is strikingly similar to the bound conformation of HIV TAR, suggesting new strategies for the development of HIV TAR binding molecules.