Conformational Analysis of Heparin-Analogue Pentasaccharides by Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulations

Synthesis of a Heparinoid Pentasaccharide Containing l-Guluronic Acid Instead of l-Iduronic Acid with Preserved Anticoagulant Activity

l-Iduronic acid is a key constituent of heparin and heparan sulfate polysaccharides due to its unique conformational plasticity, which facilitates the binding of polysaccharides to proteins. At the same time, this is the synthetically most challenging unit of heparinoid oligosaccharides; therefore, there is a high demand for its replacement with a more easily accessible sugar unit. In the case of idraparinux, an excellent anticoagulant heparinoid pentasaccharide, we demonstrated that l-iduronic acid can be replaced by an easier-to-produce l-sugar while maintaining its essential biological activity. From the inexpensive d-mannose, through a highly functionalized phenylthio mannoside, the l-gulose donor was prepared by C-5 epimerization in 10 steps with excellent yield. This unit was incorporated into the pentasaccharide by α-selective glycosylation and oxidized to l-guluronic acid. The complete synthesis required only 36 steps, with 21 steps for the longest linear route. The guluronate containing pentasaccharide inhibited coagulation factor Xa by 50% relative to the parent compound, representing an excellent anticoagulant activity. To the best of our knowledge, this is the first biologically active heparinoid anticoagulant which contains a different sugar unit instead of l-iduronic acid.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Effect of Sulfation on the Conformational Dynamics of Dermatan Sulfate Glycosaminoglycan: A Gaussian Accelerated Molecular Dynamics Study

Glycosaminoglycans (GAGs) are anionic biopolymers present on cell surfaces as a part of proteoglycans. The biological activities of GAGs depend on the sulfation pattern. In our study, we have considered three octadecasaccharide dermatan sulfate (DS) chains with increasing order of sulfation (dp6s, dp7s, and dp12s) to illuminate the role of sulfation on the GAG units and its chain conformation through 10 μs-long Gaussian accelerated molecular dynamics simulations. DS is composed of repeating disaccharide units of iduronic acid (IdoA) and N-acetylgalactosamine (N-GalNAc). Here, N-GalNAc is linked to IdoA via β(1-4), while IdoA is linked to N-GalNAc through α(1-3). With the increase in sulfation, the DS structure becomes more rigid and linear, as is evident from the distribution of root-mean-square deviations (RMSDs) and end-to-end distances. The tetrasaccharide linker region of the main chain shows a rigid conformation in terms of the glycosidic linkage. We have observed that upon sulfation (i.e., dp12s), the ring flip between two chair forms vanished for IdoA. The dynamic cross-correlation analysis reveals that the anticorrelation motions in dp12s are reduced significantly compared to dp6s or dp7s. An increase in sulfation generates relatively more stable hydrogen-bond networks, including water bridging with the neighboring monosaccharides. Despite the favorable linear structures of the GAG chains, our study also predicts few significant bendings related to the different puckering states, which may play a notable role in the function of the DS. The relation between the global conformation with the micro-level parameters such as puckering and water-mediated hydrogen bonds shapes the overall conformational space of GAGs. Overall, atomistic details of the DS chain provided in this study will help understand their functional and mechanical roles, besides developing new biomaterials.

Characterization of Heparin's Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin's solution conformational dynamics and structure-function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule's underlying conformational ensemble, which then provides important information necessary for understanding structure-function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin's conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. The anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

Investigation of the Molecular Details of the Interactions of Selenoglycosides and Human Galectin-3

Human galectin-3 (hGal-3) is involved in a variety of biological processes and is implicated in wide range of diseases. As a result, targeting hGal-3 for clinical applications has become an intense area of research. As a step towards the development of novel hGal-3 inhibitors, we describe a study of the binding of two Se-containing hGal-3 inhibitors, specifically that of di(β-D-galactopyranosyl)selenide (SeDG), in which two galactose rings are linked by one Se atom and a di(β-D-galactopyranosyl)diselenide (DSeDG) analogue with a diseleno bond between the two sugar units. The binding affinities of these derivatives to hGal-3 were determined by ¹⁵N-¹H HSQC NMR spectroscopy and fluorescence anisotropy titrations in solution, indicating a slight decrease in the strength of interaction for SeDG compared to thiodigalactoside (TDG), a well-known inhibitor of hGal-3, while DSeDG displayed a much weaker interaction strength. NMR and FA measurements showed that both seleno derivatives bind to the canonical S face site of hGal-3 and stack against the conserved W181 residue also confirmed by X-ray crystallography, revealing canonical properties of the interaction. The interaction with DSeDG revealed two distinct binding modes in the crystal structure which are in fast exchange on the NMR time scale in solution, explaining a weaker interaction with hGal-3 than SeDG. Using molecular dynamics simulations, we have found that energetic contributions to the binding enthalpies mainly differ in the electrostatic interactions and in polar solvation terms and are responsible for weaker binding of DSeDG compared to SeDG. Selenium-containing carbohydrate inhibitors of hGal-3 showing canonical binding modes offer the potential of becoming novel hydrolytically stable scaffolds for a new class of hGal-3 inhibitors.