Computational drill down on FGF1-heparin interactions through methodological evaluation

Conformational preferences of heparan sulfate to recognize the CXCL8 dimer in aqueous medium: degree of sulfation and hydrogen bonds

The sulfation pattern and epimerization of the long-chain sulfated polysaccharide heparan sulfate (HS) cause structural diversity and regulate various physiological and pathological processes when binding with proteins. In this work, we performed a series of molecular dynamics simulations of three variants of the octadecasaccharide HS with varying sulfation positions in aqueous medium in their free forms and in the presence of the chemokine CXCL8 dimer. The free energy of binding depicts the sulfation at the 6-O position of GlcNAc (HS6S), and both 3-O and 6-O positions of GlcNAc (HS3S6S) of HS variants are more likely to bind with the CXCL8 dimer than the triply sulfated HS2S3S6S, which is sulfated at the 2-O position of GlcUA additionally along with 3-O and 6-O positions of GlcNAc. Binding between HS and CXCL8 was driven by electrostatic and van der Waals interactions predominantly regardless of the sulfation pattern; however, unfavorable entropic contribution suppressed the interaction between HS and CXCL8. The contribution of different amino acid residues to the binding energetics suggested that basic amino acids line up the binding site of CXCL8. This study further acknowledges the role of interfacial water that is structured and bound with HS through hydrogen bonds, exhibiting differential hydrogen bond relaxation dynamics compared to when the HS molecules are free. Moreover, this study identifies that with the increase in sulfation, the HS-water hydrogen bond relaxation occurs faster with the complexation, while the reverse trend is followed in their free forms. Significant structural adaptation of the different sulfated HS molecules, as verified from the free energy landscapes generated from various reaction coordinates, root-mean-square-deviations, end-to-end distances, including ring pucker angles, dihedral flexibility, and the high conformational entropy cost arising from the glycosidic bonds, suggests that the different sulfated variants of HS undergo significant structural transformation to bind with CXCL8. The presence of a CXCL8 dimer imposes the bound forms of HS to adopt non-linear structures with skew-boat conformations. The atomistic details of the study would help in understanding the selectivity and conformational diversity, as well as the role of solvents in the recognition of CXCL8 by different sulfated variants of HS molecules.

Sulfation Effects of Chondroitin Sulfate to Bind a Chemokine in Aqueous Medium: Conformational Heterogeneity and Dynamics from Molecular Simulation

The sulfation patterns and degree of sulfation of chondroitin sulfate (CS), an important class of glycosaminoglycans (GAG), and their interactions with chemokines are accountable for various diseases. To realize the underlying mechanism of such complex biological phenomena at a molecular level and their application in rational drug design, a study on conformations and dynamics of CSs is necessary. To explore this, in this study, we performed a series of atomistic molecular dynamics (MD) simulations with different sulfated variants of octadecasaccharide CS, like CS-C, CS-E, and CS-T, in their free forms and when bound to the protein chemokine CXCL8 dimer in an aqueous medium. The calculated binding free energy of CSs with the CXCL8 dimer is favorable, and the degree of sulfation favors the complexation process further with prominent hydrophobic and hydrogen-bonded interactions. We find that the recognition is associated with the configurational entropy loss of the CS molecules as calculated from the Gaussian mixture approach, which supports that the degree of sulfation regulates the process. Cluster analysis through the k-means algorithm and end-to-end distance measurement revealed that although the free CS molecules adopted linear conformations, the nonlinear conformations during binding with protein were noted. Adaptation of nonlinear forms in the bound forms is noteworthy for the less-sulfated CS-C and CS-E. Apart from favorable 4C1 conformations, the occasional appearance of skew-boat forms from the free-energy map of ring pucker for the GlcUA unit was observed, which remains unaffected by the sulfation. We find that during recognition, the average relaxation time of intra-CS and inter-CS-CXCL8 hydrogen bonds (HBs) is about a magnitude lesser than that of CS-water HBs, most prominent on the involvement of higher sulfated CS-T analogues. The translational motion of surrounded water molecules in CSs exhibited sublinear diffusion, and the degree of sublinearity increases around the heavily sulfated molecules due to the hindrance created by them as well as the presence of the chemokine and exhibited markedly slow heterogeneous diffusion.

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.

Prioritizing Small Sets of Molecules for Synthesis through in-silico Tools: A Comparison of Common Ranking Methods

Prioritizing molecules for synthesis is a key role of computational methods within medicinal chemistry. Multiple tools exist for ranking molecules, from the cheap and popular molecular docking methods to more computationally expensive molecular-dynamics (MD)-based methods. It is often questioned whether the accuracy of the more rigorous methods justifies the higher computational cost and associated calculation time. Here, we compared the performance on ranking the binding of small molecules for seven scoring functions from five docking programs, one end-point method (MM/GBSA), and two MD-based free energy methods (PMX, FEP+). We investigated 16 pharmaceutically relevant targets with a total of 423 known binders. The performance of docking methods for ligand ranking was strongly system dependent. We observed that MD-based methods predominantly outperformed docking algorithms and MM/GBSA calculations. Based on our results, we recommend the application of MD-based free energy methods for prioritization of molecules for synthesis in lead optimization, whenever feasible.

Binding affinity estimation from restrained umbrella sampling simulations

The protein-ligand binding affinity quantifies the binding strength between a protein and its ligand. Computer modeling and simulations can be used to estimate the binding affinity or binding free energy using data- or physics-driven methods or a combination thereof. Here we discuss a purely physics-based sampling approach based on biased molecular dynamics simulations. Our proposed method generalizes and simplifies previously suggested stratification strategies that use umbrella sampling or other enhanced sampling simulations with additional collective-variable-based restraints. The approach presented here uses a flexible scheme that can be easily tailored for any system of interest. We estimate the binding affinity of human fibroblast growth factor 1 to heparin hexasaccharide based on the available crystal structure of the complex as the initial model and four different variations of the proposed method to compare against the experimentally determined binding affinity obtained from isothermal titration calorimetry experiments.

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.

Development of a system of heparin multilayers on titanium surfaces for dual growth factor release

The aim of the present study was to establish a modular platform of poly-L-lysine-heparin (PLL-Hep) polyelectrolyte multilayer (PEM) coatings on titanium surfaces for dual growth factor delivery of recombinant human bone morphogenic protein 2 (rhBMP2) and recombinant human vascular endothelial growth factor 165 (rhVEGF165) in clinically relevant quantities. Release characteristics for both growth factors differed significantly depending on film architecture. rhBMP2 induced activation of alkaline phosphatase in C2C12 cells and proliferation of human mesenchymal stem cells (hMSCs). rhVEGF mediated induction of von Willebrand factor (vWF) in hMSCs and proliferation of human umbilical vein endothelial cells. Osteogenic and angiogenic effects were modified by variation in cross-linking and architecture of the PEMs. By creating multilayer films with distinct zones, release characteristics and proportion of both growth factor delivery could be tuned and surface-activity modified to enhance angiogenic or osteogenic function in various ways. In summary, the system provides a modular platform for growth factor delivery that allows for individual composition and accentuation of angiogenic and osteogenic surface properties.

Aqueous Molecular Dynamics for Understanding Glycosaminoglycan Recognition by Proteins

Glycosaminoglycans (GAGs) are biopolymers that exist in most organisms. GAGs are known to bind to hundreds of proteins and partake in multiple biological processes such as growth, morphogenesis, inflammation, infection, and others. Their intrinsic structural heterogeneity and conformational variability introduce major challenges in experimental studies. On the other hand, recent advances in force field development and computational technology have yielded phenomenal opportunity to study thousands of GAG sequences simultaneously to understand recognition of target protein(s). Here, we describe experimental setup for conventional molecular dynamics simulations of GAGs to position an experimental biologist favorably in performance, analysis and interpretation of stability, specificity, and conformational properties of GAGs, while also elucidating their interactions with amino acid residues of a protein at an atomistic level in presence of water.

Mechanistic Picture for Monomeric Human Fibroblast Growth Factor 1 Stabilization by Heparin Binding

Human fibroblast growth factor (FGF) 1 or hFGF1 is a member of the FGF family that is involved in various vital processes such as cell proliferation, cell differentiation, angiogenesis, and wound healing. hFGF1, which is associated with low stability in vivo, is known to be stabilized by binding heparin sulfate, a glycosaminoglycan that aids the protein in the activation of its cell surface receptor. The poor thermal and proteolytic stability of hFGF1 and the stabilizing role of heparin have long been observed experimentally; however, the mechanistic details of these phenomena are not well understood. Here, we have used microsecond-level equilibrium molecular dynamics (MD) simulations to quantitatively characterize the structural dynamics of monomeric hFGF1 in the presence and absence of heparin hexasaccharide. We have observed a conformational change in the heparin-binding pocket of hFGF1 that occurs only in the absence of heparin. Several intramolecular interactions were also identified within the heparin-binding pocket that form only when hFGF1 interacts with heparin. The loss of both intermolecular and intramolecular interactions in the absence of heparin plausibly leads to the observed conformational change. This conformational transition results in increased flexibility of the heparin-binding pocket and provides an explanation for the susceptibility of apo hFGF1 to proteolytic degradation and thermal instability. This study provides a glimpse into mechanistic details of the heparin-mediated stabilization of hFGF1 and encourages the use of microsecond-level MD in studying the effect of binding on protein structure and dynamics. In addition, the observed differential behavior of hFGF1 in the absence and presence of heparin provides an example, where microsecond-level all-atom MD simulations are necessary to see functionally relevant biomolecular phenomena that otherwise will not be observed on sub-microsecond time scales.

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.

Further analyses of APRIL/APRIL-Receptor/Glycosaminoglycan interactions by biochemical assays linked to computational studies

A proliferation-inducing ligand (APRIL) is a member of the tumor necrosis factor superfamily. APRIL is quite unique in this superfamily for at least for two reasons: i) it binds to glycosaminoglycans (GAGs) via its positively charged N-terminus; ii) one of its signaling receptor, the transmembrane activator CAML interactor (TACI) was also reported to bind GAGs. Here, as provided by biochemical evidences with the use of an APRIL deletion mutant linked to computational studies, APRIL-GAG interaction involved other regions than the APRIL N-terminus. Preferential interaction of APRIL with heparin followed by chondroitin sulfate E were confirmed by in silico analysis. Both computational and experimental approaches did not reveal heparan sulfate binding to TACI. Together, computational results corroborated experiments contributing with atomistic details to the knowledge on this biologically relevant trimolecular system. Additionally, a high-throughput rigorous analysis of the free energy calculations data was performed to critically evaluate the applied computational methodologies.

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.

Suramin derivatives play an important role in blocking the interaction between FGF1 and FGFRD2 to inhibit cell proliferation

The inhibition of protein function by small compounds plays a critical role in controlling cell proliferation. We report on a new class of small molecule (NCTU-Alan-2026) inhibitors for cell proliferation. NCTU-Alan-2026 blocks the interaction between FGF1 and its receptor FGF1R2D2. Extensive NMR studies combined with fluorescence experiments provided the specific mechanism of how NCTU-Alan-2026 could inhibit cell proliferation. We describe an innovative therapeutic approach for anti-proliferation and demonstrate an example of inhibition of small molecules by blocking the protein-protein interaction. We found that the compound NCTU-Alan-2026 blocked the interaction between the two proteins FGF1 and FGF1R2D2 and inhibited cell proliferation. The toxicity of NCTU-Alan-2026 is lower than that of suramin. Thus, NCTU-Alan-2026 could be a better drug than suramin in the treatment of cancer.

Heparin and heparan sulfate proteoglycans promote HIV-1 p17 matrix protein oligomerization: computational, biochemical and biological implications

p17 matrix protein released by HIV+ cells interacts with leukocytes heparan sulfate proteoglycans (HSPGs), CXCR1 and CXCR2 exerting different cytokine-like activities that contribute to AIDS pathogenesis. Since the bioactive form of several cytokines is represented by dimers/oligomers and oligomerization is promoted by binding to heparin or HSPGs, here we evaluated if heparin/HSPGs also promote p17 oligomerization. Heparin favours p17 dimer, trimer and tetramer assembly, in a time- and biphasic dose-dependent way. Heparin-induced p17 oligomerization is of electrostatic nature, being it prevented by NaCl, by removing negative sulfated groups of heparin and by neutralizing positive lysine residues in the p17 N-terminus. A new computational protocol has been implemented to study heparin chains up to 24-mer accommodating a p17 dimer. Molecular dynamics show that, in the presence of heparin, two p17 molecules undergo conformational modifications creating a continuous “electropositive channel” in which heparin sulfated groups interact with p17 basic amino acids, promoting its dimerization. At the cell surface, HSPGs induce p17 oligomerization, as demonstrated by using B-lymphoblastoid Namalwa cells overexpressing the HSPG Syndecan-1. Also, HSPGs on the surface of BJAB and Raji human B-lymphoblastoid cells are required to p17 to induce ERK_1/2 activation, suggesting that HS-induced oligomerization plays a role in p17-induced lymphoid dysregulation during AIDS.

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.

Structural basis of oligosaccharide processing by glycosaminoglycan sulfotransferases

Heparan sulfate (HS) is a sulfated polysaccharide that plays a key role in morphogenesis, physiology and pathogenesis. The biosynthesis of HS takes place in the Golgi apparatus by a group of enzymes that polymerize, epimerize and sulfate the sugar chain. This biosynthetic process introduces varying degrees of sulfate substitution, which are tightly regulated and directly dictate binding specificity to different cytokines, morphogens and growth factors. Here, we report the use of molecular dynamics simulations to investigate the dynamics of substrate recognition of two glycosaminoglycan (GAG) sulfotransferases, N-deacetylase-N-sulfotransferase and 2-O-sulfotransferase to the HS chain during the biosynthetic process. We performed multiple simulations of the binding of the sulfotransferase domains to both the HS oligosaccharide substrate and sulfate donor, 3'-phosphoadenosine-5'-phosphosulfate. Analysis of extended simulations provide detailed and useful insights into the atomic interactions that are at work during oligosaccharide processing. The fast information matching method was used to detect the enzyme global dynamics and to predict the pairwise contact of residues responsible for GAG-enzyme binding and unbinding. The correlation between HS displacement and the location of the modified GAG chain were calculated, indicating a possible route for HS and heparin during sulfotransferase processing. Our data also show sulfotransferases contain a conserved interspaced positively charged amino acid residues that form a patch which controls the protein-GAG binding equilibrium. Together, our findings provide further understanding on the fine-tuned complex mechanism of GAG biosynthesis. Our findings can also be extrapolated to other systems for calculating rates of protein-GAG binding.

Modeling large protein-glycosaminoglycan complexes using a fragment-based approach

Glycosaminoglycans (GAGs), a major constituent of the extracellular matrix, participate in cell-signaling by binding specific proteins. Structural data on protein-GAG interactions are crucial to understand and modulate these signaling processes, with potential applications in regenerative medicine. However, experimental and theoretical approaches used to study GAG-protein systems are challenged by GAGs high flexibility limiting the conformational sampling above a certain size, and by the scarcity of GAG-specific docking tools compared to protein-protein or protein-drug docking approaches. We present for the first time an automated fragment-based method for docking GAGs on a protein binding site. In this approach, trimeric GAG fragments are flexibly docked to the protein, assembled based on their spacial overlap, and refined by molecular dynamics. The method appeared more successful than the classical full-ligand approach for most of 13 tested complexes with known structure. The approach is particularly promising for docking of long GAG chains, which represents a bottleneck for classical docking approaches applied to these systems. © 2019 Wiley Periodicals, Inc.

So you think computational approaches to understanding glycosaminoglycan-protein interactions are too dry and too rigid? Think again!

Glycosaminoglycans (GAGs) play key roles in virtually all biologic responses through their interaction with proteins. A major challenge in understanding these roles is their massive structural complexity. Computational approaches are extremely useful in navigating this bottleneck and, in some cases, the only avenue to gain comprehensive insight. We discuss the state-of-the-art on computational approaches and present a flowchart to help answer most basic, and some advanced, questions on GAG-protein interactions. For example, firstly, does my protein bind to GAGs?; secondly, where does the GAG bind?; thirdly, does my protein preferentially recognize a particular GAG type?; fourthly, what is the most optimal GAG chain length?; fifthly, what is the structure of the most favored GAG sequence?; and finally, is my GAG-protein system 'specific', 'non-specific', or a combination of both? Recent advances show the field is now poised to enable a non-computational researcher perform advanced experiments through the availability of various tools and online servers.

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

Glycosaminoglycans are long unbranched and complex polysaccharides that are an essential component of mammalian extracellular matrices. Characterization of their molecular structure, dynamics and interactions are essential to understand important biological phenomena in health and disease, and will lead to novel therapeutics and medical devices. However, this has proven to be a challenge experimentally and theoretical techniques are needed to develop new hypotheses, and interpret experiments. This review aims to examine the current theoretical (rather than experimental) methods used by researchers to investigate glycosaminoglycan structure, dynamics and interactions, from the monosaccharide to the macromolecular scale. It will consider techniques such as quantum mechanics, molecular mechanics, molecular dynamics, coarse graining and docking.

Characterization of structural motifs for interactions between glycosaminoglycans and proteins

Glycosaminoglycans (GAGs) are a family of linear and anionic polysaccharides that play essential roles in many biological and physiological processes. Interactions between GAGs and proteins regulate function in many proteins and are related to many human diseases and disorders. The structural motifs and mechanisms for interactions between GAGs and proteins are not fully understood. Specific bindings, including minor but unique sequences sporadically distributed along the GAG chains or variably sulfated domains interspersed by undersulfated regions, may be specifically recognized by defined domains of a variety of proteins. Understanding the molecular basis of these interactions will provide a template for developing novel glycotherapeutic agents. The present article reviews recent methodologies and progress on the characterization of structural motifs in both GAGs and proteins involved in GAG-protein interactions. The analytical approaches are categorized into three groups: affinity-based methods; molecular docking, nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography; and mass spectrometry (MS) techniques. The advantages and limitations of each category of methods are discussed and are based on examples of using these techniques to investigate binding between GAGs and proteins.