13C-Labeled Idohexopyranosyl Rings: Effects of Methyl Glycosidation and C6 Oxidation on Ring Conformational Equilibria

Conformational Preferences at the Glycosidic Linkage of Saccharides in Solution as Deduced from NMR Experiments and MD Simulations: Comparison to Crystal Structures

Glycans are central to information content and regulation in biological systems. These carbohydrate molecules are active either as oligo- or polysaccharides, often in the form of glycoconjugates. The monosaccharide entities are joined by glycosidic linkages and stereochemical arrangements are of utmost importance in determining conformation and flexibility of saccharides. The conformational preferences and population distributions at the glycosidic torsion angles φ and ψ have been investigated for O-methyl glycosides of three disaccharides where the substitution takes place at a secondary alcohol, viz., in α-l-Fucp-(1→3)-β-d-Glcp-OMe, α-l-Fucp-(1→3)-α-d-Galp-OMe and α-d-Glcp-(1→4)-α-d-Galp-OMe, corresponding to disaccharide structural elements present in bacterial polysaccharides. Stereochemical differences at or adjacent to the glycosidic linkage were explored by solution state NMR spectroscopy using one-dimensional 1 H,1 H-NOESY NMR experiments to obtain transglycosidic proton-proton distances and one- and two-dimensional heteronuclear NMR experiments to obtain 3 JCH transglycosidic coupling constants related to torsion angles φ and ψ. Computed effective proton-proton distances from molecular dynamics (MD) simulations showed excellent agreement to experimentally derived distances for the α-(1→3)-linked disaccharides and revealed that for the bimodal distribution at the ψ torsion angle for the α-(1→4)-linked disaccharide experiment and simulation were at variance with each other, calling for further force field developments. The MD simulations disclosed a highly intricate inter-residue hydrogen bonding pattern for the α-(1→4)-linked disaccharide, including a nonconventional hydrogen bond between H5' in the glucosyl residue and O3 in the galactosyl residue, supported by a large downfield 1 H NMR chemical shift displacement compared to α-d-Glcp-OMe. Comparison of population distributions of the glycosidic torsion angles φ and ψ in the disaccharide entities to those of corresponding crystal structures highlighted the potential importance of solvation on the preferred conformation.

Glycan Shape, Motions, and Interactions Explored by NMR Spectroscopy

Glycans in the form of oligosaccharides, polysaccharides, and glycoconjugates are ubiquitous in nature, and their structures range from linear assemblies to highly branched and decorated constructs. Solution state NMR spectroscopy facilitates elucidation of preferred conformations and shapes of the saccharides, motions, and dynamic aspects related to processes over time as well as the study of transient interactions with proteins. Identification of intermolecular networks at the atomic level of detail in recognition events by carbohydrate-binding proteins known as lectins, unraveling interactions with antibodies, and revealing substrate scope and action of glycosyl transferases employed for synthesis of oligo- and polysaccharides may efficiently be analyzed by NMR spectroscopy. By utilizing NMR active nuclei present in glycans and derivatives thereof, including isotopically enriched compounds, highly detailed information can be obtained by the experiments. Subsequent analysis may be aided by quantum chemical calculations of NMR parameters, machine learning-based methodologies and artificial intelligence. Interpretation of the results from NMR experiments can be complemented by extensive molecular dynamics simulations to obtain three-dimensional dynamic models, thereby clarifying molecular recognition processes involving the glycans.

Exploring Ring Conformation in Uronate Monosaccharides: Insights from Ab Initio Calculations and Classical Molecular Dynamics Simulations

The study focuses on the conformational properties of biologically relevant monosaccharides belonging to the group of uronates: α-l-iduronate, O2-sulfated-α-l-iduronate, and O2-sulfated-α-l-guluronate, either unfunctionalized or O1-methylated. We applied the previously proposed two-step methodology, combining classical MD simulations and subsequent ab initio (QM) calculations, performed on a rationally subsampled set of molecular configurations. We found that, regardless of the number of molecular configurations considered, the level of theory, and the weighting scheme applied, none of the QM approaches is capable of predicting the correct conformational equilibrium of sulfated iduronates as long as the tight counterion binding is not considered. Multicenter, ring-shape-specific binding of either Na+ or Ca2+ ions drastically shifts the conformational equilibrium of the pyranose ring in sulfated iduronates toward 1C4 but does not significantly affect the conformation of non-sulfated compounds. A similar shift is observed upon the protonation of carboxyl groups in all iduronates. In addition, we report a set of average J-coupling constant values related to vicinal protons bound to the pyranose ring of iduronates and corresponding to each of the three main groups of ring conformers, i.e., 4C1, B/S (boat/skew boat), and 1C4. In combination with the conformational energies or with the experimental data, these values allowed the relative proportions of the ring conformers to be estimated and the Karplus-type equations linking the 3JHH-coupling constants to the torsion angles within the pyranose ring to be refined.

Impact of Polarization on the Ring Puckering Dynamics of Hexose Monosaccharides

Analysis of crystal structures of hexose monosaccharides α-d-mannose (α-MAN), β-d-mannose (β-MAN), α-d-glucose (α-GLC), β-d-glucose (β-GLC), α-d-galactose (α-GAL), β-d-galactose (β-GAL), α-d-altrose (α-ALT), β-d-altrose (β-ALT), α-d-idose (α-IDO), and β-d-idose (β-IDO) reveals that the monosaccharide ring adopts multiple ring conformations. These ring conformations can be broadly classified as chair, half-chair, envelope, boat, and skew-boat conformations. The ability of the monosaccharide ring to adopt multiple conformations has been closely tied with their bioactivity. However, it has been difficult to capture the dynamic information of these conformations from experimental studies. Even from simulations, capturing these different conformations is challenging because of the energy barriers involved in the transitions between the stable ⁴C₁ and ¹C₄ chair forms. In this study, we analyze the influence of the polarizable force field on the ring dynamics of five major types of unsubstituted aldohexoses─glucose, mannose, galactose, altrose, and idose─and their anomers. We simulate microsecond trajectories to capture the influence of the CHARMM36 additive and polarizable carbohydrate force fields on the ring dynamics. The microsecond trajectories allow us to comment on the issues associated with equilibrium molecular dynamics simulations. Further, we use the extended system adaptive biasing force (eABF) method to compare the conformational sampling efficiencies of the additive and polarizable force fields. Our studies reveal that inclusion of polarization enhances the sampling of ring conformations and lowers the energy barriers between the ⁴C₁ and ¹C₄ conformations. Overall, the CHARMM36 additive force field is observed to be rigid and favor the ⁴C₁ conformations. Although the inclusion of polarizability results in enhancing ring flexibility, we observe sampling that does not agree with experimental results, warranting a revision of the polarizable Drude parameters.

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.

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.

Pyranose Ring Puckering Thermodynamics for Glycan Monosaccharides Associated with Vertebrate Proteins

The conformational properties of carbohydrates can contribute to protein structure directly through covalent conjugation in the cases of glycoproteins and proteoglycans and indirectly in the case of transmembrane proteins embedded in glycolipid-containing bilayers. However, there continue to be significant challenges associated with experimental structural biology of such carbohydrate-containing systems. All-atom explicit-solvent molecular dynamics simulations provide a direct atomic resolution view of biomolecular dynamics and thermodynamics, but the accuracy of the results depends on the quality of the force field parametrization used in the simulations. A key determinant of the conformational properties of carbohydrates is ring puckering. Here, we applied extended system adaptive biasing force (eABF) all-atom explicit-solvent molecular dynamics simulations to characterize the ring puckering thermodynamics of the ten common pyranose monosaccharides found in vertebrate biology (as represented by the CHARMM carbohydrate force field). The results, along with those for idose, demonstrate that the CHARMM force field reliably models ring puckering across this diverse set of molecules, including accurately capturing the subtle balance between 4C1 and 1C4 chair conformations in the cases of iduronate and of idose. This suggests the broad applicability of the force field for accurate modeling of carbohydrate-containing vertebrate biomolecules such as glycoproteins, proteoglycans, and glycolipids.

The conformation of the idopyranose ring revisited: How subtle O-substituent induced changes can be deduced from vicinal 1H-NMR coupling constants

The idopyranose ring plays a pivotal role in the conformational, dynamical, and intermolecular binding aspects of glycosaminoglycans like heparin and dermatan sulfate and it was early on assigned a role in the Sugar Code governing biological recognition processes. There is consensus that next to the two canonical ¹C₄ and ⁴C₁ chair conformations, the conformational space accessible to the idopyranose ring entails a ²S_O skew-boat conformation, but the equilibrium between these three ring puckers has evaded satisfactory quantification. In this study a meta-analysis of X-ray solid-state data and vicinal NMR coupling constants is presented, based on the Truncated Fourier Puckering (TFP) formalism and the generalized Karplus (CAGPLUS) equation. This approach yields a model-free, granular and consistent reckoning of 159 idopyranose solution puckering equilibria studied by NMR and allows us to reproduce the involved 636 NMR vicinal couplings with an overall residual RMS(J_obs-J_calc) of 0.184 Hz. Our analyses show that for all ring systems examined, the idopyranosyl chair conformations take up the same ring pucker irrespective of the ring substituent pattern or a vast variety in experimental conditions. Instead, it is the (skew-)boat conformation that adapts to the substitution pattern of the idopyranose ring or a specific sulfation pattern of neighboring saccharides. All idopyranose rings are involved in conformational equilibria that subsume the aforementioned conformers which turn out to differ only a few kJ/mole in conformational energy. Thus, the plasticity and flexibility of idopyranose remains intact under practically all circumstances and, as the glycosidic linkages in heparin are considered to be relatively stiff, the iduronic moiety functions as the linchpin of heparin flexibility thereby being rather a "space(r)" than a "letter" in the alleged Sugar Code alphabet.

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.

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.

Transition-Metal-Free Synthesis of C-Glycosylated Phenanthridines via K2S2O8-Mediated Oxidative Radical Decarboxylation of Uronic Acids

We have developed an efficient protocol for the synthesis of C-glycosylated phenanthridines. Tetrafuranos-4-yl and pentapyranos-5-yl radicals, generated from K₂S₂O₈-mediated oxidative decarboxylation of furan- and pyranuronic acids, undergo attack to 2-isocyanodiphenyls and ensuing homolytic aromatic substitution to provide diverse C-glycosylated phenanthridines in satisfactory yields without resort to transition metals. This reaction tolerates various functional groups, and enables ready synthesis of complex oligosaccharide-based phenanthridines. The C-glycosylated phenanthridine derived from β-cyclodextrin has been prepared, which might be potential in medicinal and biological chemistry due to its flexible conformation.