Conformational Populations of β-(1→4) O-Glycosidic Linkages Using Redundant NMR J-Couplings and Circular Statistics

Geminal 13C-1H NMR Spin-Coupling Constants in Furanose Rings: New Empirical Correlations with Conformation

Density functional theory (DFT) calculations have been used to develop a new approach to interpreting geminal (two-bond) 2 J CCH NMR spin-coupling constants in saccharides containing aldofuranosyl (five-membered) rings. In the biologically important β-d-ribofuranosyl and 2-deoxy-β-d-ribofuranosyl (2-deoxy-β-d-erythro-pentofuranosyl) rings that were used as models, many of the 2 J CCH values associated with coupling pathways involving an endocyclic C-C bond depend linearly on P/π, a measure of ring conformation. In most cases, the endocyclic C-C bond is present in the coupling pathway. In other cases, the 2 J CCH value depends linearly on either an adjacent C-C bond torsion angle or shows no systematic relationship with any endocyclic C-C bond torsion angle. In the latter case, secondary (remote) structural effects, defined as those that primarily affect C-C or C-H bond lengths in the C-C-H coupling pathway, cause the 2 J CCH value to behave with less predictability. These effects apparently cancel and lead to linearity involving an adjacent C-C bond in some cases. These findings provide a new conceptual framework to understand and exploit the dependencies of geminal 13C-1H NMR spin-couplings on furanose ring conformation. They also reveal the effect of exocyclic C-O bond torsion angles on the magnitudes and signs of 2 J CCH values in saccharides, a complication that remains to be addressed before 2 J CCH values can be used quantitatively in single- and multi-state MA'AT modeling of redundant NMR J-values in furanosyl rings.

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.

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.

Revealing elusive conformations of sucrose from hydrogen bond J-coupling in H2O: A combined NMR and quantum mechanics study

Hydrogen bonding is a crucial feature of biomolecules, but its characterization in glycans dissolved in aqueous solutions is challenging due to rapid hydrogen exchange between hydroxyl groups and H2O. In principle, the scalar (J) coupling constant can reveal the relative orientation of the atoms in the molecule. In contrast to J-coupling through H-bonds reported in proteins and nucleic acids, research on J-coupling through H-bonds in glycans dissolved in water is lacking. Here, we use sucrose as a model system for H-bonding studies; its structure, which consists of glucose (Glc) and fructose (Frc), is well-studied, and it is readily available. We apply the in-phase, antiphase-HSQC-TOCSY and quantify previously unreported through H-bond J-values for Frc-OH1-Glc-OH2 in H2O. While earlier reports of Brown and Levy indicate this H-bond as having only a single direction, our reported findings indicate the potential presence of two involving these same atoms, namely, G2OH ➔ F1O and F1OH ➔ G2O (where F and G stand for Frc and Glc, respectively). The calculated density functional theory J-values for the G2OH ➔ F1O agree with the experimental values. Additionally, we detected four other possible H-bonds in sucrose, which require different phi, psi (ϕ, ψ) torsion angles. The ϕ, ψ values are consistent with previous predictions of du Penhoat et al. and Venable et al. Our results will provide new insights into the molecular structure of sucrose and its interactions with proteins.

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.

Influence of Selective Deoxyfluorination on the Molecular Structure of Type-2 N-Acetyllactosamine

N-Acetyllactosamine is a common saccharide motif found in various biologically active glycans. This motif usually works as a backbone for additional modifications and thus significantly influences glycan conformational behavior and biological activity. In this work, we have investigated the type-2 N-acetyllactosamine scaffold using the complete series of its monodeoxyfluorinated analogs. These glycomimetics have been studied by molecular mechanics, quantum mechanics, X-ray crystallography, and various NMR techniques, which have provided a comprehensive and complete insight into the role of individual hydroxyl groups in the conformational behavior and lipophilicity of N-acetyllactosamine.

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.

Exploring long-range fluorine-carbon J-coupling for conformational analysis of deoxyfluorinated disaccharides: A combined computational and NMR study

In this study, we investigated the potential of long-range fluorine-carbon J-coupling for determining the structures of deoxyfluorinated disaccharides. Three disaccharides, previously synthesized as potential galectin inhibitors, exhibited through-space fluorine-carbon J-couplings. In our independent conformational analysis of these disaccharide derivatives, we employed a combination of density functional theory (DFT) calculations and nuclear magnetic resonance (NMR) experiments. By comparing the calculated nuclear shieldings with the experimental carbon chemical shifts, we were able to identify the most probable conformers for each compound. A model comprising fluoromethane and methane molecules was used to study the relationship between molecular arrangements and intermolecular through-space J-coupling. Our study demonstrates the important effect of internuclear distance and molecular orientation on the magnitude of fluorine-carbon coupling. The experimental values for the fluorine-carbon through-space couplings (TSCs) of the disaccharides corresponded with values calculated for the most probable conformers identified by the conformational analysis. These results unlock the broader application of fluorine-carbon TSCs as powerful tools for conformational analysis of flexible molecules, offering valuable insights for future structural investigations.

Does Inter-Residue Hydrogen Bonding in β-(1→4)-Linked Disaccharides Influence Linkage Conformation in Aqueous Solution?

Two disaccharides, methyl β-d-galactopyranosyl-(1→4)-α-d-glucopyranoside (1) and methyl β-d-galactopyranosyl-(1→4)-3-deoxy-α-d-ribo-hexopyranoside (3), were prepared with selective 13C-enrichment to allow measurement of six trans-O-glycosidic J-couplings (2JCOC, 3JCOCH, and 3JCOCC) in each compound. Density functional theory (DFT) was used to parameterize Karplus-like equations that relate these J-couplings to either ϕ or ψ. MA'AT analysis was applied to both linkages to determine mean values of ϕ and ψ in each disaccharide and their associated circular standard deviations (CSDs). Results show that deoxygenation at C3 of 1 has little effect on both the mean values and librational motions of the linkage torsion angles. This finding implies that, if inter-residue hydrogen bonding between O3H and O5' of 1 is present in aqueous solution and persistent, it plays little if any role in dictating preferred linkage conformation. Hydrogen bonding may lower the energy of the preferred linkage geometry but does not determine it to any appreciable extent. Aqueous 1-μs MD simulation supports this conclusion and also indicates greater conformational flexibility in deoxydisaccharide 3 in terms of sampling several, conformationally distinct, higher-energy conformers in solution. The populations of these latter conformers are low (3-14%) and could not be validated by MA'AT analysis. If the MD model is correct, however, C3 deoxygenation does enable conformational sampling over a wider range of ϕ/ψ values, but linkage conformation in the predominant conformer is essentially identical in both 1 and 3.

MA'AT Analysis: Unbiased Multi-State Conformational Modeling of Exocyclic Hydroxymethyl Group Conformation in Methyl Aldohexopyranosides

MA'AT analysis (J. Chem. Inf. Model. 2022, 62, 3135-3141) has been applied to model exocyclic hydroxymethyl group conformation in methyl β-D-glucopyranoside (βGlcOMe), methyl β-D-galactopyranoside (βGalOMe), and methyl β-D-mannopyranoside (βManOMe) in an unbiased manner. Using up to eight NMR J-couplings sensitive to rotation about the C5-C6 bond (torsion angle ω), two-state models of ω were obtained that are qualitatively consistent with the relative populations of the gg, gt, and tg rotamers reported previously. MA'AT analysis gave consistent unbiased gt ⇌ tg models of ω in βGalOMe, with gt more populated than tg and mean values of ω for each population similar to those obtained from aqueous 1-μs MD simulation. Using different combinations of J-couplings had little effect on the βGalOMe model in terms of the mean values of ω and circular standard deviations (CSDs). In contrast, MA'AT analysis of ω in βGlcOMe and βManOMe produced more than one two-state model independent of the ensemble of J-values used in the analyses. These models were characterized by gg ⇌ gt conformer exchange as expected, but the mean values of ω in both conformers varied significantly in the different fits, especially for the gg rotamer. Constrained (biased) MA'AT analyses in which only staggered geometries about ω were allowed gave RMSDs slightly larger than those obtained from the unbiased fits, precluding an assignment of an unbiased model. It is unclear why MA'AT analysis gives consistent and predictable unbiased models of ω in βGalOMe but not in βGlcOMe and βManOMe. One possibility is that the distribution of ω in one or both of the gg and gt conformers in the latter does not conform to a von Mises function (i.e., is not Gaussian-like), but rather to a broad and/or flat distribution that cannot be fit by the current version of MA'AT. Nevertheless, the results of this study provide new evidence of the ability of MA'AT analysis to treat multi-state conformational exchange using only experimental NMR data, extending recent MA'AT applications to furanosyl ring pseudorotation (Biochemistry 2022, 61, 239-251).

MA'AT Analysis: Probability Distributions of Molecular Torsion Angles in Solution from NMR Spectroscopy

ConspectusMonosaccharides adopt multiple conformations in solution, and this structural complexity increases significantly when they are assembled into oligosaccharides and polysaccharides. Characterization of the conformational properties of saccharides in solution by NMR spectroscopy has been hampered by several complicating factors, including difficulty interpreting spectra because of significant signal overlap, population averaging of NMR parameters, and unique properties of the spectra that make accurate measurements of NMR parameters prone to error (e.g., non-first-order effects on J-couplings). Current conformational assignments rely heavily on theoretical calculations, especially molecular dynamics (MD) simulations, to interpret the experimental NMR parameters. While these studies assert that the available experimental data fit the calculated models well, a lack of independent experimental validation of the force fields from which MD models are derived and an inability to test all possible models that might be compatible with the experimental data in an unbiased manner make the approach less than ideal.NMR spin couplings or J-couplings have been used as structure constraints in organic and other types of molecules for more than six decades. The dihedral angle dependence of vicinal (three-bond) 1H-1H spin couplings (3JHH) first described by Karplus led to an explosion of applications for a wide range of conformational problems. Other vicinal J-couplings (e.g., 3JCCOP, 3JHCOP, and 3JCOCH) have been found to exhibit similar dihedral angle dependencies. 3J values have been used to assign the preferred conformation in molecules that are conformationally homogeneous. However, many molecules, particularly those in biological systems, are conformationally flexible, which complicates structural interpretations of J values in solution. Three-state staggered models are often assumed in order to deconvolute the conformationally averaged J values into conformer populations. While widely applied, this approach assumes highly idealized models of molecular torsion angles that are likely to be poor representations of those found in solution. In addition, this treatment often gives negative populations and neglects the presence of librational averaging of molecular torsion angles.Recent work in this research group has focused on the development of a hybrid experimental-computational method, MA'AT analysis, that provides probability distributions of molecular torsion angles in solution that can be superimposed on those obtained by MD. Ensembles of redundant NMR spin couplings, including 3J (vicinal), 2J (geminal), and sometimes 1J (direct) values, are used in conjunction with circular statistics to provide single- and multistate models of these angles. MA'AT analysis provides accurate mean torsion angles and circular standard deviations (CSDs) of each mean angle that describe the librational motion about the angle. Both conformational equilibria and dynamics are revealed by the method. In this Account, the salient features of MA'AT analysis are discussed, including some applications to conformational problems involving saccharides and peptides.

Primary Structure of Glycans by NMR Spectroscopy

Glycans, carbohydrate molecules in the realm of biology, are present as biomedically important glycoconjugates and a characteristic aspect is that their structures in many instances are branched. In determining the primary structure of a glycan, the sugar components including the absolute configuration and ring form, anomeric configuration, linkage(s), sequence, and substituents should be elucidated. Solution state NMR spectroscopy offers a unique opportunity to resolve all these aspects at atomic resolution. During the last two decades, advancement of both NMR experiments and spectrometer hardware have made it possible to unravel carbohydrate structure more efficiently. These developments applicable to glycans include, inter alia, NMR experiments that reduce spectral overlap, use selective excitations, record tilted projections of multidimensional spectra, acquire spectra by multiple receivers, utilize polarization by fast-pulsing techniques, concatenate pulse-sequence modules to acquire several spectra in a single measurement, acquire pure shift correlated spectra devoid of scalar couplings, employ stable isotope labeling to efficiently obtain homo- and/or heteronuclear correlations, as well as those that rely on dipolar cross-correlated interactions for sequential information. Refined computer programs for NMR spin simulation and chemical shift prediction aid the structural elucidation of glycans, which are notorious for their limited spectral dispersion. Hardware developments include cryogenically cold probes and dynamic nuclear polarization techniques, both resulting in enhanced sensitivity as well as ultrahigh field NMR spectrometers with a 1H NMR resonance frequency higher than 1 GHz, thus improving resolution of resonances. Taken together, the developments have made and will in the future make it possible to elucidate carbohydrate structure in great detail, thereby forming the basis for understanding of how glycans interact with other molecules.

Glycosidic α-linked mannopyranose disaccharides: an NMR spectroscopy and molecular dynamics simulation study employing additive and Drude polarizable force fields

D-Mannose is a structural component in N-linked glycoproteins from viruses and mammals as well as in polysaccharides from fungi and bacteria. Structural components often consist of D-Manp residues joined via α-(1→2)-, α-(1→3)-, α-(1→4)- or α-(1→6)-linkages. As models for these oligo- and polysaccharides, a series of mannose-containing disaccharides have been investigated with respect to conformation and dynamics. Translational diffusion NMR experiments were performed to deduce rotational correlation times for the molecules, 1D ¹H,¹H-NOESY and 1D ¹H,¹H-T-ROESY NMR experiments were carried out to obtain inter-residue proton-proton distances and one-dimensional long-range and 2D J-HMBC experiments were acquired to gain information about conformationally dependent heteronuclear coupling constants across glycosidic linkages. To attain further spectroscopic data, the doubly ¹³C-isotope labeled α-D-[1,2-¹³C₂]Manp-(1→4)-α-D-Manp-OMe was synthesized thereby facilitating conformational analysis based on ¹³C,¹³C coupling constants as interpreted by Karplus-type relationships. Molecular dynamics simulations were carried out for the disaccharides with explicit water as solvent using the additive CHARMM36 and Drude polarizable force fields for carbohydrates, where the latter showed broader population distributions. Both simulations sampled conformational space in such a way that inter-glycosidic proton-proton distances were very well described whereas in some cases deviations were observed between calculated conformationally dependent NMR scalar coupling constants and those determined from experiment, with closely similar root-mean-square differences for the two force fields. However, analyses of dipole moments and radial distribution functions with water of the hydroxyl groups indicate differences in the underlying physical forces dictating the wider conformational sampling with the Drude polarizable versus additive C36 force field and indicate the improved utility of the Drude polarizable model in investigating complex carbohydrates.

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.

Using NMR to Dissect the Chemical Space and O-Sulfation Effects within the O- and S-Glycoside Analogues of Heparan Sulfate

Heparan sulfate (HS), a sulfated linear carbohydrate that decorates the cell surface and extracellular matrix, is ubiquitously distributed throughout the animal kingdom and represents a key regulator of biological processes and a largely untapped reservoir of potential therapeutic targets. The temporal and spatial variations in the HS structure underpin the concept of "heparanome" and a complex network of HS binding proteins. However, despite its widespread biological roles, the determination of direct structure-to-function correlations is impaired by HS chemical heterogeneity. Attempts to correlate substitution patterns (mostly at the level of sulfation) with a given biological activity have been made. Nonetheless, these do not generally consider higher-level conformational effects at the carbohydrate level. Here, the use of NMR chemical shift analysis, NOEs, and spin-spin coupling constants sheds new light on how different sulfation patterns affect the polysaccharide backbone geometry. Furthermore, the substitution of native O-glycosidic linkages to hydrolytically more stable S-glycosidic forms leads to observable conformational changes in model saccharides, suggesting that alternative chemical spaces can be accessed and explored using such mimetics. Employing a series of systematically modified heparin oligosaccharides (as a proxy for HS) and chemically synthesized O- and S-glycoside analogues, the chemical space occupied by such compounds is explored and described.

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.

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 H₂O/²H₂O and DMSO-d₆ solvents in selectively ¹³C-labeled compounds, to model the C1-C2-N2-C1' torsion angle (θ₁) in GlcNAc and ManNAc residues. MA'AT analysis gave mean values of θ₁ 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 θ₁ values were found to be small in a βGlcNAc-(1→4)-βGlcNAc disaccharide. MA'AT analysis gave a mean value of θ₁ of 249° for αManNAc in H₂O/²H₂O, 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.

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.

Use of Raman and Raman optical activity to extract atomistic details of saccharides in aqueous solution

Sugars are crucial components in biosystems and industrial applications. In aqueous environments, the natural state of short saccharides or charged glycosaminoglycans is floating and wiggling in solution. Therefore, tools to characterize their structure in a native aqueous environment are crucial but not always available. Here, we show that a combination of Raman/ROA and, on occasions, NMR experiments with Molecular Dynamics (MD) and Quantum Mechanics (QM) is a viable method to gain insights into structural features of sugars in solutions. Combining these methods provides information about accessible ring puckering conformers and their proportions. It also provides information about the conformation of the linkage between the sugar monomers, i.e., glycosidic bonds, allowing for identifying significantly accessible conformers and their relative abundance. For mixtures of sugar moieties, this method enables the deconvolution of the Raman/ROA spectra to find the actual amounts of its molecular constituents, serving as an effective analytical technique. For example, it allows calculating anomeric ratios for reducing sugars and analyzing more complex sugar mixtures to elucidate their real content. Altogether, we show that combining Raman/ROA spectroscopies with simulations is a versatile method applicable to saccharides. It allows for accessing many features with precision comparable to other methods routinely used for this task, making it a viable alternative. Furthermore, we prove that the proposed technique can scale up by studying the complicated raffinose trisaccharide, and therefore, we expect its wide adoption to characterize sugar structural features in solution.

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.

Computational NMR of Carbohydrates: Theoretical Background, Applications, and Perspectives

This review is written amid a marked progress in the calculation of NMR parameters of carbohydrates substantiated by a vast amount of experimental data coming from several laboratories worldwide. By no means are we trying to cover in the present compilation a huge amount of all available data. The main idea of the present review was only to outline general trends and perspectives in this dynamically developing area on the background of a marked progress in theoretical and computational NMR. Presented material is arranged in three basic sections: (1)-a brief theoretical introduction; (2)-applications and perspectives in computational NMR of monosaccharides; and (3)-calculation of NMR chemical shifts and spin-spin coupling constants of di- and polysaccharides.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Conformational Isomerism of 3-Chalcogenomethyl-N-Methyl-2-Pyrrolidinones: Insights from NMR Spectroscopy and Molecular Modeling

A conformational analysis of N-methyl-2-pyrrolidinone 3-substituted by methoxyl, thiomethoxyl, and selenomethoxyl is reported by means of ¹H nuclear magnetic resonance spectroscopy and electronic structure calculations. The five-membered ring has an envelope conformation with the α-carbonyl substituent being able to assume two positions: pseudo-axial and pseudo-equatorial. In vacuum, the calculations pointed to the pseudo-axial conformer as the most stable one, and this preference increases with the size of the substituent and a decrease in its electronegativity. Natural bond orbital analysis evidenced the importance of electron delocalization on the stability, and a principal component of analysis (PCA) plot of the hyperconjugative interactions revealed the main ones. Steric and electrostatic effects were also investigated by energy decomposition analysis.

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.

Glycosidic linkage, N-acetyl side-chain, and other structural properties of methyl 2-acetamido-2-deoxy-β-D-glucopyranosyl-(1→4)-β-D-mannopyranoside monohydrate and related compounds

The crystal structure of methyl 2-acetamido-2-deoxy-β-D-glycopyranosyl-(1→4)-β-D-mannopyranoside monohydrate, C₁₅H₂₇NO₁₁·H₂O, was determined and its structural properties compared to those in a set of mono- and disaccharides bearing N-acetyl side-chains in βGlcNAc aldohexopyranosyl rings. Valence bond angles and torsion angles in these side chains are relatively uniform, but C-N (amide) and C-O (carbonyl) bond lengths depend on the state of hydrogen bonding to the carbonyl O atom and N-H hydrogen. Relative to N-acetyl side chains devoid of hydrogen bonding, those in which the carbonyl O atom serves as a hydrogen-bond acceptor display elongated C-O and shortened C-N bonds. This behavior is reproduced by density functional theory (DFT) calculations, indicating that the relative contributions of amide resonance forms to experimental C-N and C-O bond lengths depend on the solvation state, leading to expectations that activation barriers to amide cis-trans isomerization will depend on the polarity of the environment. DFT calculations also revealed useful predictive information on the dependencies of inter-residue hydrogen bonding and some bond angles in or proximal to β-(1→4) O-glycosidic linkages on linkage torsion angles φ and ψ. Hypersurfaces correlating φ and ψ with the linkage C-O-C bond angle and total energy are sufficiently similar to render the former a proxy of the latter.

Computational 1 H NMR: Part 3. Biochemical studies

This is the third and the last part of three closely interrelated reviews dealing with computation of ¹ H nuclear magnetic resonance chemical shifts and ¹ H-¹ H spin-spin coupling constants. Present review deals with the computation of these parameters in biologically active natural products, carbohydrates, and other molecules of biological origin focusing on stereochemical applications of computational ¹ H nuclear magnetic resonance to these objects.

Glycan structures and their interactions with proteins. A NMR view

Carbohydrate molecules are essential actors in key biological events, being involved as recognition points for cell-cell and cell-matrix interactions related to health and disease. Despite outstanding advances in cryoEM, X-ray crystallography and NMR still remain the most employed techniques to unravel their conformational features and to describe the structural details of their interactions with biomolecular receptors. Given the intrinsic flexibility of saccharides, NMR methods are of paramount importance to deduce the extent of motion around their glycosidic linkages and to explore their receptor-bound conformations. We herein present our particular view on the latest advances in NMR methodologies that are permitting to magnify their applications for deducing glycan conformation and dynamics and understanding the recognition events in which there are involved.

13C-13C spin-coupling constants in crystalline 13C-labeled saccharides: conformational effects interrogated by solid-state 13C NMR spectroscopy

Solid-state ¹³C NMR spectroscopy has been used in conjunction with selectively ¹³C-labeled mono- and disaccharides to measure ¹³C-¹³C spin-couplings (J_CC) in crystalline samples. This experimental approach allows direct correlation of J_CC values with specific molecular conformations since, in crystalline samples, molecular conformation is essentially static and can be determined by X-ray crystallography. J_CC values measured in the solid-state in known molecular conformations can then be compared to corresponding J_CC values calculated in the same conformations using density functional theory (DFT). The latter comparisons provide important validation of DFT-calculated J-couplings, which is not easily obtained by other approaches and is fundamental to obtaining reliable experiment-based conformational models from redundant J-couplings by MA'AT analysis. In this study, representative ¹J_CC, ²J_CCC and ³J_COCC values were studied as either intra-residue couplings in the aldohexopyranosyl rings of monosaccharides or inter-residue (trans-glycoside) couplings in disaccharides. The results demonstrate that (a) accurate J_CC values can be measured in crystalline saccharides that have been suitably labeled with ¹³C, and (b) DFT-calculated J_CC values compare favorably with those determined by solid-state ¹³C NMR when molecular conformation is a constant in both determinations.

Solution Structure of Mannobioses Unravelled by Means of Raman Optical Activity

Structural analysis of carbohydrates is a complicated endeavour, due to the complexity and diversity of the samples at hand. Herein, we apply a combined computational and experimental approach, employing molecular dynamics (MD) and density functional theory (DFT) calculations together with NMR and Raman optical activity (ROA) measurements, in the structural study of three mannobiose disaccharides, consisting of two mannoses with varying glycosidic linkages. The disaccharide structures make up the scaffold of high mannose glycans and are therefore important targets for structural analysis. Based on the MD population analysis and NMR, the major conformers of each mannobiose were identified and used as input for DFT analysis. By systematically varying the solvent models used to describe water interacting with the molecules and applying overlap integral analysis to the resulting calculational ROA spectra, we found that a full quantum mechanical/molecular mechanical approach is required for an optimal calculation of the ROA parameters. Subsequent normal mode analysis of the predicted vibrational modes was attempted in order to identify possible marker bands for glycosidic linkages. However, the normal mode vibrations of the mannobioses are completely delocalised, presumably due to conformational flexibility in these compounds, rendering the identification of isolated marker bands unfeasible.

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.

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.

Theoretical calculations of carbon-hydrogen spin-spin coupling constants

Structural applications of theoretical calculations of carbon-hydrogen spin-spin coupling constants are reviewed covering papers published mainly during the last 10-15 years with a special emphasis on the most notable studies of hybridization, substitution and stereoelectronic effects together with the investigation of hydrogen bonding and intermolecular interactions. The wide scope of different applications of calculated carbon-hydrogen couplings in the structural elucidation of particular classes of organic and bioorganic molecules is reviewed, concentrating mainly on saturated, unsaturated, aromatic and heteroaromatic compounds and their functional derivatives, as well as on natural compounds and carbohydrates. The review is dedicated to Professor Emeritus Michael Barfield in view of his invaluable pioneering contribution to this field.