Free energy calculations of glycosaminoglycan-protein interactions

Coarse-Grained Model of Glycosaminoglycans for Biomolecular Simulations

Proteoglycans contain glycosaminoglycans (GAGs) which are negatively charged linear polymers made of repeating disaccharide units of uronic acid and hexosamine units. They play vital roles in numerous physiological and pathological processes, particularly in governing cellular communication and attachment. Depending on their sulfonation state, acetylation, and glycosidic linkages, GAGs belong to different families. The high molecular weight, heterogeneity, and flexibility of GAGs hamper their characterization at atomic resolution, but this may be circumvented via coarse-grained (CG) approaches. In this work, we report a CG model for a library of common GAG types in their isolated or proteoglycan-linked states compatible with version 2.2 (v2.2) of the widely popular CG Martini force field. The model reproduces conformational and thermodynamic properties for a wide variety of GAGs, as well as matching structural and binding data for selected proteoglycan test systems. The parameters developed here may thus be employed to study a range of GAG-containing biomolecular systems, thereby benefiting from the efficiency and broad applicability of the Martini framework.

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.

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

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

GlycoGrip: Cell Surface-Inspired Universal Sensor for Betacoronaviruses

Inspired by the role of cell-surface glycoproteins as coreceptors for pathogens, we report the development of GlycoGrip: a glycopolymer-based lateral flow assay for detecting SARS-CoV-2 and its variants. GlycoGrip utilizes glycopolymers for primary capture and antispike antibodies labeled with gold nanoparticles for signal-generating detection. A lock-step integration between experiment and computation has enabled efficient optimization of GlycoGrip test strips which can selectively, sensitively, and rapidly detect SARS-CoV-2 and its variants in biofluids. Employing the power of the glycocalyx in a diagnostic assay has distinct advantages over conventional immunoassays as glycopolymers can bind to antigens in a multivalent capacity and are highly adaptable for mutated strains. As new variants of SARS-CoV-2 are identified, GlycoGrip will serve as a highly reconfigurable biosensor for their detection. Additionally, via extensive ensemble-based docking simulations which incorporate protein and glycan motion, we have elucidated important clues as to how heparan sulfate and other glycocalyx components may bind the spike glycoprotein during SARS-CoV-2 host-cell infection. GlycoGrip is a promising and generalizable alternative to costly, labor-intensive RT-PCR, and we envision it will be broadly useful, including for rural or low-income populations that are historically undertested and under-reported in infection statistics.

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.

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

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

Combinatorial Virtual Library Screening Study of Transforming Growth Factor-β2-Chondroitin Sulfate System

Transforming growth factor-beta (TGF-β), a member of the TGF-β cytokine superfamily, is known to bind to sulfated glycosaminoglycans (GAGs), but the nature of this interaction remains unclear. In a recent study, we found that preterm human milk TGF-β2 is sequestered by chondroitin sulfate (CS) in its proteoglycan form. To understand the molecular basis of the TGF-β2-CS interaction, we utilized the computational combinatorial virtual library screening (CVLS) approach in tandem with molecular dynamics (MD) simulations. All possible CS oligosaccharides were generated in a combinatorial manner to give 24 di- (CS02), 192 tetra- (CS04), and 1536 hexa- (CS06) saccharides. This library of 1752 CS oligosaccharides was first screened against TGF-β2 using the dual filter CVLS algorithm in which the GOLDScore and root-mean-square-difference (RMSD) between the best bound poses were used as surrogate markers for in silico affinity and in silico specificity. CVLS predicted that both the chain length and level of sulfation are critical for the high affinity and high specificity recognition of TGF-β2. Interestingly, CVLS led to identification of two distinct sites of GAG binding on TGF-β2. CVLS also deduced the preferred composition of the high specificity hexasaccharides, which were further assessed in all-atom explicit solvent MD simulations. The MD results confirmed that both sites of binding form stable GAG-protein complexes. More specifically, the highly selective CS chains were found to engage the TGF-β2 monomer with high affinity. Overall, this work present key principles of recognition with regard to the TGF-β2-CS system. In the process, it led to the generation of the in silico library of all possible CS oligosaccharides, which can be used for advanced studies on other protein-CS systems. Finally, the study led to the identification of unique CS sequences that are predicted to selectively recognize TGF-β2 and may out-compete common natural CS biopolymers.

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.

Solid state structure of sodium β-1-thiophenyl glucuronate identifies 5-coordinate sodium with three independent glucoronates

Glucuronic acid is a key component of the glycosaminoglycans (GAGs) Chrondroitin Sulfate (CS), Heparin/Heparan sulfate (HS) and Hyaluronic Acid (HA), as well an important metabolite derivative. In biological systems the carboxylate of uronic acids in GAGs is involved in important H-binding interactions, and the role of metal coordination, such as sodiated systems, has indications associated with a number of biological effects, and physiological GAG-related processes. In synthetic approaches to GAG fragments, thioglycoside intermediates, or derivatives from these, are commonly employed. Of the reported examples of sodium coordination in carbohydrates, 6-coordinate systems are usually observed often with water ligands involved, Herein we report an unexpected 5-coordinate sodiated GlcA crystal structure of the parent GlcA, but as a thioglycoside derivative, whose crystal coordination differs from previous examples, with no involvement of water as a ligand and containing a distorted trigonal bypramidal sodium with each GlcA having five of 6 oxygens sodium-coordinated.

A Systems View of the Heparan Sulfate Interactome

Heparan sulfate proteoglycans consist of a small family of proteins decorated with one or more covalently attached heparan sulfate glycosaminoglycan chains. These chains have intricate structural patterns based on the position of sulfate groups and uronic acid epimers, which dictate their ability to engage a large repertoire of heparan sulfate-binding proteins, including extracellular matrix proteins, growth factors and morphogens, cytokines and chemokines, apolipoproteins and lipases, adhesion and growth factor receptors, and components of the complement and coagulation system. This review highlights recent progress in the characterization of the so-called "heparan sulfate interactome," with a major focus on systems-wide strategies as a tool for discovery and characterization of this subproteome. In addition, we compiled all heparan sulfate-binding proteins reported in the literature to date and grouped them into a few major functional classes by applying a networking approach.

Rigorous analysis of free solution glycosaminoglycan dynamics using simple, new tools

Heparin/heparan sulfates (H/HS) are ubiquitous biopolymers that interact with many proteins to induce a range of biological functions. Unfortunately, how these biopolymers recognize their preferred protein targets remain poorly understood. It is suggested that computational simulations offer attractive avenues but a number of challenges, e.g., difficulty of selecting a comprehensive force field, few simple tools to interpret data, among others, remain. This work addresses several such challenges so as to help ease the implementation and analysis of computational experiments. First, this work presents a rigorous comparison of two different recent force fields, CHARMM36 and GLYCAM06, for H/HS studies. Second, it introduces two new straightforward parameters, i.e., end-to-end distance and minimum volume enclosing ellipsoid, to understand the myriad conformational forms of oligosaccharides that evolve over time in water. Third, it presents an application to elucidate the number and nature of inter and intramolecular, nondirect bridging water molecules, which help stabilize unique forms of H/HS. The results show that nonspecialists can use either CHARMM36 or GLYCAM06 force fields because both gave comparable results, albeit with small differences. The comparative study shows that the HS hexasaccharide samples a range of conformations with nearly equivalent energies, which could be the reason for its recognition by different proteins. Finally, analysis of the nondirect water bridges across the dynamics trajectory shows their importance in stabilization of certain conformational forms, which may become important for protein recognition. Overall, the work aids nonspecialists employ computational studies for understanding the solution behavior of H/HS.

Structural Insights Into How Proteoglycans Determine Chemokine-CXCR1/CXCR2 Interactions: Progress and Challenges

Proteoglycans (PGs), present in diverse environments, such as the cell membrane surface, extracellular milieu, and intracellular granules, are fundamental to life. Sulfated glycosaminoglycans (GAGs) are covalently attached to the core protein of proteoglycans. PGs are complex structures, and are diverse in terms of amino acid sequence, size, shape, and in the nature and number of attached GAG chains, and this diversity is further compounded by the phenomenal diversity in GAG structures. Chemokines play vital roles in human pathophysiology, from combating infection and cancer to leukocyte trafficking, immune surveillance, and neurobiology. Chemokines mediate their function by activating receptors that belong to the GPCR class, and receptor interactions are regulated by how, when, and where chemokines bind GAGs. GAGs fine-tune chemokine function by regulating monomer/dimer levels and chemotactic/haptotactic gradients, which are also coupled to how they are presented to their receptors. Despite their small size and similar structures, chemokines show a range of GAG-binding geometries, affinities, and specificities, indicating that chemokines have evolved to exploit the repertoire of chemical and structural features of GAGs. In this review, we summarize the current status of research on how GAG interactions regulate ELR-chemokine activation of CXCR1 and CXCR2 receptors, and discuss knowledge gaps that must be overcome to establish causal relationships governing the impact of GAG interactions on chemokine function in human health and disease.

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.

NMR and molecular modeling reveal specificity of the interactions between CXCL14 and glycosaminoglycans

CXCL14, chemokine (C-X-C motif) ligand 14, is a novel highly conserved chemokine with unique features. Despite exhibiting the typical chemokine fold, it has a very short N-terminus of just two amino acid residues responsible for chemokine receptor activation. CXCL14 actively participates in homeostatic immune surveillance of skin and mucosae, is linked to metabolic disorders and fibrotic lung diseases and possesses strong anti-angiogenic properties in early tumor development. In this work, we investigated the interaction of CXCL14 with various glycosaminoglycans (GAGs) by nuclear magnetic resonance spectroscopy, microscale thermophoresis, analytical heparin (HE) affinity chromatography and in silico approaches to understand the molecular basis of GAG-binding. We observed different GAG-binding modes specific for the GAG type used in the study. In particular, the CXCL14 epitope for HE suggests a binding pose distinguishable from the ones of the other GAGs investigated (hyaluronic acid, chondroitin sulfate-A/C, −D, dermatan sulfate). This observation is also supported by computational methods that included molecular docking, molecular dynamics and free energy calculations. Based on our results, we suggest that distinct GAG sulfation patterns confer specificity beyond simple electrostatic interactions usually considered to represent the driving forces in protein–GAG interactions. The CXCL14–GAG system represents a promising approach to investigate the specificity of GAG–protein interactions, which represents an important topic for developing the rational approaches to novel strategies in regenerative medicine.

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.

Syntheses of defined sulfated oligohyaluronans reveal structural effects, diversity and thermodynamics of GAG-protein binding

High binding affinities of GAG toward extracellular regulatory proteins are governed by recognition diversity, sulfation pattern, length, and anomeric functionalization.

Binding of sulfated glycosaminoglycans (GAG) to a wide spectrum of extracellular regulatory proteins is crucial for physiological processes such as cell growth, migration, tissue homeostasis and repair. Thus, GAG derivatives exhibit great relevance in the development of innovative biomaterials for tissue regeneration therapies. We present a synthetic strategy for the preparation of libraries of defined sulfated oligohyaluronans as model GAG systematically varied in length, sulfation pattern and anomeric substitution in order to elucidate the effects of these parameters on GAG recognition by regulatory proteins. Through an experimental and computational approach using fluorescence polarization, ITC, docking and molecular dynamics simulations we investigate the binding of these functionalized GAG derivatives to ten representative regulatory proteins including IL-8, IL-10, BMP-2, sclerostin, TIMP-3, CXCL-12, TGF-β, FGF-1, FGF-2, and AT-III, and we establish structure–activity relationships for GAG recognition. Binding is mainly driven by enthalpy with only minor entropic contributions. In several cases binding is determined by GAG length, and in all cases by the position and number of sulfates. Affinities strongly depend on the anomeric modification of the GAG. Highest binding affinities are effected by anomeric functionalization with large fluorophores and by GAG dimerization. Our experimental and theoretical results suggest that the diversity of GAG binding sites and modes is responsible for the observed high affinities and other binding features. The presented new insights into GAG–protein recognition will be of relevance to guide the design of GAG derivatives with customized functions for the engineering of new biomaterials.

Structural basis, stoichiometry, and thermodynamics of binding of the chemokines KC and MIP2 to the glycosaminoglycan heparin

Keratinocyte-derived chemokine (KC or mCXCL1) and macrophage inflammatory protein 2 (MIP2 or mCXCL2) play nonredundant roles in trafficking blood neutrophils to sites of infection and injury. The functional responses of KC and MIP2 are intimately coupled to their interactions with glycosaminoglycans (GAGs). GAG interactions orchestrate chemokine concentration gradients and modulate receptor activity, which together regulate neutrophil trafficking. Here, using NMR, molecular dynamics (MD) simulations, and isothermal titration calorimetry (ITC), we characterized the molecular basis of KC and MIP2 binding to the GAG heparin. Both chemokines reversibly exist as monomers and dimers, and the NMR analysis indicates that the dimer binds heparin with higher affinity. The ITC experiments indicate a stoichiometry of two GAGs per KC or MIP2 dimer and that the enthalpic and entropic contributions vary significantly between the two chemokine-heparin complexes. NMR-based structural models of heparin-KC and heparin-MIP2 complexes reveal that different combinations of residues from the N-loop, 40s turn, β₃-strand, and C-terminal helix form a binding surface within a monomer and that both conserved residues and residues unique to a particular chemokine mediate the binding interactions. MD simulations indicate significant residue-specific differences in their contribution to binding and affinity for a given chemokine and between chemokines. On the basis of our observations that KC and MIP2 bind to GAG via distinct molecular interactions, we propose that the differences in these GAG interactions lead to differences in neutrophil recruitment and play nonoverlapping roles in resolution of inflammation.

Insights into the structure and dynamics of lysyl oxidase propeptide, a flexible protein with numerous partners

Lysyl oxidase (LOX) catalyzes the oxidative deamination of lysine and hydroxylysine residues in collagens and elastin, which is the first step of the cross-linking of these extracellular matrix proteins. It is secreted as a proenzyme activated by bone morphogenetic protein-1, which releases the LOX catalytic domain and its bioactive N-terminal propeptide. We characterized the recombinant human propeptide by circular dichroism, dynamic light scattering, and small-angle X-ray scattering (SAXS), and showed that it is elongated, monomeric, disordered and flexible (D_max: 11.7 nm, R_g: 3.7 nm). We generated 3D models of the propeptide by coarse-grained molecular dynamics simulations restrained by SAXS data, which were used for docking experiments. Furthermore, we have identified 17 new binding partners of the propeptide by label-free assays. They include four glycosaminoglycans (hyaluronan, chondroitin, dermatan and heparan sulfate), collagen I, cross-linking and proteolytic enzymes (lysyl oxidase-like 2, transglutaminase-2, matrix metalloproteinase-2), a proteoglycan (fibromodulin), one growth factor (Epidermal Growth Factor, EGF), and one membrane protein (tumor endothelial marker-8). This suggests new roles for the propeptide in EGF signaling pathway.

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.

Heparin-based hydrogels with tunable sulfation & degradation for anti-inflammatory small molecule delivery

Sustained release of anti-inflammatory agents remains challenging for small molecule drugs due to their low molecular weight and hydrophobicity. Therefore, the goal of this study was to control the release of a small molecule anti-inflammatory agent, crystal violet (CV), from hydrogels fabricated with heparin, a highly sulfated glycosaminoglycan capable of binding positively-charged molecules such as CV. In this system, both electrostatic interactions between heparin and CV and hydrogel degradation were tuned simultaneously by varying the level of heparin sulfation and varying the amount of dithiothreitol within hydrogels, respectively. It was found that heparin sulfation significantly affected CV release, whereby more sulfated heparin hydrogels (Hep and Hep(-N)) released CV with near zero-order release kinetics (R-squared values between 0.96-0.99). Furthermore, CV was released more quickly from fast-degrading hydrogels than slow-degrading hydrogels, providing a method to tune total CV release between 5-15 days while maintaining linear release kinetics. In particular, N-desulfated heparin hydrogels exhibited efficient CV loading (∼90% of originally included CV), near zero-order CV release kinetics, and maintenance of CV bioactivity after release, making this hydrogel formulation a promising CV delivery vehicle for a wide range of inflammatory diseases.

Computational drill down on FGF1-heparin interactions through methodological evaluation

Glycosaminoglycans (GAGs) exhibit a key role in cellular communication processes through interactions with target proteins of the extracellular matrix (ECM). The sandwich-like interaction established between Fibroblast growth factor (FGF) and heparin (HE) represents quite a peculiar protein-GAG-protein system, which has been both structurally and functionally intensively studied. The molecular recognition characteristics of this system have been exploited in various computational studies in order to deepen understanding of GAG-protein interactions. Here, we drill down on the interactions established in this peculiar macromolecular complex by analyzing the applicability of docking techniques and molecular dynamics (MD)-based approaches, and we dissect the molecular recognition properties exhibited by FGF towards a series of HE derivatives. We examine the sensitivity of MM-GBSA free energy calculations in terms of receptor conformational space sampling and changes in the ligand structures. Furthermore, we investigate its predictive power in combination with other computational methods, namely the well-established Autodock3 (AD3) and dynamic molecular docking (DMD), a targeted MD-based docking method specifically developed to account for flexibility and solvent in computer simulations of protein-GAG systems. Our results show that a site-mapping approach can be effectively combined with AD3 and DMD calculations to accurately reproduce available experimental data and, furthermore, to determine specific GAG recognition patterns. This study deepens our understanding of the applicability of available theoretical approaches to the investigation of molecular recognition in protein-GAG systems.

Structural Aspects of Heparan Sulfate Binding to Robo1-Ig1-2

Roundabout 1, or Robo1, is a cell surface signaling molecule important in axon guidance. Its interaction with heparan sulfate (HS) and members of the Slit protein family is essential to its activity, making characterization of these interactions by structural methods, such as NMR, highly desirable. However, the fact that Robo1 is a glycosylated protein prevents employment of commonly used bacterial hosts for expression of properly glycosylated forms with the uniform 15N, 13C, and 2H labeling needed for NMR studies. Here, we apply an alternative methodology, based on labeling with a single amino acid type and high structural content NMR data, to characterize a two-domain construct of glycosylated Robo1 (Robo1-Ig1-2) interacting with a synthetic HS tetramer (IdoA-GlcNS6S-IdoA2S-GlcNS6S-(CH2)5NH2). Significant chemical shift perturbations of the crosspeak from K81 on titration with the tetramer provide initial evidence for the location of a binding site and allow determination of a 255 μM disassociation constant. The binding epitopes, bound conformation, and binding site placement of the HS tetramer have been further characterized by saturation transfer difference (STD), transferred nuclear Overhauser effect (trNOE), and paramagnetic perturbation experiments. A model of the complex has been generated using constraints derived from the various NMR experiments. Postprocessing energetic analysis of this model provides a rationale for the role each glycan residue plays in the binding event, and examination of the binding site in the context of a previous Robo-Slit structure provides a rationale for modulation of Robo-Slit interactions by HS.

Estimating glycosaminoglycan-protein interaction affinity: water dominates the specific antithrombin-heparin interaction

Glycosaminoglycan (GAG)-protein interactions modulate many important biological processes. Structure-function studies on GAGs may reveal probes and drugs, but their structural complexity and highly acidic nature confound such work. Productivity will increase if we are able to identify tight-binding oligosaccharides in silico. An extension of the CHARMM force field is presented to enable modeling of polysaccharides containing sulfamate functionality, and is used to develop a reliable alchemical free-energy perturbation protocol that estimates changes in affinity for the prototypical heparin-antithrombin system to within 2.3 kcal/mol using modest simulation times. Inclusion of water is crucial during simulation as solvation energy was equal in magnitude to the sum of all other thermodynamic factors. In summary, we have identified and optimized a reliable method for estimation of GAG-protein binding affinity, and shown that solvation is a crucial component in GAG-protein interactions.

Computational analysis of interactions in structurally available protein-glycosaminoglycan complexes

Glycosaminoglycans represent a class of linear anionic periodic polysaccharides, which play a key role in a variety of biological processes in the extracellular matrix via interactions with their protein targets. Computationally, glycosaminoglycans are very challenging due to their high flexibility, periodicity and electrostatics-driven nature of the interactions with their protein counterparts. In this work, we carry out a detailed computational characterization of the interactions in protein-glycosaminoglycan complexes from the Protein Data Bank (PDB), which are split into two subsets accounting for their intrinsic nature: non-enzymatic-protein-glycosaminoglycan and enzyme-glycosaminoglycan complexes. We apply molecular dynamics to analyze the differences in these two subsets in terms of flexibility, retainment of the native interactions in the simulations, free energy components of binding and contributions of protein residue types to glycosaminoglycan binding. Furthermore, we systematically demonstrate that protein electrostatic potential calculations, previously found to be successful for glycosaminoglycan binding sites prediction for individual systems, are in general very useful for proposing protein surface regions as putative glycosaminoglycan binding sites, which can be further used for local docking calculations with these particular polysaccharides. Finally, the performance of six different docking programs (Autodock 3, Autodock Vina, MOE, eHiTS, FlexX and Glide), some of which proved to perform well for particular protein-glycosaminoglycan complexes in previous work, is evaluated on the complete protein-glycosaminoglycan data set from the PDB. This work contributes to widen our knowledge of protein-glycosaminoglycan molecular recognition and could be useful to steer a choice of the strategies to be applied in theoretical studies of these systems.

A Simple Method for Discovering Druggable, Specific Glycosaminoglycan-Protein Systems. Elucidation of Key Principles from Heparin/Heparan Sulfate-Binding Proteins

Glycosaminoglycans (GAGs) affect human physiology and pathology by modulating more than 500 proteins. GAG-protein interactions are generally assumed to be ionic and nonspecific, but specific interactions do exist. Here, we present a simple method to identify the GAG-binding site (GBS) on proteins that in turn helps predict high specific GAG-protein systems. Contrary to contemporary thinking, we found that the electrostatic potential at basic arginine and lysine residues neither identifies the GBS consistently, nor its specificity. GBSs are better identified by considering the potential at neutral hydrogen bond donors such as asparagine or glutamine sidechains. Our studies also reveal that an unusual constellation of ionic and non-ionic residues in the binding site leads to specificity. Nature engineers the local environment of Asn45 of antithrombin, Gln255 of 3-O-sulfotransferase 3, Gln163 and Asn167 of 3-O-sulfotransferase 1 and Asn27 of basic fibroblast growth factor in the respective GBSs to induce specificity. Such residues are distinct from other uncharged residues on the same protein structure in possessing a significantly higher electrostatic potential, resultant from the local topology. In contrast, uncharged residues on nonspecific GBSs such as thrombin and serum albumin possess a diffuse spread of electrostatic potential. Our findings also contradict the paradigm that GAG-binding sites are simply a collection of contiguous Arg/Lys residues. Our work demonstrates the basis for discovering specifically interacting and druggable GAG-protein systems based on the structure of protein alone, without requiring access to any structure-function relationship data.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

The Interaction of Heparin Tetrasaccharides with Chemokine CCL5 Is Modulated by Sulfation Pattern and pH

Interactions between chemokines such as CCL5 and glycosaminoglycans (GAGs) are essential for creating haptotactic gradients to guide the migration of leukocytes into inflammatory sites, and the GAGs that interact with CCL5 with the highest affinity are heparan sulfates/heparin. The interaction between CCL5 and its receptor on monocytes, CCR1, is mediated through residues Arg-17 and -47 in CCL5, which overlap with the GAG-binding (44)RKNR(47) "BBXB" motifs. Here we report that heparin and tetrasaccharide fragments of heparin are able to inhibit CCL5-CCR1 binding, with IC50 values showing strong dependence on the pattern and extent of sulfation. Modeling of the CCL5-tetrasaccharide complexes suggested that interactions between specific sulfate and carboxylate groups of heparin and residues Arg-17 and -47 of the protein are essential for strong inhibition; tetrasaccharides lacking the specific sulfation pattern were found to preferentially bind CCL5 in positions less favorable for inhibition of the interaction with CCR1. Simulations of a 12-mer heparin fragment bound to CCL5 indicated that the oligosaccharide preferred to interact simultaneously with both (44)RKNR(47) motifs in the CCL5 homodimer and engaged residues Arg-47 and -17 from both chains. Direct engagement of these residues by the longer heparin oligosaccharide provides a rationalization for its effectiveness as an inhibitor of CCL5-CCR1 interaction. In this mode, histidine (His-23) may contribute to CCL5-GAG interactions when the pH drops just below neutral, as occurs during inflammation. Additionally, an examination of the contribution of pH to modulating CCL5-heparin interactions suggested a need for careful interpretation of experimental results when experiments are performed under non-physiological conditions.

Designing "high-affinity, high-specificity" glycosaminoglycan sequences through computerized modeling

The prediction of high-affinity and/or high-specificity protein-glycosaminoglycan (GAG) interactions is an inherently difficult task, due to several factors including the shallow nature of the typical GAG-binding site and the inherent size, flexibility, diversity, and polydisperse nature of the GAG molecules. Here, we present a generally applicable methodology termed Combinatorial Library Virtual Screening (CVLS) that can identify potential high-affinity, high-specificity protein-GAG interactions from very large GAG combinatorial libraries and a suitable GAG-binding protein. We describe the CVLS approach along with the rationale behind it and provide validation for the method using the well-known antithrombin-thrombin-heparin system.

Enhanced conformational sampling of carbohydrates by Hamiltonian replica-exchange simulation

Knowledge of the structure and conformational flexibility of carbohydrates in an aqueous solvent is important to improving our understanding of how carbohydrates function in biological systems. In this study, we extend a variant of the Hamiltonian replica-exchange molecular dynamics (MD) simulation to improve the conformational sampling of saccharides in an explicit solvent. During the simulations, a biasing potential along the glycosidic-dihedral linkage between the saccharide monomer units in an oligomer is applied at various levels along the replica runs to enable effective transitions between various conformations. One reference replica runs under the control of the original force field. The method was tested on disaccharide structures and further validated on biologically relevant blood group B, Lewis X and Lewis A trisaccharides. The biasing potential-based replica-exchange molecular dynamics (BP-REMD) method provided a significantly improved sampling of relevant conformational states compared with standard continuous MD simulations, with modest computational costs. Thus, the proposed BP-REMD approach adds a new dimension to existing carbohydrate conformational sampling approaches by enhancing conformational sampling in the presence of solvent molecules explicitly at relatively low computational cost.

Fragment-based drug discovery using a multidomain, parallel MD-MM/PBSA screening protocol

We have developed a rigorous computational screening protocol to identify novel fragment-like inhibitors of N(5)-CAIR mutase (PurE), a key enzyme involved in de novo purine synthesis that represents a novel target for the design of antibacterial agents. This computational screening protocol utilizes molecular docking, graphics processing unit (GPU)-accelerated molecular dynamics, and Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) free energy estimations to investigate the binding modes and energies of fragments in the active sites of PurE. PurE is a functional octamer comprised of identical subunits. The octameric structure, with its eight active sites, provided a distinct advantage in these studies because, for a given simulation length, we were able to place eight separate fragment compounds in the active sites to increase the throughput of the MM/PBSA analysis. To validate this protocol, we have screened an in-house fragment library consisting of 352 compounds. The theoretical results were then compared with the results of two experimental fragment screens, Nuclear Magnetic Resonance (NMR) and Surface Plasmon Resonance (SPR) binding analyses. In these validation studies, the protocol was able to effectively identify the competitive binders that had been independently identified by experimental testing, suggesting the potential utility of this method for the identification of novel fragments for future development as PurE inhibitors.

Restricted dynamics of water around a protein-carbohydrate complex: computer simulation studies

Water-mediated protein-carbohydrate interaction is a complex phenomenon responsible for different biological processes in cellular environment. One of the unexplored but important issues in this area is the role played by water during the recognition process and also in controlling the microscopic properties of the complex. In this study, we have carried out atomistic molecular dynamics simulations of a protein-carbohydrate complex formed between the hyaluronan binding domain of the murine Cd44 protein and the octasaccharide hyaluronan in explicit water. Efforts have been made to explore the heterogeneous influence of the complex on the dynamic properties of water present in different regions around it. It is revealed from our analyses that the heterogeneous dynamics of water around the complex are coupled with differential time scales of formation and breaking of hydrogen bonds at the interface. Presence of a highly rigid thin layer of motionally restricted water molecules bridging the protein and the carbohydrate in the common region of the complex has been identified. Such water molecules are expected to play a crucial role in controlling properties of the complex. Importantly, it is demonstrated that the formation of the protein-carbohydrate complex affects the transverse and longitudinal degrees of freedom of the interfacial water molecules in a heterogeneous manner.

On the roles of polyvalent binding in immune recognition: perspectives in the nanoscience of immunology and the immune response to nanomedicines

Immunology often conveys the image of large molecules, either in the soluble state or in the membrane of leukocytes, forming multiple contacts with a target for actions of the immune system. Avidity names the ability of a polyvalent molecule to form multiple connections of the same kind with ligands tethered to the same surface. Polyvalent interactions are vastly stronger than their monovalent equivalent. In the present review, the functional consequences of polyvalent interactions are explored in a perspective of recent theoretical advances in understanding the thermodynamics of such binding. From insights on the structural biology of soluble pattern recognition molecules as well as adhesion molecules in the cell membranes or in their proteolytically shed form, this review documents the prominent role of polyvalent interactions in making the immune system a formidable barrier to microbial infection as well as constituting a significant challenge to the application of nanomedicines.

Biomimetic hydrogels for controlled biomolecule delivery to augment bone regeneration

The regeneration of large bone defects caused by trauma or disease remains a significant clinical problem. Although osteoinductive growth factors such as bone morphogenetic proteins have entered clinics, transplantation of autologous bone remains the gold standard to treat bone defects. The effective treatment of bone defects by protein therapeutics in humans requires quantities that exceed the physiological doses by several orders of magnitude. This not only results in very high treatment costs but also bears considerable risks for adverse side effects. These issues have motivated the development of biomaterials technologies allowing to better control biomolecule delivery from the solid phase. Here we review recent approaches to immobilize biomolecules by affinity binding or by covalent grafting to biomaterial matrices. We focus on biomaterials concepts that are inspired by extracellular matrix (ECM) biology and in particular the dynamic interaction of growth factors with the ECM. We highlight the value of synthetic, ECM-mimicking matrices for future technologies to study bone biology and develop the next generation of 'smart' implants.

Differential effects of fondaparinux and bemiparin on angiogenic and vasculogenesis-like processes

INTRODUCTION

Conventional therapy for venous thromboembolism or acute coronary syndrome involves the administration of glycoanticoagulants (heparins) or oligosaccharides (fondaparinux). We evaluated the effects of such drugs on angiogenesis and vasculogenesis-like models.

MATERIALS AND METHODS

Human umbilical vein endothelial cells or human endothelial progenitor cells were treated with bemiparin, fondaparinux or unfractionated heparin, at concentrations reflecting the doses used in clinical practice. After 24h, cell viability, proliferation, tubule formation and angiogenic molecular mechanisms, such as activation of the serine/threonine kinase AKT, were assessed. In vivo angiogenesis was studied using a Matrigel sponge assay in mice.

RESULTS

Bemiparin gave a significant decrease of in vitro angiogenesis as shown by the reduction of endothelial cell tubule network, while both fondaparinux and unfractionated heparin did not show any significant effect. In assays of Matrigel sponge invasion in mice, unfractionated heparin was able to stimulate angiogenesis and, conversely, bemiparin inhibited angiogenesis. Furthermore, both bemiparin and fondaparinux caused a significant reduction in an in vitro vasculogenesis-like model, as demonstrated by the decrease of tubule network after co-seeding of endothelial progenitor cells and human umbilical vein endothelial cells. In addition, unfractionated heparin but not bemiparin was able to increase AKT phosphorylation.

CONCLUSIONS

In in vitro experiments, bemiparin was the only drug to show an anti-angiogenic and vasculogenic-like effect, unfractionated heparin showed only a trend to increase in angiogenesis assay and fondaparinux affected only the vasculogenesis-like model. Notably, the in vivo experiments corroborated these data. Such results are important for the choice of a patient-tailored therapy.

Conformational flexibility of a protein-carbohydrate complex and the structure and ordering of surrounding water

Protein-carbohydrate non-covalent interactions are important to understand various biological processes in living organisms. One of the important issues in protein-carbohydrate binding is how the protein identifies the target carbohydrate and recognizes its conformational features. Surrounding water molecules are expected to play a critical role not only in mediating the recognition process but also in maintaining the structure of the complex. We carried out atomistic molecular dynamics (MD) simulations of an aqueous solution of the protein-carbohydrate complex formed between the hyaluronan binding domain (HABD) of the murine Cd44 protein and the octasaccharide hyaluronan (HA(8)). The conformational flexibilities of the protein and the carbohydrate, and the microscopic structure and ordering of water molecules around them in the complexed form have been explored. It is revealed that the formation of the complex is associated with significant immobilization of the monosaccharide units of the carbohydrate moiety that are involved in binding. Further, reduction in water densities around the binding residues of the two molecules in the complex with respect to their free forms clearly demonstrated that the recognition between the protein and the carbohydrate is facilitated by removal of a fraction of water molecules from regions around the binding domains.