Model-free analysis of protein backbone motion from residual dipolar couplings

On the basis of the measurement of NH residual dipolar couplings (RDCs) in 11 different alignment media, an RDC-based order parameter is derived for each residue in the protein ubiquitin. Dipolar couplings are motionally averaged in the picosecond to millisecond time range and, therefore, reflect motion slower than the inverse overall tumbling correlation time of the protein. It is found that there is considerable motion that is slower than the correlation time and could not be detected with previous NMR methodology. Amplitudes and anisotropies of the motion can be derived from the model-free analysis. The method can be applied provided that at least five sufficiently different alignment media can be found for the biomolecule under investigation.     

Spatial structure of peptides determined by residual dipolar couplings analysis

The gated decoupled (13)C NMR spectra of a dipeptide (Glu-Trp) and a tetrapeptide (NAc-Ser-Phe-Val-Gly-OMe) were recorded in D(2)O and in a lyotropic alignment medium (pentaethylene glycol monododecyl ether/n-hexanol). The residual dipolar couplings were extracted as the differences between the observed couplings for the magnetic nuclei dissolved in the latter and former media. Using a computational optimization, the spatial structures of the compounds were calculated starting from their respective low energy conformations obtained on a semiempirical basis. The uniformity of each conformation was confirmed by the solid-state (13)C NMR spectra of powder samples. Differences between the starting structures and final ones, optimized when employing residual dipolar couplings, are discussed.     

Comparison of aqueous molecular dynamics with NMR relaxation and residual dipolar couplings favors internal motion in a mannose oligosaccharide.

An investigation has been performed to assess how aqueous dynamical simulations of flexible molecules can be compared against NMR data. The methodology compares state-of-the-art NMR data (residual dipolar coupling, NOESY, and (13)C relaxation) to molecular dynamics simulations in water over several nanoseconds. In contrast to many previous applications of residual dipolar coupling in structure investigations of biomolecules, the approach described here uses molecular dynamics simulations to provide a dynamic representation of the molecule. A mannose pentasaccharide, alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->2)-D-Manp, was chosen as the model compound for this study. The presence of alpha-linked mannan is common to many glycopeptides, and therefore an understanding of the structure and the dynamics of this molecule is of both chemical and biological importance. This paper sets out to address the following questions. (1) Are the single structures which have been used to interpret residual dipolar couplings a useful representation of this molecule? (2) If dynamic flexibility is included in a representation of the molecule, can relaxation and residual dipolar coupling data then be simultaneously satisfied? (3) Do aqueous molecular dynamics simulations provide a reasonable representation of the dynamics present in the molecule and its interaction with water? In summary, two aqueous molecular dynamics simulations, each of 20 ns, were computed. They were started from two distant conformations and both converged to one flexible ensemble. The measured residual dipolar couplings were in agreement with predictions made by averaging the whole ensemble and from a specific single structure selected from the ensemble. However, the inclusion of internal motion was necessary to rationalize the relaxation data. Therefore, it is proposed that although residual dipolar couplings can be interpreted as a single-structure, this may not be a correct interpretation of molecular conformation in light of other experimental data. Second, the methodology described here shows that the ensembles from aqueous molecular dynamics can be effectively tested against experimental data sets. In the simulation, significant conformational motion was observed at each of the linkages, and no evidence for intramolecular hydrogen bonds at either alpha(1-->2) or alpha(1-->3) linkages was found. This is in contrast to simulations of other linkages, such as beta(1-->4), which are often predicted to maintain intramolecular hydrogen bonds and are coincidentally predicted to have less conformational freedom in solution.     

Complete relative stereochemistry of multiple stereocenters using only residual dipolar couplings

Residual dipolar couplings (RDCs), in combination with molecular order matrix calculations, were used to unambiguously determine the complete relative stereochemistry of an organic compound with five stereocenters. Three simple one-dimensional experiments were utilized for the measurements of (13)C-(1)H, (13)C-(19)F, (19)F-(1)H, and (1)H-(1)H RDCs. The order matrix calculation was performed on each chiral isomer independently. The fits were evaluated by the comparison of the root-mean-square deviation (rmsd) of calculated and measured RDCs. The order tensor simulations based on two different sets of RDC data collected with phage and bicelles are consistent. The resulting stereochemical assignments of the stereocenters obtained from using only RDCs are in perfect agreement with those obtained from the single-crystal X-ray structure. Six RDCs are found to be necessary to run the simulation, and seven are the minimum to get an acceptable result for the investigated compound. It was also shown that (13)C-(1)H and (1)H-(1)H RDCs, which are the easiest to measure, are also the most important and information-rich data for the order matrix calculation. The effect of each RDC on the calculation depends on the location of the corresponding vector in the structure. The direct RDC of a stereocenter is important to the configuration determination, but the configuration of stereocenters devoid of protons can also be obtained from analysis of nearby RDCs.     

Protein structural motif recognition via NMR residual dipolar couplings.

NMR residual dipolar couplings have great potential to provide rapid structural information for proteins in the solution state. This information even at low resolution may be used to advantage in proteomics projects that seek to annotate large numbers of gene products for entire genomes. In this paper, we describe a novel approach to the structural interpretation of dipolar couplings which is based on structural motif pattern recognition, where a predefined gapless structural template for a motif is used to search a set of residual dipolar couplings for good matches. We demonstrate the applicability of the method using synthetic and experimental data. We also provide an analysis of the statistical power of the method and the effects of order tensor frame orientation, motif size, and structural complexity on motif detection. Finally, we discuss remaining problems that must be overcome before the method can be used routinely to identify protein homologies.     

Application of residual dipolar couplings in organic compounds

Residual dipolar couplings (RDCs) induced by anisotropic media are a powerful tool for the structure determination of biomolecules through NMR spectroscopy. Recent advances have proven it to be a valuable tool for determination of the stereochemistry of organic molecules. By simple inspection or order matrix calculations, RDCs provide unambiguous information about the relative configurations or complete stereochemistry of organic compounds.     

Probing the diastereotopicity of methylene protons in strychnine using residual dipolar couplings

Residual dipolar couplings were successfully used to distinguish between the two diastereotopic protons on C-20 of strychnine dissolved in an organic liquid crystal (PBLG/CDCl(3)). The results presented here strongly suggest that this method will be of help in organic structure determination, making the determination of relative stereochemistry in the absence of NOE data possible.     

Calculation of residual dipolar couplings from disordered state ensembles using local alignment

Residual dipolar couplings (RDCs) have been observed in disordered states of several proteins. While their nonuniform values were initially surprising, it has been shown that reasonable approximation of experimental RDCs can be obtained using simple statistical coil models and assuming global alignment of each structure, provided that many thousands of conformers are averaged. Here we show that, by using short local alignment tensors, we can achieve good agreement between experimental and simulated RDCs with far fewer structures than required when using global alignment. This makes the possibility of using RDCs as direct restraints in structural calculations of disordered proteins much more feasible. In addition, it provides insight into the nature of RDCs in disordered states, suggesting that they are primarily reporting on local structure.     

De novo determination of protein backbone structure from residual dipolar couplings using Rosetta

As genome-sequencing projects rapidly increase the database of protein sequences, the gap between known sequences and known structures continues to grow exponentially, increasing the demand to accelerate structure determination methods. Residual dipolar couplings (RDCs) are an attractive source of experimental restraints for NMR structure determination, particularly rapid, high-throughput methods, because they yield both local and long-range orientational information and can be easily measured and assigned once the backbone resonances of a protein have been assigned. While very extensive RDC data sets have been used to determine the structure of ubiquitin, it is unclear to what extent such methods will generalize to larger proteins with less complete data sets. Here we incorporate experimental RDC restraints into Rosetta, an ab initio structure prediction method, and demonstrate that the combined algorithm provides a general method for de novo determination of a variety of protein folds from RDC data. Backbone structures for multiple proteins up to approximately 125 residues in length and spanning a range of topological complexities are rapidly and reproducibly generated using data sets that are insufficient in isolation to uniquely determine the protein fold de novo, although ambiguities and errors are observed for proteins with symmetry about an axis of the alignment tensor. The models generated are not high-resolution structures completely defined by experimental data but are sufficiently accurate to accelerate traditional high-resolution NMR structure determination and provide structure-based functional insights.     

Glycosidic torsional motions in a bicelle-associated disaccharide from residual dipolar couplings

An analysis of torsional motions about glycosidic bonds in a disaccharide is undertaken using residual dipolar coupling measurements and selective immobilization of the reducing end sugar to provide a suitable motional reference. The immobilization is accomplished by using the short chain of an alkyl glycoside to anchor the disaccharide to a bilayer medium aligned in magnetic field. Motions about the beta-(1-4) linkage of the n-butyl-4-O-beta-d-galactopyranosyl-alpha-d-mannopyranoside are shown to be substantial (+/-40 degrees ) and in good agreement with predictions of a fully solvated molecular dynamics simulation.     

Periodicity in residual dipolar couplings and nucleic acid structures

The periodicity in nucleic acid duplex structures is shown to be correlated to the periodicity in residual dipolar couplings (RDCs) in the form of an "RDC wave". This "RDC wave" is characteristic of the alignment of the duplex in the magnetic field, and hence fitting of the data allows the duplex global orientation (, Phi) to be extracted. Further, because the "RDC wave" is fit as a data set of a corresponding secondary structure element, the degeneracy problem is greatly reduced. Consequently, with the global orientation (, Phi) determined, local bond vector conformations are defined. The fit is demonstrated in the examples of the imino RDCs of the negative regulator of splicing RNA fragment (NRS23) and for the C1'H1' RDCs of the Dickerson dodecamer.     

On the origin of residual dipolar couplings from denatured proteins

Effects of steric obstruction on random flight chains are examined. Spatial probability distributions are elaborated to calculate residual dipolar couplings and residual chemical shift anisotropy, parameters that are acquired by NMR spectroscopy from solutes dissolved in dilute liquid crystals. Calculations yield chain length and residue position-dependent values in good agreement with simulations to provide understanding of recently acquired data from denatured proteins.     

Physical interpretation of residual dipolar couplings in neutral aligned media

A novel method is described for rapidly calculating alignment tensors from hydrodynamic shape, required for the prediction of residual dipolar couplings in neutral aligned media. Simulations of alignment were used to show that for steric restriction at a planar surface, the alignment process is dependent on linear hydrodynamic length. However, as discussed, previous methods are not in agreement with this observation. Therefore, the method presented here is the first to provide simple, accurate predictions of the alignment tensor for neutral and dilute media, while being consistent with simulations of alignment. It provides predictions in a fraction of the time of a simulation approach, while aiding physical intuition by providing a direct link between shape and alignment. Not only is this physically gratifying, but it also permits residual dipolar couplings to be applied in demanding situations where simulations of alignment are not desirable, such as in studies of molecular dynamics.     

Simultaneous definition of high resolution protein structure and backbone conformational dynamics using NMR residual dipolar couplings.

Despite the importance of molecular dynamics for biological activity, most approaches to protein structure determination, whether based on crystallographic or solution studies, propose three-dimensional atomic representations of a single configuration that take no account of conformational fluctuation. Non-averaged anisotropic NMR interactions, such as residual dipolar couplings, that become measurable under conditions of weak alignment, provide sensitive probes of both molecular structure and dynamics. Residual dipolar couplings are becoming increasingly powerful for the study of proteins in solution. In this minireview we present their use for the simultaneous determination of protein structure and dynamics.     

Configuration assignment in small organic molecules via residual dipolar couplings

Here we propose a new method to assign relative configurations of stereocenters in small organic molecules by using residual dipolar couplings; the main advantage of this method is that spatial proximity of the stereocenters is not required.