Multiscale modeling of glycosaminoglycan structure and dynamics: current methods and challenges

In silico approaches for better understanding cysteine cathepsin-glycosaminoglycan interactions

Cysteine cathepsins constitute the largest cathepsin family, with 11 proteases in human that are present primarily within acidic endosomal and lysosomal compartments. They are involved in the turnover of intracellular and extracellular proteins. They are synthesized as inactive procathepsins that are converted to mature active forms. Cathepsins play important roles in physiological and pathological processes and, therefore, receive increasing attention as potential therapeutic targets. Their maturation and activity can be regulated by glycosaminoglycans (GAGs), long linear negatively charged polysaccharides composed of recurring dimeric units. In this review, we summarize recent computational progress in the field of (pro)cathepsin-GAG complexes analyses.

Secondary deficiency of neuraminidase 1 contributes to CNS pathology in neurological mucopolysaccharidoses via hypersialylation of brain glycoproteins

Mucopolysaccharidoses (MPS) are lysosomal storage diseases caused by defects in catabolism of glycosaminoglycans. MPS I, II, III and VII are associated with lysosomal accumulation of heparan sulphate and manifest with neurological deterioration. Most of these neurological MPS currently lack effective treatments. Here, we report that, compared to controls, neuraminidase 1 (NEU1) activity is drastically reduced in brain tissues of neurological MPS patients and in mouse models of MPS I, II, IIIA, IIIB and IIIC, but not of other neurological lysosomal disorders not presenting with heparan sulphate storage. We further show that accumulated heparan sulphate disrupts the lysosomal multienzyme complex of NEU1 with cathepsin A (CTSA), β-galactosidase (GLB1) and glucosamine-6-sulfate sulfatase (GALNS) necessary to maintain enzyme activity, and that NEU1 deficiency is linked to partial deficiencies of GLB1 and GALNS in cortical tissues and iPSC-derived cortical neurons of neurological MPS patients. Increased sialylation of N-linked glycans in brain samples of human MPS III patients and MPS IIIC mice implicated insufficient processing of brain N-linked sialylated glycans, except for polysialic acid, which was reduced in the brains of MPS IIIC mice. Correction of NEU1 activity in MPS IIIC mice by lentiviral gene transfer ameliorated previously identified hallmarks of the disease, including memory impairment, behavioural traits, and reduced levels of the excitatory synapse markers VGLUT1 and PSD95. Overexpression of NEU1 also restored levels of VGLUT1-/PSD95-positive puncta in cortical neurons derived from iPSC of an MPS IIIA patient. Together, our data demonstrate that heparan sulphate-induced secondary NEU1 deficiency and aberrant sialylation of glycoproteins implicated in synaptogenesis, memory, and behaviour constitute a novel pathological pathway in neurological MPS spectrum crucially contributing to CNS pathology.

Graphical abstract

Differential Solvent DEEP-STD NMR and MD Simulations Enable the Determinants of the Molecular Recognition of Heparin Oligosaccharides by Antithrombin to Be Disentangled

The interaction of heparin with antithrombin (AT) involves a specific sequence corresponding to the pentasaccharide GlcNAc/NS6S-GlcA-GlcNS3S6S-IdoA2S-GlcNS6S (AGA*IA). Recent studies have revealed that two AGA*IA-containing hexasaccharides, which differ in the sulfation degree of the iduronic acid unit, exhibit similar binding to AT, albeit with different affinities. However, the lack of experimental data concerning the molecular contacts between these ligands and the amino acids within the protein-binding site prevents a detailed description of the complexes. Differential epitope mapping (DEEP)-STD NMR, in combination with MD simulations, enables the experimental observation and comparison of two heparin pentasaccharides interacting with AT, revealing slightly different bound orientations and distinct affinities of both glycans for AT. We demonstrate the effectiveness of the differential solvent DEEP-STD NMR approach in determining the presence of polar residues in the recognition sites of glycosaminoglycan-binding proteins.

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.

Computational modeling of protein-carbohydrate interactions: Current trends and future challenges

The article leads the reader through an up-to-date presentation of the concepts, developments, and main applications of computational modeling to study protein-carbohydrate interactions. It follows with the presentation of some current issues and perspectives arising from the expected evolution of generic methodological developments in deep learning, immersive analytics, and virtual reality for molecular visualization and data management. Such methodological developments for macromolecular interactions would greatly benefit a wide range of scientific endeavors in the field of carbohydrate chemistry and biochemistry, including the following interrelated efforts dealing with highly crowded media, with examples concerning glycoside transferases, the extracellular matrix, and the exploration of interactions between complex carbohydrates and intrinsically disordered 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.

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.

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.

Atomic-Resolution Experimental Structural Biology and Molecular Dynamics Simulations of Hyaluronan and Its Complexes

This review summarizes the atomic-resolution structural biology of hyaluronan and its complexes available in the Protein Data Bank, as well as published studies of atomic-resolution explicit-solvent molecular dynamics simulations on these and other hyaluronan and hyaluronan-containing systems. Advances in accurate molecular mechanics force fields, simulation methods and software, and computer hardware have supported a recent flourish in such simulations, such that the simulation publications now outnumber the structural biology publications by an order of magnitude. In addition to supplementing the experimental structural biology with computed dynamic and thermodynamic information, the molecular dynamics studies provide a wealth of atomic-resolution information on hyaluronan-containing systems for which there is no atomic-resolution structural biology either available or possible. Examples of these summarized in this review include hyaluronan pairing with other hyaluronan molecules and glycosaminoglycans, with ions, with proteins and peptides, with lipids, and with drugs and drug-like molecules. Despite limitations imposed by present-day computing resources on system size and simulation timescale, atomic-resolution explicit-solvent molecular dynamics simulations have been able to contribute significant insight into hyaluronan's flexibility and capacity for intra- and intermolecular non-covalent interactions.

Impact of calcium ions on the structural and dynamic properties of heparin oligosaccharides by computational analysis

Heparin (HP) belongs to glycosaminoglycans (GAGs), anionic linear polysaccharides composed of repetitive disaccharide units. They are key players in many biological processes occurring in the extracellular matrix and at the cell surface. GAGs are challenging molecules for computational research due to their high chemical heterogeneity, flexibility, periodicity, pseudosymmetry, predominantly electrostatics-driven nature of interactions with their protein partners and potential multipose binding. The molecular mechanisms underlying GAG interactions mediated by divalent ions, which are important for GAG binding to several proteins, are not well understood. The goal of this study was to characterize the binding of Ca²⁺ to two HP oligosaccharides of different lengths (dp10 and dp18, dp: degree of polymerization) and their impact on HP conformational space and their dynamic behavior with the use of molecular dynamics (MD)-based approaches with two Ca²⁺ parameter sets. MD data suggested that the flexibility of the monosaccharides, the glycosidic linkages and ring puckering were not affected by the presence of Ca²⁺, in contrast to H-bond propensities and the calculated R_g for a fraction of the oligosaccharide populations in both dp10 and dp18. Moreover, the essential differences in the data obtained by using two Ca²⁺ parameter sets were reported.

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.

In-Depth Molecular Dynamics Study of All Possible Chondroitin Sulfate Disaccharides Reveals Key Insight into Structural Heterogeneity and Dynamism

GAGs exhibit a high level of conformational and configurational diversity, which remains untapped in terms of the recognition and modulation of proteins. Although GAGs are suggested to bind to more than 800 biologically important proteins, very few therapeutics have been designed or discovered so far. A key challenge is the inability to identify, understand and predict distinct topologies accessed by GAGs, which may help design novel protein-binding GAG sequences. Recent studies on chondroitin sulfate (CS), a key member of the GAG family, pinpointing its role in multiple biological functions led us to study the conformational dynamism of CS building blocks using molecular dynamics (MD). In the present study, we used the all-atom GLYCAM06 force field for the first time to explore the conformational space of all possible CS building blocks. Each of the 16 disaccharides was solvated in a TIP3P water box with an appropriate number of counter ions followed by equilibration and a production run. We analyzed the MD trajectories for torsional space, inter- and intra-molecular H-bonding, bridging water, conformational spread and energy landscapes. An in-house phi and psi probability density analysis showed that 1→3-linked sequences were more flexible than 1→4-linked sequences. More specifically, phi and psi regions for 1→4-linked sequences were held within a narrower range because of intra-molecular H-bonding between the GalNAc O5 atom and GlcA O3 atom, irrespective of sulfation pattern. In contrast, no such intra-molecular interaction arose for 1→3-linked sequences. Further, the stability of 1→4-linked sequences also arose from inter-molecular interactions involving bridged water molecules. The energy landscape for both classes of CS disaccharides demonstrated increased ruggedness as the level of sulfation increased. The results show that CS building blocks present distinct conformational dynamism that offers the high possibility of unique electrostatic surfaces for protein recognition. The fundamental results presented here will support the development of algorithms that help to design longer CS chains for protein recognition.

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 Influences of Sulphation, Salt Type, and Salt Concentration on the Structural Heterogeneity of Glycosaminoglycans

The increasing recognition of the biochemical importance of glycosaminoglycans (GAGs) has in recent times made them the center of attention of recent research investigations. It became evident that subtle conformational factors play an important role in determining the relationship between the chemical composition of GAGs and their activity. Therefore, a thorough understanding of their structural flexibility is needed, which is addressed in this work by means of all-atom molecular dynamics (MD) simulations. Four major GAGs with different substitution patterns, namely hyaluronic acid as unsulphated GAG, heparan-6-sulphate, chondroitin-4-sulphate, and chondroitin-6-sulphate, were investigated to elucidate the influence of sulphation on the dynamical features of GAGs. Moreover, the effects of increasing NaCl and KCl concentrations were studied as well. Different structural parameters were determined from the MD simulations, in combination with a presentation of the free energy landscape of the GAG conformations, which allowed us to unravel the conformational fingerprints unique to each GAG. The largest effects on the GAG structures were found for sulphation at position 6, as well as binding of the metal ions in the absence of chloride ions to the carboxylate and sulphate groups, which both increase the GAG conformational flexibility.

Conformational Properties of Glycosaminoglycan Disaccharides: A Molecular Dynamics Study

The structure and conformation of glycosaminoglycans (GAGs) are of central importance to understand the mechanisms behind their functions in biological systems. Due to the inherent chemical and structural heterogeneity of GAGs, focusing on longer, naturally existing GAG chains hinders drawing conclusions on the influence of the chemical functionalization on the basic conformational degree of freedom, that is, the dynamic shape of glycosidic linkage present in the particular disaccharide repeating unit. In the present study, we have considered the complete set of 106 GAG-related disaccharides, being potential building blocks for longer GAG chains (including hyaluronan, chondroitin, keratan, dermatan, and heparan). Both the unfunctionalized units and all possible combinations of either partially or fully sulfated derivatives contribute to this number. The unbiased and enhanced sampling molecular dynamics simulations provide a link to understand the influence of sulfation on the conformational properties of GAG glycosidic linkages. Residue-residue hydrogen bonding is not significant for either the glycosidic linkage conformation or its flexibility. It was found that in the majority of cases, the dominating conformation of the linkage is weakly affected by sulfation and the main role is played by the steric and stereoelectronic effects. However, there exist numerous cases where sulfation increases the contribution of alternative conformations to a nonnegligible extent and, in some rare cases (restricted to disaccharides building heparan), leads to the reorientation of the glycosidic linkage. The identified sulfation sites, being the most important in this context, are C₆ and C₃ at the GlcNAc residue. Finally, the full set of free energy maps relying on the glycosidic dihedral angle values for diverse GAG disaccharides are provided; they may be used for further studies, focused on longer GAG chains.

Chondroitin Sulfate/Dermatan Sulfate-Protein Interactions and Their Biological Functions in Human Diseases: Implications and Analytical Tools

Chondroitin sulfate (CS) and dermatan sulfate (DS) are linear anionic polysaccharides that are widely present on the cell surface and in the cell matrix and connective tissue. CS and DS chains are usually attached to core proteins and are present in the form of proteoglycans (PGs). They not only are important structural substances but also bind to a variety of cytokines, growth factors, cell surface receptors, adhesion molecules, enzymes and fibrillary glycoproteins to execute series of important biological functions. CS and DS exhibit variable sulfation patterns and different sequence arrangements, and their molecular weights also vary within a large range, increasing the structural complexity and diversity of CS/DS. The structure-function relationship of CS/DS PGs directly and indirectly involves them in a variety of physiological and pathological processes. Accumulating evidence suggests that CS/DS serves as an important cofactor for many cell behaviors. Understanding the molecular basis of these interactions helps to elucidate the occurrence and development of various diseases and the development of new therapeutic approaches. The present article reviews the physiological and pathological processes in which CS and DS participate through their interactions with different proteins. Moreover, classic and emerging glycosaminoglycan (GAG)-protein interaction analysis tools and their applications in CS/DS-protein characterization are also discussed.

Sulfation and Calcium Favor Compact Conformations of Chondroitin in Aqueous Solutions

The effects of sulfation and calcium cations (Ca²⁺) on the atomic-resolution conformational properties of chondroitin carbohydrate polymers in aqueous solutions are not well studied owing to experimental challenges. Here, we compare all-atom explicit-solvent molecular dynamics simulations results for pairs of O-type (nonsulfated) and A-type (GlcNAc 4-O-sulfated) chondroitin 20-mers in 140 mM NaCl with and without Ca²⁺ and find that both sulfation and Ca²⁺ favor more compact polymer conformations. We also show that subtle differences in force-field parametrization can have dramatic effects on Ca²⁺ binding to chondroitin carboxylate and sulfate functional groups and thereby determine Ca²⁺-mediated intra- and interstrand association. In addition to providing an atomic-resolution picture of the interaction of Ca²⁺ with sulfated and nonsulfated chondroitin polymers, the molecular dynamics data emphasize the importance of careful force-field parametrization to balance ion-water and ion-chondroitin interactions and suggest additional parametrization efforts to tune interactions involving sulfate.

A Bittersweet Computational Journey among Glycosaminoglycans

Glycosaminoglycans (GAGs) are linear polysaccharides. In proteoglycans (PGs), they are attached to a core protein. GAGs and PGs can be found as free molecules, associated with the extracellular matrix or expressed on the cell membrane. They play a role in the regulation of a wide array of physiological and pathological processes by binding to different proteins, thus modulating their structure and function, and their concentration and availability in the microenvironment. Unfortunately, the enormous structural diversity of GAGs/PGs has hampered the development of dedicated analytical technologies and experimental models. Similarly, computational approaches (in particular, molecular modeling, docking and dynamics simulations) have not been fully exploited in glycobiology, despite their potential to demystify the complexity of GAGs/PGs at a structural and functional level. Here, we review the state-of-the art of computational approaches to studying GAGs/PGs with the aim of pointing out the “bitter” and “sweet” aspects of this field of research. Furthermore, we attempt to bridge the gap between bioinformatics and glycobiology, which have so far been kept apart by conceptual and technical differences. For this purpose, we provide computational scientists and glycobiologists with the fundamentals of these two fields of research, with the aim of creating opportunities for their combined exploitation, and thereby contributing to a substantial improvement in scientific knowledge.

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

The interaction of regulatory proteins with extracellular matrix or cell surface-anchored glycosaminoglycans (GAGs) plays important roles in molecular recognition, wound healing, growth, inflammation and many other processes. In spite of their high biological relevance, protein-GAG complexes are significantly underrepresented in structural databases because standard tools for structure determination experience difficulties in studying these complexes. Co-crystallization with subsequent X-ray analysis is hampered by the high flexibility of GAGs. NMR spectroscopy experiences difficulties related to the periodic nature of the GAGs and the sparse proton network between protein and GAG with distances that typically exceed the detection limit of nuclear Overhauser enhancement spectroscopy. In contrast, computer modeling tools have advanced over the last years delivering specific protein-GAG docking approaches successfully complemented with molecular dynamics (MD)-based analysis. Especially the combination of NMR spectroscopy in solution providing sparse structural constraints with molecular docking and MD simulations represents a useful synergy of forces to describe the structure of protein-GAG complexes. Here we review recent methodological progress in this field and bring up examples where the combination of new NMR methods along with cutting-edge modeling has yielded detailed structural information on complexes of highly relevant cytokines with GAGs.

Insights into the roles of charged residues in substrate binding and mode of action of mannuronan C-5 epimerase AlgE4

Mannuronan C-5 epimerases catalyse the epimerization of monomer residues in the polysaccharide alginate, changing the physical properties of the biopolymer. The enzymes are utilized to tailor alginate to numerous biological functions by alginate-producing organisms. The underlying molecular mechanisms that control the processive movement of the epimerase along the substrate chain is still elusive. To study this, we have used an interdisciplinary approach combining molecular dynamics simulations with experimental methods from mutant studies of AlgE4, where initial epimerase activity and product formation were addressed with NMR spectroscopy, and characteristics of enzyme-substrate interactions were obtained with isothermal titration calorimetry and optical tweezers. Positive charges lining the substrate-binding groove of AlgE4 appear to control the initial binding of poly-mannuronate, and binding also seems to be mediated by both electrostatic and hydrophobic interactions. After the catalytic reaction, negatively charged enzyme residues might facilitate dissociation of alginate from the positive residues, working like electrostatic switches, allowing the substrate to translocate in the binding groove. Molecular simulations show translocation increments of two monosaccharide units before the next productive binding event resulting in MG-block formation, with the epimerase moving with its N-terminus towards the reducing end of the alginate chain. Our results indicate that the charge pair R343-D345 might be directly involved in conformational changes of a loop that can be important for binding and dissociation. The computational and experimental approaches used in this study complement each other, allowing for a better understanding of individual residues' roles in binding and movement along the alginate chains.

Evaluation of replica exchange with repulsive scaling approach for docking glycosaminoglycans

Glycosaminoglycans (GAGs), long linear periodic anionic polysaccharides, are key molecules in the extracellular matrix (ECM). Therefore, deciphering their role in the biologically relevant context is important for fundamental understanding of the processes ongoing in ECM and for establishing new strategies in the regenerative medicine. Although GAGs represent a number of computational challenges, molecular docking is a powerful tool for analysis of their interactions. Despite the recent development of GAG-specific docking approaches, there is plenty of room for improvement. Here, replica exchange molecular dynamics with repulsive scaling (REMD-RS) recently proved to be a successful approach for protein-protein complexes, was applied to dock GAGs. In this method, effective pairwise radii are increased in different Hamiltonian replicas. REMD-RS is shown to be an attractive alternative to classical docking approaches for GAGs. This work contributes to setting up of GAG-specific computational protocols and provides new insights into the nature of these biological systems.

NMR Characterization of the Interactions Between Glycosaminoglycans and Proteins

Glycosaminoglycans (GAGs) constitute a considerable fraction of the glycoconjugates found on cellular membranes and in the extracellular matrix of virtually all mammalian tissues. The essential role of GAG-protein interactions in the regulation of physiological processes has been recognized for decades. However, the underlying molecular basis of these interactions has only emerged since 1990s. The binding specificity of GAGs is encoded in their primary structures, but ultimately depends on how their functional groups are presented to a protein in the three-dimensional space. This review focuses on the application of NMR spectroscopy on the characterization of the GAG-protein interactions. Examples of interpretation of the complex mechanism and characterization of structural motifs involved in the GAG-protein interactions are given. Selected families of GAG-binding proteins investigated using NMR are also described.

Computational insights into the role of calcium ions in protein–glycosaminoglycan systems

The prediction power of computational methodologies for studying the role of ions in protein–glycosaminoglycan interactions was critically assessed.

Role of Oligosaccharide Chain Polarity in Protein-Glycosaminoglycan Interactions

Glycosaminoglycans (GAGs) are long unbranched anionic polysaccharides made up of repetitive disaccharide units involved in biologically relevant processes in the extracellular matrix such as cell proliferation and communication. A GAG can be bound in antiparallel energetically comparable orientations on the protein surface, and these orientations are, therefore, difficult to distinguish both experimentally and computationally. In this study, for the first time we analyzed the impact of the GAG chain polarity on the interactions with Fibroblast Growth Factors-1 and -2 (FGF-1 and FGF-2). We performed a series of 1 μs molecular dynamics simulations of the FGF-1 and FGF-2 complexes with heparin (HP), a GAG representative, of different length. We analyzed the relationship between the HP orientation, energetic, and conformational space characteristics of FGF-1-HP and FGF-2-HP complexes. We concluded that HP can be bound by these proteins in the same binding sites but in different orientations, while the orientation present in the experimental structure might be favorable. Our data presented in this study provide a novel view on the impact of GAG polarity on the specificity of protein-GAG complex formation, which is an essential aspect for the proper understanding of the intermolecular interactions in these systems.

GAG-DB, the New Interface of the Three-Dimensional Landscape of Glycosaminoglycans

Glycosaminoglycans (GAGs) are complex linear polysaccharides. GAG-DB is a curated database that classifies the three-dimensional features of the six mammalian GAGs (chondroitin sulfate, dermatan sulfate, heparin, heparan sulfate, hyaluronan, and keratan sulfate) and their oligosaccharides complexed with proteins. The entries are structures of GAG and GAG-protein complexes determined by X-ray single-crystal diffraction methods, X-ray fiber diffractometry, solution NMR spectroscopy, and scattering data often associated with molecular modeling. We designed the database architecture and the navigation tools to query the database with the Protein Data Bank (PDB), UniProtKB, and GlyTouCan (universal glycan repository) identifiers. Special attention was devoted to the description of the bound glycan ligands using simple graphical representation and numerical format for cross-referencing to other databases in glycoscience and functional data. GAG-DB provides detailed information on GAGs, their bound protein ligands, and features their interactions using several open access applications. Binding covers interactions between monosaccharides and protein monosaccharide units and the evaluation of quaternary structure. GAG-DB is freely available.

Glycosaminoglycan-Protein Interactions: The First Draft of the Glycosaminoglycan Interactome

The six mammalian glycosaminoglycans (GAGs), chondroitin sulfate, dermatan sulfate, heparin, heparan sulfate, hyaluronan, and keratan sulfate, are linear polysaccharides. Except for hyaluronan, they are sulfated to various extent, and covalently attached to proteins to form proteoglycans. GAGs interact with growth factors, morphogens, chemokines, extracellular matrix proteins and their bioactive fragments, receptors, lipoproteins, and pathogens. These interactions mediate their functions, from embryonic development to extracellular matrix assembly and regulation of cell signaling in various physiological and pathological contexts such as angiogenesis, cancer, neurodegenerative diseases, and infections. We give an overview of GAG-protein interactions (i.e., specificity and chemical features of GAG- and protein-binding sequences), and review the available GAG-protein interaction networks. We also provide the first comprehensive draft of the GAG interactome composed of 832 biomolecules (827 proteins and five GAGs) and 932 protein-GAG interactions. This network is a scaffold, which in the future should integrate structures of GAG-protein complexes, quantitative data of the abundance of GAGs in tissues to build tissue-specific interactomes, and GAG interactions with metal ions such as calcium, which plays a major role in the assembly of the extracellular matrix and its interactions with cells. This contextualized interactome will be useful to identify druggable GAG-protein interactions for therapeutic purpose.

Role of Glycosaminoglycans in Procathepsin B Maturation: Molecular Mechanism Elucidated by a Computational Study

Procathepsins are an inactive, immature form of cathepsins, predominantly cysteine proteases present in the extracellular matrix (ECM) and in lysosomes that play a key role in various biological processes such as bone resorption or intracellular proteolysis. The enzymatic activity of cathepsins can be mediated by glycosaminoglycans (GAGs), long unbranched periodic negatively charged polysaccharides found in ECM that take part in many biological processes such as anticoagulation, angiogenesis, and tissue regeneration. In addition to the known effects on mature cathepsins, GAGs can mediate the maturation process of procathepsins, in particular, procathepsin B. However, the detailed mechanism of this mediation at the molecular level is still unknown. In this study, for the first time, we aimed to unravel the role of GAGs in this process using computational approaches. We rigorously analyzed procathepsin B–GAG complexes in terms of their dynamics, energetics, and potential allosteric regulation. We revealed that GAGs can stabilize the conformation of the procathepsin B structure with the active site accessible for the substrate and concluded that GAGs most probably bind to procathepsin B once the zymogen adopts the enzymatically active conformation. Our data provided a novel mechanistic view of the maturation process of procathepsin B, while the approaches elaborated here might be useful to study other procathepsins. Furthermore, our data can serve as a rational guide for experimental work on procathepsin–GAG systems that are not characterized in vivo and in vitro yet.

GAG Builder: a web-tool for modeling 3D structures of glycosaminoglycans

Here we report the launch of a web-tool (the GLYCAM-Web GAG Builder, www.glycam.org/gag) for the rapid and straightforward prediction of 3D structural models for glycosaminoglycans (GAGs). The tool provides the user with coordinate files (PDB format) for use in visualization, as well as files for performing MD simulation with the AMBER software package. Counter ions and water may also be added as desired. The tool is designed with the non-expert in mind, and as such has implemented typical default values for structural parameters, which the user may change if desired. Multiple GAG types are supported, including Heparin/Heparan Sulfate, Chondroitin Sulfate, Dermatan Sulfate, Keratan Sulfate, and Hyaluronan; however, the user may alter the default sulfation patterns to create novel sequences. The common non-natural unsaturated uronic acid (ΔUA) and its sulfated derivative are also supported.

Effect of solvent and ions on the structure and dynamics of a hyaluronan molecule

Hyaluronic acid (hyaluronan, HA) is a negatively charged polysaccharide forming highly swollen random coils in aqueous solutions. Their size decreases along with growing salt concentration, but the mechanism of this phenomenon remains unclear. We carry out molecular-dynamics simulations of a 48-monosaccharide HA oligomer in varying salt concentration and temperature. They identify the interaction points of Na⁺ ions with the HA chain and reveal their influence on the HA solvation-shell structure. The salt-dependent variation of the molecular size does not consist in the distribution of the dihedral angles of the glycosidic connections but is driven by the random flips of individual dihedral angles, which cause the formation of temporary hairpin-like structures effectively shortening the chain. They are induced by the frequency of cation-chain interactions that grow with the salt concentration, but is reduced by the simultaneous decrease of ions' activities. This leads to an anomalous random-coil shrinkage at 0.6 M salt concentration.

Analysis of Procollagen C-Proteinase Enhancer-1/Glycosaminoglycan Binding Sites and of the Potential Role of Calcium Ions in the Interaction

In this study, we characterize the interactions between the extracellular matrix protein, procollagen C-proteinase enhancer-1 (PCPE-1), and glycosaminoglycans (GAGs), which are linear anionic periodic polysaccharides. We applied molecular modeling approaches to build a structural model of full-length PCPE-1, which is not experimentally available, to predict GAG binding poses for various GAG lengths, types and sulfation patterns, and to determine the effect of calcium ions on the binding. The computational data are analyzed and discussed in the context of the experimental results previously obtained using surface plasmon resonance binding assays. We also provide experimental data on PCPE-1/GAG interactions obtained using inhibition assays with GAG oligosaccharides ranging from disaccharides to octadecasaccharides. Our results predict the localization of GAG-binding sites at the amino acid residue level onto PCPE-1 and is the first attempt to describe the effects of ions on protein-GAG binding using modeling approaches. In addition, this study allows us to get deeper insights into the in silico methodology challenges and limitations when applied to GAG-protein interactions.

The Expression, Regulation, and Biomarker Potential of Glypican-1 in Cancer

Glypican-1 (GPC-1) and other glypicans are a family of heparan sulfate proteoglycans. These proteins are highly expressed on the cell membrane and in the extracellular matrix, functioning mainly as modulators of growth factor signaling. Some of them are abnormally expressed in cancer, possibly involved in tumorigenesis, and detectable in blood as potential clinical biomarkers. GPC-1 is another glypican member that has been found to be associated with some cancers, and has increasingly interested the cancer field. Here we provide a brief review about GPC-1 in its expression, signaling and potential as a cancer biomarker.

Docking software performance in protein-glycosaminoglycan systems

We present a benchmarking study for protein-glycosaminoglycan systems with eight docking programs: Dock, rDock, ClusPro, PLANTS, HADDOCK, Hex, SwissDock and ATTRACT. We used a non-redundant representative dataset of 28 protein-glycosaminoglycan complexes with experimentally available structures, where a glycosaminoglycan ligand was longer than a trimer. Overall, the ligand binding poses could be correctly predicted in many cases by the tested docking programs, however the ranks of the docking poses are often poorly assigned. Our results suggest that Dock program performs best in terms of the pose placement, has the most suitable scoring function, and its performance did not depend on the ligand size. This suggests that the implementation of the electrostatics as well as the shape complementarity procedure in Dock are the most suitable for docking glycosaminoglycan ligands. We also analyzed how free energy patterns of the benchmarking complexes affect the performance of the evaluated docking software.

A pipeline to translate glycosaminoglycan sequences into 3D models. Application to the exploration of glycosaminoglycan conformational space

Mammalian glycosaminoglycans are linear complex polysaccharides comprising heparan sulfate, heparin, dermatan sulfate, chondroitin sulfate, keratan sulfate and hyaluronic acid. They bind to numerous proteins and these interactions mediate their biological activities. GAG-protein interaction data reported in the literature are curated mostly in MatrixDB database (http://matrixdb.univ-lyon1.fr/). However, a standard nomenclature and a machine-readable format of GAGs together with bioinformatics tools for mining these interaction data are lacking. We report here the building of an automated pipeline to (i) standardize the format of GAG sequences interacting with proteins manually curated from the literature, (ii) translate them into the machine-readable GlycoCT format and into SNFG (Symbol Nomenclature For Glycan) images and (iii) convert their sequences into a format processed by a builder generating three-dimensional structures of polysaccharides based on a repertoire of conformations experimentally validated by data extracted from crystallized GAG-protein complexes. We have developed for this purpose a converter (the CT23D converter) to automatically translate the GlycoCT code of a GAG sequence into the input file required to construct a three-dimensional model.

Design, synthesis, and biomedical applications of synthetic sulphated polysaccharides

This review summarizes the synthetic methods to sulphated polysaccharides, describes their compositional and structural diversity in regards to activity, and showcases their biomedical applications.

Exploring Structure⁻Property Relationships of GAGs to Tailor ECM-Mimicking Hydrogels

Glycosaminoglycans (GAGs) are a class of linear polysaccharides that are ubiquitous in the extracellular matrix (ECM) and on cell surfaces. Due to their key role in development, homeostasis, pathogenesis, and regeneration, GAGs are increasingly used in the design of ECM-mimicking hydrogels to stimulate tissue formation and regenerative processes via specifically orchestrated cell-instructive signals. These applications first and foremost build on the ability of GAGs to effectively bind, protect, and release morphogens. The specificity and strength of morphogen-GAG interactions are largely governed by the number and spatial distribution of negatively charged sulfate groups carried by GAGs. Herein, we summarize a mean-field approach to quantify the density of ionizable groups, GAG concentration, and cross-linking degree of GAG-containing hydrogels on the basis of microslit electrokinetic experiments. We further present and discuss a continuum model of mucosa that accounts for charge regulation by glycan-ion pairing in biological contexts and under conditions of macromolecular crowding. Finally, we discuss the modulation of the morphogen binding and transport in GAG hydrogels by selective desulfation of the GAG component.

Solution Conformation of Heparin Tetrasaccharide. DFT Analysis of Structure and Spin⁻Spin Coupling Constants

Density functional theory (DFT) has provided detailed information on the molecular structure and spin⁻spin coupling constants of heparin tetrasaccharide (GlcNS,6S-IdoA2S-GlcNS,6S-IdoA2S-OMe) representing the predominant heparin repeating-sequence. The fully optimised molecular structures of two tetrasaccharide conformations (differing from each other in the conformational form of the sulphated iduronic acid residue⁻one ¹C₄ and the other ²S₀) were obtained using the B3LYP/6-311+G(d,p) level of theory and applying explicit water molecules to simulate the presence of a solvent. The theoretical data provided insight into variations of the bond lengths, bond angles and torsion angles, formations of intra- and intermolecular hydrogen bonds and ionic interactions. Optimised molecular structures indicated the formation of a complex hydrogen bond network, including interresidue and intraresidue bonds. The ionic interactions strongly influence the first hydration shell and, together with hydrogen bonds, play an important role in shaping the 3D tetrasaccharide structure. DFT-derived indirect three⁻bond proton⁻proton coupling constants (³J_H-C-C-H) showed that the best agreement with experiment was obtained with a weighted average of 67:33 (¹C₄:²S₀) of the IdoA2S forms. Detailed analysis of Fermi-contact contributions to ³J_H-C-C-H showed that important contributions arise from the oxygen lone pairs of neighbouring oxygen atoms. The analysis also showed that the magnitude of diamagnetic spin⁻orbit contributions are sufficiently large to determine the magnitude of some proton⁻proton coupling constants. The data highlight the need to use appropriate quantum-chemical calculations for a detailed understanding of the solution properties of heparin oligosaccharides.