Sulfation and cation effects on the conformational properties of the glycan backbone of chondroitin sulfate disaccharides

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.

The role of glycosaminoglycans in blood pressure regulation

Essential hypertension (HT) is the global health problem and is a major risk factor for the development of cardiovascular and kidney disease. High salt intake has been associated with HT and impaired kidney sodium excretion is considered to be a major mechanism for the development of HT. Although kidney has a very important role in regulation of BP, this traditional view of BP regulation was challenged by recent findings suggesting that nonosmotic tissue sodium deposition is very important for BP regulation. This new paradigm indicates that sodium can be stored and deposited nonosmotically in the interstitium without water retention and without increased BP. One of the major determinants of this deposition is glycosaminoglycans (GAGs). By binding to GAGs found in the endothelial surface layer (ESL) which contains glycocalyx, sodium is osmotically inactivated and not induce concurrent water retention. Thus, GAGs has important function for homeostatic BP and sodium regulation. In the current review, we summarized the role of GAGs in ESL and BP regulation.

A Review of Chondroitin Sulfate's Preparation, Properties, Functions, and Applications

Chondroitin sulfate (CS) is a natural macromolecule polysaccharide that is extensively distributed in a wide variety of organisms. CS is of great interest to researchers due to its many in vitro and in vivo functions. CS production derives from a diverse number of sources, including but not limited to extraction from various animals or fish, bio-synthesis, and fermentation, and its purity and homogeneity can vary greatly. The structural diversity of CS with respect to sulfation and saccharide content endows this molecule with distinct complexity, allowing for functional modification. These multiple functions contribute to the application of CS in medicines, biomaterials, and functional foods. In this article, we discuss the preparation of CS from different sources, the structure of various forms of CS, and its binding to other relevant molecules. Moreover, for the creation of this article, the functions and applications of CS were reviewed, with an emphasis on drug discovery, hydrogel formation, delivery systems, and food supplements. We conclude that analyzing some perspectives on structural modifications and preparation methods could potentially influence future applications of CS in medical and biomaterial research.

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.

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.

3-O-Sulfation induces sequence-specific compact topologies in heparan sulfate that encode a dynamic sulfation code

Heparan sulfate (HS) is arguably the most diverse linear biopolymer that is known to modulate hundreds of proteins. Whereas the configurational and conformational diversity of HS is well established in terms of varying sulfation patterns and iduronic acid (IdoA) puckers, a linear helical topology resembling a cylindrical rod is the only topology thought to be occupied by the biopolymer. We reasoned that 3-O-sulfation, a rare modification in natural HS, may induce novel topologies that contribute to selective recognition of proteins. In this work, we studied a library of 24 distinct HS hexasaccharides using molecular dynamics (MD). We discovered novel compact (C) topologies that are populated significantly by a unique group of 3-O-sulfated sequences containing IdoA residues. 3-O-sulfated sequences containing glucuronic acid (GlcA) residue and sequences devoid of 3-O-sulfate groups did not exhibit high levels of the C topology and primarily exhibited only the canonical linear (L) form. The C topology arises under dynamical conditions due to rotation around an IdoA → GlcN glycosidic linkage, especially in psi (Ψ) torsion. At an atomistic level, the L → C transformation is a multi-factorial phenomenon engineered to reduce like-charge repulsion, release one or more HS-bound water molecules, and organize a bi-dentate "IdoA-cation-IdoA" interaction. These forces also drive an L → C transformation in a 3-O-sulfated octasaccharide, which has shown evidence of the unique C topology in the co-crystallized state. The 3-O-sulfate-based generation of unique, sequence-specific, compact topologies indicate that natural HS encodes a dynamic sulfation code that could be exploited for selective recognition of target proteins.

Effect of Sulfation on the Conformational Dynamics of Dermatan Sulfate Glycosaminoglycan: A Gaussian Accelerated Molecular Dynamics Study

Glycosaminoglycans (GAGs) are anionic biopolymers present on cell surfaces as a part of proteoglycans. The biological activities of GAGs depend on the sulfation pattern. In our study, we have considered three octadecasaccharide dermatan sulfate (DS) chains with increasing order of sulfation (dp6s, dp7s, and dp12s) to illuminate the role of sulfation on the GAG units and its chain conformation through 10 μs-long Gaussian accelerated molecular dynamics simulations. DS is composed of repeating disaccharide units of iduronic acid (IdoA) and N-acetylgalactosamine (N-GalNAc). Here, N-GalNAc is linked to IdoA via β(1-4), while IdoA is linked to N-GalNAc through α(1-3). With the increase in sulfation, the DS structure becomes more rigid and linear, as is evident from the distribution of root-mean-square deviations (RMSDs) and end-to-end distances. The tetrasaccharide linker region of the main chain shows a rigid conformation in terms of the glycosidic linkage. We have observed that upon sulfation (i.e., dp12s), the ring flip between two chair forms vanished for IdoA. The dynamic cross-correlation analysis reveals that the anticorrelation motions in dp12s are reduced significantly compared to dp6s or dp7s. An increase in sulfation generates relatively more stable hydrogen-bond networks, including water bridging with the neighboring monosaccharides. Despite the favorable linear structures of the GAG chains, our study also predicts few significant bendings related to the different puckering states, which may play a notable role in the function of the DS. The relation between the global conformation with the micro-level parameters such as puckering and water-mediated hydrogen bonds shapes the overall conformational space of GAGs. Overall, atomistic details of the DS chain provided in this study will help understand their functional and mechanical roles, besides developing new biomaterials.

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

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

Pyranose Ring Puckering Thermodynamics for Glycan Monosaccharides Associated with Vertebrate Proteins

The conformational properties of carbohydrates can contribute to protein structure directly through covalent conjugation in the cases of glycoproteins and proteoglycans and indirectly in the case of transmembrane proteins embedded in glycolipid-containing bilayers. However, there continue to be significant challenges associated with experimental structural biology of such carbohydrate-containing systems. All-atom explicit-solvent molecular dynamics simulations provide a direct atomic resolution view of biomolecular dynamics and thermodynamics, but the accuracy of the results depends on the quality of the force field parametrization used in the simulations. A key determinant of the conformational properties of carbohydrates is ring puckering. Here, we applied extended system adaptive biasing force (eABF) all-atom explicit-solvent molecular dynamics simulations to characterize the ring puckering thermodynamics of the ten common pyranose monosaccharides found in vertebrate biology (as represented by the CHARMM carbohydrate force field). The results, along with those for idose, demonstrate that the CHARMM force field reliably models ring puckering across this diverse set of molecules, including accurately capturing the subtle balance between 4C1 and 1C4 chair conformations in the cases of iduronate and of idose. This suggests the broad applicability of the force field for accurate modeling of carbohydrate-containing vertebrate biomolecules such as glycoproteins, proteoglycans, and glycolipids.

The Influences of Sulphation, Salt Type, and Salt Concentration on the Structural Heterogeneity of Glycosaminoglycans

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

Conformational Properties of Glycosaminoglycan Disaccharides: A Molecular Dynamics Study

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

Sulfation and Calcium Favor Compact Conformations of Chondroitin in Aqueous Solutions

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

Constructing 3-Dimensional Atomic-Resolution Models of Nonsulfated Glycosaminoglycans with Arbitrary Lengths Using Conformations from Molecular Dynamics

Glycosaminoglycans (GAGs) are the linear carbohydrate components of proteoglycans (PGs) and are key mediators in the bioactivity of PGs in animal tissue. GAGs are heterogeneous, conformationally complex, and polydisperse, containing up to 200 monosaccharide units. These complexities make studying GAG conformation a challenge for existing experimental and computational methods. We previously described an algorithm we developed that applies conformational parameters (i.e., all bond lengths, bond angles, and dihedral angles) from molecular dynamics (MD) simulations of nonsulfated chondroitin GAG 20-mers to construct 3-D atomic-resolution models of nonsulfated chondroitin GAGs of arbitrary length. In the current study, we applied our algorithm to other GAGs, including hyaluronan and nonsulfated forms of dermatan, keratan, and heparan and expanded our database of MD-generated GAG conformations. Here, we show that individual glycosidic linkages and monosaccharide rings in 10- and 20-mers of hyaluronan and nonsulfated dermatan, keratan, and heparan behave randomly and independently in MD simulation and, therefore, using a database of MD-generated 20-mer conformations, that our algorithm can construct conformational ensembles of 10- and 20-mers of various GAG types that accurately represent the backbone flexibility seen in MD simulations. Furthermore, our algorithm efficiently constructs conformational ensembles of GAG 200-mers that we would reasonably expect from MD simulations.

Efficient Construction of Atomic-Resolution Models of Non-Sulfated Chondroitin Glycosaminoglycan Using Molecular Dynamics Data

Glycosaminoglycans (GAGs) are linear, structurally diverse, conformationally complex carbohydrate polymers that may contain up to 200 monosaccharides. These characteristics present a challenge for studying GAG conformational thermodynamics at atomic resolution using existing experimental methods. Molecular dynamics (MD) simulations can overcome this challenge but are only feasible for short GAG polymers. To address this problem, we developed an algorithm that applies all conformational parameters contributing to GAG backbone flexibility (i.e., bond lengths, bond angles, and dihedral angles) from unbiased all-atom explicit-solvent MD simulations of short GAG polymers to rapidly construct models of GAGs of arbitrary length. The algorithm was used to generate non-sulfated chondroitin 10- and 20-mer ensembles which were compared to MD-generated ensembles for internal validation. End-to-end distance distributions in constructed and MD-generated ensembles have minimal differences, suggesting that our algorithm produces conformational ensembles that mimic the backbone flexibility seen in simulation. Non-sulfated chondroitin 100- and 200-mer ensembles were constructed within a day, demonstrating the efficiency of the algorithm and reduction in time and computational cost compared to simulation.

Rapid Differentiation of Chondroitin Sulfate Isomers by Gas-phase Hydrogen-deuterium Exchange

Chondroitin sulfate (CS)-glycosaminoglycans (GAGs) are linear, negatively charged polysaccharides attached to CS proteoglycans that make up a major component of biological matrices throughout both central and peripheral tissues. The position of their attached sulfate groups to the CS disaccharide is predicted to influence protein-glycan interactions and biological function. Although traditional immunohistochemical analysis of CS-GAGs in biological tissues has provided information regarding changes in GAG abundance during developmental and disease states, quantitative analysis of their specific sulfation patterns is limited due to the inherent complexity of separating CS isomers. While methods have been developed to analyze and quantify sulfation isomers using liquid phase separation, new techniques are still needed to elucidate the full biology of CS-GAGs. Here, we examine ion mobility spectrometry and gas-phase hydrogen-deuterium exchange to resolve positional sulfation isomers in the most common sulfated 4S- and 6S-CS disaccharides. The mobilities for these two isomers are highly similar and could not be resolved effectively with any drift gas tested. In contrast, gas-phase hydrogen-deuterium exchange showed very different rates of deuterium uptake with several deuterium exchange reagents, thereby presenting a promising novel and rapid approach for resolving CS isomers.

Acute body sodium depletion induces skin sodium mobilization in female Wistar rats

NEW FINDINGS

What is the central question of this study? Can Na⁺ depletion mobilize Na⁺ from the skin reservoir in ovariectomized rats? Does oestrogen replacement change the amount and the dynamics of skin Na⁺ storage? Is the reduced salt appetite after Na⁺ depletion in ovariectomized rats with oestrogen replacement related to changes in the skin Na⁺ ? What is the main finding and its importance? This work demonstrated that acute body Na⁺ depletion induced by frusemide mobilized the osmotically inactive skin Na⁺ reservoir to become osmotically active. Oestrogen treatment decreased the induced Na⁺ intake in ovariectomized rats but did not modulate the inactive Na⁺ reservoir in control conditions or its mobilization induced by Na⁺ depletion.

ABSTRACT

Oestradiol, which is an important hormone for water and electrolyte balance, also has a role in the inhibition of induced Na⁺ appetite. Sodium can be stored in the skin in osmotically active or inactive forms, and this skin Na⁺ reservoir may be involved in the control of body Na⁺ levels during physiopathological challenges. In this study, we investigated whether the effect of sodium depletion by frusemide can mobilize Na⁺ from the skin reservoir and whether oestradiol replacement changes or mobilizes the Na⁺ reserves in the skin. Ovariectomized Wistar rats were treated with vehicle or oestradiol for 7 days to evaluate the effects of oestrogen on the hydroelectrolyte balance, intake responses and skin Na⁺ and water content in basal conditions. Furthermore, the effects of oestrogen were evaluated after 24 h frusemide-induced whole-body Na⁺ depletion. Oestradiol-replaced rats exhibited reduced water intake without any significant changes in salt intake, Na⁺ excretion or water and Na⁺ skin content in basal conditions. After sodium depletion, both vehicle- and oestradiol-treated rats exhibited an increase in the osmotically active skin Na⁺ , which was associated with a decrease of the inactive skin Na⁺ reservoir. Oestrogen decreased the hypertonic saline intake induced by Na⁺ depletion, but it was not associated with any significant changes in the skin Na⁺ reservoir. Thus, sodium depletion is able to change the inactive-active skin Na⁺ reservoir balance. However, the oestrogenic modulation of sodium appetite after Na⁺ depletion is probably not related to the action of this hormone in the skin Na⁺ reservoir balance.

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.

Perspective on computational simulations of glycosaminoglycans

Glycosaminoglycans (GAGs) represent a formidable frontier for chemists, biochemists, biologists, medicinal chemists and drug delivery specialists because of massive structural complexity. GAGs are arguably the most complex, natural linear biopolymers with theoretical diversity orders of magnitude higher than proteins and nucleic acids. Yet, this diversity remains generally untapped. Computational approaches offer major routes to understand GAG structure and dynamics so as to enable novel applications of these biopolymers. In fact, computational algorithms, softwares, online tools and techniques have reached a level of sophistication that help understand atomistic details of conformational variation and protein recognition of individual GAG sequences. This review describes current approaches and challenges in computational study of GAGs. It presents a history of major findings since the earliest mention of GAGs (the 1960s), the development of parameters and force fields specific for GAGs, and the application of these tools in understanding GAG structure-function relationship. This review also presents a section on how to perform simulation of GAGs, which is directed toward researchers interested in entering this promising field with potential to impact therapy.

The dual role of the glycosaminoglycan chondroitin-6-sulfate in the development, progression and metastasis of cancer

The remarkable structural heterogeneity of chondroitin sulfate (CS) and dermatan sulfate (DS) generates biological information that can be unique to each of these glycosaminoglycans (GAGs), and changes in their composition are translated into alterations in the binding profiles of these molecules. CS/DS can bind to various cytokines and growth factors, cell surface receptors, adhesion molecules, enzymes and fibrillar glycoproteins of the extracellular matrix, thereby influencing both cell behavior and the biomechanical and biochemical properties of the matrix. In this review, we summarize the current knowledge concerning CS/DS metabolism in the human cancer stroma. The remodeling of the GAG profile in the tumor niche is manifested as a substantial increase in the CS content and a gradual decrease in the proportion between DS and CS. Furthermore, the composition of CS and DS is also affected, which results in a substantial increase in the 6‐O‐sulfated and/or unsulfated disaccharide content, which is concomitant with a decrease in the 4‐O‐sulfation level. Here, we discuss the possible impact of alterations in the CS/DS sulfation pattern on the binding capacity and specificity of these GAGs. Moreover, we propose potential consequences of the stromal accumulation of chondroitin‐6‐sulfate for the progression and metastasis of cancer.

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

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

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

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

High sodium diet converts renal proteoglycans into pro-inflammatory mediators in rats

BACKGROUND

High dietary sodium aggravates renal disease by affecting blood pressure and by its recently shown pro-inflammatory and pro-fibrotic effects. Moreover, pro-inflammatory modification of renal heparan sulfate (HS) can induce tissue remodeling. We aim to investigate if high sodium intake in normotensive rats converts renal HS into a pro-inflammatory phenotype, able to bind more sodium and orchestrate inflammation, fibrosis and lymphangiogenesis.

METHODS

Wistar rats received a normal diet for 4 weeks, or 8% NaCl diet for 2 or 4 weeks. Blood pressure was monitored, and plasma, urine and tissue collected. Tissue sodium was measured by flame spectroscopy. Renal HS and tubulo-interstitial remodeling were studied by biochemical, immunohistochemical and qRT-PCR approaches.

RESULTS

High sodium rats showed a transient increase in blood pressure (week 1; p<0.01) and increased sodium excretion (p<0.05) at 2 and 4 weeks compared to controls. Tubulo-interstitial T-cells, myofibroblasts and mRNA levels of VCAM1, TGF-β1 and collagen type III significantly increased after 4 weeks (all p<0.05). There was a trend for increased macrophage infiltration and lymphangiogenesis (both p = 0.07). Despite increased dermal sodium over time (p<0.05), renal concentrations remained stable. Renal HS of high sodium rats showed increased sulfation (p = 0.05), increased L-selectin binding to HS (p<0,05), and a reduction of sulfation-sensitive anti-HS mAbs JM403 (p<0.001) and 10E4 (p<0.01). Hyaluronan expression increased under high salt conditions (p<0.01) without significant changes in the chondroitin sulfate proteoglycan versican. Statistical analyses showed that sodium-induced tissue remodeling responses partly correlated with observed HS changes.

CONCLUSION

We show that high salt intake by healthy normotensive rats convert renal HS into high sulfated pro-inflammatory glycans involved in tissue remodeling events, but not in increased sodium storage.

Rigidity and flexibility in the tetrasaccharide linker of proteoglycans from atomic-resolution molecular simulation

Proteoglycans (PGs) are covalent conjugates between protein and carbohydrate (glycosaminoglycans). Certain classes of glycosaminoglycans such as chondroitin sulfate/dermatan sulfate and heparan sulfate utilize a specific tetrasaccharide linker for attachment to the protein component: GlcAβ1-3Galβ1-3Galβ1-4Xylβ1-O-Ser. Toward understanding the conformational preferences of this linker, the present work used all-atom explicit-solvent molecular dynamics (MD) simulations combined with Adaptive Biasing Force (ABF) sampling to determine high-resolution, high-precision conformational free energy maps ΔG(φ, ψ) for each glycosidic linkage between constituent disaccharides, including the variant where GlcA is substituted with IdoA. These linkages are characterized by single, predominant (> 97% occupancy), and broad (45° × 60° for ΔG(φ, ψ) < 1 kcal/mol) free-energy minima, while the Xyl-Ser linkage has two such minima similar in free-energy, and additional flexibility from the Ser sidechain dihedral. Conformational analysis of microsecond-scale standard MD on the complete tetrasaccharide-O-Ser conjugate is consistent with ABF data, suggesting (φ, ψ) probabilities are independent of the linker context, and that the tetrasaccharide acts as a relatively rigid unit whereas significant conformational heterogeneity exists with respect to rotation about bonds connecting Xyl to Ser. © 2017 Wiley Periodicals, Inc.

NMR-based conformation and dynamics of a tetrasaccharide-repeating sulfated fucan substituted by different counterions

The sulfated fucan from the sea urchin Lytechinus variegatus is composed of the repetitive sequence [-3)-α-l-Fucp-4( OSO3-)-(1-3)-α-l-Fucp-2,4-di( OSO3-)-(1-3)-α-l-Fucp-2( OSO3-)-(1-3)-α-l-Fucp-2( OSO3-)-(1-]n . Conformation (of rings and chains) and dynamics of this tetrasaccharide-repeating sulfated fucan substituted by Na(+) , Ca(2+) , and Li(+) as counterions have been examined through experiments of liquid-state nuclear magnetic resonance spectroscopy. Scalar coupling and nuclear Overhauser effect (NOE)-based data have confirmed that all composing units occur as (1) C4 chair conformer regardless of the cation type, unit position within the repeating sequence, and sulfation type. Chain conformation determined by NOE signal pattern assisted by molecular modeling for a theoretical octasaccharide has shown a similar linear 3D structure for the three differently substituted forms. Data derived from spin-relaxation measurements have indicated a contribution of counterion type to dynamics. The calcium-based preparation has shown the highest mobility while the sodiated one showed the lowest mobility. The set of results from this work suggests that counterion type can affect the physicochemical properties of the structurally well-defined sulfated fucan. The counterion effect seems to impact more on the structural mobility than on average conformation of the studied sulfated glycan in solution.

Force Fields for Carbohydrate-Divalent Cation Interactions

We report molecular dynamics simulations to study intermolecular interactions for carbohydrate-divalent cation complexes. We observed that common force fields from literature with standard Lorentz-Berthelot combining rules are unable to reproduce the experimental stability constants for model carbohydrate monomer (α-d-Allopyranose) and alkali earth metal cation (Mg(2+), Ca(2+), Sr(2+), or Ba(2+)) complexes. A modified combining rule with rescaled effective cross-interaction radius between cations and the hydroxyl oxygens on the carbohydrates was introduced to reproduce the experimental stability constants, which the preferential carbohydrate-cation complexing structures through the ax-eq-ax sequence of O-1, O-2, and O-3 on α-d-Allopyranose were also observed. The effective radius scaling factor obtained from (α-d-Allopyranose)-Ca(2+) complexes was directly transferrable to the similar six-membered ring (α-d-Ribopyranose)-Ca(2+) complexes; however, reparameterization for the scaling factor may be necessary for the five-membered ring (α-d-Ribofuranose)-Ca(2+) complexes.

Determinants of Glycosaminoglycan (GAG) Structure

Proteoglycans (PGs) are glycosylated proteins of biological importance at cell surfaces, in the extracellular matrix, and in the circulation. PGs are produced and modified by glycosaminoglycan (GAG) chains in the secretory pathway of animal cells. The most common GAG attachment site is a serine residue followed by a glycine (-ser-gly-), from which a linker tetrasaccharide extends and may continue as a heparan sulfate, a heparin, a chondroitin sulfate, or a dermatan sulfate GAG chain. Which type of GAG chain becomes attached to the linker tetrasaccharide is influenced by the structure of the protein core, modifications occurring to the linker tetrasaccharide itself, and the biochemical environment of the Golgi apparatus, where GAG polymerization and modification by sulfation and epimerization take place. The same cell type may produce different GAG chains that vary, depending on the extent of epimerization and sulfation. However, it is not known to what extent these differences are caused by compartmental segregation of protein cores en route through the secretory pathway or by differential recruitment of modifying enzymes during synthesis of different PGs. The topic of this review is how different aspects of protein structure, cellular biochemistry, and compartmentalization may influence GAG synthesis.