Conformational disorder in the crystal structure of methyl 2-acetamido-2-deoxy-ß-D-glucopyranosyl-(1→4)-2-acetamido-2-deoxy-ß-D-glucopyranoside (methyl ß-chitobioside) methanol monosolvate.

Methyl 2-acetamido-2-deoxy-ß-D-glucopyranosyl-(1→4)-2-acetamido-2-deoxy-ß-D-glucopyranoside (methyl ß-chitobioside), (IV), crystallizes from aqueous methanol at room temperature to give a structure (C17H30N2O22·CH3OH) containing conformational disorder in the exocyclic hydroxymethyl group of one of its ßGlcNAc residues. As observed in other X-ray structures of disaccharides containing ß-(1→4) O-glycosidic linkages, inter-residue hydrogen bonding between O3H of the ßGlcNAc bearing the OCH3 aglycone and O5 of the adjacent ßGlcNAc is observed based on the 2.79 Šinternuclear distance between the O atoms. The structure of (IV) was compared to that determined previously for 2-acetamido-2-deoxy-ß-D-glucopyranosyl-(1→4)-2-acetamido-2-deoxy-ß-D-glucopyranose (ß-chitobiose), (III). The O-glycosidic linkage torsion angles, phi (φ) and psi (ψ), in (III) and (IV) differ by 6-8°. The N-acetyl side chain conformation in (III) and (IV) shows some context dependence, with the C1-C2-N-Ccar torsion angle 10-15° smaller for the ßGlcNAc residue involved in the internal O-glycosidic linkage. In (IV), conformational disorder is observed in the exocyclic hydroxymethyl (-CH2OH) group in the ßGlcNAc residue bearing the OCH3 aglycone, and a fitting of the electron density indicates an approximate 50:50 distribution of the gauche-gauche (gg) and gauche-trans (gt) conformers in the lattice. Similar behavior is not observed in (III), presumably due to the different packing structure in the vicinity of the -CH2OH substituent that affects its ability to hydrogen bond to proximal donors/acceptors. Unlike (IV), a re-examination of the previously reported electron density of (III) revealed conformational disorder in the N-acetyl side chain attached to the reducing-end ßGlcNAc residue caused by rotation about the C2-N bond.

Methyl α-D-galactopyranosyl-(1→3)-ß-D-galactopyranoside and methyl ß-D-galactopyranosyl-(1→3)-ß-D-galactopyranoside: Glycosidic linkage conformation determined from MA'AT analysis.

MA'AT analysis has been applied to two biologically-important O-glycosidic linkages in two disaccharides, α-D-Galp-(1→3)-ß-D-GalpOMe (3) and ß-D-Galp-(1→3)-ß-D-GalpOMe (4). Using density functional theory (DFT) to obtain parameterized equations relating a group of trans-O-glycosidic NMR spin-couplings to either phi (ϕ') or psi (ψ'), and experimental 3JCOCH, 2JCOC, and 3JCOCC spin-couplings measured in aqueous solution in 13C-labeled isotopomers, probability distributions of ϕ' and ψ' in each linkage were determined and compared to those determined by aqueous 1-µs molecular dynamics (MD) simulation. Good agreement was found between the MA'AT and single-state MD conformational models of these linkages for the most part, with modest (approximately <15°) differences in the mean values of ϕ' and ψ', although the envelope of allowed angles (encoded in circular standard deviations or CSDs) is consistently larger for ϕ' determined from MA'AT analysis than from MD for both linkages. The MA'AT model of the α-Galp-(1→3)-ß-Galp linkage agrees well with those determined previously using conventional NMR methods (3JCOCH values and/or 1H-1H NOEs), but some discrepancy was observed for the ß-Galp-(1→3)-ß-Galp linkage, which may arise from errors in the conventions used to describe the linkage torsion angles. Statistical analyses of X-ray crystal structures show ranges of ϕ' and ψ' for both linkages that include the mean angles determined from MA'AT analyses, although both angles adopt a wide range of values in the crystalline state, with ϕ' in ß-Galp-(1→3)-ß-Galp linkages showing greater-than-expected conformational variability.

One-bond 13C-13C spin-coupling constants in saccharides: a comparison of experimental and calculated values by density functional theory using solid-state 13C NMR and X-ray crystallography.

Methyl aldohexopyranosides were 13C-labeled at contiguous carbons, crystallized, and studied by single-crystal X-ray crystallography and solid-state 13C nuclear magnetic resonance (NMR) spectroscopy to examine the degree to which density functional theory (DFT) can calculate one-bond 13C-13C spin-coupling constants (1JCC) in saccharides with sufficient accuracy to permit their use in MA'AT analysis, a newly-reported hybrid DFT/NMR method that provides probability distributions of molecular torsion angles in solution (Zhang et al., J. Phys. Chem. B, 2017, 121, 3042-3058; Meredith et al., J. Chem. Inf. Model., 2022, 62, 3135-3141). Experimental 1JCC values in crystalline samples of the doubly 13C-labeled compounds were measured by solid-state 13C NMR and compared to those calculated from five different DFT models: (1) 1JCC values calculated from single structures identical to those observed in crystalline samples by X-ray crystallography (all atom refinement); (2) 1JCC values calculated from the single structures in (1) but after Hirshfeld atom refinement (HAR); (3) 1JCC values calculated from the single structures in (1) after DFT-optimization of hydrogen atoms only; and (4 and 5) 1JCC values calculated in rotamers of torsion angle θ2 (C1-C2-O2-O2H) or ω (C4-C5-C6-O6) from which either specific or generalized parameterized equations were obtained and used to calculate 1JCC values in the specific θ2 or ω rotamers observed in crystalline samples. Good qualitative agreement was observed between calculated 1JCC values and those measured by solid-state 13C NMR regardless of the DFT model, but in no cases were calculated 1JCC values quantitative, differing (over-estimated) on average by 4-5% from experimental values. These findings, and those reported recently from solution NMR studies (Tetrault et al., J. Phys. Chem. B 2022, 126, 9506-9515), indicate that improvements in DFT calculations are needed before calculated 1JCC values can be used directly as reliable constraints in MA'AT analyses of saccharides in solution.

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 JHH, JCH, and 1JCC values for use in MA'AT analysis of saccharides. In contrast, the accuracy of DFT-calculated one-bond 1JCH and 1JCC 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 [1JCHcalc - 1JCHexp] (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, [1JCHMA'AT - 1JCHexp] (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 1JCC values, with an average % difference of 5.1 ± 0.2%. Calculated 1JCC values were consistently larger than experimental values. Strategy 2 also gave calculated 1JCC 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 1JCH values in saccharides can be calculated nearly quantitatively, but 1JCC values appear to be consistently overestimated by ∼5% using current DFT methods.

N-Acetyl Side-Chain Conformation in Saccharides: Solution Models Obtained from MA'AT Analysis.

MA'AT analysis has been applied to model the conformational properties of N-acetyl side-chains in biologically important GlcNAc and ManNAc monosaccharides and in a ßGlcNAc-(1→4)-ßGlcNAc disaccharide. Density functional theory calculations were conducted to obtain parameterized equations that relate the magnitudes and signs of 10 spin-coupling constants to conformations of the C2-N2 bonds of GlcNAc and ManNAc. Six of these equations were used with experimental J-couplings, measured in H2O/2H2O and DMSO-d6 solvents in selectively 13C-labeled compounds, to model the C1-C2-N2-C1' torsion angle (θ1) in GlcNAc and ManNAc residues. MA'AT analysis gave mean values of θ1 of 106° for αGlcNAc and ∼116° for ßGlcNAc residues, with circular standard deviations (CSDs) of 21-22° in aqueous solution, in excellent agreement with those obtained by aqueous molecular dynamics (MD) simulation. Parameter space plots revealed unique MA'AT fits of the data, and root mean squared deviations (<0.2 Hz) were twofold smaller than those back-calculated from MD models, indicating that the MA'AT models better fit the experimental J-couplings. Context effects on both θ1 values were found to be small in a ßGlcNAc-(1→4)-ßGlcNAc disaccharide. MA'AT analysis gave a mean value of θ1 of 249° for αManNAc in H2O/2H2O, with a CSD of ∼19°, with both values in good agreement with MD. MA'AT models of N-acetyl side-chains are similar to those obtained previously for O-acetyl side-chains (J. Phys. Chem. B 2017, 121, 66-77). Both O- and N-acetylation conformationally constrain the C-O or C-N bonds relative to the same bonds in unsubstituted compounds. The present work confirms the ability of MA'AT analysis to reveal relatively small changes in mean molecular torsion angles in solution and provides additional evidence of the method as an experimental tool complementary to MD simulation.