Structural characterization of a mannose-binding protein-trimannoside complex using residual dipolar couplings

The ligand-binding properties of a 53 kDa homomultimeric trimer from mannose-binding protein (MBP) have been investigated using residual dipolar couplings (RDCs) that are easily measured from NMR spectra of the ligand and isotopically labeled protein. Using a limited set of 1H-15N backbone amide NMR assignments for MBP and orientational information derived from the RDC measurements in aligned media, an order tensor for MBP has been determined that is consistent with symmetry-based predictions of an axially symmetric system. 13C-1H couplings for a bound trisaccharide ligand, methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside (trimannoside) have been determined at natural abundance and used as orientational constraints. The bound ligand geometry and orientational constraints allowed docking of the trimannoside ligand in the binding site of MBP to produce a structural model for MBP-oligosaccharide interactions.     

A new approach for applying residual dipolar couplings as restraints in structure elucidation

Residual dipolar couplings are useful global structural restraints. The dipolar couplings define the orientation of a vector with respect to the alignment tensor. Although the size of the alignment tensor can be derived from the distribution of the experimental dipolar couplings, its orientation with respect to the coordinate system of the molecule is unknown at the beginning of structure determination. This causes convergence problems in the simulated annealing process. We therefore propose a protocol that translates dipolar couplings into intervector projection angles, which are independent of the orientation of the alignment tensor with respect to the molecule. These restraints can be used during the whole simulated annealing protocol.     

Conformational studies of blood group A and blood group B oligosaccharides using NMR residual dipolar couplings

The conformations of two synthetic trisaccharides of blood group A and B (alpha-L-Fucp-(1-->2)-[alpha-D-GalpNAc-(1-->3)]-alpha-D-Galp and alpha-L-Fucp-(1-->2)-[alpha-D-Galp-(1-->3)]-alpha-D-Galp, respectively) and of a type A tetrasaccharide alditol, Fucp-(1-->2)-[alpha-D-GalpNAc-(1-->3)]-beta-D-Galp-(1-->3)-GalNAc-ol, were studied by NMR measurements of one-bond C-H residual dipolar couplings in partially oriented liquid crystal solutions. The conformations of the three oligosaccharides were analyzed by generating thousands of structures using a Monte-Carlo method. Two different strategies were applied to calculate theoretical dipolar couplings for these structures. In the first method, the orientation of the molecule was calculated from the optimal fit of the molecular model to the experimental data, while in the second method the orientation tensor was calculated directly from the moment of inertia of the molecular model. Both methods of analysis give similar results but with slightly better agreement with experiment for the former one. The analysis of the results implies a single unique conformation for both blood group epitopes in solution in disagreement with theoretical models suggesting the existence of two conformers in solution.     

Determination of residual dipolar couplings in homonuclear MOCCA-SIAM experiments

In solutions with partial molecular alignment, anisotropic magnetic interactions such as the chemical shift anisotropy, the electric quadrupole interaction, and the magnetic dipole-dipole interaction are no longer averaged out to zero in contrast to isotropic solutions. The resulting residual anisotropic magnetic interactions are increasingly used in biological NMR studies for the determination of 3D structures of proteins and other biomolecules. In the present paper we propose a new approach allowing the measurement of residual H^N-H dipolar couplings of non-isotope enriched proteins based on the application of the MOCCA-SIAM experiment. This experiment allows the measurement of homonuclear coupling constants with an accuracy of ca. ±0.2 Hz and is therefore particularly well suited to determine residual dipolar couplings at relatively low degrees of molecular orientation. The agreement between experimentally determined residual H^N-H couplings and calculated values is demonstrated for BPTI.     

Theoretical framework for NMR residual dipolar couplings in unfolded proteins

A theoretical framework for the prediction of nuclear magnetic resonance (NMR) residual dipolar couplings (RDCs) in unfolded proteins under weakly aligning conditions is presented. The unfolded polypeptide chain is modeled as a random flight chain while the alignment medium is represented by a set of regularly arranged obstacles. For the case of bicelles oriented perpendicular to the magnetic field, a closed-form analytical result is derived. With the obtained analytical expression the RDCs are readily accessible for any locus along the chain, for chains of differing length, and for varying bicelle concentrations. The two general features predicted by the model are (i) RDCs in the center segments of a polypeptide chain are larger than RDCs in the end segments, resulting in a bell-shaped sequential distribution of RDCs, and (ii) couplings are larger for shorter chains than for longer chains at a given bicelle concentration. Experimental data available from the literature confirm the first prediction of the model, providing a tool for recognizing fully unfolded polypeptide chains. With less certainty experimental data appear to support the second prediction as well. However, more systematic experimental studies are needed in order to validate or disprove the predictions of the model. The presented framework is an important step towards a solid theoretical foundation for the analysis of experimentally measured RDCs in unfolded proteins in the case of alignment media such as polyacrylamide gels and neutral bicelle systems which align biomacromolecules by a steric mechanism. Various improvements and generalizations are possible within the suggested approach.     

Refinement of d(GCGAAGC) hairpin structure using one- and two-bond residual dipolar couplings

The structure of the ¹³C,¹⁵N-labeled d(GCGAAGC) hairpin, as determined by NMR spectroscopy and refined using molecular dynamics with NOE-derived distances, torsion angles, and residual dipolar couplings (RDCs), is presented. Although the studied molecule is of small size, it is demonstrated that the incorporation of diminutive RDCs can significantly improve local structure determination of regions undefined by the conventional restraints. Very good correlation between the experimental and back-calculated small one- and two-bond ¹H-¹³C, ¹H-¹⁵N, ¹³C-¹³C and ¹³C-¹⁵N coupling constants has been attained. The final structures clearly show typical features of the miniloop architecture. The structure is discussed in context of the extraordinary stability of the d(GCGAAGC) hairpin, which originates from a complex interplay between the aromatic base stacking and hydrogen bonding interactions.     

Conformational studies of Lewis X and Lewis A trisaccharides using NMR residual dipolar couplings.

The conformations of the histo-blood group carbohydrate antigens Lewis X (Le(x)) and Lewis A (Le(a)) were studied by NMR measurements of one-bond C-H residual dipolar couplings in partially oriented liquid crystal solutions. A strategy for rapid calculation of the difference between theoretical and experimental dipolar couplings of a large number of model structures generated by computer simulations was developed, resulting in an accurate model structure for the compounds. Monte Carlo simulations were used to generate models for the trisaccharides, and orientations of each model were sought that could reproduce the experimental residual dipolar coupling values. For both, Le(a) and Le(x), single low energy models giving excellent agreement with experiment were found, implying a compact rigidly folded conformation for both trisaccharides. The new approach was also applied to the pentasaccharides lacto-N-fucopentaose 2 (LNF-2) and lacto-N-fucopentaose 3 (LNF-3) proving its consistency and robustness. For describing the conformation of tightly folded oligosaccharides, a definition for characterization of ring planes in pyranoside chairs is proposed and applied to the analysis of the relation between the fucose and galactose residues in the epitopes, revealing the structural similarity between them.     

Assessment of molecular structure using frame-independent orientational restraints derived from residual dipolar couplings

Residual dipolar couplings measured in weakly aligning liquid-crystalline solvent contain valuable information on the structure of biomolecules in solution. Here we demonstrate that dipolar couplings (DCs) can be used to derive a comprehensive set of pairwise angular restraints that do not depend on the orientation of the alignment tensor principal axes. These restraints can be used to assess the agreement between a trial protein structure and a set of experimental dipolar couplings by means of a graphic representation termed a `DC consistency map'. Importantly, these maps can be used to recognize structural elements consistent with the experimental DC data and to identify structural parameters that require further refinement, which could prove important for the success of DC-based structure calculations. This approach is illustrated for the 42 kDa maltodextrin-binding protein.     

Evaluation of the influence of anisotropic indirect nuclear spin-spin coupling tensors on effective residual dipolar couplings for model peptides

Residual dipolar couplings (RDCs) observed between nuclear spins in molecules in partially oriented media have become a valuable source of information for NMR spectroscopists seeking to structurally characterize biological macromolecules. Examination of the form of the direct (D) and indirect (J) nuclear spin-spin coupling Hamiltonians indicates that all observed RDCs contain an unknown contribution from the anisotropic part of J (J) in addition to the direct dipolar contribution, D _PQ. Here, we evaluate the influence of J on RDCs through a series of DFT calculations on model peptides. Very small corrections to one-bond RDCs measured between heavy atoms in peptides and proteins are recommended: +0.51% for N-C spin pairs, and +0.45% for C-C spin pairs. The corrections to RDCs involving at least one proton are negligible. This latter point is likely to be equally applicable to nucleic acids and oligosaccharides in addition to peptides and proteins. Finally, the orientations of the J(N, C) and J(C, C) tensors in the molecular framework are reported for glycylglycine.     

A novel interactive tool for rigid-body modeling of multi-domain macromolecules using residual dipolar couplings

Residual dipolar couplings (RDC), measured by dissolving proteins in dilute liquid crystal media, or by studying naturally paramagnetic molecules, have rapidly become established as routine measurements in the investigation of the structure of macromolecules by NMR. One of the most obvious applications of the previously inaccessible long-range angular information afforded by RDC is the accurate definition of domain orientation in multi-module macromolecules or complexes. In this paper we describe a novel program developed to allow the determination of alignment tensor parameters for individual or multiple domains in macromolecules from residual dipolar couplings and to facilitate their manipulation to construct low-resolution models of macromolecular structure. For multi-domain systems the program determines the relative orientation of individual structured domains, and provides graphical user-driven rigid-body modeling of the different modules relative to the common tensorial frame. Translational freedom in the common frame, and equivalent rotations about the diagonalized (x,y,z) axes are used to position the different modules in the common frame to find a model in best agreement with experimentally measured couplings alone or in combination with additional experimental or covalent information.     

Exact solutions for chemical bond orientations from residual dipolar couplings

New methods for determining chemical structures from residual dipolar couplings are presented. The fundamental dipolar coupling equation is converted to an elliptical equation in the principal alignment frame. This elliptical equation is then combined with other angular or dipolar coupling constraints to form simple polynomial equations that define discrete solutions for the unit vector(s). The methods are illustrated with residual dipolar coupling data on ubiquitin taken in a single anisotropic medium. The protein backbone is divided into its rigid groups (namely, its peptide planes and C frames), which may be solved for independently. A simple procedure for recombining these independent solutions results in backbone dihedral angles and that resemble those of the known native structure. Subsequent refinement of these - angles by the ROSETTA program produces a structure of ubiquitin that agrees with the known native structure to 1.1 Å C rmsd.     

Refinement of the solution structure of the heparin-binding domain of vascular endothelial growth factor using residual dipolar couplings

Previous NMR structural studies of the heparin-binding domain of vascular endothelial growth factor (VEGF₁₆₅) revealed a novel fold comprising two subdomains, each containing two disulfide bridges and a short two-stranded antiparallel -sheet. The mutual orientation of the two subdomains was poorly defined by the NMR data. Heteronuclear relaxation data suggested that this disorder resulted from a relative lack of experimental restraints due to the limited size of the interface, rather than inherent high-frequency flexibility. Refinement of the structure using ¹H^N-¹⁵N residual dipolar coupling restraints results in significantly improved definition of the relative subdomain orientations.     

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.     

Determination of protein global folds using backbone residual dipolar coupling and long-range NOE restraints

We report the determination of the global fold of human ubiquitin using protein backbone NMR residual dipolar coupling and long-range nuclear Overhauser effect (NOE) data as conformational restraints. Specifically, by use of a maximum of three backbone residual dipolar couplings per residue (N_i-H^N _i, N_i-C_i–1, H^N _i - C_i–1) in two tensor frames and only backbone H^N-H^N NOEs, a global fold of ubiquitin can be derived with a backbone root-mean-square deviation of 1.4 Å with respect to the crystal structure. This degree of accuracy is more than adequate for use in databases of structural motifs, and suggests a general approach for the determination of protein global folds using conformational restraints derived only from backbone atoms.     

NMR experiments for the measurement of proton-proton and carbon-carbon residual dipolar couplings in uniformly labelled oligosaccharides

A 2D-HSQC-carbon selective/proton selective-constant time COSY, 2D-HSQC-(sel C, sel H)-CT COSY experiment, which is applicable to uniformly ¹³C isotopically enriched samples (U-¹³C) of oligosaccharides or oligonucleotides is proposed for the measurement of proton–proton RDC in crowded regions of 2D-spectra. In addition, a heteronuclear constant time-COSY experiment, ¹³C-¹³C CT-COSY, is proposed for the measurement of one bond carbon–carbon RDC in these molecules. These two methods provide an extension, to U-¹³C molecules, of the original homonuclear constant time-COSY experiment proposed by Tian et al. (1999) for saccharides. The combination of a number of these RDC with NOE data may provide the method of choice to study oligosaccharide conformation in the free and receptor-bound state.     

Characterization of magnetically oriented phospholipid micelles for measurement of dipolar couplings in macromolecules

Weak alignment of solute molecules with the magnetic field can be achieved in a dilute liquid crystalline medium, consisting of an aqueous mixture of dimyristoyl-phosphatidylcholine (DMPC) and dihexanoyl-phosphatidylcholine (DHPC). For a certain range of molar ratios, DMPC and DHPC can form large, disc-shaped particles, commonly referred to as bicelles (Sanders and Schwonek, 1992), which cooperatively align in the magnetic field and induce a small degree of alignment on asymmetrically shaped solute molecules. As a result, dipolar couplings between pairs of ¹H, ¹³C or ¹⁵N nuclei are no longer averaged to zero by rotational diffusion and they can be readily measured, providing valuable structural information. The stability of these liquid crystals and the degree of alignment of the solute molecules depend strongly on experimental variables such as the DMPC:DHPC ratio and concentration, the preparation protocol of the DMPC/DHPC mixtures, as well as salt, temperature, and pH. The lower temperature limit for which the liquid crystalline phase is stable can be reduced to 20 °C by using a ternary mixture of DHPC, DMPC, and 1-myristoyl-2-myristoleoyl-sn-glycero-3-phosphocholine, or a binary mixture of DHPC and ditridecanoyl-phosphatidylcholine. These issues are discussed, with an emphasis on the use of the medium for obtaining weak alignment of biological macromolecules.     

Measurements of residual dipolar couplings in peptide inhibitors weakly aligned by transient binding to peptide amyloid fibrils**

In this communication, we suggest that transferred residual dipolar couplings (trRDCs) can be employed to restrain the structure of peptide inhibitors transiently binding to beta-amyloid fibrils. The effect is based on the spontaneous alignment of amyloid fibrils with the fibril axis parallel to the magnetic field. This alignment is transferred to the transiently binding peptide inhibitor and is reflected in the size of the trRDCs. We find that the peptide inhibitor adopts a beta-sheet conformation with the backbone N-H and C-H dipolar vectors aligned preferentially parallel and perpendicular, respectively, to the fibril axis.     

Structure and orientation of a G protein fragment in the receptor bound state from residual dipolar couplings

Residual dipolar couplings for a ligand that is in fast exchange between a free state and a state where it is bound to a macroscopically ordered membrane protein carry precise information on the structure and orientation of the bound ligand. The couplings originate in the bound state but can be detected on the free ligand using standard high resolution NMR. This approach is used to study an analog of the C-terminal undecapeptide of the alpha-subunit of the heterotrimeric G protein transducin when bound to photo-activated rhodopsin. Rhodopsin is the major constituent of disk-shaped membrane vesicles from rod outer segments of bovine retinas, which align spontaneously in the NMR magnet. Photo-activation of rhodopsin triggers transient binding of the peptide, resulting in measurable dipolar contributions to 1J(NH) and 1J(CH) splittings. These dipolar couplings report on the time-averaged orientation of bond vectors in the bound peptide relative to the magnetic field, i.e. relative to the membrane normal. Approximate distance restraints of the bound conformation were derived from transferred NOEs, as measured from the difference of NOESY spectra recorded prior to and after photo-activation. The N-terminal eight residues of the bound undecapeptide adopt a near-ideal alpha-helical conformation. The helix is terminated by an alpha(L) type C-cap, with Gly9 at the C' position in the center of the reverse turn. The angle between the helix axis and the membrane normal is 40 degrees (+/-4) degrees. Peptide protons that make close contact with the receptor are identified by analysis of the NOESY cross-relaxation pattern and include the hydrophobic C terminus of the peptide.     

Quantitative conformational analysis of the core region of N-glycans using residual dipolar couplings, aqueous molecular dynamics, and steric alignment

A method is described for quantitatively investigating the dynamic conformation of small oligosaccharides containing an (16) linkage. It was applied to the oligosaccharide Man-(13) {Man- (16)}Man--O-Me, which is a core region frequently observed in N-linked glycans. The approach tests an aqueous molecular dynamics simulation, capable of predicting microscopic dynamics, against experimental residual dipolar couplings, by assuming that alignment is caused purely by steric hindrance. The experimental constraints were heteronuclear and homonuclear residual dipolar couplings, and in particular those within the (16) linkage itself. Powerful spin-state-selective pulse sequences and editing schemes were used to obtain the most relevant couplings for testing the model. Molecular dynamics simulations in water over a period of 50 ns were not able to predict the correct rotamer population at the (16) linkage to agree with the experimental data. However, this sampling problem could be corrected using a simple maximum likelihood optimisation, indicating that the simulation was modelling local dynamics correctly. The maximum likelihood prediction of the residual dipolar couplings was found to be an almost equal population of the gg and gt rotamer conformations at the (16) linkage, and the tg conformation was predicted to be unstable and unpopulated in aqueous solution. In this case all twelve measured residual dipolar couplings could be satisfied. This conformer population could also be used to make predictions of scalar couplings with the use of a previously derived empirical equation, and is qualitatively in agreement with previous predictions based on NMR, X-ray crystallography and optical data.     

Transverse relaxation optimised spin-state selective NMR experiments for measurement of residual dipolar couplings

Three transverse relaxation optimised NMR experiments (TROSY) for the measurement of scalar and dipolar couplings suitable for proteins dissolved in aqueous iso- and anisotropic solutions are described. The triple-spin-state-selective experiments yield couplings between ¹H^N-¹³C, ¹⁵N-¹³C, ¹H^N-¹³C _i–1, ¹⁵N-¹³C _i–1, ¹H^N-¹³C_i–1, ¹⁵N-¹³C_i–1, and ¹³C_i–1-¹³C _i–1 without introducing nonessential spectral crowding compared with an ordinary two-dimensional ¹⁵N-¹H correlation spectrum and without requiring explicit knowledge of carbon assignments. This set of /-J-TROSY experiments is most useful for perdeuterated proteins in studies of structure–activity relationships by NMR to observe, in addition to epitopes for ligands, also conformational changes induced by binding of ligands.