Investigation of the structure of regulatory proteins interacting with glycosaminoglycans by combining NMR spectroscopy and molecular modeling - the beginning of a wonderful friendship

A biological guide to glycosaminoglycans: current perspectives and pending questions

Mammalian glycosaminoglycans (GAGs), except hyaluronan (HA), are sulfated polysaccharides that are covalently attached to core proteins to form proteoglycans (PGs). This article summarizes key biological findings for the most widespread GAGs, namely HA, chondroitin sulfate/dermatan sulfate (CS/DS), keratan sulfate (KS), and heparan sulfate (HS). It focuses on the major processes that remain to be deciphered to get a comprehensive view of the mechanisms mediating GAG biological functions. They include the regulation of GAG biosynthesis and postsynthetic modifications in heparin (HP) and HS, the composition, heterogeneity, and function of the tetrasaccharide linkage region and its role in disease, the functional characterization of the new PGs recently identified by glycoproteomics, the selectivity of interactions mediated by GAG chains, the display of GAG chains and PGs at the cell surface and their impact on the availability and activity of soluble ligands, and on their move through the glycocalyx layer to reach their receptors, the human GAG profile in health and disease, the roles of GAGs and particular PGs (syndecans, decorin, and biglycan) involved in cancer, inflammation, and fibrosis, the possible use of GAGs and PGs as disease biomarkers, and the design of inhibitors targeting GAG biosynthetic enzymes and GAG-protein interactions to develop novel therapeutic approaches.

Sweet but Challenging: Tackling the Complexity of GAGs with Engineered Tailor-Made Biomaterials

Glycosaminoglycans (GAGs) play a crucial role in tissue homeostasis by regulating the activity and diffusion of bioactive molecules. Incorporating GAGs into biomaterials has emerged as a widely adopted strategy in medical applications, owing to their biocompatibility and ability to control the release of bioactive molecules. Nevertheless, immobilized GAGs on biomaterials can elicit distinct cellular responses compared to their soluble forms, underscoring the need to understand the interactions between GAG and bioactive molecules within engineered functional biomaterials. By controlling critical parameters such as GAG type, density, and sulfation, it becomes possible to precisely delineate GAG functions within a biomaterial context and to better mimic specific tissue properties, enabling tailored design of GAG-based biomaterials for specific medical applications. However, this requires access to pure and well-characterized GAG compounds, which remains challenging. This review focuses on different strategies for producing well-defined GAGs and explores high-throughput approaches employed to investigate GAG-growth factor interactions and to quantify cellular responses on GAG-based biomaterials. These automated methods hold considerable promise for improving the understanding of the diverse functions of GAGs. In perspective, the scientific community is encouraged to adopt a rational approach in designing GAG-based biomaterials, taking into account the in vivo properties of the targeted tissue for medical applications.

Molecular Dynamics Simulation-Based Prediction of Glycosaminoglycan Interactions with Drug Molecules

Glycosaminoglycans (GAGs) are a class of long linear anionic periodic polysaccharides. Their biological activities are very broad including tissue remodeling, regulation of cell proliferation, cell migration, cell differentiation, participation in bacterial/viral infections, and immune response. They can interact with many important biomolecular partners in the extracellular matrix of the cell including small drug molecules. Recently, several GAG-bioactive small molecule complexes have been experimentally and theoretically studied. Some of these compounds in complexes with GAGs may potentially interfere with protein-GAG or peptide-GAG multimolecular systems affecting the processes of cellular differentiation or have anti-inflammatory, antiviral as well as antithrombotic effects. Although many studies have been conducted on GAG-drug complexes, the molecular mechanisms of the formation of such complexes are still poorly understood. At the same time, the complexity of their physicochemical properties renders the use of both experimental and computational methods to study these molecular systems challenging. Here, we present the molecular dynamics-based protocols successfully employed to in silico analyze GAG-small molecule interactions.

Studying specificity in protein-glycosaminoglycan recognition with umbrella sampling

In the past few decades, glycosaminoglycan (GAG) research has been crucial for gaining insights into various physiological, pathological, and therapeutic aspects mediated by the direct interactions between the GAG molecules and diverse proteins. The structural and functional heterogeneities of GAGs as well as their ability to bind specific proteins are determined by the sugar composition of the GAG, the size of the GAG chains, and the degree and pattern of sulfation. A deep understanding of the interactions in protein-GAG complexes is essential to explain their biological functions. In this study, the umbrella sampling (US) approach is used to pull away a GAG ligand from the binding site and then pull it back in. We analyze the binding interactions between GAGs of three types (heparin, desulfated heparan sulfate, and chondroitin sulfate) with three different proteins (basic fibroblast growth factor, acidic fibroblast growth factor, and cathepsin K). The main focus of our study was to evaluate whether the US approach is able to reproduce experimentally obtained structures, and how useful it can be for getting a deeper understanding of GAG properties, especially protein recognition specificity and multipose binding. We found that the binding free energy landscape in the proximity of the GAG native binding pose is complex and implies the co-existence of several binding poses. The sliding of a GAG chain along a protein surface could be a potential mechanism of GAG particular sequence recognition by proteins.

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.

Extension of the SUGRES-1P Coarse-Grained Model of Polysaccharides to Heparin

Heparin is an unbranched periodic polysaccharide composed of negatively charged monomers and involved in key biological processes, including anticoagulation, angiogenesis, and inflammation. Its structure and dynamics have been studied extensively using experimental as well as theoretical approaches. The conventional approach of computational chemistry applied to the analysis of biomolecules is all-atom molecular dynamics, which captures the interactions of individual atoms by solving Newton's equation of motion. An alternative is molecular dynamics simulations using coarse-grained models of biomacromolecules, which offer a reduction of the representation and consequently enable us to extend the time and size scale of simulations by orders of magnitude. In this work, we extend the UNIfied COarse-gRaiNed (UNICORN) model of biological macromolecules developed in our laboratory to heparin. We carried out extensive tests to estimate the optimal weights of energy terms of the effective energy function as well as the optimal Debye-Hückel screening factor for electrostatic interactions. We applied the model to study unbound heparin molecules of polymerization degree ranging from 6 to 68 residues. We compare the obtained coarse-grained heparin conformations with models obtained from X-ray diffraction studies of heparin. The SUGRES-1P force field was able to accurately predict the general shape and global characteristics of heparin molecules.

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.

Understanding Conformational Properties and Role of Hydrogen Bonds in Glycosaminoglycans-Interleukin8 Complexes in Aqueous Medium by Molecular Dynamics Simulation

Atomistic molecular dynamics simulations were performed under ambient conditions to explore the conformational features and binding affinities of hexameric glycosaminoglycans (GAGs) with chemokine Interleukin8 (IL8) in an aqueous medium. We tried to understand the role of hydrogen bonds (HBs) involving conserved water in mediating the interactions. The Luzar-Chandler model was adopted to study the kinetics of HB breaking and formation concerning different water-mediated HBs. The conformational flexibilities of bound GAGs are due to the flexible glycosidic linkages than the occasional/rare ring pucker conformation. The free energy landscape constructed with ϕ, and ψ, depicted that different conformational minima associated with the glycosidic linkage flexibility of the GAGs in bound states are separated by energy barriers. The binding affinities of IL8 towards GAGs are favored through the electrostatic and non-polar solvation interactions. 4-different types of conserved water were explored in the solvent-mediated binding of GAGs with IL8. The average lifetime of the IL8-GAG direct HB pairs was ∼ten times less than the IL8-GAG-shared water HBs. This is due to the rapid establishment of HB breaking and reformation kinetics involving water of a shared layer. We find that despite the highly negatively charged surface of GAGs, the IL8 surface populated by non-cationic amino acids could serve as a promising binding site in addition to the cationic surface of the protein.

Molecular Dynamics-Based Comparative Analysis of Chondroitin and Dermatan Sulfates

Glycosaminoglycans (GAGs) are a class of linear anionic periodic polysaccharides containing disaccharide repetitive units. These molecules interact with a variety of proteins in the extracellular matrix and so participate in biochemically crucial processes such as cell signalling affecting tissue regeneration as well as the onset of cancer, Alzheimer's or Parkinson's diseases. Due to their flexibility, periodicity and chemical heterogeneity, often termed "sulfation code", GAGs are challenging molecules both for experiments and computation. One of the key questions in the GAG research is the specificity of their intermolecular interactions. In this study, we make a step forward to deciphering the "sulfation code" of chondroitin sulfates-4,6 (CS4, CS6, where the numbers correspond to the position of sulfation in NAcGal residue) and dermatan sulfate (DS), which is different from CSs by the presence of IdoA acid instead of GlcA. We rigorously investigate two sets of these GAGs in dimeric, tetrameric and hexameric forms with molecular dynamics-based descriptors. Our data clearly suggest that CS4, CS6 and DS are substantially different in terms of their structural, conformational and dynamic properties, which contributes to the understanding of how these molecules can be different when they bind proteins, which could have practical implications for the GAG-based drug design strategies in the regenerative medicine.

Molecular Dynamics Approaches Dissect Molecular Mechanisms Underlying Methylene Blue-Glycosaminoglycan Interactions

Glycosaminoglycans (GAGs) are a class of periodic anionic linear polysaccharides involved in a number of biologically relevant processes in the extracellular matrix via interactions with various types of molecules including proteins, peptides and small organic molecules. The metachromatic dye methylene blue (MB) is a GAG binding agent. This molecule possesses a tricyclic, monocationic phenothiazine ring system, while the terminal methyl groups attached to the nitrogen atoms bear the most positive charges of the cation and, therefore, represent potential binding sites for negatively charged GAGs. In this study, we rigorously explored molecular mechanisms underlying these interactions for several GAG types: heparin, heparan and chondroitin sulfates. We found that GAG-MB interactions are predominantly electrostatically driven, with the particularly important role of sulfate groups. MB oligomeric stack formation was favored in the presence of GAGs. Furthermore, the impact of MB binding on the conformation of GAGs was also evaluated. The novel results allow for better quantitative analytics of GAG composition in the studied biochemical systems using MB dye as a GAG-specific marker. Our data add to the knowledge on small molecule-GAG interactions and could be potentially useful for novel developments in drug design and putative disease therapies in which GAGs are involved.

Pentosan Polysulphate (PPS), a Semi-Synthetic Heparinoid DMOAD With Roles in Intervertebral Disc Repair Biology emulating The Stem Cell Instructive and Tissue Reparative Properties of Heparan Sulphate

This review highlights the attributes of pentosan polysulphate (PPS) in the promotion of intervertebral disc (IVD) repair processes. PPS has been classified as a disease modifying osteoarthritic drug (DMOAD) and many studies have demonstrated its positive attributes in the countering of degenerative changes occurring in cartilaginous tissues during the development of osteoarthritis (OA). Degenerative changes in the IVD also involve inflammatory cytokines, degradative proteases and cell signalling pathways similar to those operative in the development of OA in articular cartilage. PPS acts as a heparan sulphate (HS) mimetic to effect its beneficial effects in cartilage. The IVD contains small cell membrane HS-proteoglycans (HSPGs) such as syndecan, and glypican and a large multifunctional HS/chondroitin sulphate (CS) hybrid proteoglycan (HSPG2/perlecan) that have important matrix stabilising properties and sequester, control and present growth factors from the FGF, VEGF, PDGF and BMP families to cellular receptors to promote cell proliferation, differentiation and matrix synthesis. HSPG2 also has chondrogenic properties and stimulates the synthesis of extracellular matrix (ECM) components, expansion of cartilaginous rudiments and has roles in matrix stabilisation and repair. Perlecan is a perinuclear and nuclear proteoglycan in IVD cells with roles in chromatin organisation and control of transcription factor activity, immunolocalises to stem cell niches in cartilage, promotes escape of stem cells from quiescent recycling, differentiation and attainment of pluripotency and migratory properties. These participate in tissue development and morphogenesis, ECM remodelling and repair. PPS also localises in the nucleus of stromal stem cells, promotes development of chondroprogenitor cell lineages, ECM synthesis and repair and discal repair by resident disc cells. The availability of recombinant perlecan and PPS offer new opportunities in repair biology. These multifunctional agents offer welcome new developments in repair strategies for the IVD.