NMR studies of RNA dynamics and structural plasticity using NMR residual dipolar couplings

Magnetic control of self-assembly and disassembly in organic materials

Because organic molecules and materials are generally insensitive or weakly sensitive to magnetic fields, a certain means to enhance their magnetic responsiveness needs to be exploited. Here we show a strategy to amplify the magnetic responsiveness of self-assembled peptide nanostructures by synergistically combining the concepts of perfect α-helix and rod-coil supramolecular building blocks. Firstly, we develop a monomeric, nonpolar, and perfect α-helix (MNP-helix). Then, we employ the MNP-helix as the rod block of rod-coil amphiphiles (rod-coils) because rod-coils are well-suited for fabricating responsive assemblies. We show that the self-assembly processes of the designed rod-coils and disassembly of rod-coil/DNA complexes can be controlled in a magnetically responsive manner using the relatively weak magnetic field provided by the ordinary neodymium magnet [0.07 ~ 0.25 Tesla (T)]. These results demonstrate that magnetically responsive organic assemblies usable under practical conditions can be realized by using rod-coil supramolecular building blocks containing constructively organized diamagnetic moieties.

Weakly aligned Ti3C2Tx MXene liquid crystals: measuring residual dipolar coupling in multiple co-solvent systems

Residual Dipolar Coupling (RDC), acquired relying on weakly alignment media, is highly valuable for the structural elucidation of organic molecules. Arising from the striking features of no background signals and low critical concentrations, two-dimensional (2D) liquid crystals (LCs) show the clear advantages of acting as alignment media to measure RDCs. So far, creating multisolvent compatible 2D LC media through a simple and versatile method is still formidably challenging. Herein, we report the rapid creation of aligned media based on the Ti₃C₂T_x MXene, which self-aligned in multiple co-solvents including CH₃OH-H₂O, DMSO-H₂O, DMF-H₂O, and acetone-H₂O. We demonstrated the applicability of these aligned media for the RDC measurement of small organic molecules with different polarities and solubilities. Notably, Ti₃C₂T_x MXene LCs without chemical modification enabled RDC measurements on aromatic molecules. The straightforward preparation of Ti₃C₂T_x media and its compatibility with multiple solvents will push RDC measurement as a routine methodology for structural elucidation. It may also facilitate the investigation of solvation effects on conformational dynamics.

Advances in NMR Spectroscopy of Weakly Aligned Biomolecular Systems

The measurement and application of residual dipolar couplings (RDCs) in solution NMR studies of biological macromolecules has become well established over the past quarter of a century. Numerous methods for generating the requisite anisotropic orientational molecular distribution have been demonstrated, each with its specific strengths and weaknesses. In parallel, an enormous number of pulse schemes have been introduced to measure the many different types of RDCs, ranging from the most widely measured backbone amide ¹⁵N-¹H RDCs, to ¹H-¹H RDCs and couplings between low-γ nuclei. Applications of RDCs range from structure validation and refinement to the determination of relative domain orientations, the measurement of backbone and domain motions, and de novo structure determination. Nevertheless, it appears that the power of the RDC methodology remains underutilized. This review aims to highlight the practical aspects of sample preparation and RDC measurement while describing some of the most straightforward applications that take advantage of the exceptionally precise information contained in such data. Some emphasis will be placed on more recent developments that enable the accurate measurement of RDCs in larger systems, which is key to the ongoing shift in focus of biological NMR spectroscopy from structure determination toward gaining improved understanding of how molecular flexibility drives protein function.

Different views of the dynamic landscape covered by the 5'-hairpin of the 7SK small nuclear RNA

The 7SK small nuclear RNA (7SKsnRNA) plays a key role in the regulation of RNA polymerase II by sequestrating and inhibiting the positive transcription elongation factor b (P-TEFb) in the 7SK ribonucleoprotein complex (7SKsnRNP), a process mediated by interaction with the protein HEXIM. P-TEFb is also an essential cellular factor recruited by the viral protein Tat to ensure the replication of the viral RNA in the infection cycle of the human immunodeficiency virus (HIV-1). Tat promotes the release of P-TEFb from the 7SKsnRNP and subsequent activation of transcription, by displacing HEXIM from the 5'-hairpin of the 7SKsnRNA. This hairpin (HP1), comprising the signature sequence of the 7SKsnRNA, has been the subject of three independent structural studies aimed at identifying the structural features that could drive the recognition by the two proteins, both depending on arginine-rich motifs (ARM). Interestingly, four distinct structures were determined. In an attempt to provide a comprehensive view of the structure-function relationship of this versatile RNA, we present here a structural analysis of the models, highlighting how HP1 is able to adopt distinct conformations with significant impact on the compactness of the molecule. Since these models are solved under different conditions by nuclear magnetic resonance (NMR) and crystallography, the impact of the buffer composition on the conformational variation was investigated by complementary biophysical approaches. Finally, using isothermal titration calorimetry, we determined the thermodynamic signatures of the Tat-ARM and HEXIM-ARM peptide interactions with the RNA, showing that they are associated with distinct binding mechanisms.

Reduction in Dynamics of Base pair Opening upon Ligand Binding by the Cocaine-Binding Aptamer

We have used magnetization transfer NMR experiments to measure the exchange rate constant (k_ex) of the imino protons in the unbound, cocaine-bound, and quinine-bound forms of the cocaine-binding DNA aptamer. Both long-stem 1 (MN4) and short-stem 1 (MN19) variants were analyzed, corresponding to structures with a prefolded secondary structure and ligand-induced-folding versions of this aptamer, respectively. The k_ex values were measured as a function of temperature from 5 to 45°C to determine the thermodynamics of the base pair opening for MN4. We find that the base pairs close to the ligand-binding site become stronger upon ligand binding, whereas those located away from the binding site do not strengthen. With the buffer conditions used in this study, we observe imino ¹H signals in MN19 not previously seen, which leads us to conclude that in the free form, both stem 2 and parts of stem 3 are formed and that the base pairs in stem 1 become structured or more rigid upon binding. This is consistent with the k_ex values for MN19 decreasing in both stem 1 and at the ligand-binding site. Based on the temperature dependence of the k_ex values, we find that MN19 is more dynamic than MN4 in the free and both ligand-bound forms. For MN4, ligand-binding results in the reduction of dynamics that are localized to the binding site. These results demonstrate that an aptamer in which the base pairs are preformed also experiences a reduction in dynamics with ligand binding.

Using magnetic coupling to improve multiple resonance NMR probe circuits

Multiple resonant NMR circuits are an especially important aspect of modern NMR detectors. Today NMR researchers can carry out direct observe, indirect observe and double or triple resonance experiments on many different nuclides without having to change the NMR probe. The price for this convenience is loss of efficiency. In this paper it is explained how by using magnetic coupling and symmetric multiple resonance circuits the efficiency of these experiments can be improved over capacitively coupled circuits to 50 Ω ports.

Site-Specific Spin Labeling of RNA for NMR and EPR Structural Studies

Many RNA architectures were discovered to be involved in essential biological pathways acting as catalysts and/or regulators of gene expression, transcription, translation, splicing, or viral infection. The key to understand their diverse biological functions is to investigate their structure and dynamic. Nuclear Magnetic Resonance (NMR) is a powerful method to gain insight into these properties. However, the study of high-molecular-weight RNAs by NMR remains challenging. Advances in biochemical and NMR methods over the recent years allow to overcome the limitation of NMR. In particular, the incorporation of paramagnetic probes, coupled to the measurement of the induced effects on nuclear spins, has become an efficient tool providing long-range distance restraints and information on dynamic in solution. At the same time, the use of spin label enabled the application of Electron Paramagnetic Resonance (EPR) to study biological macromolecules. Combining NMR and EPR is emerging as a new approach to investigate the architecture of biological systems.Here, we describe an efficient protocol to introduce a paramagnetic probe into a RNA at a specific position. This method enables various combinations of isotopic labeling for NMR and is also of interest for EPR studies.

Accuracy of MD solvent models in RNA structure refinement assessed via liquid-crystal NMR and spin relaxation data

Molecular dynamics (MD) simulations play an important role in characterizing Ribonucleic Acid (RNA) structure, augmenting information from experimental techniques such as Nuclear Magnetic Resonance (NMR). In this work, we examine the accuracy of structural representation resulting from application of a number of explicit and implicit solvent models and refinement protocols against experimental data ranging from high density of residual dipolar coupling (RDC) restraints to completely unrestrained simulations. For a prototype A-form RNA helix, our results indicate that AMBER RNA force field with either implicit or explicit solvent can produce a realistic dynamic representation of RNA helical structure, accurately cross-validating with respect to a diverse array of NMR observables. In refinement against NMR distance restraints, modern MD force fields are found to be equally adequate, with high fidelity cross-validation to the residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs), while slightly over-estimating structural order as monitored via NMR relaxation data. With restraints trimmed to encode only for base pairing information, cross-validation quality significantly deteriorates, now exhibiting a pronounced dependence on the choice of the solvent model. This deterioration is found to be partially reversible by increasing planarity restraints on the nucleobase geometry. For completely unrestrained MD simulations, the choice of water model becomes very important, with the best-performing TIP4P-Ew accurately reproducing both the RDC and RCSA data, while closely matching the NMR-derived order parameters. The information provided here will serve as a foundation for MD-based refinement of solution state NMR structures of RNA.

Maximizing accuracy of RNA structure in refinement against residual dipolar couplings

Structural information about ribonucleic acid (RNA) is lagging behind that of proteins, in part due to its high charge and conformational variability. Molecular dynamics (MD) has played an important role in describing RNA structure, complementing information from both nuclear magnetic resonance (NMR), or X-ray crystallography. We examine the impact of the choice of the empirical force field for RNA structure refinement using cross-validation against residual dipolar couplings (RDCs) as structural accuracy reporter. Four force fields, representing both the state-of-the art in RNA simulation and the most popular selections in NMR structure determination, are compared for a prototypical A-RNA helix. RNA structural accuracy is also evaluated as a function of both density and nature of input NMR data including RDCs, anisotropic chemical shifts, and distance restraints. Our results show a complex interplay between the experimental restraints and the force fields indicating two best-performing choices: high-fidelity refinement in explicit solvent, and the conformational database-derived potentials. Accuracy of RNA models closely tracks the density of 1-bond C-H RDCs, with other data types having beneficial, but smaller effects. At lower RDC density, or when refining against NOEs only, the two selected force fields are capable of accurately describing RNA helices with little or no experimental RDC data, making them available for the higher order structure assembly or better quantification of the intramolecular dynamics. Unrestrained simulations of simple RNA motifs with state-of-the art MD force fields appear to capture the flexibility inherent in nucleic acids while also maintaining a good agreement with the experimental observables.

Interactions between cyclic nucleotides and common cations: an ab initio molecular dynamics study

We present the first, to the best of our knowledge, ab initio molecular dynamics (AIMD) investigation on three aqueous solutions where an abasic cyclic nucleotide model is solvated in the presence of distinct cations (i.e., Na⁺, K⁺ and Mg²⁺). We elucidate the typical modalities of interaction between those ionic species and the nucleotide moiety by first-principles numerical simulations, starting from an inner-shell binding configuration on a time scale of 100 ps (total simulation time of ∼600 ps). Whereas the strong "structure-maker" Mg²⁺ is permanently bound to one of the two oxygen atoms of the phosphate group of the nucleotide model, Na⁺ and K⁺ show binding times τ_b of 65 ps and 10-15 ps, respectively, thus reflecting their chemical nature in aqueous solutions. Furthermore, we qualitatively relate these findings to approximate free-energy barriers of the cations' unbinding obtained by means of exploratory well-tempered metadynamics. With the aim of shedding light on the features of commonly employed force-fields (FFs), classical MD simulations (almost 200 trajectories with a total simulation time of ∼18 μs) using the biomolecular AMBER FF are also reported. By choosing several combinations of the parametrization for the water environment (i.e., TIP3P, SPC/E and OPC) and cations (i.e., Joung-Cheatham, Li-Merz 12-6 and Li-Merz 12-6-4), we found significant differences in the radial distribution functions and residence times compared to the ab initio results. The Na⁺ and K⁺ ions wrongly show quasi-identical radial distribution functions and the Li & Merz 12-6-4 Lennard-Jones parameters for Mg²⁺ were found to be essential in quickly reaching the binding state consistent with AIMD.

RNA Structural Differentiation: Opportunities with Pattern Recognition

Our awareness and appreciation of the many regulatory roles of RNA have dramatically increased in the past decade. This understanding, in addition to the impact of RNA in many disease states, has renewed interest in developing selective RNA-targeted small molecule probes. However, the fundamental guiding principles in RNA molecular recognition that could accelerate these efforts remain elusive. While high-resolution structural characterization can provide invaluable insight, examples of well-characterized RNA structures, not to mention small molecule:RNA complexes, remain limited. This Perspective provides an overview of the current techniques used to understand RNA molecular recognition when high-resolution structural information is unavailable. We will place particular emphasis on a new method, pattern recognition of RNA with small molecules (PRRSM), that provides rapid insight into critical components of RNA recognition and differentiation by small molecules as well as into RNA structural features.

Comprehensive review of chemical strategies for the preparation of new aminoglycosides and their biological activities

A systematic analysis of all synthetic and chemoenzymatic methodologies for the preparation of aminoglycosides for a variety of applications (therapeutic and agricultural) reported in the scientific literature up to 2017 is presented. This comprehensive analysis of derivatization/generation of novel aminoglycosides and their conjugates is divided based on the types of modifications used to make the new derivatives. Both the chemical strategies utilized and the biological results observed are covered. Structure-activity relationships based on different synthetic modifications along with their implications for activity and ability to avoid resistance against different microorganisms are also presented.

Determination of the conformational ensemble of the TAR RNA by X-ray scattering interferometry

The conformational ensembles of structured RNA's are crucial for biological function, but they remain difficult to elucidate experimentally. We demonstrate with HIV-1 TAR RNA that X-ray scattering interferometry (XSI) can be used to determine RNA conformational ensembles. X-ray scattering interferometry (XSI) is based on site-specifically labeling RNA with pairs of heavy atom probes, and precisely measuring the distribution of inter-probe distances that arise from a heterogeneous mixture of RNA solution structures. We show that the XSI-based model of the TAR RNA ensemble closely resembles an independent model derived from NMR-RDC data. Further, we show how the TAR RNA ensemble changes shape at different salt concentrations. Finally, we demonstrate that a single hybrid model of the TAR RNA ensemble simultaneously fits both the XSI and NMR-RDC data set and show that XSI can be combined with NMR-RDC to further improve the quality of the determined ensemble. The results suggest that XSI-RNA will be a powerful approach for characterizing the solution conformational ensembles of RNAs and RNA-protein complexes under diverse solution conditions.

Development and Implementation of an HTS-Compatible Assay for the Discovery of Selective Small-Molecule Ligands for Pre-microRNAs

microRNAs (miRNAs) are small gene regulatory RNAs, and their expression has been found to be dysregulated in a number of human diseases. To facilitate the discovery of small molecules capable of selectively modulating the activity of a specific miRNA, we have utilized new high-throughput screening technology targeting Dicer-mediated pre-miRNA maturation. Pilot screening of ~50,000 small molecules and ~33,000 natural product extract libraries against pre-miR-21 processing indicated the potential of our assay for this goal, yielding a campaign Z' factor of 0.52 and an average plate signal-to-background (S/B) ratio of 13. Using two-dimensional screening against a second pre-miRNA, pre-let-7d, we evaluated the selectivity of confirmed hits. The results presented demonstrate how high-throughput screening can be used to identify selective small molecules for a target RNA.

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

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

Structure of RNA Stem Loop B from the Picornavirus Replication Platform

The presumptive RNA cloverleaf at the start of the 5'-untranslated region of the picornavirus genome is an essential element in replication. Stem loop B (SLB) of the cloverleaf is a recognition site for the host polyC-binding protein, which initiates a switch from translation to replication. Here we present the solution structure of human rhinovirus isotype 14 SLB using nuclear magnetic resonance spectroscopy. SLB adopts a predominantly A-form helical structure. The stem contains five Watson-Crick base pairs and one wobble base pair and is capped by an eight-nucleotide loop. The wobble base pair introduces perturbations into the helical parameters but does not appear to introduce flexibility. However, the helix major groove appears to be accessible. Flexibility is seen throughout the loop and in the terminal nucleotides. The pyrimidine-rich region of the loop, the apparent recognition site for the polyC-binding protein, is the most disordered region of the structure.

Structure Calculation and Reconstruction of Discrete-State Dynamics from Residual Dipolar Couplings

Residual dipolar couplings (RDCs) acquired by nuclear magnetic resonance (NMR) spectroscopy are an indispensable source of information in investigation of molecular structures and dynamics. Here, we present a comprehensive strategy for structure calculation and reconstruction of discrete-state dynamics from RDC data that is based on the singular value decomposition (SVD) method of order tensor estimation. In addition to structure determination, we provide a mechanism of producing an ensemble of conformations for the dynamical regions of a protein from RDC data. The developed methodology has been tested on simulated RDC data with ±1 Hz of error from an 83 residue α protein (PDB ID 1A1Z ) and a 213 residue α/β protein DGCR8 (PDB ID 2YT4 ). In nearly all instances, our method reproduced the structure of the protein including the conformational ensemble to within less than 2 Å. On the basis of our investigations, arc motions with more than 30° of rotation are identified as internal dynamics and are reconstructed with sufficient accuracy. Furthermore, states with relative occupancies above 20% are consistently recognized and reconstructed successfully. Arc motions with a magnitude of 15° or relative occupancy of less than 10% are consistently unrecognizable as dynamical regions within the context of ±1 Hz of error.

NMR Methods for Characterization of RNA Secondary Structure

Knowledge of RNA secondary structure is often sufficient to identify relationships between the structure of RNA and processing pathways, and the design of therapeutics. Nuclear magnetic resonance (NMR) can identify types of nucleotide base pairs and the sequence, thus limiting possible secondary structures. Because NMR experiments, like chemical mapping, are performed in solution, not in single crystals, experiments can be initiated as soon as the biomolecule is expressed and purified. This chapter summarizes NMR methods that permit rapid identification of RNA secondary structure, information that can be used as supplements to chemical mapping, and/or as preliminary steps required for 3D structure determination. The primary aim is to provide guidelines to enable a researcher with minimal knowledge of NMR to quickly extract secondary structure information from basic datasets. Instrumental and sample considerations that can maximize data quality are discussed along with some details for optimal data acquisition and processing parameters. Approaches for identifying base pair types in both unlabeled and isotopically labeled RNA are covered. Common problems, such as missing signals and overlaps, and approaches to address them are considered. Programs under development for merging NMR data with structure prediction algorithms are briefly discussed.

The "Speedy" Synthesis of Atom-Specific (15)N Imino/Amido-Labeled RNA

Although numerous reports on the synthesis of atom-specific (15)N-labeled nucleosides exist, fast and facile access to the corresponding phosphoramidites for RNA solid-phase synthesis is still lacking. This situation represents a severe bottleneck for NMR spectroscopic investigations on functional RNAs. Here, we present optimized procedures to speed up the synthesis of (15)N(1) adenosine and (15)N(1) guanosine amidites, which are the much needed counterparts of the more straightforward-to-achieve (15)N(3) uridine and (15)N(3) cytidine amidites in order to tap full potential of (1)H/(15)N/(15)N-COSY experiments for directly monitoring individual Watson-Crick base pairs in RNA. Demonstrated for two preQ1 riboswitch systems, we exemplify a versatile concept for individual base-pair labeling in the analysis of conformationally flexible RNAs when competing structures and conformational dynamics are encountered.

Dynamic motions of the HIV-1 frameshift site RNA

The HIV-1 frameshift site (FS) plays a critical role in viral replication. During translation, the HIV-1 FS transitions from a 3-helix to a 2-helix junction RNA secondary structure. The 2-helix junction structure contains a GGA bulge, and purine-rich bulges are common motifs in RNA secondary structure. Here, we investigate the dynamics of the HIV-1 FS 2-helix junction RNA. Interhelical motions were studied under different ionic conditions using NMR order tensor analysis of residual dipolar couplings. In 150 mM potassium, the RNA adopts a 43°(±4°) interhelical bend angle (β) and displays large amplitude, anisotropic interhelical motions characterized by a 0.52(±0.04) internal generalized degree of order (GDOint) and distinct order tensor asymmetries for its two helices (η = 0.26(±0.04) and 0.5(±0.1)). These motions are effectively quenched by addition of 2 mM magnesium (GDOint = 0.87(±0.06)), which promotes a near-coaxial conformation (β = 15°(±6°)) of the two helices. Base stacking in the bulge was investigated using the fluorescent purine analog 2-aminopurine. These results indicate that magnesium stabilizes extrahelical conformations of the bulge nucleotides, thereby promoting coaxial stacking of helices. These results are highly similar to previous studies of the HIV transactivation response RNA, despite a complete lack of sequence similarity between the two RNAs. Thus, the conformational space of these RNAs is largely determined by the topology of their interhelical junctions.

A fully enzymatic method for site-directed spin labeling of long RNA

Site-directed spin labeling is emerging as an essential tool to investigate the structural and dynamical features of RNA. We propose here an enzymatic method, which allows the insertion of a paramagnetic center at a specific position in an RNA molecule. The technique is based on a segmental approach using a ligation protocol with T4 RNA ligase 2. One transcribed acceptor RNA is ligated to a donor RNA in which a thio-modified nucleotide is introduced at its 5′-end by in vitro transcription with T7 RNA polymerase. The paramagnetic thiol-specific reagent is subsequently attached to the RNA ligation product. This novel strategy is demonstrated by introducing a paramagnetic probe into the 55 nucleotides long RNA corresponding to K-turn and Specifier Loop domains from the Bacillus subtilis tyrS T-Box leader RNA. The efficiency of the coupling reaction and the quality of the resulting spin-labeled RNA were assessed by Mass Spectrometry, Electron Paramagnetic Resonance (EPR) and Nuclear Magnetic Resonance (NMR). This method enables various combinations of isotopic segmental labeling and spin labeling schemes, a strategy that will be of particular interest to investigate the structural and dynamical properties of large RNA complexes by NMR and EPR spectroscopies.

Structural dynamics of a single-stranded RNA-helix junction using NMR

Many regulatory RNAs contain long single strands (ssRNA) that adjoin secondary structural elements. Here, we use NMR spectroscopy to study the dynamic properties of a 12-nucleotide (nt) ssRNA tail derived from the prequeuosine riboswitch linked to the 3' end of a 48-nt hairpin. Analysis of chemical shifts, NOE connectivity, (13)C spin relaxation, and residual dipolar coupling data suggests that the first two residues (A25 and U26) in the ssRNA tail stack onto the adjacent helix and assume an ordered conformation. The following U26-A27 step marks the beginning of an A6-tract and forms an acute pivot point for substantial motions within the tail, which increase toward the terminal end. Despite substantial internal motions, the ssRNA tail adopts, on average, an A-form helical conformation that is coaxial with the helix. Our results reveal a surprising degree of structural and dynamic complexity at the ssRNA-helix junction, which involves a fine balance between order and disorder that may facilitate efficient pseudoknot formation on ligand recognition.

Structural determinants for ligand capture by a class II preQ1 riboswitch

Prequeuosine (preQ1) riboswitches are RNA regulatory elements located in the 5' UTR of genes involved in the biosynthesis and transport of preQ1, a precursor of the modified base queuosine universally found in four tRNAs. The preQ1 class II (preQ1-II) riboswitch regulates preQ1 biosynthesis at the translational level. We present the solution NMR structure and conformational dynamics of the 59 nucleotide Streptococcus pneumoniae preQ1-II riboswitch bound to preQ1. Unlike in the preQ1 class I (preQ1-I) riboswitch, divalent cations are required for high-affinity binding. The solution structure is an unusual H-type pseudoknot featuring a P4 hairpin embedded in loop 3, which forms a three-way junction with the other two stems. (13)C relaxation and residual dipolar coupling experiments revealed interhelical flexibility of P4. We found that the P4 helix and flanking adenine residues play crucial and unexpected roles in controlling pseudoknot formation and, in turn, sequestering the Shine-Dalgarno sequence. Aided by divalent cations, P4 is poised to act as a "screw cap" on preQ1 recognition to block ligand exit and stabilize the binding pocket. Comparison of preQ1-I and preQ1-II riboswitch structures reveals that whereas both form H-type pseudoknots and recognize preQ1 using one A, C, or U nucleotide from each of three loops, these nucleotides interact with preQ1 differently, with preQ1 inserting into different grooves. Our studies show that the preQ1-II riboswitch uses an unusual mechanism to harness exquisite control over queuosine metabolism.

A method of determining RNA conformational ensembles using structure-based calculations of residual dipolar couplings

We describe a method of determining the conformational fluctuations of RNA based on the incorporation of nuclear magnetic resonance (NMR) residual dipolar couplings (RDCs) as replica-averaged structural restraints in molecular dynamics simulations. In this approach, the alignment tensor required to calculate the RDCs corresponding to a given conformation is estimated from its shape, and multiple replicas of the RNA molecule are simulated simultaneously to reproduce in silico the ensemble-averaging procedure performed in the NMR measurements. We provide initial evidence that with this approach it is possible to determine accurately structural ensembles representing the conformational fluctuations of RNA by applying the reference ensemble test to the trans-activation response element of the human immunodeficiency virus type 1.

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.

Advances in the REDCAT software package

Background

Residual Dipolar Couplings (RDCs) have emerged in the past two decades as an informative source of experimental restraints for the study of structure and dynamics of biological macromolecules and complexes. The REDCAT software package was previously introduced for the analysis of molecular structures using RDC data. Here we report additional features that have been included in this software package in order to expand the scope of its analyses. We first discuss the features that enhance REDCATs user-friendly nature, such as the integration of a number of analyses into one single operation and enabling convenient examination of a structural ensemble in order to identify the most suitable structure. We then describe the new features which expand the scope of RDC analyses, performing exercises that utilize both synthetic and experimental data to illustrate and evaluate different features with regard to structure refinement and structure validation.

Results

We establish the seamless interaction that takes place between REDCAT, VMD, and Xplor-NIH in demonstrations that utilize our newly developed REDCAT-VMD and XplorGUI interfaces. These modules enable visualization of RDC analysis results on the molecular structure displayed in VMD and refinement of structures with Xplor-NIH, respectively. We also highlight REDCAT’s Error-Analysis feature in reporting the localized fitness of a structure to RDC data, which provides a more effective means of recognizing local structural anomalies. This allows for structurally sound regions of a molecule to be identified, and for any refinement efforts to be focused solely on locally distorted regions.

Conclusions

The newly engineered REDCAT software package, which is available for download via the WWW from http://ifestos.cse.sc.edu, has been developed in the Object Oriented C++ environment. Our most recent enhancements to REDCAT serve to provide a more complete RDC analysis suite, while also accommodating a more user-friendly experience, and will be of great interest to the community of researchers and developers since it hides the complications of software development.

Probing protein ensemble rigidity and hydrogen-deuterium exchange

Protein rigidity and flexibility can be analyzed accurately and efficiently using the program floppy inclusion and rigid substructure topography (FIRST). Previous studies using FIRST were designed to analyze the rigidity and flexibility of proteins using a single static (snapshot) structure. It is however well known that proteins can undergo spontaneous sub-molecular unfolding and refolding, or conformational dynamics, even under conditions that strongly favor a well-defined native structure. These (local) unfolding events result in a large number of conformers that differ from each other very slightly. In this context, proteins are better represented as a thermodynamic ensemble of 'native-like' structures, and not just as a single static low-energy structure. Working with this notion, we introduce a novel FIRST-based approach for predicting rigidity/flexibility of the protein ensemble by (i) averaging the hydrogen bonding strengths from the entire ensemble and (ii) by refining the mathematical model of hydrogen bonds. Furthermore, we combine our FIRST-ensemble rigidity predictions with the ensemble solvent accessibility data of the backbone amides and propose a novel computational method which uses both rigidity and solvent accessibility for predicting hydrogen-deuterium exchange (HDX). To validate our predictions, we report a novel site specific HDX experiment which characterizes the native structural ensemble of Acylphosphatase from hyperthermophile Sulfolobus solfataricus (Sso AcP). The sub-structural conformational dynamics that is observed by HDX data, is closely matched with the FIRST-ensemble rigidity predictions, which could not be attained using the traditional single 'snapshot' rigidity analysis. Moreover, the computational predictions of regions that are protected from HDX and those that undergo exchange are in very good agreement with the experimental HDX profile of Sso AcP.

A general method for constructing atomic-resolution RNA ensembles using NMR residual dipolar couplings: the basis for interhelical motions revealed

The ability to modulate alignment and measure multiple independent sets of NMR residual dipolar couplings (RDCs) has made it possible to characterize internal motions in proteins at atomic resolution and with time scale sensitivity ranging from picoseconds up to milliseconds. The application of such methods to the study of RNA dynamics, however, remains fundamentally limited by the inability to modulate alignment and by strong couplings between internal and overall motions that complicate the quantitative interpretation of RDCs. Here, we address this problem by showing that RNA alignment can be generally modulated, in a controlled manner, by variable elongation of A-form helices and that the information contained within the measured RDCs can be extracted even in the presence of strong couplings between motions and overall alignment via structure-based prediction of alignment. Using this approach, four RDC data sets, and a broad conformational pool obtained from a 8.2 μs molecular dynamics simulation, we successfully construct and validate an atomic resolution ensemble of human immunodeficiency virus type I transactivation response element RNA. This ensemble reveals local motions in and around the bulge involving changes in stacking and hydrogen-bonding interactions, which are undetectable by traditional spin relaxation and drive global changes in interhelical orientation. This new approach broadens the scope of using RDCs in characterizing the dynamics of nucleic acids.

The internal dynamics of mini c TAR DNA probed by electron paramagnetic resonance of nitroxide spin-labels at the lower stem, the loop, and the bulge

Electron paramagnetic resonance (EPR) at 236.6 and 9.5 GHz probed the tumbling of nitroxide spin probes in the lower stem, in the upper loop, and near the bulge of mini c TAR DNA. High-frequency 236.6 GHz EPR, not previously applied to spin-labeled oligonucleotides, was notably sensitive to fast, anisotropic, hindered local rotational motion of the spin probe, occurring approximately about the NO nitroxide axis. Labels attached to the 2'-aminocytidine sugar in the mini c TAR DNA showed such anisotropic motion, which was faster in the lower stem, a region previously thought to be partially melted. More flexible labels attached to phosphorothioates at the end of the lower stem tumbled isotropically in mini c TAR DNA, mini TAR RNA, and ψ(3) RNA, but at 5 °C, the motion became more anisotropic for the labeled RNAs, implying more order within the RNA lower stems. As observed by 9.5 GHz EPR, the slowing of nanosecond motions of large segments of the oligonucleotide was enhanced by increasing the ratio of the nucleocapsid protein NCp7 to mini c TAR DNA from 0 to 2. The slowing was most significant at labels in the loop and near the bulge. At a 4:1 ratio of NCp7 to mini c TAR DNA, all labels reported tumbling times of >5 ns, indicating a condensation of NCp7 and TAR DNA. At the 4:1 ratio, pulse dipolar EPR spectroscopy of bilabels attached near the 3' and 5' termini showed evidence of an NCp7-induced increase in the 3'-5' end-to-end distance distribution and a partially melted stem.

An RNA Mapping DataBase for curating RNA structure mapping experiments

SUMMARY

We have established an RNA mapping database (RMDB) to enable structural, thermodynamic and kinetic comparisons across single-nucleotide-resolution RNA structure mapping experiments. The volume of structure mapping data has greatly increased since the development of high-throughput sequencing techniques, accelerated software pipelines and large-scale mutagenesis. For scientists wishing to infer relationships between RNA sequence/structure and these mapping data, there is a need for a database that is curated, tagged with error estimates and interfaced with tools for sharing, visualization, search and meta-analysis. Through its on-line front-end, the RMDB allows users to explore single-nucleotide-resolution mapping data in heat-map, bar-graph and colored secondary structure graphics; to leverage these data to generate secondary structure hypotheses; and to download the data in standardized and computer-friendly files, including the RDAT and community-consensus SNRNASM formats. At the time of writing, the database houses 53 entries, describing more than 2848 experiments of 1098 RNA constructs in several solution conditions and is growing rapidly.

AVAILABILITY

Freely available on the web at http://rmdb.stanford.edu.

CONTACT

rhiju@stanford.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics Online.

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.

Synthesis of (6-(13)C)pyrimidine nucleotides as spin-labels for RNA dynamics

We present a (13)C-based isotope labeling protocol for RNA. Using (6-(13)C)pyrimidine phosphoramidite building blocks, site-specific labels can be incorporated into a target RNA via chemical oligonucleotide solid-phase synthesis. This labeling scheme is particularly useful for studying milli- to microsecond dynamics via NMR spectroscopy, as an isolated spin system is a crucial prerequisite to apply Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion type experiments. We demonstrate the applicability for the characterization and detection of functional dynamics on various time scales by incorporating the (6-(13)C)uridine and -cytidine labels into biologically relevant RNAs. The refolding kinetics of a bistable terminator antiterminator segment involved in the gene regulation process controlled by the preQ(1) riboswitch class I was investigated. Using (13)C CPMG relaxation dispersion NMR spectroscopy, the milli- to microsecond dynamics of the HIV-1 transactivation response element RNA and the Varkud satellite stem loop V motif was addressed.

Flexibility analysis of biomacromolecules with application to computer-aided drug design

Flexibility characteristics of biomacromolecules can be efficiently determined down to the atomic level by a graph-theoretical technique as implemented in the FIRST (Floppy Inclusion and Rigid Substructure Topology) and ProFlex software packages. The method has been successfully applied to a series of protein and nucleic acid structures. Here, we describe practical guidelines for setting up and performing a flexibility analysis, discuss current bottlenecks of the approach, and provide sample applications as to how this technique can support computer-aided drug design approaches.

Mapping the landscape of RNA dynamics with NMR spectroscopy

Among the three major classes of biomacromolecules (DNA, RNA, and proteins) RNA's pronounced dynamics are the most explicitly linked to its wide variety of functions, which include catalysis and the regulation of transcription, translation, and splicing. These functions are mediated by a range of RNA biomachinery, including such varied examples as macromolecular noncoding RNAs, microRNAs, small interfering RNAs, riboswitch RNAs, and RNA thermometers. In each case, the functional dynamics of an interconversion is characterized by an associated rate constant. In this Account, we provide an introduction to NMR spectroscopic characterization of the landscape of RNA dynamics. We introduce strategies for measuring NMR parameters at various time scales as well as the underlying models for describing the corresponding rate constants. RNA exhibits significant dynamic motion, which can be modulated by (i) intermolecular interactions, including specific and nonspecific binding of ions (such as Mg(2+) and tertiary amines), (ii) metabolites in riboswitches or RNA aptamers, and (iii) macromolecular interactions within ribonucleic protein particles, including the ribosome and the spliceosome. Our understanding of the nature of these dynamic changes in RNA targets is now being incorporated into RNA-specific approaches in the design of RNA inhibitors. Interactions of RNA with proteins, other RNAs, or small molecules often occur through binding mechanisms that follow an induced fit mechanism or a conformational selection mechanism, in which one of several populated RNA conformations is selected through ligand binding. The extent of functional dynamics, including the kinetic formation of a specific RNA tertiary fold, is dependent on the messenger RNA (mRNA) chain length. Thus, during de novo synthesis of mRNA, both in prokaryotes and eukaryotes, nascent mRNA of various lengths will adopt different secondary and tertiary structures. The speed of transcription has a critical influence on the functional dynamics of the RNA being synthesized. In addition to modulating the local dynamics of a conformational RNA ensemble, a given RNA sequence may adopt more than one global, three-dimensional structure. RNA modification is one way to select among these alternative structures, which are often characterized by nearly equal stability, but with high energy barriers for conformational interconversion. The refolding of different secondary and tertiary structures has been found to be a major regulatory mechanism for transcription and translation. These conformational transitions can be characterized with NMR spectroscopy, for any given RNA sequence, in response to external stimuli.

Accommodating variety in iron-responsive elements: Crystal structure of transferrin receptor 1 B IRE bound to iron regulatory protein 1

Iron responsive elements (IREs) are short stem-loop structures found in several mRNAs encoding proteins involved in cellular iron metabolism. Iron regulatory proteins (IRPs) control iron homeostasis through differential binding to the IREs, accommodating any sequence or structural variations that the IREs may present. Here we report the structure of IRP1 in complex with transferrin receptor 1 B (TfR B) IRE, and compare it to the complex with ferritin H (Ftn H) IRE. The two IREs are bound to IRP1 through nearly identical protein-RNA contacts, although their stem conformations are significantly different. These results support the view that binding of different IREs with IRP1 depends both on protein and RNA conformational plasticity, adapting to RNA variation while retaining conserved protein-RNA contacts.

Characterizing RNA dynamics at atomic resolution using solution-state NMR spectroscopy

Many recently discovered noncoding RNAs do not fold into a single native conformation but sample many different conformations along their free-energy landscape to carry out their biological function. Here we review solution-state NMR techniques that measure the structural, kinetic and thermodynamic characteristics of RNA motions spanning picosecond to second timescales at atomic resolution, allowing unprecedented insights into the RNA dynamic structure landscape. From these studies a basic description of the RNA dynamic structure landscape is emerging, bringing new insights into how RNA structures change to carry out their function as well as applications in RNA-targeted drug discovery and RNA bioengineering.

Unraveling the structural complexity in a single-stranded RNA tail: implications for efficient ligand binding in the prequeuosine riboswitch

Single-stranded RNAs (ssRNAs) are ubiquitous RNA elements that serve diverse functional roles. Much of our understanding of ssRNA conformational behavior is limited to structures in which ssRNA directly engages in tertiary interactions or is recognized by proteins. Little is known about the structural and dynamic behavior of free ssRNAs at atomic resolution. Here, we report the collaborative application of nuclear magnetic resonance (NMR) and replica exchange molecular dynamics (REMD) simulations to characterize the 12 nt ssRNA tail derived from the prequeuosine riboswitch. NMR carbon spin relaxation data and residual dipolar coupling measurements reveal a flexible yet stacked core adopting an A-form-like conformation, with the level of order decreasing toward the terminal ends. An A-to-C mutation within the polyadenine tract alters the observed dynamics consistent with the introduction of a dynamic kink. Pre-ordering of the tail may increase the efficacy of ligand binding above that achieved by a random-coil ssRNA. The REMD simulations recapitulate important trends in the NMR data, but suggest more internal motions than inferred from the NMR analysis. Our study unmasks a previously unappreciated level of complexity in ssRNA, which we believe will also serve as an excellent model system for testing and developing computational force fields.

3D maps of RNA interhelical junctions

More than 50% of RNA secondary structure is estimated to be A-form helices, which are linked together by various junctions. Here we describe a protocol for computing three interhelical Euler angles describing the relative orientation of helices across RNA junctions. 5' and 3' helices, H1 and H2, respectively, are assigned based on the junction topology. A reference canonical helix is constructed using an appropriate molecular builder software consisting of two continuous idealized A-form helices (iH1 and iH2) with helix axis oriented along the molecular Z-direction running toward the positive direction from iH1 to iH2. The phosphate groups and the carbon and oxygen atoms of the sugars are used to superimpose helix H1 of a target interhelical junction onto the corresponding iH1 of the reference helix. A copy of iH2 is then superimposed onto the resulting H2 helix to generate iH2'. A rotation matrix R is computed, which rotates iH2' into iH2 and expresses the rotation parameters in terms of three Euler angles α(h), β(h) and γ(h). The angles are processed to resolve a twofold degeneracy and to select an overall rotation around the axis of the reference helix. The three interhelical Euler angles define clockwise rotations around the 5' (-γ(h)) and 3' (α(h)) helices and an interhelical bend angle (β(h)). The angles can be depicted graphically to provide a 'Ramachandran'-type view of RNA global structure that can be used to identify unusual conformations as well as to understand variations due to changes in sequence, junction topology and other parameters.

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.

Structurally conserved five nucleotide bulge determines the overall topology of the core domain of human telomerase RNA

Telomerase is a unique ribonucleoprotein complex that catalyzes the addition of telomeric DNA repeats onto the 3' ends of linear chromosomes. All vertebrate telomerase RNAs contain a catalytically essential core domain that includes the template and a pseudoknot with extended helical subdomains. Within these helical regions is an asymmetric 5-nt internal bulge loop (J2a/b) flanked by helices (P2a and P2b) that is highly conserved in its location but not sequence. NMR structure determination reveals that J2a/b forms a defined S-shape and creates an ∼90 ° bend with a surprisingly low twist (∼10 °) between the flanking helices. A search of RNA structures revealed only one other example of a 5-nt bulge, from hepatitis C virus internal ribosome entry site, with a different sequence but the same structure. J2a/b is intrinsically flexible but the interhelical motions across the loop are remarkably restricted. Nucleotide substitutions in J2a/b that affect the bend angle, direction, and interhelical dynamics are correlated with telomerase activity. Based on the structures of P2ab (J2a/b and flanking helices), the conserved region of the pseudoknot (P2b/P3, previously determined) and the remaining helical segment (P2a.1-J2a.1 refined using residual dipolar couplings and the modeling program MC-Sym) we have calculated an NMR-based model of the full-length pseudoknot. The model and dynamics analysis show that J2a/b serves as a dominant structural and dynamical element in defining the overall topology of the core domain, and suggest that interhelical motions in P2ab facilitate nucleotide addition along the template and template translocation.

Molecular recognition of RNA: challenges for modelling interactions and plasticity

There is growing interest in molecular recognition processes of RNA because of RNA's widespread involvement in biological processes. Computational approaches are increasingly used for analysing and predicting binding to RNA, fuelled by encouraging progress in developing simulation, free energy and docking methods for nucleic acids. These developments take into account challenges regarding the energetics of RNA-ligand binding, RNA plasticity, and the presence of water molecules and ions in the binding interface. Accordingly, we will detail advances in force field and scoring function development for molecular dynamics (MD) simulations, free energy computations and docking calculations of nucleic acid complexes. Furthermore, we present methods that can detect moving parts within RNA structures based on graph-theoretical approaches or normal mode analysis (NMA). As an example of the successful use of these developments, we will discuss recent structure-based drug design approaches that focus on the bacterial ribosomal A-site RNA as a drug target.