Characterization of the interaction of interleukin-8 with hyaluronan, chondroitin sulfate, dermatan sulfate and their sulfated derivatives by spectroscopy and molecular modeling

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.

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.

Heparan sulfates and heparan sulfate proteoglycans in hematopoiesis

ABSTRACT

From signaling mediators in stem cells to markers of differentiation and lineage commitment to facilitators for the entry of viruses, such as HIV-1, cell surface heparan sulfate (HS) glycans with distinct modification patterns play important roles in hematopoietic biology. In this review, we provide an overview of the importance of HS and the proteoglycans (HSPGs) to which they are attached within the major cellular subtypes of the hematopoietic system. We summarize the roles of HSPGs, HS, and HS modifications within each main hematopoietic cell lineage of both myeloid and lymphoid arms. Lastly, we discuss the biological advances in the detection of HS modifications and their potential to further discriminate cell types within hematopoietic tissue.

Structural Elucidation of Glycosaminoglycans in the Tissue of Flounder and Isolation of Chondroitin Sulfate C

Glycosaminoglycans (GAGs) are valuable bioactive polysaccharides with promising biomedical and pharmaceutical applications. In this study, we analyzed GAGs using HPLC-MS/MS from the bone (B), muscle (M), skin (S), and viscera (V) of Scophthalmus maximus (SM), Paralichthysi (P), Limanda ferruginea (LF), Cleisthenes herzensteini (G), Platichthys bicoloratus (PB), Pleuronichthys cornutus (PC), and Cleisthenes herzensteini (CH). Unsaturated disaccharide products were obtained by enzymatic hydrolysis of the GAGs and subjected to compositional analysis of chondroitin sulfate (CS), heparin sulfate (HS), and hyaluronic acid (HA), including the sulfation degree of CS and HS, as well as the content of each GAG. The contents of GAGs in the tissues and the sulfation degree differed significantly among the fish. The bone of S. maximus contained more than 12 μg of CS per mg of dry tissue. Although the fish typically contained high levels of CSA (CS-4S), some fish bone tissue exhibited elevated levels of CSC (CS-6S). The HS content was found to range from 10-150 ug/g, primarily distributed in viscera, with a predominant non-sulfated structure (HS-0S). The structure of HA is well-defined without sulfation modification. These analytical results are independent of biological classification. We provide a high-throughput rapid detection method for tissue samples using HPLC-MS/MS to rapidly screen ideal sources of GAG. On this basis, four kinds of CS were prepared and purified from flounder bone, and their molecular weight was determined to be 23-28 kDa by HPGPC-MALLS, and the disaccharide component unit was dominated by CS-6S, which is a potential substitute for CSC derived from shark cartilage.

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.

Versican Associates with Tumor Immune Phenotype and Limits T-cell Trafficking via Chondroitin Sulfate

Immunotherapies for cancers of epithelial origin have limited efficacy, and a growing body of evidence links the composition of extracellular matrix (ECM) with the likelihood of a favorable response to treatment. The ECM may be considered an immunologic barrier, restricting the localization of cytotoxic immune cells to stromal areas and inhibiting their contact with tumor cells. Identifying ECM components of this immunologic barrier could provide targets that whether degraded in situ may support antitumor immunity and improve immunotherapy response. Using a library of primary triple-negative breast cancer tissues, we correlated CD8+ T-cell tumor contact with ECM composition and identified a proteoglycan, versican (VCAN), as a putative member of the immunologic barrier. Our analysis reveals that CD8+ T-cell contact with tumor associates with the location of VCAN expression, the specific glycovariant of VCAN [defined through the pattern of posttranslational attachments of glycosaminoglycans (GAG)], and the cell types that produce the variant. In functional studies, the isomers of chondroitin sulfate presented on VCAN have opposing roles being either supportive or inhibiting of T-cell trafficking, and removal of the GAGs ameliorates these effects on T-cell trafficking. Overall, we conclude that VCAN can either support or inhibit T-cell trafficking within the tumor microenvironment depending on the pattern of GAGs present, and that VCAN is a major component of the ECM immunologic barrier that defines the type of response to immunotherapy.

SIGNIFICANCE

The response to immunotherapy has been poor toward solid tumors despite immune cells infiltrating into the tumor. The ECM has been associated with impacting T-cell infiltration toward the tumor and in this article we have identified VCAN and its structural modification, chondroitin sulfate as having a key role in T-cell invasion.

Dimerization and Crowding in the Binding of Interleukin 8 to Dendritic Glycosaminoglycans as Artificial Proteoglycans

The interactions of glycosaminoglycans (GAG) with proteins of the extracellular matrix govern and regulate complex physiological functions including cellular growth, immune response, and inflammation. Repetitive presentation of GAG binding motifs, as found in native proteoglycans, might enhance GAG-protein binding through multivalent interactions. Here, we report the chemical synthesis of dendritic GAG oligomers constructed of nonasulfated hyaluronan tetrasaccharides for investigating the binding of the protein chemokine interleukin 8 (IL-8) to artificial, well-defined proteoglycan architectures. Binding of mutant monomeric and native dimerizable IL-8 was investigated by NMR spectroscopy and isothermal titration calorimetry. Dendritic oligomerization of GAG increased the binding affinity of both monomeric and dimeric IL-8. Monomeric IL-8 bound to monomeric and dimeric GAG with KD values of 7.3 and 0.108 μM, respectively. The effect was less pronounced for dimerizable wild-type IL-8, for which GAG dimerization improved the affinity from 34 to 5 nM. Binding of dimeric IL-8 to oligomeric GAG was limited by steric crowding effects, strongly reducing the affinity of subsequent binding events. In conclusion, the strongest effect of GAG oligomerization was the amplified binding of IL-8 monomers, which might concentrate monomeric protein in the extracellular matrix and thus promote protein dimerization under physiological conditions.

Assessing Genetic Algorithm-Based Docking Protocols for Prediction of Heparin Oligosaccharide Binding Geometries onto Proteins

Although molecular docking has evolved dramatically over the years, its application to glycosaminoglycans (GAGs) has remained challenging because of their intrinsic flexibility, highly anionic character and rather ill-defined site of binding on proteins. GAGs have been treated as either fully "rigid" or fully "flexible" in molecular docking. We reasoned that an intermediate semi-rigid docking (SRD) protocol may be better for the recapitulation of native heparin/heparan sulfate (Hp/HS) topologies. Herein, we study 18 Hp/HS-protein co-complexes containing chains from disaccharide to decasaccharide using genetic algorithm-based docking with rigid, semi-rigid, and flexible docking protocols. Our work reveals that rigid and semi-rigid protocols recapitulate native poses for longer chains (5→10 mers) significantly better than the flexible protocol, while 2→4-mer poses are better predicted using the semi-rigid approach. More importantly, the semi-rigid docking protocol is likely to perform better when no crystal structure information is available. We also present a new parameter for parsing selective versus non-selective GAG-protein systems, which relies on two computational parameters including consistency of binding (i.e., RMSD) and docking score (i.e., GOLD Score). The new semi-rigid protocol in combination with the new computational parameter is expected to be particularly useful in high-throughput screening of GAG sequences for identifying promising druggable targets as well as drug-like Hp/HS sequences.

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.

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.

In silico insights into procathepsin S maturation mediated by glycosaminoglycans

Procathepsins, inactive precursors of cathepsins are present in the extracellular matrix (ECM) and in lysosomes. Their active forms are involved in a number of biologically relevant processes, including bone resorption, intracellular proteolysis and regulation of programmed cell death. These processes might be mediated by glycosaminoglycans (GAGs), long unbranched periodic negatively charged polysaccharides. GAGs are also present in ECM and play important role in anticoagulation, angiogenesis and tissue regeneration. GAGs not only mediate the enzymatic activity of cathepsins but can also regulate the process of procathepsin maturation, as it was shown for procathepsin B and S. In this study, we propose the molecular mechanism underlying the biological role of GAGs in procathepsin S maturation and compare our findings with computational data obtained for procathepsin B. We rigorously analyse procathepsin S-GAG complexes in terms of their dynamics, free energy and potential allosteric regulation. We conclude that the GAG binding region might have an effect on the dynamics of procathepsin S structure and so affect its maturation by two different mechanisms.

Sulfotransferase 1C2 promotes hepatocellular carcinoma progression by enhancing glycolysis and fatty acid metabolism

BACKGROUND

Hepatocellular carcinoma (HCC) is aggressive liver cancer. Despite advanced imaging and other diagnostic measures, HCC in a significant portion of patients had reached the advanced stage at the first diagnosis. Unfortunately, there is no cure for advanced HCC. As a result, HCC is still a leading cause of cancer death, and there is a pressing need for new diagnostic markers and therapeutic targets.

METHODS

We investigated sulfotransferase 1C2 (SUTL1C2), which we recently showed was overexpressed in human HCC cancerous tissues. Specifically, we analyzed the effects of SULT1C2 knockdown on the growth, survival, migration, and invasiveness of two HCC cell lines, i.e., HepG2 and Huh7 cells. We also studied the transcriptomes and metabolomes in the two HCC cell lines before and after SULT1C2 knockdown. Based on the transcriptome and metabolome data, we further investigated the SULT1C2 knockdown-mediated shared changes, i.e., glycolysis and fatty acid metabolism, in the two HCC cell lines. Finally, we performed rescue experiments to determine whether the inhibitory effects of SULT1C2 knockdown could be rescued via overexpression.

RESULTS

We showed that SULT1C2 overexpression promoted the growth, survival, migration, and invasiveness of HCC cells. In addition, SULT1C2 knockdown resulted in a wide range of gene expression and metabolome changes in HCC cells. Moreover, analysis of shared alterations showed that SULT1C2 knockdown significantly suppressed glycolysis and fatty acid metabolism, which could be rescued via SULT1C2 overexpression.

CONCLUSIONS

Our data suggest that SULT1C2 is a potential diagnostic marker and therapeutic target for human HCC.

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

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

Molecular Dynamics-Based Comparative Analysis of Chondroitin and Dermatan Sulfates

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

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.

Creating a Functional Biomimetic Cartilage Implant Using Hydrogels Based on Methacrylated Chondroitin Sulfate and Hyaluronic Acid

The load-bearing function of articular cartilage tissue contrasts with the poor load-bearing capacity of most soft hydrogels used for its regeneration. The present study explores whether a hydrogel based on the methacrylated natural polymers chondroitin sulfate (CSMA) and hyaluronic acid (HAMA), injected into warp-knitted spacer fabrics, could be used to create a biomimetic construct with cartilage-like mechanical properties. The swelling ratio of the combined CSMA/HAMA hydrogels in the first 20 days was higher for hydrogels with a higher CSMA concentration, and these hydrogels also degraded quicker, whereas those with a 1.33 wt% of HAMA were stable for more than 120 days. When confined by a polyamide 6 (PA6) spacer fabric, the volumetric swelling of the combined CSMA/HAMA gels (10 wt%, 6.5 × CSMA:HAMA ratio) was reduced by ~53%. Both the apparent peak and the equilibrium modulus significantly increased in the PA6-restricted constructs compared to the free-swelling hydrogels after 28 days of swelling, and no significant differences in the moduli and time constant compared to native bovine cartilage were observed. Moreover, the cell viability in the CSMA/HAMA PA6 constructs was comparable to that in gelatin–methacrylamide (GelMA) PA6 constructs at one day after polymerization. These results suggest that using a HydroSpacer construct with an extracellular matrix (ECM)-like biopolymer-based hydrogel is a promising approach for mimicking the load-bearing properties of native cartilage.

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.

Sulfation Pattern Dependent Iron(III) Mediated Interleukin-8 Glycan Binding

Cytokines such as interleukin-8 activate the immune system during infection and interact with sulfated glycosaminoglycans with specific sulfation patterns. In some cases, these interactions are mediated by metal ion binding which can be used to tune surface-based glycan-protein interactions. We evaluated the effect of both hyaluronan sulfation degree and Fe³⁺ on interleukin-8 binding by electrochemical impedance spectroscopy and surface characterizations. Our results show that sulfation degree and metal ion interactions have a synergistic effect in tuning the electrochemical response of the glycated surfaces to the cytokine.

Hydrolytic (In)stability of Methacrylate Esters in Covalently Cross-Linked Hydrogels Based on Chondroitin Sulfate and Hyaluronic Acid Methacrylate

Chondroitin sulfate (CS) and hyaluronic acid (HA) methacrylate (MA) hydrogels are under investigation for biomedical applications. Here, the hydrolytic (in)stability of the MA esters in these polysaccharides and hydrogels is investigated. Hydrogels made with glycidyl methacrylate-derivatized CS (CSGMA) or methacrylic anhydride (CSMA) degraded after 2-25 days in a cross-linking density-dependent manner (pH 7.4, 37 °C). HA methacrylate (HAMA) hydrogels were stable over 50 days under the same conditions. CS(G)MA hydrogel degradation rates increased with pH, due to hydroxide-driven ester hydrolysis. Desulfated chondroitin MA hydrogels also degrade, indicating that sulfate groups are not responsible for CS(G)MA's hydrolytic sensitivity (pH 7.0-8.0, 37 °C). This sensitivity is likely because CS(G)MA's N-acetyl-galactosamines do not form hydrogen bonds with adjacent glucuronic acid oxygens, whereas HAMA's N-acetyl-glucosamines do. This bond absence allows CS(G)MA higher chain flexibility and hydration and could increase ester hydrolysis sensitivity in CS(G)MA networks. This report helps in biodegradable hydrogel development based on endogenous polysaccharides for clinical applications.

Stabilizing Role of Water Solvation on Anion-π Interactions in Proteins

In this work, anion-π interactions between sulfate groups (SO₄ ^2-) and protein aromatic amino acids (AAs) (histidine protonated (HisP), histidine neutral (HisN), tyrosine (Tyr), tryptophan (Trp), and phenylalanine (Phe)) in an aqueous environment have been analyzed using quantum chemical (QC) calculations and molecular dynamics (MD) simulations. Sulfates can occur naturally in solution and can be contained in biomolecules playing relevant roles in their biological function. In particular, the presence of sulfate groups in glycosaminoglycans such as heparin and heparan sulfate has been shown to be relevant for protein and cellular communication and, consequently, for tissue regeneration. Therefore, anion-π interactions between sulfate groups and aromatic residues represent a relevant aspect to investigate. QC results show that such an anion-π mode of interaction between SO₄ ^2- and aromatic AAs is only possible in the presence of water molecules, in the absence of any other cooperative non-covalent interactions. Protonated histidine stands out in terms of its enhancement in the magnitude of interaction strength on solvation. Other AAs such as non-protonated histidine, tyrosine, and phenylalanine can stabilize anion-π interactions on solvation, albeit with weak interaction energy. Tryptophan does not exhibit any anion-π mode of interaction with SO₄ ^2-. The order of magnitude of the interaction of aromatic AAs with SO₄ ^2- on microsolvation is HisP > HisN > Tyr > Trp > Phe. Atoms in molecules (AIM) analysis illustrates the significance of water molecules in stabilizing the divalent SO₄ ^2- anion over the π surface of the aromatic AAs. MD simulation analysis shows that the order of magnitude of the interaction of SO₄ ^2- with aromatic AAs in macroscopic solvation is HisP > HisN, Tyr, Trp > Phe, which is very much in line with the QC results. Spatial distribution function analysis illustrates that protonated histidine alone is capable of establishing the anion-π interaction with SO₄ ^2- in the solution phase. This study sheds light on the understanding of anion-π interactions between SO₄ ^2- and aromatic AAs such as His and Tyr observed in protein crystal structures and the significance of water molecules in stabilizing such interactions, which is not feasible otherwise.

Modeling Protein-Glycosaminoglycan Complexes: Does the Size Matter?

Docking glycosaminoglycans (GAGs) has been challenging because of the complex nature of these long periodic linear and negatively charged polysaccharides. Although standard docking tools like Autodock3 are successful when docking GAGs up to hexameric length, they experience challenges to properly dock longer GAGs. Similar limitations concern other docking approaches typically developed for docking ligands of limited size to proteins. At the same time, most of more advanced docking approaches are challenging for a user who is inexperienced with complex in silico methodologies. In this work, we evaluate the binding energies of complexes with different lengths of GAGs using all-atom molecular dynamics simulations. Based on this analysis, we propose a new docking protocol for long GAGs that consists of conventional docking of short GAGs and further elongation with the use of a coarse-grained representation of the GAG parts not being in direct contact with its protein receptor. This method automated by a simple script is straightforward to use within the Autodock3 framework but also useful in combination with other standard docking tools. We believe that this method with some minor case-specific modifications could also be used for docking other linear charged polymers.

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.

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.

Immunologic Roles of Hyaluronan in Dermal Wound Healing

Hyaluronic acid (HA), a glycosaminoglycan ubiquitous in the skin, has come into the limelight in recent years for its role in facilitating dermal wound healing. Specifically, HA's length of linearly repeating disaccharides-in other words, its molecular weight (MW)-determines its effects. High molecular weight (HMW)-HA serves an immunosuppressive and anti-inflammatory role, whereas low molecular weight (LMW)-HA contributes to immunostimulation and thus inflammation. During the inflammatory stage of tissue repair, direct and indirect interactions between HA and the innate and adaptive immune systems are of particular interest for their long-lasting impact on wound repair. This review seeks to synthesize the literature on wound healing with a focus on HA's involvement in the immune subsystems.

The potential role of glycosaminoglycans in serum amyloid A fibril formation by in silico approaches

•

SAA dimer binds GAGs stronger than its monomer.

•

Due to its net negative charge SAA prefers to bind short GAG sequences.

•

GAG binding by SAA is electrostatics-driven and rather unspecific.

•

GAG binding site is constituted predominantly by the N-terminal helix residues.

•

GAG binding could potenitally attenuate unfolding of the N-terminal helix.

Serum amyloid A (SAA) is actively involved in such pathological processes as atherosclerosis, rheumatoid arthritis, cancer and Alzheimer's disease by its aggregation. One of the factors that can attenuate its aggregation and so affects its physiological role is its interactions with glycosminoglycans (GAGs), linear anionic periodic polysaccharides. These molecules located in the extracellular matrix of the cell are highly variable in their chemical composition and sulfation patterns. Despite the available experimental evidence of SAA-GAG interactions, no mechanistic details at atomic level have been reported for these systems so far. In our work we aimed to apply diverse computational tools to characterize SAA-GAG complexes formation and to answer questions about their potential specificity, energetic patterns, particular SAA residues involved in these interactions, favourable oligomeric state of the protein and the potential influence of GAGs on SAA aggregation. Molecular docking, conventional and replica exchange molecular dynamics approaches were applied to corroborate the experimental knowledge and to propose the corresponding molecular models. SAA-GAG complex formation was found to be electrostatics-driven and rather unspecific of a GAG sulfation pattern, more favorable for the dimer than for the monomer when binding to a short GAG oligosaccharide through its N-terminal helix, potentially contributing to the unfolding of this helix, which could lead to the promotion of the protein aggregation. The data obtained add to the specific knowledge on SAA-GAG systems and deepen the general understanding of protein-GAG interactions that is of a considerable value for the development of GAG-based approaches in a broad theurapeutic context.

Adipocyte-Derived Versican and Macrophage-Derived Biglycan Control Adipose Tissue Inflammation in Obesity

Obesity is characterized by adipose tissue inflammation. Because proteoglycans regulate inflammation, here we investigate their role in adipose tissue inflammation in obesity. We find that adipose tissue versican and biglycan increase in obesity. Versican is produced mainly by adipocytes and biglycan by adipose tissue macrophages. Both proteoglycans are also present in adipose tissue from obese human subjects undergoing gastric bypass surgery. Deletion of adipocyte-specific versican or macrophage-specific biglycan in mice reduces macrophage accumulation and chemokine and cytokine expression, although only adipocyte-specific versican deletion leads to sustained improvement in glucose tolerance. Macrophage-derived biglycan activates inflammatory genes in adipocytes. Versican expression increases in cultured adipocytes exposed to excess glucose, and adipocyte-conditioned medium stimulates inflammation in resident peritoneal macrophages, in part because of a versican breakdown product, versikine. These findings provide insights into the role of adipocyte- and macrophage-derived proteoglycans in adipose tissue inflammation in obesity.

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.

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.

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.

Characterization of defined sulfated heparin-like oligosaccharides by electrospray ionization ion trap mass spectrometry

Glycosaminoglycans (GAG) as long, unbranched polysaccharides are major components of the extracellular matrix. Many studies provided additional evidence of a specific binding between mediators and sulfated GAG, at which the sulfation code-which means the number and positions of sulfate groups along the polysaccharide chain-plays an important role. GAG from natural sources are very inhomogeneous regarding their sulfation patterns and molecular weight. Additionally, there is a high risk of contamination. This results in a growing interest in the careful characterization of native GAG and the synthesis of artificial GAG. Additionally, chemically oversulfated GAG analogues show many favorable properties. However, the structural characterization of these carbohydrates by mass spectrometry remains challenging. One significant problem is the sulfate loss during the ionization, which increases with the number of sulfate residues. We used the sulfated pentasaccharide fondaparinux as model substance to optimize sample preparation and measurement conditions, compared different established desalination methods and already existing protocols for sulfated oligosaccharides, and investigated their impact on the quality of the mass spectra. After optimization of the measurement conditions, we could establish a gentle and fast protocol for the mass spectrometry characterization of (fully) sulfated, artificial GAG-like oligosaccharides with minimized sulfate loss in the positive and negative ion mode. Here, the negative ion mode was more sensitive in comparison with the positive one, and fondaparinux species with sulfate loss were not detectable under the optimized conditions in the positive ion mode.

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.

In vitro Studies of Transendothelial Migration for Biological and Drug Discovery

Inflammatory diseases and cancer metastases lack concrete pharmaceuticals for their effective treatment despite great strides in advancing our understanding of disease progression. One feature of these disease pathogeneses that remains to be fully explored, both biologically and pharmaceutically, is the passage of cancer and immune cells from the blood to the underlying tissue in the process of extravasation. Regardless of migratory cell type, all steps in extravasation involve molecular interactions that serve as a rich landscape of targets for pharmaceutical inhibition or promotion. Transendothelial migration (TEM), or the migration of the cell through the vascular endothelium, is a particularly promising area of interest as it constitutes the final and most involved step in the extravasation cascade. While in vivo models of cancer metastasis and inflammatory diseases have contributed to our current understanding of TEM, the knowledge surrounding this phenomenon would be significantly lacking without the use of in vitro platforms. In addition to the ease of use, low cost, and high controllability, in vitro platforms permit the use of human cell lines to represent certain features of disease pathology better, as seen in the clinic. These benefits over traditional pre-clinical models for efficacy and toxicity testing are especially important in the modern pursuit of novel drug candidates. Here, we review the cellular and molecular events involved in leukocyte and cancer cell extravasation, with a keen focus on TEM, as discovered by seminal and progressive in vitro platforms. In vitro studies of TEM, specifically, showcase the great experimental progress at the lab bench and highlight the historical success of in vitro platforms for biological discovery. This success shows the potential for applying these platforms for pharmaceutical compound screening. In addition to immune and cancer cell TEM, we discuss the promise of hepatocyte transplantation, a process in which systemically delivered hepatocytes must transmigrate across the liver sinusoidal endothelium to successfully engraft and restore liver function. Lastly, we concisely summarize the evolving field of porous membranes for the study of TEM.

The abnormal accumulation of heparan sulfate in patients with mucopolysaccharidosis prevents the elastolytic activity of cathepsin V

Mucopolysaccharidosis (MPS) are rare inherited diseases characterized by accumulation of lysosomal glycosaminoglycans, including heparan sulfate (HS). Patients exhibit progressive multi-visceral dysfunction and shortened lifespan mainly due to a severe cardiac/respiratory decline. Cathepsin V (CatV) is a potent elastolytic protease implicated in extracellular matrix (ECM) remodeling. Whether CatV is inactivated by HS in lungs from MPS patients remained unknown. Herein, CatV colocalized with HS in MPS bronchial epithelial cells. HS level correlated positively with the severity of respiratory symptoms and negatively to the overall endopeptidase activity of cysteine cathepsins. HS bound tightly to CatV and impaired its activity. Withdrawal of HS by glycosidases preserved exogenous CatV activity, while addition of Surfen, a HS antagonist, restored elastolytic CatV-like activity in MPS samples. Our data suggest that the pathophysiological accumulation of HS may be deleterious for CatV-mediated ECM remodeling and for lung tissue homeostasis, thus contributing to respiratory disorders associated to MPS diseases.

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.

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.

Targeting Chemokine-Glycosaminoglycan Interactions to Inhibit Inflammation

Leukocyte migration into tissues depends on the activity of chemokines that form concentration gradients to guide leukocytes to a specific site. Interaction of chemokines with their specific G protein-coupled receptors (GPCRs) on leukocytes induces leukocyte adhesion to the endothelial cells, followed by extravasation of the leukocytes and subsequent directed migration along the chemotactic gradient. Interaction of chemokines with glycosaminoglycans (GAGs) is crucial for extravasation in vivo. Chemokines need to interact with GAGs on endothelial cells and in the extracellular matrix in tissues in order to be presented on the endothelium of blood vessels and to create a concentration gradient. Local chemokine retention establishes a chemokine gradient and prevents diffusion and degradation. During the last two decades, research aiming at reducing chemokine activity mainly focused on the identification of inhibitors of the interaction between chemokines and their cognate GPCRs. This approach only resulted in limited success. However, an alternative strategy, targeting chemokine-GAG interactions, may be a promising approach to inhibit chemokine activity and inflammation. On this line, proteins derived from viruses and parasites that bind chemokines or GAGs may have the potential to interfere with chemokine-GAG interactions. Alternatively, chemokine mimetics, including truncated chemokines and mutant chemokines, can compete with chemokines for binding to GAGs. Such truncated or mutated chemokines are characterized by a strong binding affinity for GAGs and abrogated binding to their chemokine receptors. Finally, Spiegelmers that mask the GAG-binding site on chemokines, thereby preventing chemokine-GAG interactions, were developed. In this review, the importance of GAGs for chemokine activity in vivo and strategies that could be employed to target chemokine-GAG interactions will be discussed in the context of inflammation.

Rat cathepsin K: Enzymatic specificity and regulation of its collagenolytic activity

Human cathepsin K (hCatK), which is highly expressed in osteoclasts, has the noteworthy ability to cleave type I and II collagens in their helical domain. Its collagenase potency depends strictly on the formation of an oligomeric complex with chondroitin 4-sulfate (C4-S). Accordingly, hCatK is a pivotal protease involved in bone resorption and is an attractive target for the treatment of osteoporosis. As rat is a common animal model for the evaluation of hCatK inhibitors, we conducted a comparative analysis of rat CatK (rCatK) and hCatK, which share a high degree of identity (88%) and similarity (93%). The pH activity profile of both enzymes displayed a similar bell-shaped curve (optimal pH: 6.4). Presence of Ser134 and Val160 in the S2 pocket of rCatK instead of Ala and Leu residues, respectively, in hCatK, led to a weaker peptidase activity, as observed for mouse CatK. Also, regardless of the presence of C4-S, rCatK cleaved in the nonhelical telopeptide regions of both type I (tail) and type II (articular joint) rat collagens. Structure-based computational analyses (electrostatic potential, molecular docking, molecular dynamics, free energy calculations) sustained that the C4-S mediated collagenolytic activity of rCatK obeys distinct molecular interactions from those of hCatK. Additionally, T-kininogen (a.k.a. thiostatin), a unique rat serum acute phase molecule, acted as a tight-binding inhibitor of hCatK (Ki = 0.11 ± 0.05 nM). Taken into account the increase of T-Kininogen level in inflamed rat sera, this may raise the question of the appropriateness to evaluate pharmacological hCatK inhibitors in this peculiar animal model.

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.

Structural basis of chemokine interactions with heparan sulfate, chondroitin sulfate, and dermatan sulfate

Chemokines play diverse roles in human pathophysiology, ranging from trafficking leukocytes and immunosurveillance to the regulation of metabolism and neural function. Chemokine function is intimately coupled to binding tissue glycosaminoglycans (GAGs), heparan sulfate (HS), chondroitin sulfate (CS), and dermatan sulfate (DS). Currently, very little is known about how the structural features and sequences of a given chemokine, the structure and sulfation pattern of a given GAG, and structural differences among GAGs and among chemokines impact binding interactions. In this study, we used solution NMR spectroscopy to characterize the binding interactions of two related neutrophil-activating chemokines, CXCL1 and CXCL5, with HS, CS, and DS. For both chemokines, the dimer bound all three GAGs with higher affinity than did the monomer, and affinities of the chemokines for CS and DS were lower than for HS. NMR-based structural models reveal diverse binding geometries and show that the binding surfaces for each of the three GAGs were different between the two chemokines. However, a given chemokine had similar binding interactions with CS and DS that were different from HS. Considering the fact that CXCL1 and CXCL5 activate the same CXCR2 receptor, we conclude that GAG interactions play a role in determining the nature of chemokine gradients, levels of free chemokine available for receptor activation, how chemokines bind their receptors, and that differences in these interactions determine chemokine-specific function.

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.

Tick saliva protein Evasin-3 modulates chemotaxis by disrupting CXCL8 interactions with glycosaminoglycans and CXCR2

Chemokines are a group of chemotaxis proteins that regulate cell trafficking and play important roles in immune responses and inflammation. Ticks are blood-sucking parasites that secrete numerous immune-modulatory agents in their saliva to evade host immune responses. Evasin-3 is a small salivary protein that belongs to a class of chemokine-binding proteins isolated from the brown dog tick, Rhipicephalus sanguineus. Evasin-3 has been shown to have a high affinity for chemokines CXCL1 and CXCL8 and to diminish inflammation in mice. In the present study, solution NMR spectroscopy was used to investigate the structure of Evasin-3 and its CXCL8–Evasin-3 complex. Evasin-3 is found to disrupt the glycosaminoglycan-binding site of CXCL8 and inhibit the interaction of CXCL8 with CXCR2. Structural data were used to design two novel CXCL8-binding peptides. The linear tEv3 17–56 and cyclic tcEv3 16–56 dPG Evasin-3 variants were chemically synthesized by solid-phase peptide synthesis. The affinity of these newly synthesized variants to CXCL8 was measured by surface plasmon resonance biosensor analysis. The K_d values of tEv3 17–56 and tcEv3 16–56 dPG were 27 and 13 nm, respectively. Both compounds effectively inhibited CXCL8-induced migration of polymorphonuclear neutrophils. The present results suggest utility of synthetic Evasin-3 variants as scaffolds for designing and fine-tuning new chemokine-binding agents that suppress immune responses and inflammation.

Protein Interaction with Charged Macromolecules: From Model Polymers to Unfolded Proteins and Post-Translational Modifications

Interaction of proteins with charged macromolecules is involved in many processes in cells. Firstly, there are many naturally occurred charged polymers such as DNA and RNA, polyphosphates, sulfated glycosaminoglycans, etc., as well as pronouncedly charged proteins such as histones or actin. Electrostatic interactions are also important for "generic" proteins, which are not generally considered as polyanions or polycations. Finally, protein behavior can be altered due to post-translational modifications such as phosphorylation, sulfation, and glycation, which change a local charge of the protein region. Herein we review molecular modeling for the investigation of such interactions, from model polyanions and polycations to unfolded proteins. We will show that electrostatic interactions are ubiquitous, and molecular dynamics simulations provide an outstanding opportunity to look inside binding and reveal the contribution of electrostatic interactions. Since a molecular dynamics simulation is only a model, we will comprehensively consider its relationship with the experimental data.