Review: Use of residual dipolar couplings to determine the structure of carbohydrates

Exploring the impact of alignment media on RDC analysis of phosphorus-containing compounds: a molecular docking approach

Residual dipolar couplings (RDCs) are employed in NMR analysis when conventional methods, such as J-couplings and nuclear Overhauser effects (NOEs) fail. Low-energy (optimized) conformers are often used as input structures in RDC analysis programs. However, these low-energy structures do not necessarily resemble conformations found in anisotropic environments due to interactions with the alignment medium, especially if the analyte molecules are flexible. Considering interactions with alignment media in RDC analysis, we developed and evaluated a molecular docking-based approach to generate more accurate conformer ensembles for compounds in the presence of the poly-γ-benzyl-L-glutamate alignment medium. We designed chiral phosphorus-containing compounds that enabled us to utilize 31P NMR parameters for the stereochemical analysis. Using P3D/PALES software to evaluate diastereomer discrimination, we found that our conformer ensembles outperform moderately the standard, low-energy conformers in RDC analysis. To further improve our results, we (i) averaged the experimental values of the molecular docking-based conformers by applying the Boltzmann distribution and (ii) optimized the structures through normal mode relaxation, thereby enhancing the Pearson correlation factor R and even diastereomer discrimination in some cases. Nevertheless, we presume that significant differences between J-couplings in isotropic and in anisotropic environments may preclude RDC measurements for flexible molecules. Therefore, generating conformer ensembles based on molecular docking enhances RDC analysis for mildly flexible systems while flexible molecules may require applying more advanced approaches, in particular approaches including dynamical effects.

Residual Dipolar Coupling Based Conformational Analysis Allows the Configurational Assessment of Steroids with up to Eight Stereocenters

Residual dipolar couplings (RDCs) induced by anisotropic media have been proved as a powerful tool for the structure elucidation of organic molecules in solution in nuclear magnetic resonance (NMR) based analysis. The value of dipolar couplings to solve complex conformational and configurational problems represents indeed an appealing analytical tool for the pharmaceutical industry particularly focusing on the stereochemistry characterization of NCEs since the early phase of the drug development process. In our work, RDCs were used for the conformational and configurational study of synthetic steroids with multiple stereocenters - prednisone and beclomethasone dipropionate (BDP) -. For both molecules the correct relative configuration was identified among all the possible diastereoisomers (32 and 128 respectively) arising from the compounds stereogenic carbons. Only for prednisone the use of additional experimental data (i. e. rOes) was necessary to resolve the right stereochemical structure.

The Flexibility of Oligosaccharides Unveiled Through Residual Dipolar Coupling Analysis

The intrinsic flexibility of glycans complicates the study of their structures and dynamics, which are often important for their biological function. NMR has provided insights into the conformational, dynamic and recognition features of glycans, but suffers from severe chemical shift degeneracy. We employed labelled glycans to explore the conformational behaviour of a β(1-6)-Glc hexasaccharide model through residual dipolar couplings (RDCs). RDC delivered information on the relative orientation of specific residues along the glycan chain and provided experimental clues for the existence of certain geometries. The use of two different aligning media demonstrated the adaptability of flexible oligosaccharide structures to different environments.

Extending the applicability of P3D for structure determination of small molecules

The application of anisotropic nuclear magnetic resonance (NMR) parameters for the correct structural assignment of small molecules requires the use of partially ordered media. Previously we demonstrated that the use of P3D simulations using poly(γ -benzyl-L-glutamate) (PBLG) as an alignment medium allows for the determination of the correct diastereomer from extremely sparse NMR data. Through the analysis of the structural characteristics of small molecules in different alignment media, here we show that when steric or electrostatic factors dominate the alignment, P3D-PBLG retains its diastereomer discrimination power. We also demonstrate that P3D simulations can define the different conformations of a flexible small molecule from sparse NMR data.

Synthesis and O-Glycosidic Linkage Conformational Analysis of 13C-Labeled Oligosaccharide Fragments of an Antifreeze Glycolipid

NMR studies of two ¹³C-labeled disaccharides and a tetrasaccharide were undertaken that comprise the backbone of a novel thermal hysteresis glycolipid containing a linear glycan sequence of alternating [βXyl p-(1→4)-βMan p-(1→4)] _n dimers. Experimental trans-glycoside NMR J-couplings, parameterized equations obtained from density functional theory (DFT) calculations, and an in-house circular statistics package ( MA'AT) were used to derive conformational models of linkage torsion angles ϕ and ψ in solution, which were compared to those obtained from molecular dynamics simulations. Modeling using different probability distribution functions showed that MA'AT models of ϕ in βMan(1→4)βXyl and βXyl(1→4)βMan linkages are very similar in the disaccharide building blocks, whereas MA'AT models of ψ differ. This pattern is conserved in the tetrasaccharide, showing that linkage context does not influence linkage geometry in this linear system. Good agreement was observed between the MA'AT and MD models of ψ with respect to mean values and circular standard deviations. Significant differences were observed for ϕ, indicating that revision of the force-field employed by GLYCAM is probably needed. Incorporation of the experimental models of ϕ and ψ into the backbone of an octasaccharide fragment leads to a helical amphipathic topography that may affect the thermal hysteresis properties of the glycolipid.

Combining pure shift and J-edited spectroscopies: A strategy for extracting chemical shifts and scalar couplings from highly crowded proton spectra of oligomeric saccharides

We report the application of pure shift and J-edited nuclear magnetic resonance spectroscopies to the structural analysis of a protected maltotrioside synthetic intermediate whose crowded ¹ H spectrum displays highly crowded regions. The analytical strategy is based on the implementation of J-edited and TOCSY experiments whose resolution is optimized by the use of broadband homonuclear decoupling and selective refocusing techniques, to assign and measure chemical shifts and homonuclear scalar couplings with high accuracy. The resulting data show a high level of complementarity, providing a detailed insight into each subunit of this oligomeric saccharide, even for proton sites whose nuclear magnetic resonance signals strongly overlap. This approach allowed for fully assigning proton chemical shifts and extracting 80% of the ³ J_HH couplings that are in excellent agreement with those expected for D-gluco-pyranosyl units in ⁴ C₁ conformations.

Residual Dipolar Couplings in Structure Determination of Natural Products

The determination of natural products stereochemistry remains a formidable task. Residual dipolar couplings (RDCs) induced by anisotropic media are a powerful tool for determination of the stereochemistry of organic molecule in solution. This review will provide a short introduction on RDCs-based methodology for the structural elucidation of natural products. Special attention is given to the current availability of alignment media in organic solvents. The applications of RDCs for structural analysis of some examples of natural products were discussed and summarized. This review provides a short introduction on RDCs-based methodology for the structural elucidation of natural products. Special attention is given to the current availability of alignment media in organic solvents. The applications of RDCs for structural analysis of some examples of natural products were discussed and summarized.

The β-reducing end in α(2-8)-polysialic acid constitutes a unique structural motif

Over the years, structural characterizations of α(2-8)-polysialic acid (polySia) in solution have produced inconclusive results. Efforts for obtaining detailed information in this important antigen have focused primarily on the α-linked residues and not on the distinctive characteristics of the terminal ones. The thermodynamically preferred anomeric configuration for the reducing end of sialic acids is β, which has the [I]CO2- group equatorial and the OH ([I]OH2) axial, while for all other residues the CO2- group is axial. We show that this purportedly minor difference has distinct consequences for the structure of α(2-8)-polySia near the reducing end, as the β configuration places the [I]OH2 in a favorable position for the formation of a hydrogen bond with the carboxylate group of the following residue ([II]CO2-). Molecular dynamics (MD) simulations predicted the hydrogen bond, which we subsequently directly detected by NMR. The combination of MD and residual dipolar couplings shows that the net result for the structure of Sia2-βOH is a stable conformation with well-defined hydration and charge patterns, and consistent with experimental NOE-based hydroxyl and aliphatic inter-proton distances. Moreover, we provide evidence that this distinct conformation is preserved on Sia oligosaccharides, thus constituting a motif that determines the structure and dynamics of α(2-8)-polySia for at least the first two residues of the polymer. We suggest the hypothesis that this structural motif sheds light on a longtime puzzling observation for the requirement of 10 residues of α(2-8)-polySia in order to bind effectively to specific antibodies, about four units more than for analogous cases.

Direct NOE simulation from long MD trajectories

A software package, MD2NOE, is presented which calculates Nuclear Overhauser Effect (NOE) build-up curves directly from molecular dynamics (MD) trajectories. It differs from traditional approaches in that it calculates correlation functions directly from the trajectory instead of extracting inverse sixth power distance terms as an intermediate step in calculating NOEs. This is particularly important for molecules that sample conformational states on a timescale similar to molecular reorientation. The package is tested on sucrose and results are shown to differ in small but significant ways from those calculated using an inverse sixth power assumption. Results are also compared to experiment and found to be in reasonable agreement despite an expected underestimation of water viscosity by the water model selected.

Paramagnetic NMR probes for characterization of the dynamic conformations and interactions of oligosaccharides

Paramagnetism-assisted nuclear magnetic resonance (NMR) techniques have recently been applied to a wide variety of biomolecular systems, using sophisticated immobilization methods to attach paramagnetic probes, such as spin labels and lanthanide-chelating groups, at specific sites of the target biomolecules. This is also true in the field of carbohydrate NMR spectroscopy. NMR analysis of oligosaccharides is often precluded by peak overlap resulting from the lack of variability of local chemical structures, by the insufficiency of conformational restraints from nuclear Overhauser effect (NOE) data due to low proton density, and moreover, by the inherently flexible nature of carbohydrate chains. Paramagnetic probes attached to the reducing ends of oligosaccharides cause paramagnetic relaxation enhancements (PREs) and/or pseudocontact shifts (PCSs) resolve the peak overlap problem. These spectral perturbations can be sources of long-range atomic distance information, which complements the local conformational information derived from J couplings and NOEs. Furthermore, paramagnetic NMR approaches, in conjunction with computational methods, have opened up possibilities for the description of dynamic conformational ensembles of oligosaccharides in solution. Several applications of paramagnetic NMR techniques are presented to demonstrate their utility for characterizing the conformational dynamics of oligosaccharides and for probing the carbohydrate-recognition modes of proteins. These techniques can be applied to the characterization of transient, non-stoichiometric interactions and will contribute to the visualization of dynamic biomolecular processes involving sugar chains.

A tensor-free method for the structural and dynamical refinement of proteins using residual dipolar couplings

Residual dipolar couplings (RDCs) are parameters measured in nuclear magnetic resonance spectroscopy that can provide exquisitely detailed information about the structure and dynamics of biological macromolecules. We describe here a method of using RDCs for the structural and dynamical refinement of proteins that is based on the observation that the RDC between two atomic nuclei depends directly on the angle ϑ between the internuclear vector and the external magnetic field. For every pair of nuclei for which an RDC is available experimentally, we introduce a structural restraint to minimize the deviation from the value of the angle ϑ derived from the measured RDC and that calculated in the refinement protocol. As each restraint involves only the calculation of the angle ϑ of the corresponding internuclear vector, the method does not require the definition of an overall alignment tensor to describe the preferred orientation of the protein with respect to the alignment medium. Application to the case of ubiquitin demonstrates that this method enables an accurate refinement of the structure and dynamics of this protein to be obtained.

Detecting low-level flexibility using residual dipolar couplings: a study of the conformation of cellobiose

We have developed novel NMR methods for the measurement of heteronuclear residual dipolar couplings (RDCs) in molecules with severely overlapping NMR resonances. These and other methods enabled us to obtain 31 RDCs for α-D-cellobiose and 24 RDCs for β-D-cellobiose. The interpretation of the data in the approximation of a rigid disaccharide structure, using RDCs and interglycosidic (3)J coupling constants, yielded conformation that is very close to that determined using X-ray crystallography. However, depending on which ring was used to calculate the order parameters, the dihedral angle ψH varied up to 30° or 40°, while the φH angle was always the same. This indicates residual flexibility of the glycosidic linkage between the two monosaccharide rings and was observed for both α- and β-D-cellobiose. The RDC analysis using rigid fragments rather than a complete molecule has thus shown that the glycosidic bond of cellobiose is not completely rigid and exhibits low-level flexibility. The sources of this flexibility are discussed and evidence presented to support a hypothesis that it is associated with the ψ more than the φ angle.

NMR of glycans: shedding new light on old problems

The diversity in molecular arrangements and dynamics displayed by glycans renders traditional NMR strategies, employed for proteins and nucleic acids, insufficient. Because of the unique properties of glycans, structural studies often require the adoption of a different repertoire of tailor-made experiments and protocols. We present an account of recent developments in NMR techniques that will deepen our understanding of structure-function relations in glycans. We open with a survey and comparison of methods utilized to determine the structure of proteins, nucleic acids and carbohydrates. Next, we discuss the structural information obtained from traditional NMR techniques like chemical shifts, NOEs/ROEs, and coupling-constants, along with the limitations imposed by the unique intrinsic characteristics of glycan structure on these approaches: flexibility, range of conformers, signal overlap, and non-first-order scalar (strong) coupling. Novel experiments taking advantage of isotopic labeling are presented as an option for overcoming spectral overlap and raising sensitivity. Computational tools used to explore conformational averaging in conjunction with NMR parameters are described. In addition, recent developments in hydroxyl detection and hydrogen bond detection in protonated solvents, in contrast to traditional sample preparations in D2O for carbohydrates, further increase the tools available for both structure information and chemical shift assignments. We also include previously unpublished data in this context. Accurate determination of couplings in carbohydrates has been historically challenging due to the common presence of strong-couplings. We present new strategies proposed for dealing with their influence on NMR signals. We close with a discussion of residual dipolar couplings (RDCs) and the advantages of using (13)C isotope labeling that allows gathering one-bond (13)C-(13)C couplings with a recently improved constant-time COSY technique, in addition to the commonly measured (1)H-(13)C RDCs.

Advances in the REDCAT software package

Background

Residual Dipolar Couplings (RDCs) have emerged in the past two decades as an informative source of experimental restraints for the study of structure and dynamics of biological macromolecules and complexes. The REDCAT software package was previously introduced for the analysis of molecular structures using RDC data. Here we report additional features that have been included in this software package in order to expand the scope of its analyses. We first discuss the features that enhance REDCATs user-friendly nature, such as the integration of a number of analyses into one single operation and enabling convenient examination of a structural ensemble in order to identify the most suitable structure. We then describe the new features which expand the scope of RDC analyses, performing exercises that utilize both synthetic and experimental data to illustrate and evaluate different features with regard to structure refinement and structure validation.

Results

We establish the seamless interaction that takes place between REDCAT, VMD, and Xplor-NIH in demonstrations that utilize our newly developed REDCAT-VMD and XplorGUI interfaces. These modules enable visualization of RDC analysis results on the molecular structure displayed in VMD and refinement of structures with Xplor-NIH, respectively. We also highlight REDCAT’s Error-Analysis feature in reporting the localized fitness of a structure to RDC data, which provides a more effective means of recognizing local structural anomalies. This allows for structurally sound regions of a molecule to be identified, and for any refinement efforts to be focused solely on locally distorted regions.

Conclusions

The newly engineered REDCAT software package, which is available for download via the WWW from http://ifestos.cse.sc.edu, has been developed in the Object Oriented C++ environment. Our most recent enhancements to REDCAT serve to provide a more complete RDC analysis suite, while also accommodating a more user-friendly experience, and will be of great interest to the community of researchers and developers since it hides the complications of software development.