Conformational analysis of alpha-D-Fuc-(1-->4)-beta-D-GlcNAc-OMe. One-dimensional transient NOE experiments and Metropolis Monte Carlo simulations

Context Effects on Human Milk Oligosaccharide Linkage Conformation and Dynamics Revealed by MA'AT Analysis

An emerging NMR method, MA'AT analysis, has been applied to investigate context effects on the conformational properties of several human milk oligosaccharides (HMOs). The MA'AT model of the β-(1→4) linkage in the disaccharide, methyl β-lactoside (MeL), was compared to those obtained for the same linkage in the HMO trisaccharides, methyl 2'-fucosyllactoside (Me2'FL) and methyl 3-fucosyllactoside (Me3FL), and in the tetrasaccharide, methyl 2',3-difucosyllactoside (Me2',3DFL). MA'AT analysis revealed significant context effects on the mean values and circular standard deviations (CSDs) of the psi (ψ) torsion angles in these linkages. α-Fucosylation at both O2'Gal and O3Glc of MeL to give Me2',3DFL significantly constrained librational motion about ψ (70% reduction in the CSD) and shifted its mean value by ∼18°. α-Fucosylation at the O3Glc of MeL to give Me3FL constrained ψ more than α-fucosylation at the O2Gal to give Me2'FL. These effects can be explained by the expected solution conformation of Me3FL, which closely resembles the Lewisx trisaccharide. Comparisons of MA'AT models of ψ to those obtained by 1 μs aqueous molecular dynamics simulation (GLYCAM06) revealed identical trends, that is, MA'AT analysis was able to recapitulate molecular behavior in solution that was heretofore only available from MD simulation. The results highlight the capabilities of MA'AT analysis to determine probability distributions of molecular torsion angles in solution as well as degrees of librational averaging of these angles.

MA'AT analysis of the O-glycosidic linkages of oligosaccharides using nonconventional NMR J-couplings: MA'AT and MD models of phi

MA'AT analysis (Meredith et al., J. Chem. Inf. Model. 2022, 62, 3135-3141) is a new NMR-based method to treat ensembles of redundant NMR spin-coupling constants (J-couplings) to obtain experiment-based probability distributions of molecular torsion angles in solution. Work reported to date on modeling the conformations of O-glycosidic linkages of oligosaccharides using three conventional J-coupling constraints (2 J COC, 3 J COCH, 3 J COCC) has shown that the method gives mean torsion angles and circular standard deviations (CSDs) for psi in very good agreement with those obtained by MD simulation. On the other hand, CSDs for phi determined by MA'AT analysis have consistently been much larger than those determined by MD, calling into question either the reliability of MA'AT analysis or MD to accurately predict this behavior. Prior work has shown that this discrepancy does not stem from the limitations of DFT-based J-coupling equation parameterization where secondary conformational dependencies can introduce uncertainties. The present work re-visits this problem by incorporating a new nonconventional J-coupling constraint into MA'AT analyses of phi, namely, a geminal (two-bond) 2 J CCH J-value that exhibits a strong primary dependence on phi. The latter property pertains explicitly to linkages contributed by GlcNAc pyranosyl rings and pyranosyl rings devoid of substituents at C2 (i.e., deoxy residues) where known secondary contributions to 2 J CCH magnitude caused by C-O bond rotation involving the coupled carbon are negligible or absent. The results show that when 2 J CCH values are added to the analysis, phi CSDs reduce considerably, bringing them into better alignment with those obtained by MD simulation. The cause of the discrepancy when only three conventional J-couplings are used to treat phi appears to be associated with the two-bond 2 J COC, which has properties that make it less effective than the non-conventional 2 J CCH as a discriminator of different conformational models of phi.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Synthesis of novel ammonium and selenonium ions and their evaluation as inhibitors of UDP-galactopyranose mutase

The syntheses of two ammonium salts of 1,4-dideoxy-1,4-imino-d-galactitol containing erythritol sulfate side chains are described. The parent compound is a known inhibitor of the enzyme UDP-galactopyranose mutase (UDP-galactopyranose furanomutase, E.C. 5.4.99.9), which is responsible for the conversion of UDP-galactopyranose into UDP-galactofuranose and is presumably protonated in its active form. The side chain was chosen because it is present in a known sulfonium ion alpha-glucosidase inhibitor, salacinol. The syntheses of the selenonium analogues derived from 1,4-dideoxy-1,4-seleno-d-galactitol are also described. The synthetic strategy in the syntheses of all four salts involved the nucleophilic attack of a protected derivative of the alditol at the least hindered carbon of 2,4-O-benzylidene d- or l-erythritol-1,3-cyclic sulfate to give adducts that were subsequently deprotected. The importance of different protecting groups used in the various syntheses is also highlighted. Enzyme inhibition assays carried out on these compounds, and the corresponding sulfonium ions synthesized previously, show that they are poor inhibitors of UDP-galactopyranose mutase.

Complexes of glucoamylase with maltoside heteroanalogues: bound ligand conformations by use of transferred NOE experiments and molecular modeling

Transferred nuclear Overhauser effect (trNOE) experiments have been performed to investigate the conformations of the competitive inhibitors, methyl 5'-thio-4-N-alpha-maltoside 3a and methyl 5'-thio-4-S-alpha-maltoside 4 when bound to the catalytic subunit of the enzyme glucoamylase. These NMR data suggest that, although each of the free ligands populates two conformational families, both heteroanalogues are bound by the enzyme in conformations in the area of the global energy minimum. These conformations have been used as initial points for docking into the active site of the enzyme taken from a X-ray crystal structure of the related glucoamylase-D-gluco-dihydroacarbose 2 complex. Minimization of the resulting complexes has yielded structures for the bound complexes. Corroboration of the structures is provided by fast T(1)(rho)-relaxation effects for certain ligand protons as a result of close contacts with protons in the enzyme active site. The results auger well for the combined use of transferred NOE spectroscopy and molecular modeling based on X-ray crystal structures of complexes of suitable congeners for the rapid analysis of ligand-receptor interactions.

Computational carbohydrate chemistry: what theoretical methods can tell us

Computational methods have had a long history of application to carbohydrate systems and their development in this regard is discussed. The conformational analysis of carbohydrates differs in several ways from that of other biomolecules. Many glycans appear to exhibit numerous conformations coexisting in solution at room temperature and a conformational analysis of a carbohydrate must address both spatial and temporal properties. When solution nuclear magnetic resonance data are used for comparison, the simulation must give rise to ensemble-averaged properties. In contrast, when comparing to experimental data obtained from crystal structures a simulation of a crystal lattice, rather than of an isolated molecule, is appropriate. Molecular dynamics simulations are well suited for such condensed phase modeling. Interactions between carbohydrates and other biological macromolecules are also amenable to computational approaches. Having obtained a three-dimensional structure of the receptor protein, it is possible to model with accuracy the conformation of the carbohydrate in the complex. An example of the application of free energy perturbation simulations to the prediction of carbohydrate-protein binding energies is presented.

Modeling polysaccharides: present status and challenges

The most recent tools that have been developed for modeling the three-dimensional features of polysaccharides and carbohydrate polymers are presented. The presentation starts with a description of the conformations of the monosaccharides, and of the flexible rings such as in the case of five-membered rings, and a thorough description of the conformational space that is available for a disaccharide unit, either in vacuo or in an aqueous phase. The extension to the modeling of the parent polysaccharides is addressed, based on the assumption that owing to the size and relative rigidity of the intervening monosaccharides units, the rotations at a particular linkage can be, under some conditions, considered as independent of nearest neighbor interactions. Appropriate modeling techniques are described that can provide insights into the dimensions of the chain in a solution which is best described as a random coil accompanied by the occurrence of local "helical" regions. With the help of such descriptors such as helical parameters, the ordered state of polysaccharide strands can be readily characterized. The generation of double or triple helices can be then attempted in order to explore the occurrence of such multi-stranded arrangements that may be energetically stable. The final step in the determination of the structure of polysaccharides in the ordered state, is the investigation of the interactions of different helices. This may lead to either the best arrangement(s) between two polymeric chains, or to the prediction of the dimensions, and the symmetry of a three-dimensional lattice. Some of the tools which have been developed should allow automatic scarches for meaningful correlations between structures and functions, through exploratory data analysis. Structure-function or structure-property correlation could be then used to model changes arising from structural alterations. This would open the field of polysaccharide engineering.

Conformational analysis of blood group A trisaccharide in solution and in the binding site of Dolichos biflorus lectin using transient and transferred nuclear Overhauser enhancement (NOE) and rotating-frame NOE experiments

The present study is concerned with the elucidation of the conformation of the blood group A trisaccharide (alpha-D-GalNAc(1-3)[alpha-L-Fuc(1-->2)] beta-D-Gal-O-R) in the combining site of Dolichos biflorus seed lectin by use of 400-MHz and 600-MHz NMR spectroscopy. D. biflorus lectin displays a unique specificity for GalNAc residues. It occurs in solution as a tetrameric assembly having a molecular mass of 110 kDa, with two carbohydrate-binding sites per molecule. First, NOE build-up curves were obtained for the free blood group A trisaccharide from one-dimensional transient NOE experiments. Simulated NOE build-up curves were constructed from an ensemble of low-energy conformers derived from previous investigations. The comparison of theoretical and experimental data indicates that an equilibrium between two families of low-energy conformers most likely reflects the solution behavior of the trisaccharide in solution. Two-dimensional transferred NOE and rotating-frame enhancements (ROE) were subsequently measured for the trisaccharide complexed with the D. biflorus seed lectin. In addition to the NOEs observed for the free trisaccharide, the transferred NOESY spectrum showed several new NOEs that were identified as spin diffusion using a rotating-frame NOESY (ROESY) experiment. Experimental interglycosidic transferred nuclear Overhauser effect (TRNOE) build-up curves were compared to theoretical curves calculated for both low-energy conformers located in the D. biflorus lectin-binding site. Calculations of theoretical TRNOE were performed using a combination of the full relaxation matrix and the protein-ligand exchange matrix. Comparison between experimental and simulated TRNOE volumes leads to the conclusion that one conformation of blood group A trisaccharide is selected upon binding by D. biflorus lectin.

Three-dimensional structures of oligosaccharides

Oligosaccharides represent a particularly challenging class of molecules for conformational analysis. Recent advances in experimental and theoretical methods have begun to yield further insight into their conformational behavior; however, general rules governing their conformational preferences have not yet emerged. X-ray and NMR techniques may provide vital insights into protein-bound oligosaccharide conformations, but these do not necessarily represent highly populated solution conformations. Moreover, an oligosaccharide's inherent flexibility and lack of strong intermolecular interactions places extreme demands on theoretical methods.