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.Copyright 2001 John Wiley & Sons, Inc.     

A unifying probabilistic framework for analyzing residual dipolar couplings

Residual dipolar couplings provide complementary information to the nuclear Overhauser effect measurements that are traditionally used in biomolecular structure determination by NMR. In a de novo structure determination, however, lack of knowledge about the degree and orientation of molecular alignment complicates the analysis of dipolar coupling data. We present a probabilistic framework for analyzing residual dipolar couplings and demonstrate that it is possible to estimate the atomic coordinates, the complete molecular alignment tensor, and the error of the couplings simultaneously. As a by-product, we also obtain estimates of the uncertainty in the coordinates and the alignment tensor. We show that our approach encompasses existing methods for determining the alignment tensor as special cases, including least squares estimation, histogram fitting, and elimination of an explicit alignment tensor in the restraint energy.     

15N-1H Residual dipolar coupling analysis of native and alkaline-K79A Saccharomyces cerevisiae cytochrome c

Residual dipolar couplings (RDCs) and pseudocontact shifts are experimentally accessible properties in nuclear magnetic resonance that are related to structural parameters and to the magnetic susceptibility anisotropy. We have determined RDCs due to field-induced orientation of oxidized-K79A and reduced cytochrome c at pH 7.0 and oxidized-K79A cytochrome c at pH 11.1 through measurements of amide (15)N-(1)H (1)J couplings at 800 and 500 MHz. The pH 7.0 RDCs for Fe(III)- and Fe(II)-cytochrome c together with available nuclear Overhauser effects were used to recalculate solution structures that were consistent with both sets of constraints. Molecular magnetic susceptibility anisotropy values were calculated for both redox states of the protein. By subtracting the residual dipolar couplings (RDCs) of the reduced form from those of the oxidized form measured at the same magnetic field (800 MHz), we found the RDC contribution of the paramagnetic metal ion in the oxidized protein. The magnetic susceptibility anisotropy, which was calculated from the structure, was found to be the same as that of the paramagnetic metal ion obtained independently from pseudocontact shifts, thereby indicating that the elements of secondary structure either are rigid or display the same mobility in both oxidation states. The residual dipolar coupling values of the alkaline-K79A form are small with respect to those of oxidized native cytochrome, whereas the pseudocontact shifts are essentially of the same magnitude, indicating local mobility. Importantly, this is the first time that mobility has been found through comparison of RDCs with pseudocontact shifts.     

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.     

The utility of residual dipolar couplings in detecting motion in carbohydrates: application to sucrose

The solution structure and dynamics of sucrose are examined using a combination of NMR residual dipolar coupling and molecular mechanics force fields. It is found that the alignment tensors of the individual rings are different, and that fitting 35 measured residual dipolar couplings to structures with specific phi, psi values indicates the presence of three major conformations: phi, psi=(120 degrees ,270 degrees), (45 degrees, 300 degrees) and (90 degrees ,180 degrees). Furthermore, fitting two structures simultaneously to the 35 residual dipolar couplings results in a substantial improvement in the fits. The existence of multiple conformations having similar stabilities is a strong indication of motion, due to the interconversion among these states. Results from four molecular mechanics force fields are in general agreement with the experimental results. However, there are major disagreements between force fields. Because fits of residual dipolar couplings to structures are dependent on the force field used to calculate the structures, multiple force fields were used to interpret NMR data. It is demonstrated that the pucker of the fructofuranosyl ring affects the calculated potential energy surface, and the fit to the residual dipolar couplings data. Previously published 13C nuclear relaxation results suggesting that sucrose is rigid are not inconsistent with the present results when motional timescales are considered.     

Structure determination of membrane proteins by NMR spectroscopy.

Current strategies for determining the structures of membrane proteins in lipid environments by NMR spectroscopy rely on the anisotropy of nuclear spin interactions, which are experimentally accessible through experiments performed on weakly and completely aligned samples. Importantly, the anisotropy of nuclear spin interactions results in a mapping of structure to the resonance frequencies and splittings observed in NMR spectra. Distinctive wheel-like patterns are observed in two-dimensional 1H-15N heteronuclear dipolar/15N chemical shift PISEMA (polarization inversion spin-exchange at the magic angle) spectra of helical membrane proteins in highly aligned lipid bilayer samples. One-dimensional dipolar waves are an extension of two-dimensional PISA (polarity index slant angle) wheels that map protein structures in NMR spectra of both weakly and completely aligned samples. Dipolar waves describe the periodic wave-like variations of the magnitudes of the heteronuclear dipolar couplings as a function of residue number in the absence of chemical shift effects. Since weakly aligned samples of proteins display these same effects, primarily as residual dipolar couplings, in solution NMR spectra, this represents a convergence of solid-state and solution NMR approaches to structure determination.     

Recognition of protein folds via dipolar couplings

Alignment of proteins in dilute liquid crystalline medium gives rise to residual dipolar couplings which provide orientational information of vectors connecting the interacting nuclei. Considering that proteins are mainly composed of regular secondary structures in a finite number of different mutual orientations, main chain dipolar couplings appear sufficient to reveal structural resemblance. Similarity between dipolar couplings measured from a protein and corresponding values computed from a known structure imply homologous structures. For dissimilar structures the agreement between experimental and calculated dipolar couplings remains poor. In this way protein folds can be readily recognized prior to a comprehensive structure determination. This approach has been demonstrated by showing the similarity in fold between the hitherto unknown structure of calerythrin and sarcoplasmic calcium-binding proteins from Nereis diversicolor and Branchiostoma lanceolatum with known crystal structures.     

High-resolution X-ray and NMR structures of the SMN Tudor domain: conformational variation in the binding site for symmetrically dimethylated arginine residues

The SMN protein, which is linked to spinal muscular atrophy (SMA), plays an important role in the assembly of the spliceosomal small nuclear ribonucleoprotein complexes. This function requires binding of SMN to the arginine-glycine (RG) rich C-terminal tails of the Sm proteins, which contain symmetrically dimethylated arginine residues (sDMA) in vivo. Using NMR titrations, we show that the SMN Tudor domain recognizes these sDMAs in the methylated RG repeats. Upon complex formation a cluster of conserved aromatic residues in the SMN Tudor domain interacts with the sDMA methyl groups. We present two high resolution structures of the uncomplexed SMN Tudor domain, a 1.8A crystal structure and an NMR structure that has been refined against a large number of backbone and side-chain residual dipolar couplings. The backbone conformation of both structures is very similar, however, differences are observed for the cluster of conserved aromatic side-chains in the sDMA binding pocket. In order to validate these variations we introduce a novel application of residual dipolar couplings for aromatic rings. We show that structural information can be derived from aromatic ring residual dipolar couplings, even in the presence of internal motions such as ring flipping. These residual dipolar couplings and ring current shifts independently confirm that the SMN Tudor domain adopts two different conformations in the sDMA binding pocket. The observed structural variations may play a role for the recognition of sDMAs.     

Solution structure of a LewisX analogue by off-resonance 1H NMR spectroscopy without use of an internal distance reference

Summary A recent ¹H NMR method has been applied to the determination of the solution structure and internal dynamics of a synthetic mixed C/O trisaccharide related to sialyl Lewis^x. Varying the rf field offset in ROESY-type experiments enabled the measurement of longitudinal and transverse dipolar cross-relaxation rates with high accuracy. Assuming that for each proton pair the motion could be represented by a single exponential autocorrelation function, it was possible to derive geometrical parameters (r) and dynamic parameters _cp. With this assumption, 224 cross-relaxation rates have been transformed into 30 interproton distance constraints and 30 dipolar correlation times. The distance constraints have been used in a simulated-annealing procedure. This trisaccharide exhibits a structure close to the O-glycosidic analogue, but its flexibility seems highly reduced. On the basis of the determined structure and dynamics, it is shown that no conformational exchange occurs, the molecule existing in the form of a unique family in aqueous solution. In order to assess the quality of the resulting structures and to validate this new experimental procedure of distance extraction, we finally compare these solution structures to the ones obtained using three different sets of distances deduced from three choices of internal reference. It appears that this procedure allows the determination of the most precise and accurate solution.Abbreviations COSY correlation spectroscopy - NOE nuclear Overhauser enhancement - NOESY nuclear Overhauser enhancement spectroscopy; rmsd, root-mean-square deviation - ROESY rotating frame Overhauser enhancement spectroscopy - SLe^x sialyl Lewis^x - TOCSY total correlation spectroscopy     

Quantitative internuclear distancesvia two-dimensional nuclear magnetic resonance spectra: A test case and a DNA octamer duplex

Two-dimensional proton nuclear magnetic resonance nuclear Overhauser effect experiments have been performed at a series of mixing times on proflavine and on a DNA octamer duplex [d-(GGAATTCC)]₂ in solution. Using the complete matrix approach recently explored theoretically (Keepers and James, 1984), proton-proton internuclear distances were determined quantitatively for proflavine from the two-dimensional nuclear Overhauser effect results. Since proflavine is a rigid molecule with X-ray crystal structure determined, interproton distances obtained from the two-dimensional nuclear Overhauser effect experiments in solution can be compared with those for the crystalline compound agreement is better than 10 %. Experimental two-dimensional nuclear Overhauser effect spectral data for [d-(GGAATTCC)]₂ were analyzed by comparison with theoretical two-dimensional nuclear Overhauser effect spectra at each mixing time calculated using the complete 70 × 70 relaxation matrix. The theoretical spectra were calculated using two structures: a standard B-form DNA structure and an energy-minimized structure based on similarity of the octamer's six internal residues with those of [d-(CGCGAATTCGCG)]₂, for which the crystal structure has been determined. Neither the standard B-DNA nor the energy-minimized structure yield theoretical two-dimensional nuclear Overhauser effect spectra which accurately reproduce all experimental peak intensities. But many aspects of the experimental spectra can be represented by both the B-DNA and the energy-minimized structure. In general, the energy-minimized structure yields theoretical two-dimensional nuclear Overhauser effect spectra which mimic many, if not all, features of the experimental, spectra including structural characteristics at the purine-pyrimidine junction.     

Structure determination of a new protein from backbone-centered NMR data and NMR-assisted structure prediction

Targeting of proteins for structure determination in structural genomic programs often includes the use of threading and fold recognition methods to exclude proteins belonging to well-populated fold families, but such methods can still fail to recognize preexisting folds. The authors illustrate here a method in which limited amounts of structural data are used to improve an initial homology search and the data are subsequently used to produce a structure by data-constrained refinement of an identified structural template. The data used are primarily NMR-based residual dipolar couplings, but they also include additional chemical shift and backbone-nuclear Overhauser effect data. Using this methodology, a backbone structure was efficiently produced for a 10 kDa protein (PF1455) from Pyrococcus furiosus. Its relationship to existing structures and its probable function are discussed.     

Structure and disorder in the ribonuclease S-peptide probed by NMR residual dipolar couplings

NMR residual dipolar couplings for the S-peptide of ribonuclease A aligned in C8E5/n-octanol liquid crystals are consistent with the presence of a native-like alpha-helix structure undergoing dynamic fraying. Residues 3-13, which correspond to the first alpha-helix of ribonuclease A, show couplings that become more negative at low temperature and in the presence of salt, conditions which stabilize alpha-helical structure in the S-peptide. By contrast, dipolar couplings from the N and C termini of the peptide are close to zero and remain nearly invariant with changes in solution conditions. Torsion angle dynamics simulations using a gradient of dihedral restraint bounds that increase from the center to the ends of the peptide reproduce the experimentally observed sequence dependence of dipolar couplings. The magnitudes of residual dipolar couplings depend on the anisotropy of the solute. Native proteins often achieve nearly spherical shapes due to the hydrophobic effect. Embryonic partially folded structures such as the S-peptide alpha-helix have an intrinsically greater potential for anisotropy that can result in sizable residual dipolar couplings in the absence of long-range structure.     

High-resolution NMR structure of an RNA model system: the 14-mer cUUCGg tetraloop hairpin RNA

We present a high-resolution nuclear magnetic resonance (NMR) solution structure of a 14-mer RNA hairpin capped by cUUCGg tetraloop. This short and very stable RNA presents an important model system for the study of RNA structure and dynamics using NMR spectroscopy, molecular dynamics (MD) simulations and RNA force-field development. The extraordinary high precision of the structure (root mean square deviation of 0.3 Å) could be achieved by measuring and incorporating all currently accessible NMR parameters, including distances derived from nuclear Overhauser effect (NOE) intensities, torsion-angle dependent homonuclear and heteronuclear scalar coupling constants, projection-angle-dependent cross-correlated relaxation rates and residual dipolar couplings. The structure calculations were performed with the program CNS using the ARIA setup and protocols. The structure quality was further improved by a final refinement in explicit water using OPLS force field parameters for non-bonded interactions and charges. In addition, the 2′-hydroxyl groups have been assigned and their conformation has been analyzed based on NOE contacts. The structure currently defines a benchmark for the precision and accuracy amenable to RNA structure determination by NMR spectroscopy. Here, we discuss the impact of various NMR restraints on structure quality and discuss in detail the dynamics of this system as previously determined.     

Evaluation of uncertainty in alignment tensors obtained from dipolar couplings

Residual dipolar couplings and their corresponding alignment tensors are useful for structural analysis of macromolecules. The error in an alignment tensor, derived from residual dipolar couplings on the basis of a known structure, is determined not only by the accuracy of the measured couplings but also by the uncertainty in the structure (structural noise). This dependence is evaluated quantitatively on the basis of simulated structures using Monte-Carlo type analyses. When large numbers of dipolar couplings are available, structural noise is found to result in a systematic underestimate of the magnitude of the alignment tensor. Particularly in cases where only few dipolar couplings are available, structural noise can cause significant errors in best-fitted alignment tensor values, making determination of the relative orientation of small fragments and evaluation of local backbone mobility from dipolar couplings difficult. An example for the protein ubiquitin demonstrates the inherent limitations in characterizing motions on the basis of local alignment tensor magnitudes.     

High-resolution structure determination of the CylR2 homodimer using paramagnetic relaxation enhancement and structure-based prediction of molecular alignment

Structure determination of homooligomeric proteins by NMR spectroscopy is difficult due to the lack of chemical shift perturbation data, which is very effective in restricting the binding interface in heterooligomeric systems, and the difficulty of obtaining a sufficient number of intermonomer distance restraints. Here we solved the high-resolution solution structure of the 15.4 kDa homodimer CylR2, the regulator of cytolysin production from Enterococcus faecalis, which deviates by 1.1 angstroms from the previously determined X-ray structure. We studied the influence of different experimental information such as long-range distances derived from paramagnetic relaxation enhancement, residual dipolar couplings, symmetry restraints and intermonomer Nuclear Overhauser Effect restraints on the accuracy of the derived structure. In addition, we show that it is useful to combine experimental information with methods of ab initio docking when the available experimental data are not sufficient to obtain convergence to the correct homodimeric structure. In particular, intermonomer distances may not be required when residual dipolar couplings are compared to values predicted on the basis of the charge distribution and the shape of ab initio docking solutions.     

Validating the use of database potentials in protein structure determination by NMR

The refinement of protein structures determined by nuclear magnetic resonance spectroscopy against database potentials of mean force allows for the exclusion of unfavourable conformations of the protein backbone during a structure calculation, resulting in protein structures with a marked improvement in Ramachandran statistics. In this communication, we use multiple sets of residual dipolar couplings as quality assessment criteria for several proteins and show that not only do the Ramachandran and structural quality statistics improve, but a significant improvement in the accuracy of structures is achieved upon refinement.     

What is the average conformation of bacteriophage T4 lysozyme in solution? A domain orientation study using dipolar couplings measured by solution NMR

Lysozyme from T4 bacteriophage is comprised of two domains that are both involved in binding substrate. Although wild-type lysozyme has been exclusively crystallized in a closed form that is similar to the peptidoglycan-bound conformation, a more open structure is thought to be required for ligand binding. To determine the relative arrangement of domains within T4 lysozyme in the solution state, dipolar couplings were measured in several different dilute liquid crystalline media by solution NMR methods. The dipolar coupling data were analyzed with a domain orientation procedure described previously that utilizes high- resolution X-ray structures. The cleft between the domains is significantly larger in the average solution structure than what is observed in the X-ray structure of the ligand-free form of the protein (approximately 17 degrees closure from solution to X-ray structures). A comparison of the solution domain orientation with X-ray-derived structures in the protein data base shows that the solution structure resembles a crystal structure obtained for the M6I mutant. Dipolar couplings were also measured on the lysozyme mutant T21C/T142C, which was oxidized to form an inter-domain disulfide bond (T4SS). In this case, the inter-domain solution structure was found to be more closed than was observed in the crystal (approximately 11 degrees). Direct refinement of lysozyme crystal structures with the measured dipolar couplings using the program CNS, establishes that this degree of closure can be accommodated whilst maintaining the inter-domain cystine bond. The differences between the average solution conformations obtained using dipolar couplings and the crystal conformations for both forms of lysozyme investigated in this study illustrate the impact that crystal packing interactions can have on the arrangement of domains within proteins and the importance of alternative methods to X-ray crystallography for evaluating inter-domain structure.     

Improvement of hydrogen bond geometry in protein NMR structures by residual dipolar couplings – an assessment of the interrelation of NMR restraints

We have examined how the hydrogen bond geometry in three different proteins is affected when structural restraints based on measurements of residual dipolar couplings are included in the structure calculations. The study shows, that including restraints based solely on (1)H(N)-(15)N residual dipolar couplings has pronounced impact on the backbone rmsd and Ramachandran plot but does not improve the hydrogen bond geometry. In the case of chymotrypsin inhibitor 2 the addition of (13)CO-(13)C(alpha) and (15)N-(13)CO one bond dipolar couplings as restraints in the structure calculations improved the hydrogen bond geometry to a quality comparable to that obtained in the 1.8 A resolution X-ray structure of this protein. A systematic restraint study was performed, in which four types of restraints, residual dipolar couplings, hydrogen bonds, TALOS angles and NOEs, were allowed in two states. This study revealed the importance of using several types of residual dipolar couplings to get good hydrogen bond geometry. The study also showed that using a small set of NOEs derived only from the amide protons, together with a full set of residual dipolar couplings resulted in structures of very high quality. When reducing the NOE set, it is mainly the side-chain to side-chain NOEs that are removed. Despite of this the effect on the side-chain packing is very small when a reduced NOE set is used, which implies that the over all fold of a protein structure is mainly determined by correct folding of the backbone.