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

NMR Coupling Constants, Karplus Equations, and Adjusted MD Statistics: Detecting Diagnostic Torsion Angles for the Solution Geometry of 6-[α-d-Mannopyranosyl]-d-Mannopyranose (Mannobiose)

The quantitative solution conformations of 2-(hydroxymethyl)-tetrahydropyran, α-methyl-d-mannopyranoside, and 6-[α-d-mannopyranosyl]-d-mannopyranose (mannobiose) are described. Parametrized Karplus equations for redundant spin pairs across the terminal ω-torsion and the glycosidic ω-torsion for mannobiose are developed, including ω/θ-hypersurfaces for the terminal hydroxymethylene group. Experimental NMR data, algorithmic spectral simulation (clustered Hamiltonian method), molecular dynamics (MD) simulations (GLYCAM06), energy minimizations by DFT, and adjusted torsion angle populations weighted over the Karplus-type equations are used. We demonstrate that spectral simulation is a powerful tool in the refinement of initial J values obtained from static GAIO DFT calculations. We also show that only as few as one of multiple redundant torsions can be diagnostic for conformational analysis of the disaccharide.

Expeditious chemical synthesis of xylomannans disproves the proposed antifreeze activities

Cold-adapted species are able to generate cryoprotective proteins and glycoproteins to prevent freezing damage. The [→4)-β-D-Manp-(1→4)-β-D-Xylp-(1→] n xylomannan from the Alaska beetle Upis ceramboides was disclosed by Walters and co-workers in 2009 as the first glycan-based antifreeze agent, which was later reported to be found in diverse taxa. Here, we report the rapid synthesis of four types of xylomannans, including the proposed antifreeze xylomannan up to a 64-mer (Type I), the regioisomeric [→3)-β-D-Manp-(1→4)-β-D-Xylp-(1→] n 16-mer (Type II), the diastereomeric [→4)-β-L-Manp-(1→4)-β-D-Xylp-(1→] n 16-mer (Type III) and the block-wise [→4)-β-D-Manp-(1→] m [→4)-β-D-Xylp-(1→] n 32-mer (Type IV), by employing a strategic iterative exponential glycan growth (IEGG) process. The nuclear magnetic resonance spectral data of the alleged natural xylomannan are in accordance only to those of the block-wise Type IV glycan and none of these synthetic xylomannans has been found to be capable of inducing thermal hysteresis. These results disprove the previous reports about the natural occurrence of antifreeze xylomannans.

Atom-economical synthesis of 1,2-bis(phosphine oxide)ethanes from calcium carbide with straightforward access to deuterium- and 13C-labeled bidentate phosphorus ligands and metal complexes

Straightforward access to bidentate phosphorus ligands and bis(phosphineoxide)ethanes is described based on atom-economic addition reaction. A practical approach was developed to incorporate 2H and 13C labels using easily available reagents.

Synthetic Approaches to Break the Chemical Shift Degeneracy of Glycans

NMR spectroscopy is the leading technique for determining glycans' three-dimensional structure and dynamic in solution as well as a fundamental tool to study protein-glycan interactions. To overcome the severe chemical shift degeneracy of these compounds, synthetic probes carrying NMR-active nuclei (e. g., 13 C or 19 F) or lanthanide tags have been proposed. These elegant strategies permitted to simplify the complex NMR analysis of unlabeled analogues, shining light on glycans' conformational aspects and interaction with proteins. Here, we highlight some key achievements in the synthesis of specifically labeled glycan probes and their contribution towards the fundamental understanding of glycans.

One-Bond 13C-1H and 13C-13C Spin-Coupling Constants as Constraints in MA'AT Analysis of Saccharide Conformation

MA'AT analysis uses ensembles of redundant experimental NMR spin-coupling constants, parametrized J-coupling equations obtained from density functional theory (DFT) calculations, and circular statistics to produce probability distributions of molecular torsion angles in solution and information on librational motions about these angles (Meredith et al., J. Chem. Info. Model. 2022, 62, 3135-3141). Current DFT methods give nearly quantitative two- and three-bond J_HH, J_CH, and ¹J_CC values for use in MA'AT analysis of saccharides. In contrast, the accuracy of DFT-calculated one-bond ¹J_CH and ¹J_CC values is more difficult to determine, preventing their use in MA'AT modeling. This report describes experimental and computational studies that address this problem using two approaches (Strategies 1 and 2). Differences [¹J_CH^calc - ¹J_CH^exp] (Strategy 1) ranged from -1.2 to 2.5 Hz, giving an average difference of 0.8 ± 1.7 Hz. Percent differences ranged from -0.8% to 1.6%, giving an average % difference of 0.5 ± 1.1%. In comparison, [¹J_CH^MA'AT - ¹J_CH^exp] (Strategy 2) ranged from -1.8 to 0.2 Hz, giving an average difference of -1.2 ± 0.7 Hz. Percent differences ranged from -1.2% to 0.1%, giving an average % difference of -0.8 ± 0.5%. Strategy 1 gave an average difference of 2.1 Hz between calculated and experimental ¹J_CC values, with an average % difference of 5.1 ± 0.2%. Calculated ¹J_CC values were consistently larger than experimental values. Strategy 2 also gave calculated ¹J_CC values that were larger than the experimental values, with an average difference of 2.3 ± 0.6 Hz, and an average % difference of 5.6 ± 1.6%. The findings of both strategies are similar and indicate that ¹J_CH values in saccharides can be calculated nearly quantitatively, but ¹J_CC values appear to be consistently overestimated by ∼5% using current DFT methods.

Computational 1 H and 13 C NMR in structural and stereochemical studies

Present review outlines the advances and perspectives of computational ¹ H and ¹³ C NMR applied to the stereochemical studies of inorganic, organic, and bioorganic compounds, involving in particular natural products, carbohydrates, and carbonium ions. The first part of the review briefly outlines theoretical background of the modern computational methods applied to the calculation of chemical shifts and spin-spin coupling constants at the DFT and the non-empirical levels. The second part of the review deals with the achievements of the computational ¹ H and ¹³ C NMR in the stereochemical investigation of a variety of inorganic, organic, and bioorganic compounds, providing in an abridged form the material partly discussed by the author in a series of parent reviews. Major attention is focused herewith on the publications of the recent years, which were not reviewed elsewhere.

Nonconventional NMR Spin-Coupling Constants in Oligosaccharide Conformational Modeling: Structural Dependencies Determined from Density Functional Theory Calculations

Nonconventional NMR spin-coupling constants were investigated to determine their potential as conformational constraints in MA'AT modeling of the O-glycosidic linkages of oligosaccharides. Four (¹ J _C1',H1', ¹ J _C1',C2', ² J _C1',H2', and ² J _C2',H1') and eight (¹ J _C4,H4, ¹ J _C3,C4, ¹ J _C4,C5, ² J _C3,H4, ² J _C4,H3, ² J _C5,H4, ² J _C4,H5, and ² J _C3,C5) spin-couplings in methyl β-d-galactopyranosyl-(1→4)-β-d-glucopyranoside (methyl β-lactoside) were calculated using density functional theory (DFT) to determine their dependencies on O-glycosidic linkage C-O torsion angles, ϕ and ψ, respectively. Long-range ⁴ J _H1',H4 was also examined as a potential conformational constraint of either ϕ or ψ. Secondary effects of exocyclic (hydroxyl) C-O bond rotation within or proximal to these coupling pathways were investigated. Based on the findings of methyl β-lactoside, analogous J-couplings were studied in five additional two-bond O-glycosidic linkages [βGlcNAc-(1→4)-βMan, 2-deoxy-βGlc-(1→4)-βGlc, αMan-(1→3)-βMan, αMan-(1→2)-αMan, and βGlcNAc(1→2)-αMan] to determine whether the coupling behaviors observed in methyl β-lactoside were more broadly observed. Of the 13 nonconventional J-couplings studied, 7 exhibit properties that may be useful in future MA'AT modeling of O-glycosidic linkages, none of which involve coupling pathways that include the linkage C-O bonds. The findings also provide new insights into the general effects of exocyclic C-O bond conformation on the magnitude of experimental spin-couplings in saccharides and other hydroxyl-containing molecules.

MA'AT: A Web-Based Application to Determine Rotamer Population Distributions in Solution from Nuclear Magnetic Resonance Spin-Coupling Constants

A hybrid experimental-computational method to determine conformational equilibria of molecules in solution has been developed based on the use of redundant nuclear magnetic resonance (NMR) spin-spin coupling constants (spin-couplings; J-couplings), density functional theory (DFT) calculations, and circular statistics. The mathematics that underpins the method, known as MA'AT analysis, is presented, and key components of a computer program that applies this algorithm are discussed. The method was tested using single-state and multi-state models to identify the factors required to obtain reliable results, to establish the limitations of the method, and to highlight techniques to evaluate the uniqueness of solution.

MA'AT Analysis of Aldofuranosyl Rings: Unbiased Modeling of Conformational Equilibria and Dynamics in Solution

MA'AT analysis has been applied to methyl β-d-ribofuranoside (3) and methyl 2-deoxy-β-d-erythro-pentofuranoside (4) to demonstrate the ability of this new experimental method to determine multi-state conformational equilibria in solution. Density functional theory (DFT) was used to obtain parameterized equations for >20 NMR spin-coupling constants sensitive to furanose ring conformation in 3 and 4, and these equations were used in conjunction with experimental spin-couplings to produce unbiased MA'AT models of ring pseudorotation. These models describe two-state north-south conformational exchange consistent with results obtained from traditional treatments of more limited sets of NMR spin-couplings (e.g., PSEUROT). While PSEUROT, MA'AT, and aqueous molecular dynamics models yielded similar two-state models, MA'AT analysis gives more reliable results since significantly more experimental observables are employed compared to PSEUROT, and no assumptions are needed to render the fitting tractable. MA'AT models indicate a roughly equal distribution of north and south ring conformers of 4 in aqueous (²H₂O) solution compared to ∼80% north forms for 3. Librational motion about the mean pseudorotation phase angles P of the preferred north and south conformers of 3 in solution is more constrained than that for 4. The greater rigidity of the β-ribo ring may be caused by synergistic stereoelectronic effects and/or noncovalent (e.g., hydrogen-bonding) interactions in solution that preferentially stabilize north forms of 3. MA'AT analysis of oligonucleotides and other furanose ring-containing biomolecules promises to improve current experimental models of sugar ring behavior in solution and help reveal context effects on ring conformation in more complex biologically important systems.

Two-bond 13C-13C spin-coupling constants in saccharides: dependencies on exocyclic hydroxyl group conformation

Seven doubly ¹³C-labeled isotopomers of methyl β-D-glucopyranoside, methyl β-D-xylopyranoside, methyl β-D-galactopyranoside, methyl β-D-galactopyranosyl-(1→4)-β-D-glucopyranoside and methyl β-D-galactopyranosyl-(1→4)-β-D-xylopyranoside were prepared, crystallized, and studied by single-crystal X-ray crystallography and solid-state ¹³C NMR spectroscopy to determine experimentally the dependence of ²J_C1,C3 values in aldopyranosyl rings on the C1-C2-O2-H torsion angle, θ₂, involving the C2 carbon of the C1-C2-C3 coupling pathway. Using X-ray crystal structures to determine θ₂ in crystalline samples and by selecting compounds that exhibit a relatively wide range of θ₂ values in the crystalline state, ²J_C1,C3 values measured in crystalline samples were plotted against θ₂ and the resulting plot compared to that obtained from density functional theory (DFT) calculations. For θ₂ values ranging from ∼90° to ∼240°, very good agreement was observed between the experimental and theoretical plots, providing strong validation of DFT-calculated spin-coupling dependencies on exocyclic C-O bond conformation involving the central carbon of geminal C-C-C coupling pathways. These findings provide new experimental evidence supporting the use of ²J_CCC values as non-conventional spin-coupling constraints in MA'AT conformational modeling of saccharides in solution, and the use of NMR spin-couplings not involving coupled hydroxyl hydrogens as indirect probes of C-O bond conformation. Solvomorphism was observed in crystalline βGal-(1→4)-βGlcOCH₃ wherein the previously-reported methanol solvate form was found to spontaneously convert to a monohydrate upon air-drying, leading to small but discernible conformational changes in, and a new crystalline form of, this disaccharide.

Isopropyl 3-deoxy-α-D-ribo-hexopyranoside (isopropyl 3-deoxy-α-D-glucopyranoside): evaluating trends in structural parameters

Isopropyl 3-deoxy-α-D-ribo-hexopyranoside (isopropyl 3-deoxy-α-D-glucopyranoside), C₉H₁₈O₅, (I), crystallizes from a methanol-ethyl acetate solvent mixture at room temperature in a ⁴C₁ chair conformation that is slightly distorted towards the ^C5S_C1 twist-boat form. A comparison of the structural parameters in (I), methyl α-D-glucopyranoside, (II), α-D-glucopyranosyl-(1→4)-D-glucitol (maltitol), (III), and 3-deoxy-α-D-ribo-hexopyranose (3-deoxy-α-D-glucopyranose), (IV), shows that most endocyclic and exocyclic bond lengths, valence bond angles and torsion angles in the aldohexopyranosyl rings are more affected by anomeric configuration, aglycone structure and/or the conformation of exocyclic substituents, such as hydroxymethyl groups, than by monodeoxygenation at C3. The structural effects observed in the crystal structures of (I)-(IV) were confirmed though density functional theory (DFT) calculations in computed structures (I)^c-(IV)^c. Exocyclic hydroxymethyl groups adopt the gauche-gauche (gg) conformation (H5 anti to O6) in (I) and (III), and the gauche-trans (gt) conformation (C4 anti to O6) in (II) and (IV). The O-glycoside linkage conformations in (I) and (III) resemble those observed in disaccharides containing β-(1→4) linkages.

A combined NMR, MD and DFT conformational analysis of 9-O-acetyl sialic acid-containing GM3 ganglioside glycan and its 9-N-acetyl mimic

O-Acetylation of carbohydrates such as sialic acids is common in nature, but its role is not clearly understood due to the lability of O-acetyl groups. We demonstrated previously that 9-acetamido-9-deoxy-N-acetylneuraminic acid (Neu5Ac9NAc) is a chemically and biologically stable mimic of the 9-O-acetyl-N-acetylneuraminic acid (Neu5,9Ac2) of the corresponding sialoglycans. Here, a systematic nuclear magnetic resonance (NMR) spectroscopic and molecular dynamics (MD) simulation study was undertaken for Neu5,9Ac2-containing GM3 ganglioside glycan (GM3-glycan) and its Neu5Ac9NAc analog. GM3-glycan with Neu5Ac as the non-O-acetyl form of Neu5,9Ac2 was used as a control. Complete 1H and 13C NMR chemical shift assignments, three-bond 1H-13C trans-glycosidic coupling constants (3JCH), accurate 1H-1H coupling constants (3JHH), nuclear Overhauser effects and hydrogen bonding detection were carried out. Results show that structural modification (O- or N-acetylation) on the C-9 of Neu5Ac in GM3 glycan does not cause significant conformational changes on either its glycosidic dihedral angles or its secondary structure. All structural differences are confined to the Neu5Ac glycerol chain, and minor temperature-dependent changes are seen in the aglycone portion. We also used Density Functional Theory (DFT) quantum mechanical calculations to improve currently used 3JHH Karplus relations. Furthermore, OH chemical shifts were assigned at -10°C and no evidence of an intramolecular hydrogen bond was observed. The results provide additional evidence regarding structural similarities between sialosides containing 9-N-acetylated and 9-O-acetylated Neu5Ac and support the opportunity of using 9-N-acetylated Neu5Ac as a stable mimic to study the biochemical role of 9-O-acetylated Neu5Ac.

Reconciling MA'AT and molecular dynamics models of linkage conformation in oligosaccharides

MA'AT conformational models of the phi torsion angles of O-glycosidic linkages differ from those obtained from MD simulation. To determine the source of the discrepancy, MA'AT analyses were performed using DFT-derived equations obtained with and without psi constraints. The resulting phi models were essentially the same, indicating a force-field problem. Circular standard deviations (CSDs) were found to provide reliable estimates of torsional averaging.

Use of Circular Statistics To Model αMan-(1→2)-αMan and αMan-(1→3)-α/βMan O-Glycosidic Linkage Conformation in 13C-Labeled Disaccharides and High-Mannose Oligosaccharides

A new experimental method, MA' AT analysis, has been applied to investigate the conformational properties of O-glycosidic linkages in several biologically important mannose-containing di- and oligosaccharides. Methyl α-d-mannopyranosyl-(1→2)-α-d-mannopyranoside (2), methyl α-d-mannopyranosyl-(1→3)-α-d-mannopyranoside (3), and methyl α-d-mannopyranosyl-(1→3)-β-d-mannopyranoside (4) were prepared with selective ¹³C-enrichment to enable the measurement of NMR scalar couplings across their internal O-glycosidic linkages. Density functional theory (DFT) was used to parameterize equations for J_CH and J_CC values in 2-4 that are sensitive to phi (ϕ) and psi (ψ). The experimental J-couplings and parameterized equations were treated using a circular statistics algorithm encoded in the MA' AT program. Conformations about ϕ and ψ treated using single-state von Mises models gave excellent fits to the ensembles of redundant J-couplings. Mean values and circular standard deviations (CSDs) for each linkage torsion angle ϕ (CSD) and ψ (CSD) in 2, -29° (25°) and 20° (22°); in 3, -36° (36°) and 8° (27°); in 4, -37° (34°) and 10° (26°); ϕ = H1'-C1'-O1'-CX and ψ = C1'-O1'-CX-HX (CX = aglycone carbon) were compared to histograms obtained from 1 μs aqueous molecular dynamics (MD) simulations and X-ray database statistical analysis. MA' AT-derived models of ψ were in very good agreement with the MD and X-ray data, but not those of ϕ, suggesting a need for force field revision. The effect of structural context on linkage conformation was also investigated in four selectively ¹³C-labeled homomannose tri- and tetrasaccharides using the MA' AT method. In the cases examined, context effects were found to be small.