Limitations in the description of conformational preferences of C-disaccharides: The (1 → 3)-C-mannobiose case

On the origin of the electronic and magnetic circular dichroism of naphthyl C-glycosides: Anomeric configuration

Aryl C-glycosides, in which the glycosidic bond is changed to a carbon-carbon bond, are an important family of biologically-active compounds. They often serve as secondary metabolites or exhibit antibiotic and cytostatic activities. Their stability to hydrolysis has made them attractive targets for new drugs. Their conformational behavior often strongly influences the resulting function. Their detailed structural and conformational description is thus highly desirable. This work studies the structure of three different naphthyl C-glycosides using UV-vis absorption as well as electronic and magnetic circular dichroism. It also describes their conformational preferences using a combination of molecular dynamics and DFT calculations. The reliability of these preferences has been verified by simulations of spectral properties and a comparison with their measured spectra. In particular, ECD spectroscopy has been shown to distinguish easily between α- and β-pseudoanomers of aryl C-glycosides. Computer simulations and spectral decomposition have revealed how the resulting ECD patterns of the naphthyl glycosides studied are influenced by different conformer populations. In conclusion, reliable ECD patterns cannot be calculated by separating the naphthyl rotation from other conformational motions. MCD patterns have been similar for all the naphthyl C-glycosides studied. No clear diagnostic features have been found for either the pseudoanomeric configuration or the preferred hydroxymethyl rotamer. Nevertheless, the work has demonstrated the potential of MCD for the study of aryl glycosides interacting with proteins.

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.

Efficiently Computing NMR 1H and 13C Chemical Shifts of Saccharides in Aqueous Environment

Determining the structure of saccharides in their native environment is crucial to understanding their function and more accurately targeting their utilization. Nuclear magnetic resonance observables such as the nuclear Overhauser effect or spin-spin coupling constants are routinely utilized to study saccharides in their native water environment. However, while highly sensitive to the local environment, chemical shifts are mostly overlooked, despite being commonly measured for compounds identification. Although chemical shifts carry considerable structural information, their direct association with structure is notoriously difficult. This is mostly due to the similarity in the chemical nature of most saccharides causing similar physicochemical environments close to sugar C and H atoms, resulting in comparable chemical shifts. The rise of computational power allows one to compute reliable chemical shifts and use them to determine atomistic details of these sugars in solution. However, any prediction is severely limited by the computational protocol used and its accuracy. In this work, we studied a set of 31 saccharides on which we evaluated various computational protocols to calculate the total number of 375 ¹H and 327 ¹³C chemical shifts of sugars in an aqueous environment. Our study proposes two cost-effective protocols for simulating ¹H and ¹³C chemical shifts that we recommend for further use. These protocols can help with the interpretation of experimental spectra, but we also show that they are also capable of structure prediction independently. This is possible because of the low mean absolute deviations of calculated shifts from the experiment (0.06 ppm for ¹H and 1.09 ppm for ¹³C). We explore different solvation methods, basis sets, and optimization schemes to reach such accuracy. A correct sampling of the conformation phase space of flexible sugar molecules is also key to obtaining accurately converged theoretical chemical shifts. The linear regression method was applied to convert the calculated isotropic nuclear magnetic shielding constants to simulated chemical shifts comparable with the experiment. The achieved level of accuracy can help in utilizing chemical shifts for elucidating the 3D atomistic structure of saccharides in aqueous solutions. All linear regression parameters obtained on our extensive set of sugars for all the tested protocols can be reutilized in future works.

Unlocking the Hydrolytic Mechanism of GH92 α‐1,2‐Mannosidases: Computation Inspires the use of C‐Glycosides as Michaelis Complex Mimics

The conformational changes in a sugar moiety along the hydrolytic pathway are key to understand the mechanism of glycoside hydrolases (GHs) and to design new inhibitors. The two predominant itineraries for mannosidases go via ^OS₂→B_2,5→¹S₅ and ³S₁→³H₄→¹C₄. For the CAZy family 92, the conformational itinerary was unknown. Published complexes of Bacteroides thetaiotaomicron GH92 catalyst with a S‐glycoside and mannoimidazole indicate a ⁴C₁→⁴H₅/¹S₅→¹S₅ mechanism. However, as observed with the GH125 family, S‐glycosides may not act always as good mimics of GH's natural substrate. Here we present a cooperative study between computations and experiments where our results predict the E₅→B_2,5/¹S₅→¹S₅ pathway for GH92 enzymes. Furthermore, we demonstrate the Michaelis complex mimicry of a new kind of C‐disaccharides, whose biochemical applicability was still a chimera.

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.