Independent alignment of RNA for dynamic studies using residual dipolar couplings

More than Proton Detection-New Avenues for NMR Spectroscopy of RNA

Ribonucleic acid oligonucleotides (RNAs) play pivotal roles in cellular function (riboswitches), chemical biology applications (SELEX-derived aptamers), cell biology and biomedical applications (transcriptomics). Furthermore, a growing number of RNA forms (long non-coding RNAs, circular RNAs) but also RNA modifications are identified, showing the ever increasing functional diversity of RNAs. To describe and understand this functional diversity, structural studies of RNA are increasingly important. However, they are often more challenging than protein structural studies as RNAs are substantially more dynamic and their function is often linked to their structural transitions between alternative conformations. NMR is a prime technique to characterize these structural dynamics with atomic resolution. To extend the NMR size limitation and to characterize large RNAs and their complexes above 200 nucleotides, new NMR techniques have been developed. This Minireview reports on the development of NMR methods that utilize detection on low-γ nuclei (heteronuclei like 13 C or 15 N with lower gyromagnetic ratio than 1 H) to obtain unique structural and dynamic information for large RNA molecules in solution. Experiments involve through-bond correlations of nucleobases and the phosphodiester backbone of RNA for chemical shift assignment and make information on hydrogen bonding uniquely accessible. Previously unobservable NMR resonances of amino groups in RNA nucleobases are now detected in experiments involving conformational exchange-resistant double-quantum 1 H coherences, detected by 13 C NMR spectroscopy. Furthermore, 13 C and 15 N chemical shifts provide valuable information on conformations. All the covered aspects point to the advantages of low-γ nuclei detection experiments in RNA.

RNA Dynamics by NMR Spectroscopy

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

Applications of NMR to structure determination of RNAs large and small

Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool to investigate the structure and dynamics of RNA, because many biologically important RNAs have conformationally flexible structures, which makes them difficult to crystallize. Functional, independently folded RNA domains, range in size between simple stem-loops of as few as 10-20 nucleotides, to 50-70 nucleotides, the size of tRNA and many small ribozymes, to a few hundred nucleotides, the size of more complex RNA enzymes and of the functional domains of non-coding transcripts. In this review, we discuss new methods for sample preparation, assignment strategies and structure determination for independently folded RNA domains of up to 100 kDa in molecular weight.

Structure of a low-population binding intermediate in protein-RNA recognition

The interaction of the HIV-1 protein transactivator of transcription (Tat) and its cognate transactivation response element (TAR) RNA transactivates viral transcription and represents a paradigm for the widespread occurrence of conformational rearrangements in protein-RNA recognition. Although the structures of free and bound forms of TAR are well characterized, the conformations of the intermediates in the binding process are still unknown. By determining the free energy landscape of the complex using NMR residual dipolar couplings in replica-averaged metadynamics simulations, we observe two low-population intermediates. We then rationally design two mutants, one in the protein and another in the RNA, that weaken specific nonnative interactions that stabilize one of the intermediates. By using surface plasmon resonance, we show that these mutations lower the release rate of Tat, as predicted. These results identify the structure of an intermediate for RNA-protein binding and illustrate a general strategy to achieve this goal with high resolution.

Structural Analysis of Multi-Helical RNAs by NMR-SAXS/WAXS: Application to the U4/U6 di-snRNA

NMR and SAXS (small-angle X-ray scattering)/WAXS (wide-angle X-ray scattering) are highly complementary approaches for the analysis of RNA structure in solution. Here we describe an efficient NMR-SAXS/WAXS approach for structural investigation of multi-helical RNAs. We illustrate this approach by determining the overall fold of a 92-nt 3-helix junction from the U4/U6 di-snRNA. The U4/U6 di-snRNA is conserved in eukaryotes and is part of the U4/U6.U5 tri-snRNP, a large ribonucleoprotein complex that comprises a major subunit of the assembled spliceosome. Helical orientations can be determined by X-ray scattering data alone, but the addition of NMR RDC (residual dipolar coupling) restraints improves the structure models. RDCs were measured in two different external alignment media and also by magnetic susceptibility anisotropy. The resulting alignment tensors are collinear, which is a previously noted problem for nucleic acids. Including WAXS data in the calculations produces models with significantly better fits to the scattering data. In solution, the U4/U6 di-snRNA forms a 3-helix junction with a planar Y-shaped structure and has no detectable tertiary interactions. Single-molecule Förster resonance energy transfer data support the observed topology. A comparison with the recently determined cryo-electron microscopy structure of the U4/U6.U5 tri-snRNP illustrates how proteins scaffold the RNA and dramatically alter the geometry of the U4/U6 3-helix junction.

Modulating RNA Alignment Using Directional Dynamic Kinks: Application in Determining an Atomic-Resolution Ensemble for a Hairpin using NMR Residual Dipolar Couplings

Approaches that combine experimental data and computational molecular dynamics (MD) to determine atomic resolution ensembles of biomolecules require the measurement of abundant experimental data. NMR residual dipolar couplings (RDCs) carry rich dynamics information, however, difficulties in modulating overall alignment of nucleic acids have limited the ability to fully extract this information. We present a strategy for modulating RNA alignment that is based on introducing variable dynamic kinks in terminal helices. With this strategy, we measured seven sets of RDCs in a cUUCGg apical loop and used this rich data set to test the accuracy of an 0.8 μs MD simulation computed using the Amber ff10 force field as well as to determine an atomic resolution ensemble. The MD-generated ensemble quantitatively reproduces the measured RDCs, but selection of a sub-ensemble was required to satisfy the RDCs within error. The largest discrepancies between the RDC-selected and MD-generated ensembles are observed for the most flexible loop residues and backbone angles connecting the loop to the helix, with the RDC-selected ensemble resulting in more uniform dynamics. Comparison of the RDC-selected ensemble with NMR spin relaxation data suggests that the dynamics occurs on the ps-ns time scales as verified by measurements of R(1ρ) relaxation-dispersion data. The RDC-satisfying ensemble samples many conformations adopted by the hairpin in crystal structures indicating that intrinsic plasticity may play important roles in conformational adaptation. The approach presented here can be applied to test nucleic acid force fields and to characterize dynamics in diverse RNA motifs at atomic resolution.

Insights into domain-domain motions in proteins and RNA from solution NMR

Many multidomain proteins and ribonucleic acids consist of domains that autonomously fold and that are linked together by flexible junctions. This architectural design allows domains to sample a wide range of positions with respect to one another, yet do so in a way that retains structural specificity, since the number of sampled conformations remains extremely small compared to the total conformations that would be sampled if the domains were connected by an infinitely long linker. This “tuned” flexibility in interdomain conformation is in turn used in many biochemical processes.

There is great interest in characterizing the dynamic properties of multidomain systems, and moving beyond conventional descriptions in terms of static structures, toward the characterization of population-weighted ensembles describing a distribution of many conformations sampled in solution. There is also great interest in understanding the design principles and underlying physical and chemical interactions that specify the nature of interdomain flexibility. NMR spectroscopy is one of the most powerful techniques for characterizing motions in complex biomolecules and has contributed greatly toward our basic understanding of dynamics in proteins and nucleic acids and its role in folding, recognition, and signaling.

Here, we review methods that have been developed in our laboratories to address these challenges. Our approaches are based on the ability of one domain of the molecule to self-align in a magnetic field, or to dominate the overall orientation of the molecule, so that the conformational freedom of other domains can be assessed by their degree of alignment induced by the aligned part. In turn, this self-alignment ability can be intrinsic or can be caused by tagging appropriate constructs to the molecule of interest. In general, self-alignment is due to magnetic susceptibility anisotropy. Nucleic acids with elongated helices have this feature, as well as several paramagnetic metal centers that can be found in, or attached to, a protein domain.

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.

Advances in the determination of nucleic acid conformational ensembles

Conformational changes in nucleic acids play a key role in the way genetic information is stored, transferred, and processed in living cells. Here, we describe new approaches that employ a broad range of experimental data, including NMR-derived chemical shifts and residual dipolar couplings, small-angle X-ray scattering, and computational approaches such as molecular dynamics simulations to determine ensembles of DNA and RNA at atomic resolution. We review the complementary information that can be obtained from diverse sets of data and the various methods that have been developed to combine these data with computational methods to construct ensembles and assess their uncertainty. We conclude by surveying RNA and DNA ensembles determined using these methods, highlighting the unique physical and functional insights obtained so far.

NMR studies of nucleic acid dynamics

Nucleic acid structures have to satisfy two diametrically opposite requirements; on one hand they have to adopt well-defined 3D structures that can be specifically recognized by proteins; on the other hand, their structures must be sufficiently flexible to undergo very large conformational changes that are required during key biochemical processes, including replication, transcription, and translation. How do nucleic acids introduce flexibility into their 3D structure without losing biological specificity? Here, I describe the development and application of NMR spectroscopic techniques in my laboratory for characterizing the dynamic properties of nucleic acids that tightly integrate a broad set of NMR measurements, including residual dipolar couplings, spin relaxation, and relaxation dispersion with sample engineering and computational approaches. This approach allowed us to obtain fundamental new insights into directional flexibility in nucleic acids that enable their structures to change in a very specific functional manner.