Insight into the CSA tensors of nucleobase carbons in RNA polynucleotides from solution measurements of residual CSA: towards new long-range orientational constraints

Probing Molecular Packing of Lipid Nanoparticles from 31P Solution and Solid-State NMR

Lipid nanoparticles (LNPs) are intricate multicomponent systems widely recognized for their efficient delivery of oligonucleotide cargo to host cells. Gaining insights into the molecular properties of LNPs is crucial for their effective design and characterization. However, analysis of their internal structure at the molecular level presents a significant challenge. This study introduces 31P nuclear magnetic resonance (NMR) methods to acquire structural and dynamic information about the phospholipid envelope of LNPs. Specifically, we demonstrate that the 31P chemical shift anisotropy (CSA) parameters serve as a sensitive indicator of the molecular assembly of distearoylphosphatidylcholine (DSPC) lipids within the particles. An analytical protocol for measuring 31P CSA is developed, which can be implemented using either solution NMR or solid-state NMR, offering wide accessibility and adaptability. The capability of this method is demonstrated using both model DSPC liposomes and real-world pharmaceutical LNP formulations. Furthermore, our method can be employed to investigate the impact of formulation processes and composition on the assembly of specifically LNP particles or, more generally, phospholipid-based delivery systems. This makes it an indispensable tool for evaluating critical pharmaceutical properties such as structural homogeneity, batch-to-batch reproducibility, and the stability of the particles.

Kinetic Resolution of the Atomic 3D Structures Formed by Ground and Excited Conformational States in an RNA Dynamic Ensemble

Knowing the 3D structures formed by the various conformations populating the RNA free-energy landscape, their relative abundance, and kinetic interconversion rates is required to obtain a quantitative and predictive understanding of how RNAs fold and function at the atomic level. While methods integrating ensemble-averaged experimental data with computational modeling are helping define the most abundant conformations in RNA ensembles, elucidating their kinetic rates of interconversion and determining the 3D structures of sparsely populated short-lived RNA excited conformational states (ESs) remains challenging. Here, we developed an approach integrating Rosetta-FARFAR RNA structure prediction with NMR residual dipolar couplings and relaxation dispersion that simultaneously determines the 3D structures formed by the ground-state (GS) and ES subensembles, their relative abundance, and kinetic rates of interconversion. The approach is demonstrated on HIV-1 TAR, whose six-nucleotide apical loop was previously shown to form a sparsely populated (∼13%) short-lived (lifetime ∼ 45 μs) ES. In the GS, the apical loop forms a broad distribution of open conformations interconverting on the pico-to-nanosecond time scale. Most residues are unpaired and preorganized to bind the Tat-superelongation protein complex. The apical loop zips up in the ES, forming a narrow distribution of closed conformations, which sequester critical residues required for protein recognition. Our work introduces an approach for determining the 3D ensemble models formed by sparsely populated RNA conformational states, provides a rare atomic view of an RNA ES, and kinetically resolves the atomic 3D structures of RNA conformational substates, interchanging on time scales spanning 6 orders of magnitude, from picoseconds to microseconds.

DNA base order parameter determination without influence of chemical exchange

Nuclear magnetic resonance (NMR) spectroscopy is a versatile tool used to investigate the dynamic properties of biological macromolecules and their complexes. NMR relaxation data can provide order parameters S2, which represent the mobility of bond vectors reorienting within a molecular frame. Determination of S2 parameters typically involves the use of transverse NMR relaxation rates. However, the accuracy in S2 determination can be diminished by elevation of the transverse relaxation rates through conformational or chemical exchange involving protonation/deprotonation or non-Watson-Crick base-pair states of nucleic acids. Here, we propose an approach for determination of S2 parameters without the influence of exchange processes. This approach utilizes transverse and longitudinal 13C chemical shift anisotropy (CSA) - dipole-dipole (DD) cross-correlation rates instead of 13C transverse relaxation rates. Anisotropy in rotational diffusion is taken into consideration. An application of this approach to nucleotide base CH groups of a uniformly 13C/15N-labeled DNA duplex is demonstrated.

Polymer-Nanodiscs as a Novel Alignment Medium for High-Resolution NMR-Based Structural Studies of Nucleic Acids

Residual dipolar couplings (RDCs) are increasingly used for high-throughput NMR-based structural studies and to provide long-range angular constraints to validate and refine structures of various molecules determined by X-ray crystallography and NMR spectroscopy. RDCs of a given molecule can be measured in an anisotropic environment that aligns in an external magnetic field. Here, we demonstrate the first application of polymer-based nanodiscs for the measurement of RDCs from nucleic acids. Polymer-based nanodiscs prepared using negatively charged SMA-EA polymer and zwitterionic DMPC lipids were characterized by size-exclusion chromatography, 1H NMR, dynamic light-scattering, and 2H NMR. The magnetically aligned polymer-nanodiscs were used as an alignment medium to measure RDCs from a 13C/15N-labeled fluoride riboswitch aptamer using 2D ARTSY-HSQC NMR experiments. The results showed that the alignment of nanodiscs is stable for nucleic acids and nanodisc-induced RDCs fit well with the previously determined solution structure of the riboswitch. These results demonstrate that SMA-EA-based lipid-nanodiscs can be used as a stable alignment medium for high-resolution structural and dynamical studies of nucleic acids, and they can also be applicable to study various other biomolecules and small molecules in general.

Advances that facilitate the study of large RNA structure and dynamics by nuclear magnetic resonance spectroscopy

The characterization of functional yet nonprotein coding (nc) RNAs has expanded the role of RNA in the cell from a passive player in the central dogma of molecular biology to an active regulator of gene expression. The misregulation of ncRNA function has been linked with a variety of diseases and disorders ranging from cancers to neurodegeneration. However, a detailed molecular understanding of how ncRNAs function has been limited; due, in part, to the difficulties associated with obtaining high-resolution structures of large RNAs. Tertiary structure determination of RNA as a whole is hampered by various technical challenges, all of which are exacerbated as the size of the RNA increases. Namely, RNAs tend to be highly flexible and dynamic molecules, which are difficult to crystallize. Biomolecular nuclear magnetic resonance (NMR) spectroscopy offers a viable alternative to determining the structure of large RNA molecules that do not readily crystallize, but is itself hindered by some technical limitations. Recently, a series of advancements have allowed the biomolecular NMR field to overcome, at least in part, some of these limitations. These advances include improvements in sample preparation strategies as well as methodological improvements. Together, these innovations pave the way for the study of ever larger RNA molecules that have important biological function. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

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.

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

High-performance virtual screening by targeting a high-resolution RNA dynamic ensemble

Dynamic ensembles hold great promise in advancing RNA-targeted drug discovery. Here we subjected the transactivation response element (TAR) RNA from human immunodeficiency virus type-1 to experimental high-throughput screening against ~100,000 drug-like small molecules. Results were augmented with 170 known TAR-binding molecules and used to generate sublibraries optimized for evaluating enrichment when virtually screening a dynamic ensemble of TAR determined by combining NMR spectroscopy data and molecular dynamics simulations. Ensemble-based virtual screening scores molecules with an area under the receiver operator characteristic curve of ~0.85-0.94 and with ~40-75% of all hits falling within the top 2% of scored molecules. The enrichment decreased significantly for ensembles generated from the same molecular dynamics simulations without input NMR data and for other control ensembles. The results demonstrate that experimentally determined RNA ensembles can significantly enrich libraries with true hits and that the degree of enrichment is dependent on the accuracy of the ensemble.

Combining asymmetric 13C-labeling and isotopic filter/edit NOESY: a novel strategy for rapid and logical RNA resonance assignment

Although ∼98% of the human genomic output is transcribed as non-protein coding RNA, <2% of the protein data bank structures comprise RNA. This huge structural disparity stems from combined difficulties of crystallizing RNA for X-ray crystallography along with extensive chemical shift overlap and broadened linewidths associated with NMR of RNA. While half of the deposited RNA structures in the PDB were solved by NMR methods, the usefulness of NMR is still limited by the high cost of sample preparation and challenges of resonance assignment. Here we propose a novel strategy for resonance assignment that combines new strategic 13C labeling technologies with filter/edit type NOESY experiments to greatly reduce spectral complexity and crowding. This new strategy allowed us to assign important non-exchangeable resonances of proton and carbon (1', 2', 2, 5, 6 and 8) nuclei using only one sample and <24 h of NMR instrument time for a 27 nt model RNA. The method was further extended to assigning a 6 nt bulge from a 61 nt viral RNA element justifying its use for a wide range RNA chemical shift resonance assignment problems.

Structural conservation in the template/pseudoknot domain of vertebrate telomerase RNA from teleost fish to human

Telomerase is an RNA-protein complex that includes a unique reverse transcriptase that catalyzes the addition of single-stranded telomere DNA repeats onto the 3' ends of linear chromosomes using an integral telomerase RNA (TR) template. Vertebrate TR contains the template/pseudoknot (t/PK) and CR4/5 domains required for telomerase activity in vitro. All vertebrate pseudoknots include two subdomains: P2ab (helices P2a and P2b with a 5/6-nt internal loop) and the minimal pseudoknot (P2b-P3 and associated loops). A helical extension of P2a, P2a.1, is specific to mammalian TR. Using NMR, we investigated the structures of the full-length TR pseudoknot and isolated subdomains in Oryzias latipes (Japanese medaka fish), which has the smallest vertebrate TR identified to date. We determined the solution NMR structure and studied the dynamics of medaka P2ab, and identified all base pairs and tertiary interactions in the minimal pseudoknot. Despite differences in length and sequence, the structure of medaka P2ab is more similar to human P2ab than predicted, and the medaka minimal pseudoknot has the same tertiary interactions as the human pseudoknot. Significantly, although P2a.1 is not predicted to form in teleost fish, we find that it forms in the full-length pseudoknot via an unexpected hairpin. Model structures of the subdomains are combined to generate a model of t/PK. These results provide evidence that the architecture for the vertebrate t/PK is conserved from teleost fish to human. The organization of the t/PK on telomerase reverse transcriptase for medaka and human is modeled based on the cryoEM structure of Tetrahymena telomerase, providing insight into function.

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

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

Rapid Exchange Between Free and Bound States in RNA-Dendrimer Polyplexes: Implications on the Mechanism of Delivery and Release

A combination of solution NMR, dynamic light scattering (DLS), and fluorescence quenching assays were employed to obtain insights into the dynamics and structural features of a polyplex system consisting of HIV-1 transactivation response element (TAR) and PEGylated generation 5 poly(amidoamine) dendrimer (G5-PEG). NMR chemical shift mapping and (13)C spin relaxation based dynamics measurements depict the polyplex system as a highly dynamic assembly where the RNA, with its local structure and dynamics preserved, rapidly exchanges ( Collapse

Key Words

• therapeutic rna delivery
• pamam
• pegylation
• tar
• dls
• nmr spin relaxation

Collapse

MESH Headings

• Dendrimers/chemistry
• HIV-1/chemistry
• Magnetic Resonance Spectroscopy/methods
• Polyamines/chemistry
• Polyethylene Glycols/chemistry
• RNA/chemistry
• Response Elements/genetics
• Transfection/methods

Collapse

Grants

• R01 AI066975 NIAID NIH HHS
• P50 GM103297 NIGMS NIH HHS
• R01 AI087463-01 NIAID NIH HHS
• R01 EB005028 NIBIB NIH HHS
• R01EB005028 NIBIB NIH HHS
• R01 GM089846 NIGMS NIH HHS

Collapse

Affiliation(s)

Anisha Shakya Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Casey A. Dougherty Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA
Yi Xue Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Hashim M. Al-Hashimi Department of Biochemistry and Chemistry, Duke University Medical Center, Durham, NC 27710, USA
Mark M. Banaszak Holl Department of Chemistry, University of Michigan, Ann Arbor, MI 48109-1055, USA

Collapse

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.

The structure of the SOLE element of oskar mRNA

mRNA localization by active transport is a regulated process that requires association of mRNPs with protein motors for transport along either the microtubule or the actin cytoskeleton. oskar mRNA localization at the posterior pole of the Drosophila oocyte requires a specific mRNA sequence, termed the SOLE, which comprises nucleotides of both exon 1 and exon 2 and is assembled upon splicing. The SOLE folds into a stem-loop structure. Both SOLE RNA and the exon junction complex (EJC) are required for oskar mRNA transport along the microtubules by kinesin. The SOLE RNA likely constitutes a recognition element for a yet unknown protein, which either belongs to the EJC or functions as a bridge between the EJC and the mRNA. Here, we determine the solution structure of the SOLE RNA by Nuclear Magnetic Resonance spectroscopy. We show that the SOLE forms a continuous helical structure, including a few noncanonical base pairs, capped by a pentanucleotide loop. The helix displays a widened major groove, which could accommodate a protein partner. In addition, the apical helical segment undergoes complex dynamics, with potential functional significance.

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.

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.

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.

Structure and sequence elements of the CR4/5 domain of medaka telomerase RNA important for telomerase function

Telomerase is a unique reverse transcriptase that maintains the 3' ends of eukaryotic chromosomes by adding tandem telomeric repeats. The RNA subunit (TR) of vertebrate telomerase provides a template for reverse transcription, contained within the conserved template/pseudoknot domain, and a conserved regions 4 and 5 (CR4/5) domain, all essential for catalytic activity. We report the nuclear magnetic resonance (NMR) solution structure of the full-length CR4/5 domain from the teleost fish medaka (Oryzias latipes). Three helices emanate from a structured internal loop, forming a Y-shaped structure, where helix P6 stacks on P5 and helix P6.1 points away from P6. The relative orientations of the three helices are Mg2+ dependent and dynamic. Although the three-way junction is structured and has unexpected base pairs, telomerase activity assays with nucleotide substitutions and deletions in CR4/5 indicate that none of these are essential for activity. The results suggest that the junction is likely to change conformation in complex with telomerase reverse transcriptase and that it provides a flexible scaffold that allows P6 and P6.1 to correctly fold and interact with telomerase reverse transcriptase.

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

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

Molecular mechanism of GTPase activation at the signal recognition particle (SRP) RNA distal end

The signal recognition particle (SRP) RNA is a universally conserved and essential component of the SRP that mediates the co-translational targeting of proteins to the correct cellular membrane. During the targeting reaction, two functional ends in the SRP RNA mediate distinct functions. Whereas the RNA tetraloop facilitates initial assembly of two GTPases between the SRP and SRP receptor, this GTPase complex subsequently relocalizes ∼100 Å to the 5',3'-distal end of the RNA, a conformation crucial for GTPase activation and cargo handover. Here we combined biochemical, single molecule, and NMR studies to investigate the molecular mechanism of this large scale conformational change. We show that two independent sites contribute to the interaction of the GTPase complex with the SRP RNA distal end. Loop E plays a crucial role in the precise positioning of the GTPase complex on these two sites by inducing a defined bend in the RNA helix and thus generating a preorganized recognition surface. GTPase docking can be uncoupled from its subsequent activation, which is mediated by conserved bases in the next internal loop. These results, combined with recent structural work, elucidate how the SRP RNA induces GTPase relocalization and activation at the end of the protein targeting reaction.

Parameterization of solvent-protein interaction and its use on NMR protein structure determination

NMR structure determination is frequently hindered by an insufficient amount of distance information for determining the correct fold of the protein in its early stages. In response we introduce a simple and general structure-based metric that can be used to incorporate NMR-based restraints on protein surface accessibility. This metric is inversely proportional to the sum of the inverse square distances to neighboring heavy atoms. We demonstrate the use of this restraint using a dataset from the water to protein magnetization transfer experiment on the protein Bax and the solvent paramagnetic relaxation enhancement experiment on the protein ubiquitin and Qua1 homodimer. The calculated solvent accessibility values using the new empirical function are well correlated with the experimental data. By incorporating an associated energy term into Xplor-NIH, we show that structure calculation with a limited number of additional experimental restraints, improves both the precision and accuracy of the resulting structures. This new empirical energy term will have general applicability to other types of solvent accessibility data.

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.

Structural NMR of protein oligomers using hybrid methods

Solving structures of native oligomeric protein complexes using traditional high-resolution NMR techniques remains challenging. However, increased utilization of computational platforms, and integration of information from less traditional NMR techniques with data from other complementary biophysical methods, promises to extend the boundary of NMR-applicable targets. This article reviews several of the techniques capable of providing less traditional and complementary structural information. In particular, the use of orientational constraints coming from residual dipolar couplings and residual chemical shift anisotropy offsets are shown to simplify the construction of models for oligomeric complexes, especially in cases of weak homo-dimers. Combining this orientational information with interaction site information supplied by computation, chemical shift perturbation, paramagnetic surface perturbation, cross-saturation and mass spectrometry allows high resolution models of the complexes to be constructed with relative ease. Non-NMR techniques, such as mass spectrometry, EPR and small angle X-ray scattering, are also expected to play increasingly important roles by offering alternative methods of probing the overall shape of the complex. Computational platforms capable of integrating information from multiple sources in the modeling process are also discussed in the article. And finally a new, detailed example on the determination of a chemokine tetramer structure will be used to illustrate how a non-traditional approach to oligomeric structure determination works in practice.

A device for the measurement of residual chemical shift anisotropy and residual dipolar coupling in soluble and membrane-associated proteins

Residual dipolar coupling (RDC) and residual chemical shift anisotropy (RCSA) report on orientational properties of a dipolar bond vector and a chemical shift anisotropy principal axis system, respectively. They can be highly complementary in the analysis of backbone structure and dynamics in proteins as RCSAs generally include a report on vectors out of a peptide plane while RDCs usually report on in-plane vectors. Both RDC and RCSA average to zero in isotropic solutions and require partial orientation in a magnetic field to become observable. While the alignment and measurement of RDC has become routine, that of RCSA is less common. This is partly due to difficulties in providing a suitable isotopic reference spectrum for the measurement of the small chemical shift offsets coming from RCSA. Here we introduce a device (modified NMR tube) specifically designed for accurate measurement of reference and aligned spectra for RCSA measurements, but with a capacity for RDC measurements as well. Applications to both soluble and membrane anchored proteins are illustrated.

Major groove width variations in RNA structures determined by NMR and impact of 13C residual chemical shift anisotropy and 1H-13C residual dipolar coupling on refinement

Ribonucleic acid structure determination by NMR spectroscopy relies primarily on local structural restraints provided by (1)H- (1)H NOEs and J-couplings. When employed loosely, these restraints are broadly compatible with A- and B-like helical geometries and give rise to calculated structures that are highly sensitive to the force fields employed during refinement. A survey of recently reported NMR structures reveals significant variations in helical parameters, particularly the major groove width. Although helical parameters observed in high-resolution X-ray crystal structures of isolated A-form RNA helices are sensitive to crystal packing effects, variations among the published X-ray structures are significantly smaller than those observed in NMR structures. Here we show that restraints derived from aromatic (1)H- (13)C residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs) can overcome NMR restraint and force field deficiencies and afford structures with helical properties similar to those observed in high-resolution X-ray structures.

Site-specific backbone amide (15)N chemical shift anisotropy tensors in a small protein from liquid crystal and cross-correlated relaxation measurements

Site-specific (15)N chemical shift anisotropy (CSA) tensors have been derived for the well-ordered backbone amide (15)N nuclei in the B3 domain of protein G (GB3) from residual chemical shift anisotropy (RCSA) measured in six different mutants that retain the native structure but align differently relative to the static magnetic field when dissolved in a liquid crystalline Pf1 suspension. This information is complemented by measurement of cross-correlated relaxation rates between the (15)N CSA tensor and either the (15)N-(1)H or (15)N-(13)C' dipolar interaction. In agreement with recent solid state NMR measurements, the (15)N CSA tensors exhibit only a moderate degree of variation from averaged values, but have larger magnitudes in alpha-helical (-173 +/- 7 ppm) than in beta-sheet (-162 +/- 6 ppm) residues, a finding also confirmed by quantum computations. The orientations of the least shielded tensor component cluster tightly around an in-peptide-plane vector that makes an angle of 19.6 +/- 2.5 degrees with the N-H bond, with the asymmetry of the (15)N CSA tensor being slightly smaller in alpha-helix (eta = 0.23 +/- 0.17) than in beta-sheet (eta = 0.31 +/- 0.11). The residue-specific (15)N CSA values are validated by improved agreement between computed and experimental (15)N R(1rho) relaxation rates measured for (15)N-{(2)H} sites in GB3, which are dominated by the CSA mechanism. Use of residue-specific (15)N CSA values also results in more uniform generalized order parameters, S(2), and predicts considerable residue-by-residue variations in the magnetic field strengths where TROSY line narrowing is most effective.

Measurement and interpretation of 15N-1H residual dipolar couplings in larger proteins

A decade ago, Dr. L.E. Kay and co-workers described an ingenious HNCO-based triple-resonance experiment from which several protein backbone RDCs can be measured simultaneously (Yang et al. (1999) [1]). They implemented a J-scaling technique in the (15)N dimension of the 3D experiment to obtain the NH RDCs. We have used this idea to carry out J-scaling in a 2D (15)N-(1)H-TROSY experiment and have found it to be an excellent method to obtain NH RDCs for larger proteins upto 70 kDa, far superior to commonly used HSQC in-phase/anti-phase and HSQC/TROSY comparisons. Here, this method, dubbed "RDC-TROSY" is discussed in detail and the limits of its utility are assessed by simulations. Prominent in the latter analysis is the evaluation of the effect of amide proton flips on the "RDC-TROSY" linewidths. The details of the technical and computational implementations of these methods for the determination of domain orientations in 45-60 kDa Hsp70 chaperone protein constructs are described.

Slight mistuning of a cryogenic probe significantly perturbs the water 1H precession frequency

A shift of the water proton precession frequency is described that can introduce errors in chemical shifts derived using the water signal as the chemical shift reference. This shift, f(s), arises as a consequence of radiation damping when the water proton and detector circuit resonance frequencies differ. Herein it is shown that experimental values of f(s), measured as a function of detector circuit tuning offset for 500 and 900 MHz cryogenic probes, are in good agreement with theory. Of importance is the fact that even a small degree of mistuning, which does not significantly impact the performance of a pulse sequence, introduces chemical shift errors of +/-0.03 ppm, that negatively impact many types of experiments. A simple remedy that attenuates the frequency shift is presented.

Chemical shift anisotropy of imino 15N nuclei in Watson-Crick base pairs from magic angle spinning liquid crystal NMR and nuclear spin relaxation

Knowledge of (15)N chemical shift anisotropy is prerequisite both for quantitative interpretation of nuclear spin relaxation rates in terms of local dynamics and for the use of residual chemical shift anisotropy (RCSA) as a constraint in structure determination. Accurate measurement of the very small RCSA from the difference in (15)N chemical shift under isotropic and weakly aligning liquid crystalline conditions is very sensitive to minute differences in sample conditions, such as pH or ionic strength. For this reason, chemical shifts were measured for the same solution, under static liquid crystalline alignment, and under magic angle spinning conditions where alignment relative to the magnetic field is removed. Measurements were made for 14 well-resolved G-N(1) and 6 U-N(3) (15)N nuclei in a sample of tRNA(Val). Fitting these RCSA data together with (15)N-(1)H dipole-CSA cross-correlated relaxation measurements to the recently refined structural model of tRNA(Val) yields the magnitude, asymmetry, and orientation of the (15)N CSA tensors.

RNA phosphodiester backbone dynamics of a perdeuterated cUUCGg tetraloop RNA from phosphorus-31 NMR relaxation analysis

We have analyzed the relaxation properties of all (31)P nuclei in an RNA cUUCGg tetraloop model hairpin at proton magnetic field strengths of 300, 600 and 900 MHz in solution. Significant H, P dipolar contributions to R (1) and R (2) relaxation are observed in a protonated RNA sample at 600 MHz. These contributions can be suppressed using a perdeuterated RNA sample. In order to interpret the (31)P relaxation data (R (1), R (2)), we measured the (31)P chemical shift anisotropy (CSA) by solid-state NMR spectroscopy under various salt and hydration conditions. A value of 178.5 ppm for the (31)P CSA in the static state (S (2) = 1) could be determined. In order to obtain information about fast time scale dynamics we performed a modelfree analysis on the basis of our relaxation data. The results show that subnanosecond dynamics detected around the phosphodiester backbone are more pronounced than the dynamics detected for the ribofuranosyl and nucleobase moieties of the individual nucleotides (Duchardt and Schwalbe, J Biomol NMR 32:295-308, 2005; Ferner et al., Nucleic Acids Res 36:1928-1940, 2008). Furthermore, the dynamics of the individual phosphate groups seem to be correlated to the 5' neighbouring nucleobases.

Constructing atomic-resolution RNA structural ensembles using MD and motionally decoupled NMR RDCs

A broad structural landscape often needs to be characterized in order to fully understand how regulatory RNAs perform their biological functions at the atomic level. We present a protocol for visualizing thermally accessible RNA conformations at atomic-resolution and with timescales extending up to milliseconds. The protocol combines molecular dynamics (MD) simulations with experimental residual dipolar couplings (RDCs) measured in partially aligned (13)C/(15)N isotopically enriched elongated RNA samples. The structural ensembles generated in this manner provide insights into RNA dynamics and its role in functionally important transitions.

Effect of local sugar and base geometry on 13C and 15N magnetic shielding anisotropy in DNA nucleosides

Density functional theory was employed to study the dependence of 13C and 15N magnetic shielding tensors on the glycosidic torsion angle (chi) and conformation of the sugar ring in 2'-deoxyadenosine, 2'-deoxyguanosine, 2'-deoxycytidine, and 2'-deoxythymidine. In general, the magnetic shielding of the glycosidic nitrogens and the sugar carbons was found to depend on both the conformation of the sugar ring and chi. Our calculations indicate that the magnetic shielding anisotropy of the C6 atom in pyrimidine and the C8 atom in purine bases depends strongly on chi. The remaining base carbons were found to be insensitive to both sugar pucker and chi re-orientation. These results call into question the underlying assumptions of currently established methods for interpreting residual chemical shift anisotropies and 13C and 15N auto- and cross-correlated relaxation rates and highlight possible limitations of DNA applications of these methods.

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

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

Dynamics of large elongated RNA by NMR carbon relaxation

We present an NMR strategy for characterizing picosecond-to-nanosecond internal motions in uniformly 13C/15N-labeled RNAs that combines measurements of R1, R1rho, and heteronuclear 13C{1H} NOEs for protonated base (C2, C5, C6, and C8) and sugar (C1') carbons with a domain elongation strategy for decoupling internal from overall motions and residual dipolar coupling (RDC) measurements for determining the average RNA global conformation and orientation of the principal axis of the axially symmetric rotational diffusion. TROSY-detected pulse sequences are presented for the accurate measurement of nucleobase carbon R1 and R1rho rates in large RNAs. The relaxation data is analyzed using a model free formalism which takes into account the very high anisotropy of overall rotational diffusion (Dratio approximately 4.7), asymmetry of the nucleobase CSAs and noncollinearity of C-C, C-H dipolar and CSA interactions under the assumption that all interaction tensors for a given carbon experience identical isotropic internal motions. The approach is demonstrated and validated on an elongated HIV-1 TAR RNA (taum approximately 18 ns) both in free form and bound to the ligand argininamide (ARG). Results show that, while ARG binding reduces the amplitude of collective helix motions and local mobility at the binding pocket, it leads to a drastic increase in the local mobility of "spacer" bulge residues linking the two helices which undergo virtually unrestricted internal motions (S2 approximately 0.2) in the ARG bound state. Our results establish the ability to quantitatively study the dynamics of RNAs which are significantly larger and more anisotropic than customarily studied by NMR carbon relaxation.

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.

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

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

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.

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

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

Resolving fast and slow motions in the internal loop containing stem-loop 1 of HIV-1 that are modulated by Mg2+ binding: role in the kissing-duplex structural transition

Stem loop 1 (SL1) is a highly conserved hairpin in the 5'-leader of the human immunodeficiency virus type I that forms a metastable kissing dimer that is converted during viral maturation into a stable duplex with the aid of the nucleocapsid (NC) protein. SL1 contains a highly conserved internal loop that promotes the kissing-duplex transition by a mechanism that remains poorly understood. Using NMR, we characterized internal motions induced by the internal loop in an SL1 monomer that may promote the kissing-duplex transition. This includes micro-to-millisecond secondary structural transitions that cause partial melting of three base-pairs above the internal loop making them key nucleation sites for exchanging strands and nanosecond rigid-body stem motions that can help bring strands into spatial register. We show that while Mg2+ binds to the internal loop and arrests these internal motions, it preserves and/or activates local mobility at internal loop residues G272 and G273 which are implicated in NC binding. By stabilizing SL1 without compromising the accessibility of G272 and G273 for NC binding, Mg2+ may increase the dependence of the kissing-duplex transition on NC binding thus preventing spontaneous transitions from taking place and ensuring that viral RNA and protein maturation occur in concert.

Structural plasticity and Mg2+ binding properties of RNase P P4 from combined analysis of NMR residual dipolar couplings and motionally decoupled spin relaxation

The P4 helix is an essential element of ribonuclease P (RNase P) that is believed to bind catalytically important metals. Here, we applied a combination of NMR residual dipolar couplings (RDCs) and a recently introduced domain-elongation strategy for measuring "motionally decoupled" relaxation data to characterize the structural dynamics of the P4 helix from Bacillus subtilis RNase P. In the absence of divalent ions, the two P4 helical domains undergo small amplitude (approximately 13 degrees) collective motions about an average interhelical angle of 10 degrees. The highly conserved U7 bulge and helical residue C8, which are proposed to be important for substrate recognition and metal binding, are locally mobile at pico- to nanosecond timescales and together form the pivot point for the collective domain motions. Chemical shift mapping reveals significant association of Mg2+ ions at the P4 major groove near the flexible pivot point at residues (A5, G22, G23) previously identified to bind catalytically important metals. The Mg2+ ions do not, however, significantly alter the structure or dynamics of P4. Analysis of results in the context of available X-ray structures of the RNA component of RNase P and structural models that include the pre-tRNA substrate suggest that the internal motions observed in P4 likely facilitate adaptive changes in conformation that take place during folding and substrate recognition, possibly aided by interactions with Mg2+ ions. Our results add to a growing view supporting the existence of functionally important internal motions in RNA occurring at nanosecond timescales.

Residue-Specific 13C‘ CSA Tensor Principal Components for Ubiquitin: Correlation between Tensor Components and Hydrogen Bonding

The residue-specific 13C' CSA tensor principal components, sigma(11), sigma(22), sigma(33), and the tensor orientation defined by the rotation angles beta and gamma have been determined by solution NMR for uniformly labeled ubiquitin partially aligned in four different media. Spurious chemical shift deviations due to solvent effects were corrected with an offset calculated by linear regression of the residual dipolar couplings and chemical shifts at increasing alignment strengths. Analysis of this effect revealed no obvious correlation to solvent exposure. Data obtained in solution from a protein offer a better sampling of 13C' CSA for different amino acid types in a complex heterogeneous environment, thereby allowing for the evaluation of structural variables that would be challenging to achieve by other methods. The 13C' CSA principal components cluster about the average values previously determined, and experimental correlations observed between sigma(11), sigma(22) tensorial components and C'O...H(N) hydrogen bonding are discussed. The inverse association of sigma(11) and sigma(22) exemplify the calculated and solid-state NMR observed effect on the tensor components by hydrogen bonding. We also show that 13C' CSA tensors are sensitive to hydrogen-bond length but not hydrogen-bond angle. This differentiation was previously unavailable. Similarly, hydrogen bonding to the conjugated NH of the same peptide plane has no detectable effect. Importantly, the observed weak correlations signify the presence of confounding influences such as nearest-neighbor effects, side-chain conformation, electrostatics, and other long-range factors to the 13C' CSA tensor. These analyses hold future potential for exploration provided that more accurate data from a larger number of proteins and alignments become available.

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

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

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.

Pseudo-CSA Restraints for NMR Refinement of Nucleic Acid Structure

Upon alignment of oligonucleotides in a magnetic field, the downfield TROSY component of the 13C-{1H} doublet changes its resonance frequency as a result of residual 13C-1H dipolar coupling (RDC) and residual 13C chemical shift anisotropy (RCSA), and the sum of these two second rank tensors is referred to as the pseudo-CSA. The experimentally measured difference in the resonance frequency of the 13C TROSY component in the aligned and isotropic samples is referred to as residual pseudo-CSA (RPCSA), and it can be used directly as a restraint during structure calculation. Because measurement of the RPCSA involves detection of the narrow TROSY 13C doublet component, it is applicable to systems with larger rotational correlation times than RDC measurement. The method is demonstrated for structure refinement of the helical region of a 24-nt stem-loop segment or ribosomal helix-35, uniformly enriched in 13C and 15N, with RPCSA values measured at 5 and 25 degrees C. Substantial cross-validated improvements in structural accuracy are obtained upon incorporation of RPCSA restraints.