A Method for Direct Determination of Helical Parameters in Nucleic Acids Using Residual Dipolar Couplings

Guanine bases in DNA G-quadruplex adopt nonplanar geometries owing to solvation and base pairing

The effect of base pairing and solvation on pyramidalization of the glycosidic nitrogen found in the residues of parallel G-quadruplex with NDB ID UDF062 is analyzed and explained with theoretical calculations. The extent of the pyramidalization depends on the local geometry of the 2'-deoxyguanosine residues, namely on their glycosidic torsion and sugar pucker, which are predetermined by the 3D-architecture of G-quadruplex. Pyramidal inversion of the glycosidic nitrogen found in 2'-deoxyguanosines of G-quadruplex is induced owing to site-specifically coordinated solvent. Different adiabatic structural constraints used for fixing the base-to-sugar orientation of 2'-deoxyguanosine in geometry optimizations result in different extents of pyramidalization and induce pyramidal inversion of the glycosidic nitrogen. These model geometry constraints helped us analyze the effect of real constraints represented by explicit molecular environment of selected residues of the G-quadruplex. The maximal extent of the glycosidic nitrogen pyramidalization found in the high-resolution crystal structure corresponds to the calculation to deformation energy of only 1 kcal mol(-1). The out-of-plane deformations of nucleobases thus provide a way for compensating the site-specific external environmental stress on the G-quadruplex.

Structure of a conserved retroviral RNA packaging element by NMR spectroscopy and cryo-electron tomography

The 5'-untranslated regions of all gammaretroviruses contain a conserved "double-hairpin motif" (Ψ(CD)) that is required for genome packaging. Both hairpins (SL-C and SL-D) contain GACG tetraloops that, in isolated RNAs, are capable of forming "kissing" interactions stabilized by two intermolecular G-C base pairs. We have determined the three-dimensional structure of the double hairpin from the Moloney murine leukemia virus ([Ψ(CD)](2), 132 nt, 42.8 kDa) using a (2)H-edited NMR-spectroscopy-based approach. This approach enabled the detection of (1)H-(1)H dipolar interactions that were not observed in previous studies of isolated SL-C and SL-D hairpin RNAs using traditional (1)H-(1)H correlated and (1)H-(13)C-edited NMR methods. The hairpins participate in intermolecular cross-kissing interactions (SL-C to SL-D' and SLC' to SL-D) and stack in an end-to-end manner (SL-C to SL-D and SL-C' to SL-D') that gives rise to an elongated overall shape (ca 95 Å×45 Å×25 Å). The global structure was confirmed by cryo-electron tomography (cryo-ET), making [Ψ(CD)](2) simultaneously the smallest RNA to be structurally characterized to date by cryo-ET and among the largest to be determined by NMR. Our findings suggest that, in addition to promoting dimerization, [Ψ(CD)](2) functions as a scaffold that helps initiate virus assembly by exposing a cluster of conserved UCUG elements for binding to the cognate nucleocapsid domains of assembling viral Gag proteins.

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.

Revisiting the planarity of nucleic acid bases: Pyramidilization at glycosidic nitrogen in purine bases is modulated by orientation of glycosidic torsion

We describe a novel, fundamental property of nucleobase structure, namely, pyramidilization at the N1/9 sites of purine and pyrimidine bases. Through a combined analyses of ultra-high-resolution X-ray structures of both oligonucleotides extracted from the Nucleic Acid Database and isolated nucleotides and nucleosides from the Cambridge Structural Database, together with a series of quantum chemical calculations, molecular dynamics (MD) simulations, and published solution nuclear magnetic resonance (NMR) data, we show that pyramidilization at the glycosidic nitrogen is an intrinsic property. This property is common to isolated nucleosides and nucleotides as well as oligonucleotides—it is also common to both RNA and DNA. Our analysis suggests that pyramidilization at N1/9 sites depends in a systematic way on the local structure of the nucleoside. Of note, the pyramidilization undergoes stereo-inversion upon reorientation of the glycosidic bond. The extent of the pyramidilization is further modulated by the conformation of the sugar ring. The observed pyramidilization is more pronounced for purine bases, while for pyrimidines it is negligible. We discuss how the assumption of nucleic acid base planarity can lead to systematic errors in determining the conformation of nucleotides from experimental data and from unconstrained MD simulations.

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.

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.

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

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

NMR structure of the 101-nucleotide core encapsidation signal of the Moloney murine leukemia virus

The full length, positive-strand genome of the Moloney Murine Leukemia Virus contains a "core encapsidation signal" that is essential for efficient genome packaging during virus assembly. We have determined the structure of a 101-nucleotide RNA that contains this signal (called mPsi) using a novel isotope-edited NMR approach. The method is robust and should be generally applicable to larger RNAs. mPsi folds into three stem loops, two of which (SL-C and SL-D) co-stack to form an extended helix. The third stem loop (SL-B) is connected to SL-C by a flexible, four-nucleotide linker. The structure contains five mismatched base-pairs, an unusual C.CG base-triple platform, and a novel "A-minor K-turn," in which unpaired adenosine bases A340 and A341 of a GGAA bulge pack in the minor groove of a proximal stem, and a bulged distal uridine (U319) forms a hydrogen bond with the phosphodiester of A341. Phylogenetic analyses indicate that these essential structural elements are conserved among the murine C-type retroviruses.

Residual dipolar couplings in nucleic acid structure determination

Solution NMR spectroscopy of nucleic acids has been limited by the short-range nature of the nuclear Overhauser effect and scalar coupling restraints normally used in structure determination. The addition of residual dipolar couplings, obtained from slightly oriented mixtures, provides bond vector angles relative to a universal alignment tensor. The accurate determination of helix curvature, domain orientation and the stoichiometry of homomultimeric nucleic acid complexes is now possible.

Refinement of local and long-range structural order in theophylline-binding RNA using (13)C-(1)H residual dipolar couplings and restrained molecular dynamics

13C-(1)H residual dipolar couplings (RDC) have been measured for the bases and sugars in the theophylline-binding RNA aptamer, dissolved in filamentous phage medium, and used to investigate the long-range structural and dynamic behavior of the molecule in the solution state. The orientation dependent RDC provide additional restraints to further refine the overall structure of the RNA-theophylline complex, whose long-range order was poorly defined in the NOE-based structural ensemble. Structure refinement using RDC normally assumes that molecular alignment can be characterized by a single tensor and that the molecule is essentially rigid. To address the validity of this assumption for the complex of interest, we have analyzed distinct domains of the RNA molecule separately, so that local structure and alignment tensors experienced by each region are independently determined. Alignment tensors for the stem regions of the molecule were allowed to float freely during a restrained molecular dynamics structure refinement protocol and found to converge to similar magnitudes. During the second stage of the calculation, a single alignment tensor was thus applied for the whole molecule and an average molecular conformation satisfying all experimental data was determined. Semirigid-body molecular dynamics calculations were used to reorient the refined helical regions to a relative orientation consistent with this alignment tensor, allowing determination of the global conformation of the molecule. Simultaneously, the local structure of the theophylline-binding core of the molecule was refined under the influence of this common tensor. The final ensemble has an average pairwise root mean square deviation of 1.50 +/- 0.19 A taken over all heavy atoms, compared to 3.5 +/- 1.1 A for the ensemble determined without residual dipolar coupling. This study illustrates the importance of considering both the local and long-range nature of RDC when applying these restraints to structure refinements of nucleic acids.

Dipolar couplings as a probe of molecular dynamics and structure in solution

The introduction of residual dipolar coupling methodology has increased the scope of structural biological problems that can be addressed by NMR spectroscopy. Conformational changes, the relative orientation of domains, and intermolecular complexes can now be characterized accurately and rapidly using NMR. The development of residual dipolar coupling methodology for the rapid recognition of homologous protein folds and for studies of submillisecond timescale dynamics has also seen considerable progress.

Incorporating residual dipolar couplings into the NMR solution structure determination of nucleic acids

NMR solution structures of nucleic acids are generally less well defined than similar-sized proteins. Most NMR structures of nucleic acids are defined only by short-range interactions, such as intrabase-pair or sequential nuclear Overhauser effects (NOEs), and J-coupling constants, and there are no long-range structural data on the tertiary structure. Residual dipolar couplings represent an extremely valuable source of distance and angle information for macromolecules but they average to zero in isotropic solutions. With the recent advent of general methods for partial alignment of macromolecules in solution, residual dipolar couplings are rapidly becoming indispensable constraints for solution NMR structural studies. These residual dipolar couplings give long-range global structural information and thus complement the strictly local structural data obtained from standard NOE and torsion angle constraints. Such global structural data are especially important in nucleic acids due to the more elongated, less-globular structure of many DNAs and RNAs. Here we review recent progress in application of residual dipolar couplings to structural studies of nucleic acids. We also present results showing how refinement procedures affect the final solution structures of nucleic acids.

NMR methodology for the study of nucleic acids

During the past few years, NMR methodology for the study of nucleic acids has benefited from new developments that greatly improved state-of-the-art technology for the precise determination of three-dimensional structures. Substantial progress has been made in designing experimental protocols for the measurement of residual dipolar couplings, in sensitivity optimization of triple-resonance experiments and in detection of hydrogen bonds and in developing computational methods for structure refinement using NMR restraints.

Improving the accuracy of NMR structures of DNA by means of a database potential of mean force describing base-base positional interactions

NMR structure determination of nucleic acids presents an intrinsically difficult problem since the density of short interproton distance contacts is relatively low and limited to adjacent base pairs. Although residual dipolar couplings provide orientational information that is clearly helpful, they do not provide translational information of either a short-range (with the exception of proton-proton dipolar couplings) or long-range nature. As a consequence, the description of the nonbonded contacts has a major impact on the structures of nucleic acids generated from NMR data. In this paper, we describe the derivation of a potential of mean force derived from all high-resolution (2 A or better) DNA crystal structures available in the Nucleic Acid Database (NDB) as of May 2000 that provides a statistical description, in simple geometric terms, of the relative positions of pairs of neighboring bases (both intra- and interstrand) in Cartesian space. The purpose of this pseudopotential, which we term a DELPHIC base-base positioning potential, is to bias sampling during simulated annealing refinement to physically reasonable regions of conformational space within the range of possibilities that are consistent with the experimental NMR restraints. We illustrate the application of the DELPHIC base-base positioning potential to the structure refinement of a DNA dodecamer, d(CGCGAATTCGCG)(2), for which NOE and dipolar coupling data have been measured in solution and for which crystal structures have been determined. We demonstrate by cross-validation against independent NMR observables (that is, both residual dipolar couplings and NOE-derived intereproton distance restraints) that the DELPHIC base-base positioning potential results in a significant increase in accuracy and obviates artifactual distortions in the structures arising from the limitations of conventional descriptions of the nonbonded contacts in terms of either Lennard-Jones van der Waals and electrostatic potentials or a simple van der Waals repulsion potential. We also demonstrate, using experimental NMR data for a complex of the male sex determining factor SRY with a duplex DNA 14mer, which includes a region of highly unusual and distorted DNA, that the DELPHIC base-base positioning potential does not in any way hinder unusual interactions and conformations from being satisfactorily sampled and reproduced. We expect that the methodology described in this paper for DNA can be equally applied to RNA, as well as side chain-side chain interactions in proteins and protein-protein complexes, and side chain-nucleic acid interactions in protein-nucleic acid complexes. Further, this approach should be useful not only for NMR structure determination but also for refinement of low-resolution (3-3.5 A) X-ray data.