Resolution-optimized NMR measurement of (1)D(CH), (1)D(CC) and (2)D(CH) residual dipolar couplings in nucleic acid bases

Unequivocal structure confirmation of a breitfussin analog by anisotropic NMR measurements

Structural features of proton-deficient heteroaromatic natural products, such as the breitfussins, can severely complicate their characterization by NMR spectroscopy. For the breitfussins in particular, the constitution of the five-membered oxazole central ring cannot be unequivocally established via conventional NMR methods when the 4′-position is halogenated. The level of difficulty is exacerbated by 4′-iodination, as the accuracy with which theoretical NMR parameters are determined relies extensively on computational treatment of the relativistic effects of the iodine atom. It is demonstrated in the present study, that the structure of a 4′-iodo breitfussin analog can be unequivocally established by anisotropic NMR methods, by adopting a reduced singular value decomposition (SVD) protocol that leverages the planar structures exhibited by its conformers.

Structural features of proton-deficient heteroaromatic natural products, such as the breitfussins, can severely complicate their characterization by NMR spectroscopy.

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.

Restricted puckering of mineralized RNA-like riboses

The pseudorotational motions of highly hydroxylated pentafuranose sugars in the free state and tethered to hydroxyapatite have been compared. The conformation pentafuranose ring remains restricted at the North region of the pseudorotational wheel, which is the one typically observed for nucleosides and nucleotides in the double helix A-RNA, when the phosphate-bearing sugar is anchored to the mineral surface. Results indicate that the severe restrictions imposed by the mineral are responsible of the double helix preservation when DNA and RNA are encapsulated in crystalline nanorods.

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.

Nuclear magnetic resonance analysis of protein-DNA interactions

Recent methodological and instrumental advances in solution-state nuclear magnetic resonance have opened up the way to investigating challenging problems in structural biology such as large macromolecular complexes. This review focuses on the experimental strategies currently employed to solve structures of protein-DNA complexes and to analyse their dynamics. It highlights how these approaches can help in understanding detailed molecular mechanisms of target recognition.

NMR structure of the A730 loop of the Neurospora VS ribozyme: insights into the formation of the active site

The Neurospora VS ribozyme is a small nucleolytic ribozyme with unique primary, secondary and global tertiary structures, which displays mechanistic similarities to the hairpin ribozyme. Here, we determined the high-resolution NMR structure of a stem–loop VI fragment containing the A730 internal loop, which forms part of the active site. In the presence of magnesium ions, the A730 loop adopts a structure that is consistent with existing biochemical data and most likely reflects its conformation in the VS ribozyme prior to docking with the cleavage site internal loop. Interestingly, the A730 loop adopts an S-turn motif that is also present in loop B within the hairpin ribozyme active site. The S-turn appears necessary to expose the Watson–Crick edge of a catalytically important residue (A756) so that it can fulfill its role in catalysis. The A730 loop and the cleavage site loop of the VS ribozyme display structural similarities to internal loops found in the active site of the hairpin ribozyme. These similarities provided a rationale to build a model of the VS ribozyme active site based on the crystal structure of the hairpin ribozyme.

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.

Liquid-crystal NMR structure of HIV TAR RNA bound to its SELEX RNA aptamer reveals the origins of the high stability of the complex

Transactivation-response element (TAR) is a stable stem-loop structure of HIV RNA, which plays a crucial role during the life cycle of the virus. The apical loop of TAR acts as a binding site for essential cellular cofactors required for the replication of HIV. High-affinity aptamers directed against the apical loop of TAR have been identified by the SELEX approach. The RNA aptamers with the highest affinity for TAR fold as hairpins and form kissing complexes with the targeted RNA through loop-loop interactions. The aptamers with the strongest binding properties all possess a GA base pair combination at the loop-closing position. Using liquid-crystal NMR methodology, we have obtained a structural model in solution of a TAR-aptamer kissing complex with an unprecedented accuracy. This high-resolution structure reveals that the GA base pair is unilaterally shifted toward the 5' strand and is stabilized by a network of intersugar hydrogen bonds. This specific conformation of the GA base pair allows for the formation of two supplementary stable base-pair interactions. By systematic permutations of the loop-closing base pair, we establish that the identified atomic interactions, which form the basis for the high stability of the complex, are maintained in several other kissing complexes. This study rationalizes the stabilizing role of the loop-closing GA base pairs in kissing complexes and may help the development or improvement of drugs against RNA loops of viruses or pathogens as well as the conception of biochemical tools targeting RNA hairpins involved in important biological functions.

13C-detection in RNA bases: revealing structure-chemical shift relationships

The chemical shifts of the unprotonated carbons in the proton-deficient nucleobases of RNA are rarely reported, despite the valuable information that they contain about base-pairing and base-stacking. We have developed 13C-detected 2D-experiments to identify the unprotonated 13C in the RNA bases and have assigned all the base nuclei of uniformly 13C,15N-labeled HIV-2 TAR-RNA. The 13C chemical shift distributions revealed perturbations correlated with the base-pairing and base-stacking properties of all four base-types. From this work, we conclude that the information contained in the chemical shift perturbations within the base rings can provide valuable restraint information for solving RNA structures, especially in conformational averaged regions, where NOE-based information is not available.

13C-detected NMR experiments for measuring chemical shifts and coupling constants in nucleic acid bases

The paper presents a set of two-dimensional experiments that utilize direct (13)C detection to provide proton-carbon, carbon-carbon and carbon-nitrogen correlations in the bases of nucleic acids. The set includes a (13)C-detected proton-carbon correlation experiment for the measurement of (13)C-(13)C couplings, the CaCb experiment for correlating two quaternary carbons, the HCaCb experiment for the (13)C-(13)C correlations in cases where one of the carbons has a proton attached, the HCC-TOCSY experiment for correlating a proton with a network of coupled carbons, and a (13)C-detected (13)C-(15)N correlation experiment for detecting the nitrogen nuclei that cannot be detected via protons. The IPAP procedure is used for extracting the carbon-carbon couplings and/or carbon decoupling in the direct dimension, while the S(3)E procedure is preferred in the indirect dimension of the carbon-nitrogen experiment to obtain the value of the coupling constant. The experiments supply accurate values of (13)C and (15)N chemical shifts and carbon-carbon and carbon-nitrogen coupling constants. These values can help to reveal structural features of nucleic acids either directly or via induced changes when the sample is dissolved in oriented media.

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

An increasing number of RNAs are being discovered that perform their functions by undergoing large changes in conformation in response to a variety of cellular signals, including recognition of proteins and small molecular targets, changes in temperature, and RNA synthesis itself. The measurement of NMR residual dipolar couplings (RDCs) in partially aligned systems is providing new insights into the structural plasticity of RNA through combined characterization of large-amplitude collective helix motions and local flexibility in noncanonical regions over a wide window of biologically relevant timescales (<milliseconds). Here, we review RDC methodology for studying RNA structural dynamics and survey what has been learnt thus far from application of these methods. Future methodological challenges are also identified.

Characterizing the relative orientation and dynamics of RNA A-form helices using NMR residual dipolar couplings

We present a protocol for determining the relative orientation and dynamics of A-form helices in 13C/15N isotopically enriched RNA samples using NMR residual dipolar couplings (RDCs). Non-terminal Watson-Crick base pairs in helical stems are experimentally identified using NOE and trans-hydrogen bond connectivity and modeled using the idealized A-form helix geometry. RDCs measured in the partially aligned RNA are used to compute order tensors describing average alignment of each helix relative to the applied magnetic field. The order tensors are translated into Euler angles defining the average relative orientation of helices and order parameters describing the amplitude and asymmetry of interhelix motions. The protocol does not require complete resonance assignments and therefore can be implemented rapidly to RNAs much larger than those for which complete high-resolution NMR structure determination is feasible. The protocol is particularly valuable for exploring adaptive changes in RNA conformation that occur in response to biologically relevant signals. Following resonance assignments, the procedure is expected to take no more than 2 weeks of acquisition and data analysis time.

Chemical shift tensors of protonated base carbons in helical RNA and DNA from NMR relaxation and liquid crystal measurements

Knowledge of (13)C chemical shift anisotropy (CSA) tensors in nucleotide bases is important for interpretation of NMR relaxation data in terms of local dynamic properties of nucleic acids and for analysis of residual chemical shift anisotropy (RCSA) resulting from weak alignment. CSA tensors for protonated nucleic acid base carbons have been derived from measurements on a uniformly (13)C-enriched helical A-form RNA segment and a helical B-form DNA dodecamer at natural (13)C abundance. The magnitudes of the derived CSA principal values are tightly restricted by the magnetic field dependencies of the (13)C transverse relaxation rates, whereas the tensor orientation and asymmetry follow from quantitative measurements of interference between (13)C-{(1)H} dipolar and (13)C CSA relaxation mechanisms. Changes in the chemical shift between the isotropic and aligned states, Deltadelta, complement these measurements and permit cross-validation. The CSA tensors are determined from the experimental Deltadelta values and relaxation rates, under the assumption that the CSA tensor of any specific carbon in a given type of base is independent of the base position in either the RNA or DNA helix. However, the experimental data indicate that for pyrimidine C(6) carbons in A-form RNA the CSA magnitude is considerably larger than in B-form DNA. This result is supported by quantum chemical calculations and is attributed in part to the close proximity between intranucleotide C(6)H and O(5)' atoms in RNA. The magnitudes of the measured CSA tensors, on average, agree better with previous solid-state NMR results obtained on powdered nucleosides than with prior results from quantum chemical calculations on isolated bases, which depend rather strongly on the level of theory at which the calculations are carried out. In contrast, previously computed orientations of the chemical shift tensors agree well with the present experimental results and exhibit less dependence on the level of theory at which the computations are performed.

Structure and thermodynamics of a conserved U2 snRNA domain from yeast and human

The spliceosome is a dynamic ribonucleoprotein complex responsible for the removal of intron sequences from pre-messenger RNA. The highly conserved 5' end of the U2 small nuclear RNA (snRNA) makes key base-pairing interactions with the intron branch point sequence and U6 snRNA. U2 stem I, a stem-loop located in the 5' region of U2, has been implicated in spliceosome assembly and may modulate the folding of the U2 and U6 snRNAs in the spliceosome active site. Here we present the NMR structures of U2 stem I from human and Saccharomyces cerevisiae. These sequences represent the two major classes of U2 stem I, distinguished by the identity of tandem wobble pairs (UU/UU in yeast and CA/GU in human) and the presence of post-transcriptional modifications (four 2'-O-methyl groups and two pseudouracils in human). The structures reveal that the UU/UU and CA/GU tandem wobble pairs are nearly isosteric. The tandem wobble pairs separate two thermodynamically distinct regions of Watson-Crick base pairs, with the modified nucleotides in human stem I conferring a significant increase in stability. We hypothesize that the separate thermodynamic stabilities of U2 stem I exist to allow the structure to transition through different folded conformations during spliceosome assembly and catalysis.

Impact of static and dynamic A-form heterogeneity on the determination of RNA global structural dynamics using NMR residual dipolar couplings

We examined how static and dynamic deviations from the idealized A-form helix propagate into errors in the principal order tensor parameters determined using residual dipolar couplings (rdcs). A 20-ns molecular dynamics (MD) simulation of the HIV-1 transactivation response element (TAR) RNA together with a survey of spin relaxation studies of RNA dynamics reveals that pico-to-nanosecond local motions in non-terminal Watson-Crick base-pairs will uniformly attenuate base and sugar one bond rdcs by approximately 7%. Gaussian distributions were generated for base and sugar torsion angles through statistical comparison of 40 RNA X-ray structures solved to <3.0 A resolution. For a typical number (>or=11) of one bond C-H base and sugar rdcs, these structural deviations together with rdc uncertainty (1.5 Hz) lead to average errors in the magnitude and orientation of the principal axis of order that are <9% and <4 degrees, respectively. The errors decrease to <5% and <4 degrees for >or=17 rdcs. A protocol that allows for estimation of error in A-form order tensors due to both angular deviations and rdc uncertainty (Aform-RDC) is validated using theoretical simulations and used to analyze rdcs measured previously in TAR in the free state and bound to four distinct ligands. Results confirm earlier findings that the two TAR helices undergo large changes in both their mean relative orientation and dynamics upon binding to different targets.

Weak alignment NMR: a hawk-eyed view of biomolecular structure

Imposing a very slight deviation from the isotropic random distribution of macromolecules in solution in an NMR sample tube permits the measurement of residual internuclear dipolar couplings (RDCs). Such interactions are very sensitive functions of the time-averaged orientation of the corresponding internuclear vectors and thereby offer highly precise structural information. In recent years, advances have been made both in the technology to measure RDCs and in the computational procedures that integrate this information in the structure determination process. The exceptional precision with which RDCs can be measured under weakly aligned conditions is also starting to reveal the mostly, but not universally, subtle effects of internal protein dynamics. Importantly, RDCs potentially can reveal motions taking place on a timescale slower than rotational diffusion and analysis is uniquely sensitive to the direction of motion, not just its amplitude.

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.

Accurate measurement of 15N-13C residual dipolar couplings in nucleic acids

New 3D HCN quantitative J (QJ) pulse schemes are presented for the precise and accurate measurement of one-bond 15N1/9-13C1', 15N1/9-13C6/8, and 15N1/9-13C2/4 residual dipolar couplings (RDCs) in weakly aligned nucleic acids. The methods employ 1H-13C multiple quantum (MQ) coherence or TROSY-type pulse sequences for optimal resolution and sensitivity. RDCs are obtained from the intensity ratio of H1'-C1'-N1/9 (MQ-HCN-QJ) or H6/8-C6/8-N1/9 (TROSY-HCN-QJ) correlations in two interleaved 3D NMR spectra, with dephasing intervals of zero (reference spectrum) and approximately 1/(2J(NC)) (attenuated spectrum). The different types of 15N-13C couplings can be obtained by using either the 3D MQ-HCN-QJ or TROSY-HCN-QJ pulse scheme, with the appropriate setting of the duration of the constant-time 15N evolution period and the offset of two frequency-selective 13C pulses. The methods are demonstrated for a uniformly 13C, 15N-enriched 24-nucleotide stem-loop RNA sequence, helix-35psi, aligned in the magnetic field using phage Pf1. For measurements of RDCs systematic errors are found to be negligible, and experiments performed on a 1.5 mM helix-35psi sample result in an estimated precision of ca. 0.07 Hz for 1D(NC), indicating the utility of the measured RDCs in structure validation and refinement. Indeed, for a complete set of 15N1/9-13C1', 15N1/9-13C6/8, and 15N1/9-13C2/4 dipolar couplings obtained for the stem nucleotides, the measured RDCs are in excellent agreement with those predicted for an NMR structure of helix-35psi, refined using independently measured observables, including 13C-1H, 13C-13C and 1H-1H dipolar couplings.