A Bittersweet Computational Journey among Glycosaminoglycans

Developing and Benchmarking Sulfate and Sulfamate Force Field Parameters via Ab Initio Molecular Dynamics Simulations To Accurately Model Glycosaminoglycan Electrostatic Interactions

Glycosaminoglycans (GAGs) are negatively charged polysaccharides found on cell surfaces, where they regulate transport pathways of foreign molecules toward the cell. The structural and functional diversity of GAGs is largely attributed to varied sulfation patterns along the polymer chains, which makes understanding their molecular recognition mechanisms crucial. Molecular dynamics (MD) simulations, thanks to their unmatched microscopic resolution, have the potential to be a reference tool for exploring the patterns responsible for biologically relevant interactions. However, the capability of molecular dynamics force fields used in biosimulations to accurately capture sulfation-specific interactions is not well established, partly due to the intrinsic properties of GAGs that pose challenges for most experimental techniques. In this work, we evaluate the performance of molecular dynamics force fields for sulfated GAGs by studying ion pairing of Ca2+ to sulfated moieties─N-methylsulfamate and methylsulfate─that resemble N- and O-sulfation found in GAGs, respectively. We tested available nonpolarizable (CHARMM36 and GLYCAM06) and explicitly polarizable (Drude and AMOEBA) force fields, and derived new implicitly polarizable models through charge scaling (prosECCo75 and GLYCAM-ECC75) that are consistent with our developed "charge-scaling" framework. The calcium-sulfamate/sulfate interaction free energy profiles obtained with the tested force fields were compared against reference ab initio molecular dynamics (AIMD) simulations, which serve as a robust alternative to experiments. AIMD simulations indicate that the preferential Ca2+ binding mode to sulfated GAG groups is solvent-shared pairing. Only our scaled-charge models agree satisfactorily with the AIMD data, while all other force fields exhibit poorer agreement, sometimes even qualitatively. Surprisingly, even explicitly polarizable force fields display a notable disagreement with the AIMD data, likely attributed to difficulties in their optimization and possible inherent limitations in depicting high-charge-density ion interactions accurately. Finally, the underperforming force fields lead to unrealistic aggregation of sulfated saccharides, which qualitatively disagrees with our understanding of the soft glycocalyx environment. Our results highlight the importance of accurately treating electronic polarization in MD simulations of sulfated GAGs and caution against over-reliance on currently available models without thorough validation and optimization.

Structural Insights into Endostatin-Heparan Sulfate Interactions Using Modeling Approaches

Glycosaminoglycans (GAGs) play a key role in a variety of biological processes in the extracellular matrix (ECM) via interactions with their protein targets. Due to their high flexibility, periodicity and electrostatics-driven interactions, GAG-containing complexes are very challenging to characterize both experimentally and in silico. In this study, we, for the first time, systematically analyzed the interactions of endostatin, a proteolytic fragment of collagen XVIII known to be anti-angiogenic and anti-tumoral, with heparin (HP) and representative heparan sulfate (HS) oligosaccharides of various lengths, sequences and sulfation patterns. We first used conventional molecular docking and a docking approach based on a repulsive scaling-replica exchange molecular dynamics technique, as well as unbiased molecular dynamic simulations, to obtain dynamically stable GAG binding poses. Then, the corresponding free energies of binding were calculated and the amino acid residues that contribute the most to GAG binding were identified. We also investigated the potential influence of Zn2+ on endostatin-HP complexes using computational approaches. These data provide new atomistic details of the molecular mechanism of HP's binding to endostatin, which will contribute to a better understanding of its interplay with proteoglycans at the cell surface and in the extracellular matrix.

Ligand binding of interleukin-8: a comparison of glycosaminoglycans and acidic peptides

Recognition and binding of regulatory proteins to glycosaminoglycans (GAGs) from the extracellular matrix is a process of high biological importance. The interaction between negatively charged sulfate or carboxyl groups of the GAGs and clusters of basic amino acids on the protein is crucial in this binding process and it is believed that electrostatics represent the key factor for this interaction. However, given the rather undirected nature of electrostatics, it is important to achieve a clear understanding of its role in protein-GAG interactions and how specificity and selectivity in these systems can be achieved, when the classical key-lock binding motif is not applicable. Here, we compare protein binding of a highly charged heparin (HP) hexasaccharide with four de novo designed decapeptides of varying negative net charge. The charge density of these peptides was comparable to typical GAGs of the extracellular matrix. We used the regulatory protein interleukin-8 (IL-8) because its interactions with GAGs are well described. All four peptide ligands bind to the same epitope of IL-8 but show much weaker binding affinity as revealed in 1H-15N HSQC NMR titration experiments. Complementary molecular docking and molecular dynamics simulations revealed further atomistic details of the interaction mode of GAG versus peptide ligands. Overall, similar contributions to the binding energy and hydrogen bond formation are determined for HP and the highly charged peptides, suggesting that the entropic loss of the peptides upon binding likely account for the remarkably different affinity of GAG versus peptide ligands to IL-8.

Glycosaminoglycans: What Remains To Be Deciphered?

Glycosaminoglycans (GAGs) are complex polysaccharides exhibiting a vast structural diversity and fulfilling various functions mediated by thousands of interactions in the extracellular matrix, at the cell surface, and within the cells where they have been detected in the nucleus. It is known that the chemical groups attached to GAGs and GAG conformations comprise "glycocodes" that are not yet fully deciphered. The molecular context also matters for GAG structures and functions, and the influence of the structure and functions of the proteoglycan core proteins on sulfated GAGs and vice versa warrants further investigation. The lack of dedicated bioinformatic tools for mining GAG data sets contributes to a partial characterization of the structural and functional landscape and interactions of GAGs. These pending issues will benefit from the development of new approaches reviewed here, namely (i) the synthesis of GAG oligosaccharides to build large and diverse GAG libraries, (ii) GAG analysis and sequencing by mass spectrometry (e.g., ion mobility-mass spectrometry), gas-phase infrared spectroscopy, recognition tunnelling nanopores, and molecular modeling to identify bioactive GAG sequences, biophysical methods to investigate binding interfaces, and to expand our knowledge and understanding of glycocodes governing GAG molecular recognition, and (iii) artificial intelligence for in-depth investigation of GAGomic data sets and their integration with proteomics.

Pharmacology of Heparin and Related Drugs: An Update

Heparin has been used extensively as an antithrombotic and anticoagulant for close to 100 years. This anticoagulant activity is attributed mainly to the pentasaccharide sequence, which potentiates the inhibitory action of antithrombin, a major inhibitor of the coagulation cascade. More recently it has been elucidated that heparin exhibits anti-inflammatory effect via interference of the formation of neutrophil extracellular traps and this may also contribute to heparin's antithrombotic activity. This illustrates that heparin interacts with a broad range of biomolecules, exerting both anticoagulant and nonanticoagulant actions. Since our previous review, there has been an increased interest in these nonanticoagulant effects of heparin, with the beneficial role in patients infected with SARS2-coronavirus a highly topical example. This article provides an update on our previous review with more recent developments and observations made for these novel uses of heparin and an overview of the development status of heparin-based drugs. SIGNIFICANCE STATEMENT: This state-of-the-art review covers recent developments in the use of heparin and heparin-like materials as anticoagulant, now including immunothrombosis observations, and as nonanticoagulant including a role in the treatment of SARS-coronavirus and inflammatory conditions.

HS, an Ancient Molecular Recognition and Information Storage Glycosaminoglycan, Equips HS-Proteoglycans with Diverse Matrix and Cell-Interactive Properties Operative in Tissue Development and Tissue Function in Health and Disease

Heparan sulfate is a ubiquitous, variably sulfated interactive glycosaminoglycan that consists of repeating disaccharides of glucuronic acid and glucosamine that are subject to a number of modifications (acetylation, de-acetylation, epimerization, sulfation). Variable heparan sulfate chain lengths and sequences within the heparan sulfate chains provide structural diversity generating interactive oligosaccharide binding motifs with a diverse range of extracellular ligands and cellular receptors providing instructional cues over cellular behaviour and tissue homeostasis through the regulation of essential physiological processes in development, health, and disease. heparan sulfate and heparan sulfate-PGs are integral components of the specialized glycocalyx surrounding cells. Heparan sulfate is the most heterogeneous glycosaminoglycan, in terms of its sequence and biosynthetic modifications making it a difficult molecule to fully characterize, multiple ligands also make an elucidation of heparan sulfate functional properties complicated. Spatio-temporal presentation of heparan sulfate sulfate groups is an important functional determinant in tissue development and in cellular control of wound healing and extracellular remodelling in pathological tissues. The regulatory properties of heparan sulfate are mediated via interactions with chemokines, chemokine receptors, growth factors and morphogens in cell proliferation, differentiation, development, tissue remodelling, wound healing, immune regulation, inflammation, and tumour development. A greater understanding of these HS interactive processes will improve therapeutic procedures and prognoses. Advances in glycosaminoglycan synthesis and sequencing, computational analytical carbohydrate algorithms and advanced software for the evaluation of molecular docking of heparan sulfate with its molecular partners are now available. These advanced analytic techniques and artificial intelligence offer predictive capability in the elucidation of heparan sulfate conformational effects on heparan sulfate-ligand interactions significantly aiding heparan sulfate therapeutics development.

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.

HIV-1 Tat and Heparan Sulfate Proteoglycans Orchestrate the Setup of in Cis and in Trans Cell-Surface Interactions Functional to Lymphocyte Trans-Endothelial Migration

HIV-1 transactivating factor Tat is released by infected cells. Extracellular Tat homodimerizes and engages several receptors, including integrins, vascular endothelial growth factor receptor 2 (VEGFR2) and heparan sulfate proteoglycan (HSPG) syndecan-1 expressed on various cells. By means of experimental cell models recapitulating the processes of lymphocyte trans-endothelial migration, here, we demonstrate that upon association with syndecan-1 expressed on lymphocytes, Tat triggers simultaneously the in cis activation of lymphocytes themselves and the in trans activation of endothelial cells (ECs). This "two-way" activation eventually induces lymphocyte adhesion and spreading onto the substrate and vascular endothelial (VE)-cadherin reorganization at the EC junctions, with consequent endothelial permeabilization, leading to an increased extravasation of Tat-presenting lymphocytes. By means of a panel of biochemical activation assays and specific synthetic inhibitors, we demonstrate that during the above-mentioned processes, syndecan-1, integrins, FAK, src and ERK_1/2 engagement and activation are needed in the lymphocytes, while VEGFR2, integrin, src and ERK_1/2 are needed in the endothelium. In conclusion, the Tat/syndecan-1 complex plays a central role in orchestrating the setup of the various in cis and in trans multimeric complexes at the EC/lymphocyte interface. Thus, by means of computational molecular modelling, docking and dynamics, we also provide a characterization at an atomic level of the binding modes of the Tat/heparin interaction, with heparin herein used as a structural analogue of the heparan sulfate chains of syndecan-1.

Advanced Molecular Dynamics Approaches to Model a Tertiary Complex APRIL/TACI with Long Glycosaminoglycans

Glycosaminoglycans (GAGs) are linear anionic periodic polysaccharides participating in a number of biologically relevant processes in the extracellular matrix via interactions with their protein targets. Due to their periodicity, conformational flexibility, pseudo-symmetry of the sulfation pattern, and the key role of electrostatics, these molecules are challenging for both experimental and theoretical approaches. In particular, conventional molecular docking applied for GAGs longer than 10-mer experiences severe difficulties. In this work, for the first time, 24- and 48-meric GAGs were docked using all-atomic repulsive-scaling Hamiltonian replica exchange molecular dynamics (RS-REMD), a novel methodology based on replicas with van der Waals radii of interacting molecules being scaled. This approach performed well for proteins complexed with oligomeric GAGs and is independent of their length, which distinguishes it from other molecular docking approaches. We built a model of long GAGs in complex with a proliferation-inducing ligand (APRIL) prebound to its receptors, the B cell maturation antigen and the transmembrane activator and calcium modulator and cyclophilin ligand interactor (TACI). Furthermore, the prediction power of the RS-REMD for this tertiary complex was evaluated. We conclude that the TACI–GAG interaction could be potentially amplified by TACI’s binding to APRIL. RS-REMD outperformed Autodock3, the docking program previously proven the best for short GAGs.