Conformer populations of L-iduronic acid residues in glycosaminoglycan sequences

Cationic Polysaccharides Bind to the Endothelial Cell Surface Extracellular Matrix Involving Heparan Sulfate

Cationic polysaccharides have been extensively studied for drug delivery via the bloodstream, yet few have progressed to clinical use. Endothelial cells lining the blood vessel wall are coated in an anionic extracellular matrix called the glycocalyx. However, we do not fully comprehend the charged polysaccharide interactions with the glycocalyx. We reveal that the cationic polysaccharide poly(acetyl, arginyl) glucosamine (PAAG) exhibits the highest association with the endothelial glycocalyx, followed by dextran (neutral) and hyaluronan (anionic). Furthermore, we demonstrate that PAAG binds heparan sulfate (HS) within the glycocalyx, leading to intracellular accumulation. Using an in vitro glycocalyx model, we demonstrate a charge-based extent of association of polysaccharides with HS. Mechanistically, we observe that PAAG binding to HS occurs via a condensation reaction and functionally protects HS from degradation. Together, this study reveals the interplay between polysaccharide charge properties and interactions with the endothelial cell glycocalyx toward improved delivery system design and application.

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

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

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.

Synthesis and ab initio conformational investigation of a series of model sulfated α-L-iduronopyranosides

Due to the all-axial orientation of the OH-groups in the 1C4 chair conformation considered standard for L-hexapyranosides, including l-iduronopyranoside - a component of many biologically and medically significant sulfated glycans, these monosaccharides can be anticipated to display unusual conformations upon the introduction of bulky and charged substituents. Herein we describe the synthesis of a series of iduronopyranoside derivatives with varying sulfation patterns, which were studied computationally using the DLPNO-MP2 approach and by means of analyzing their chemical shifts to ascertain the effects sulfation has on the conformation of the iduronopyranoside ring.

Epoxide-Mediated Trans-Thioglycosylation and Application to the Synthesis of Oligosaccharides Related to the Capsular Polysaccharides of C. jejuni HS:4

The enzyme-resistant thioglycosides are highly valuable immunogens because of their enhanced metabolic stability. We report the first synthesis of a family of thiooligosaccharides related to the capsular polysaccharides (CPS) of Campylobacter jejuni HS:4 for potential use in conjugate vaccines. The native CPS structures of the pathogen consist of a challenging repeating disaccharide formed with β(1→4)-linked 6-deoxy-β-D-ido-heptopyranoside and N-acetyl-D-glucosamine; the rare 6-deoxy-ido-heptopyranosyl backbone and β-anomeric configuration of the former monosaccharide makes the synthesis of this family of antigens very challenging. So far, no synthesis of the thioanalogs of the CPS antigens have been reported. The unprecedented synthesis presented in this work is built on an elegant approach by using β-glycosylthiolate as a glycosyl donor to open the 2,3-epoxide functionality of pre-designed 6-deoxy-β-D-talo-heptopyranosides. Our results illustrated that this key trans-thioglycosylation can be designed in a modular and regio and stereo-selective manner. Built on the success of this novel approach, we succeeded the synthesis of a family of thiooligosaccharides including a thiohexasaccharide which is considered to be the desired antigen length and complexity for immunizations. We also report the first direct conversion of base-stable but acid-labile 2-trimethylsilylethyl glycosides to glycosyl-1-thioacetates in a one-pot manner.

Synthesis of 6-deoxy-d-ido-heptopyranose-containing fragments of the Campylobacter jejuni strain CG8486 capsular polysaccharide

Campylobacters are important causes of gastrointestinal illness and the capsular polysaccharides (CPS) they produce are key virulence factors and targets for vaccine development. We report here the synthesis of two fragments of the Campylobacter jejuni CG8486 strain CPS that contain a rare 6-deoxy-d-ido-heptopyranose residue and, in one target, two O-methyl phosphoramidate (MeOPN) motifs. The synthetic approach features the stereoselective construction of the β-d-ido-heptopyranoside linkage via glycosylation with a β-d-galacto-heptopyranoside donor followed by a one-pot sequential C-2 and C-3 inversion. During the syntheses, we uncovered a number of interesting conformational effects with regard to the 6-deoxy-ido-heptopyranose ring, the glycosidic linkage connecting the two monosaccharides, and the MeOPN groups.

Solvent Model Benchmark for Molecular Dynamics of Glycosaminoglycans

In computational studies of glycosaminoglycans (GAGs), a group of anionic, periodic linear polysaccharides, so far there has been very little discussion about the role of solvent models in the molecular dynamics simulations of these molecules. Predominantly, the TIP3P water model is commonly used as one of the most popular explicit water models in general. However, there are numerous alternative explicit and implicit water models that are neglected in the computational research of GAGs. Since solvent-mediated interactions are particularly important for GAG dynamic and structural properties, it would be of great interest for the GAG community to establish the solvent model that is suited the best in terms of the quality of theoretically obtained GAG parameters and, at the same time, would be reasonably demanding in terms of computational resources required. In this study, heparin (HP) was simulated using five implicit and six explicit solvent models with the aim to find out how different solvent models influence HP's molecular descriptors in the molecular dynamics simulations. Here, we initiate the search for the most appropriate solvent representation for GAG systems and we hope to encourage other groups to contribute to this highly relevant subject.

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.

Evidence for Multiple Binding Modes in the Initial Contact Between SARS-CoV-2 Spike S1 Protein and Cell Surface Glycans

Infection of host cells by SARS-CoV-2 begins with recognition by the virus S (spike) protein of cell surface heparan sulfate (HS), tethering the virus to the extracellular matrix environment, and causing the subunit S1-RBD to undergo a conformational change into the 'open' conformation. These two events promote the binding of S1-RBD to the angiotensin converting enzyme 2 (ACE2) receptor, a preliminary step toward viral-cell membrane fusion. Combining ligand-based NMR spectroscopy with molecular dynamics, oligosaccharide analogues were used to explore the interactions between S1-RBD of SARS CoV-2 and HS, revealing several low-specificity binding modes and previously unidentified potential sites for the binding of extended HS polysaccharide chains. The evidence for multiple binding modes also suggest that highly specific inhibitors will not be optimal against protein S but, rather, diverse HS-based structures, characterized by high affinity and including multi-valent compounds, may be required.

Efficiently Computing NMR 1H and 13C Chemical Shifts of Saccharides in Aqueous Environment

Determining the structure of saccharides in their native environment is crucial to understanding their function and more accurately targeting their utilization. Nuclear magnetic resonance observables such as the nuclear Overhauser effect or spin-spin coupling constants are routinely utilized to study saccharides in their native water environment. However, while highly sensitive to the local environment, chemical shifts are mostly overlooked, despite being commonly measured for compounds identification. Although chemical shifts carry considerable structural information, their direct association with structure is notoriously difficult. This is mostly due to the similarity in the chemical nature of most saccharides causing similar physicochemical environments close to sugar C and H atoms, resulting in comparable chemical shifts. The rise of computational power allows one to compute reliable chemical shifts and use them to determine atomistic details of these sugars in solution. However, any prediction is severely limited by the computational protocol used and its accuracy. In this work, we studied a set of 31 saccharides on which we evaluated various computational protocols to calculate the total number of 375 ¹H and 327 ¹³C chemical shifts of sugars in an aqueous environment. Our study proposes two cost-effective protocols for simulating ¹H and ¹³C chemical shifts that we recommend for further use. These protocols can help with the interpretation of experimental spectra, but we also show that they are also capable of structure prediction independently. This is possible because of the low mean absolute deviations of calculated shifts from the experiment (0.06 ppm for ¹H and 1.09 ppm for ¹³C). We explore different solvation methods, basis sets, and optimization schemes to reach such accuracy. A correct sampling of the conformation phase space of flexible sugar molecules is also key to obtaining accurately converged theoretical chemical shifts. The linear regression method was applied to convert the calculated isotropic nuclear magnetic shielding constants to simulated chemical shifts comparable with the experiment. The achieved level of accuracy can help in utilizing chemical shifts for elucidating the 3D atomistic structure of saccharides in aqueous solutions. All linear regression parameters obtained on our extensive set of sugars for all the tested protocols can be reutilized in future works.

Enhanced Antiviral Function of Magnesium Chloride-Modified Heparin on a Broad Spectrum of Viruses

Previous studies reported on the broad-spectrum antiviral function of heparin. Here we investigated the antiviral function of magnesium-modified heparin and found that modified heparin displayed a significantly enhanced antiviral function against human adenovirus (HAdV) in immortalized and primary cells. Nuclear magnetic resonance analyses revealed a conformational change of heparin when complexed with magnesium. To broadly explore this discovery, we tested the antiviral function of modified heparin against herpes simplex virus type 1 (HSV-1) and found that the replication of HSV-1 was even further decreased compared to aciclovir. Moreover, we investigated the antiviral effect against the new severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) and measured a 55-fold decreased viral load in the supernatant of infected cells associated with a 38-fold decrease in virus growth. The advantage of our modified heparin is an increased antiviral effect compared to regular heparin.

Neural Tissue Homeostasis and Repair Is Regulated via CS and DS Proteoglycan Motifs

Chondroitin sulfate (CS) is the most abundant and widely distributed glycosaminoglycan (GAG) in the human body. As a component of proteoglycans (PGs) it has numerous roles in matrix stabilization and cellular regulation. This chapter highlights the roles of CS and CS-PGs in the central and peripheral nervous systems (CNS/PNS). CS has specific cell regulatory roles that control tissue function and homeostasis. The CNS/PNS contains a diverse range of CS-PGs which direct the development of embryonic neural axonal networks, and the responses of neural cell populations in mature tissues to traumatic injury. Following brain trauma and spinal cord injury, a stabilizing CS-PG-rich scar tissue is laid down at the defect site to protect neural tissues, which are amongst the softest tissues of the human body. Unfortunately, the CS concentrated in gliotic scars also inhibits neural outgrowth and functional recovery. CS has well known inhibitory properties over neural behavior, and animal models of CNS/PNS injury have demonstrated that selective degradation of CS using chondroitinase improves neuronal functional recovery. CS-PGs are present diffusely in the CNS but also form denser regions of extracellular matrix termed perineuronal nets which surround neurons. Hyaluronan is immobilized in hyalectan CS-PG aggregates in these perineural structures, which provide neural protection, synapse, and neural plasticity, and have roles in memory and cognitive learning. Despite the generally inhibitory cues delivered by CS-A and CS-C, some CS-PGs containing highly charged CS disaccharides (CS-D, CS-E) or dermatan sulfate (DS) disaccharides that promote neural outgrowth and functional recovery. CS/DS thus has varied cell regulatory properties and structural ECM supportive roles in the CNS/PNS depending on the glycoform present and its location in tissue niches and specific cellular contexts. Studies on the fruit fly, Drosophila melanogaster and the nematode Caenorhabditis elegans have provided insightful information on neural interconnectivity and the role of the ECM and its PGs in neural development and in tissue morphogenesis in a whole organism environment.

Midkine Interaction with Chondroitin Sulfate Model Synthetic Tetrasaccharides and Their Mimetics: The Role of Aromatic Interactions

Midkine (MK) is a neurotrophic factor that participates in the embryonic central nervous system (CNS) development and neural stem cell regulation, interacting with sulfated glycosaminoglycans (GAGs). Chondroitin sulfate (CS) is the natural ligand in the CNS. In this work, we describe the interactions between a library of synthetic models of CS‐types and mimics. We did a structural study of this library by NMR and MD (Molecular Dynamics), concluding that the basic shape is controlled by similar geometry of the glycosidic linkages. Their 3D structures are a helix with four residues per turn, almost linear. We have studied the tetrasaccharide‐midkine complexes by ligand observed NMR techniques and concluded that the shape of the ligands does not change upon binding. The ligand orientation into the complex is very variable. It is placed inside the central cavity of MK formed by the two structured beta‐sheets domains linked by an intrinsically disordered region (IDR). Docking analysis confirmed the participation of aromatics residues from MK completed with electrostatic interactions. Finally, we test the biological activity by increasing the MK expression using CS tetrasaccharides and their capacity in enhancing the growth stimulation effect of MK in NIH3T3 cells.

From Cancer to COVID-19: A Perspective on Targeting Heparan Sulfate-Protein Interactions

Heparan sulfate (HS) is a complex, polyanionic polysaccharide ubiquitously expressed on cell surfaces and in the extracellular matrix. HS interacts with numerous proteins to mediate a vast array of biological and pathological processes. Inhibition of HS‐protein interactions is thus an attractive approach for new therapeutic development for cancer and infectious diseases, including COVID‐19; however, synthesis of well‐defined native HS oligosaccharides remains challenging. This has aroused significant interest in the development of HS mimetics which are more synthetically tractable and have fewer side effects, such as undesired anticoagulant activity. This account provides a perspective on the design and synthesis of different classes of HS mimetics with useful properties, and the development of various assays and molecular modelling tools to progress our understanding of their interactions with HS‐binding proteins.

Conformational Analysis of Heparin-Analogue Pentasaccharides by Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulations

Elucidation and improvement of the blood coagulant properties of heparin are the focus of intense research. In this study, we performed conformational analysis using nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations on the heparin pentasaccharide analogue idraparinux, its disulfonatomethyl analogue, which features a slightly improved blood coagulation property, and a trisulfonatomethyl analogue, in which the activity has been totally abolished. As the ring conformation of the G subunit has been suggested as a major determinant of the biological properties, we analyzed the sugar ring conformations and dynamics of the interglycosidic linkages. We found that the conformation of the G ring is dominated by the ²S_O skewed boat next to the ¹C₄ chair in all three derivatives. Both the thermodynamics and the kinetics of the conformational states were found to be highly similar in the three derivatives. Molecular kinetic analysis showed that the ²S_O skewed boat state of the G ring is equally favorable in the three analogues, resulting in similar ²S_O populations. Also, the transition kinetics from the ¹C₄ chair to the ²S_O skewed boat was found to be comparable in the derivatives, which indicates a similar energy barrier between the two states of the G subunit. We also identified a slower conformational transition between the dominant ⁴C₁ chair and the boat conformations on the E subunit. Both G and E ring flips are also accompanied by changes along the interglycosidic linkages, which take place highly synchronously with the ring flips. These findings indicate that conformational plasticity of the G ring and the dominance of the ²S_O skewed boat populations do not necessarily warrant the biological activity of the derivatives and hence the impact of other factors also needs to be considered.

Conformational Modulation of Iduronic Acid-Containing Sulfated Glycosaminoglycans by a Polynuclear Platinum Compound and Implications for Development of Antimetastatic Platinum Drugs

1 H NMR spectroscopic studies on the 1:1 adduct of the pentasaccharide Fondaparinux (FPX) and the substitution-inert polynuclear platinum complex TriplatinNC show significant modulation of geometry around the glycosidic linkages of the FPX constituent monosaccharides. FPX is a valid model for the highly sulfated cell signalling molecule heparan sulfate (HS). The conformational ratio of the 1 C4 :2 S0 forms of the FPX residue IdoA(2S) is altered from ca. 35:65 (free FPX) to ca. 75:25 in the adduct; the first demonstration of a small molecule affecting conformational changes on a HS oligosaccharide. Functional consequences of such binding are suggested to be inhibition of HS cleavage in MDA-MB-231 triple-negative breast cancer (TNBC) cells. We further describe inhibition of metastasis by TriplatinNC in the TNBC 4T1 syngeneic tumour model. Our work provides insight into a novel approach for design of platinum drugs (and coordination compounds in general) with intrinsic anti-metastatic potential.

BMP6 binding to heparin and heparan sulfate is mediated by N-terminal and C-terminal clustered basic residues

BACKGROUND

The bone morphogenetic protein 6 (BMP6) is a crucial inducer of hepcidin, the peptide hormone that regulates the iron availability in our body. Hepcidin expression is influenced by hepatic heparan sulfate (HS) and by heparin administration, suggesting BMP6 interaction with heparin/HS. The BMP2/4 subfamily has been deeply characterized to have a N-terminal heparin/HS binding domain (HBD), whose basic residues contact the sulfate groups on heparin and HS. Such detailed characterization is still required for other, structurally different BMPs, including BMP6.

METHODS

BMP6 peptides encompassing potential HBDs were analysed on heparin-functionalized plates and microcantilevers, and on membrane HS expressing CHO-K1 cells. Monomeric wild-type BMP6 and mutants were produced, substituting the basic residues with non-charged ones, and their affinity to the heparin-column was measured. The BMP6-heparin interaction was also predicted at atomic level by in silico molecular dynamics.

RESULTS

N-terminal and C-terminal BMP6 peptides showed high heparin affinity in solid-phase assays. The mutation of the two sites (R5L, R6S, R7L and K126N, K127N, R129S) abolished the heparin-binding activity of the recombinant monomeric BMP6. Monomeric BMP6 and peptides specifically bound to membrane HS of CHO-K1 cells through the same domains. Molecular dynamic studies supported the role of the two HBDs, suggesting a cooperative behaviour.

CONCLUSIONS

In BMP6, N-terminal (R5, R6, R7) and C-terminal (K126, K127, R129) domains mediate the interaction with heparin and HS.

GENERAL SIGNIFICANCE

This study provides the molecular mechanism supporting the use of heparin to sequester BMP6 and inhibit hepcidin expression, a novel clinical approach for high-hepcidin iron disorders.

Comprehensive analysis of carbohydrate-protein recognition in the Protein Data Bank

Carbohydrate-protein interactions underpin wide-ranging aspects of biology. However, such interactions remain relatively unexplored in pharmaceutical and biotechnological applications, in part due to the challenges associated with their structural characterisation, both experimentally and computationally. Knowledge-based approaches have shown great success in the prediction of drug-protein and protein-protein interactions, although have not been comprehensively investigated for carbohydrate-protein interactions. In this work, carbohydrate-protein complexes from the Protein Data Bank were comprehensively obtained and analysed to identify patterns in how carbohydrate-protein interactions are mediated.

A click-flipped enzyme substrate boosts the performance of the diagnostic screening for Hunter syndrome

We report on the unexpected finding that click modification of iduronyl azides results in a conformational flip of the pyranose ring, which led to the development of a new strategy for the design of superior enzyme substrates for the diagnostic assaying of iduronate-2-sulfatase (I2S), a lysosomal enzyme related to Hunter syndrome. Synthetic substrates are essential in testing newborns for metabolic disorders to enable early initiation of therapy. Our click-flipped iduronyl triazole showed a remarkably better performance with I2S than commonly used O-iduronates. We found that both O- and triazole-linked substrates are accepted by the enzyme, irrespective of their different conformations, but only the O-linked product inhibits the activity of I2S. Thus, in the long reaction times required for clinical assays, the triazole substrate substantially outperforms the O-iduronate. Applying our click-flipped substrate to assay I2S in dried blood spots sampled from affected patients and random newborns significantly increased the confidence in discriminating between these groups, clearly indicating the potential of the click-flip strategy to control the biomolecular function of carbohydrates.

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.

Degeneracy of the Antithrombin Binding Sequence in Heparin: 2-O-Sulfated Iduronic Acid Can Replace the Critical Glucuronic Acid

Heparin binds to and activates antithrombin (AT) through a specific pentasaccharide sequence, in which a trisaccharide subsite, containing glucuronic acid (GlcA), has been considered as the initiator in the recognition of the polysaccharide by the protein. Recently it was suggested that sulfated iduronic acid (IdoA2S) could replace this "canonical" GlcA. Indeed, a heparin octasaccharidic sequence obtained by chemoenzymatic synthesis, in which GlcA is replaced with IdoA2S, has been found to similarly bind to and activate antithrombin. By using saturation-transfer-difference (STD) NMR, NOEs, transferred NOEs (tr-NOEs) NMR and molecular dynamics, we show that, upon binding to AT, this IdoA2S unit develops comparable interactions with AT as GlcA. Interestingly, two IdoA2S units, both present in a ¹ C₄ -² S₀ equilibrium in the unbound saccharide, shift to full ² S₀ and full ¹ C₄ upon binding to antithrombin, providing the best illustration of the critical role of iduronic acid conformational flexibility in biological systems.

The conformation of the idopyranose ring revisited: How subtle O-substituent induced changes can be deduced from vicinal 1H-NMR coupling constants

The idopyranose ring plays a pivotal role in the conformational, dynamical, and intermolecular binding aspects of glycosaminoglycans like heparin and dermatan sulfate and it was early on assigned a role in the Sugar Code governing biological recognition processes. There is consensus that next to the two canonical ¹C₄ and ⁴C₁ chair conformations, the conformational space accessible to the idopyranose ring entails a ²S_O skew-boat conformation, but the equilibrium between these three ring puckers has evaded satisfactory quantification. In this study a meta-analysis of X-ray solid-state data and vicinal NMR coupling constants is presented, based on the Truncated Fourier Puckering (TFP) formalism and the generalized Karplus (CAGPLUS) equation. This approach yields a model-free, granular and consistent reckoning of 159 idopyranose solution puckering equilibria studied by NMR and allows us to reproduce the involved 636 NMR vicinal couplings with an overall residual RMS(J_obs-J_calc) of 0.184 Hz. Our analyses show that for all ring systems examined, the idopyranosyl chair conformations take up the same ring pucker irrespective of the ring substituent pattern or a vast variety in experimental conditions. Instead, it is the (skew-)boat conformation that adapts to the substitution pattern of the idopyranose ring or a specific sulfation pattern of neighboring saccharides. All idopyranose rings are involved in conformational equilibria that subsume the aforementioned conformers which turn out to differ only a few kJ/mole in conformational energy. Thus, the plasticity and flexibility of idopyranose remains intact under practically all circumstances and, as the glycosidic linkages in heparin are considered to be relatively stiff, the iduronic moiety functions as the linchpin of heparin flexibility thereby being rather a "space(r)" than a "letter" in the alleged Sugar Code alphabet.

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.

Tools for the Quality Control of Pharmaceutical Heparin

Heparin is a vital pharmaceutical anticoagulant drug and remains one of the few naturally sourced pharmaceutical agents used clinically. Heparin possesses a structural order with up to four levels of complexity. These levels are subject to change based on the animal or even tissue sources that they are extracted from, while higher levels are believed to be entirely dynamic and a product of their surrounding environments, including bound proteins and associated cations. In 2008, heparin sources were subject to a major contamination with a deadly compound-an over-sulphated chondroitin sulphate polysaccharide-that resulted in excess of 100 deaths within North America alone. In consideration of this, an arsenal of methods to screen for heparin contamination have been applied, based primarily on the detection of over-sulphated chondroitin sulphate. The targeted nature of these screening methods, for this specific contaminant, may leave contamination by other entities poorly protected against, but novel approaches, including library-based chemometric analysis in concert with a variety of spectroscopic methods, could be of great importance in combating future, potential threats.

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.

Binding and structure-kinetic relationship analysis of selective TLR4-targeted immunosuppressive self-assembling heparin nanoparticles

Self-assembling aliphatic heparin derivatives were shown to inhibit the immune system by antagonizing Toll-like receptor 4/myeloid differentiation protein 2 (TLR4/MD2). In the present study, glycol split heparin-d-erythro-sphingosine conjugates (NAHNP) and its regioselectively desulfated derivatives with shortened aliphatic chains were investigated regarding their biophysical properties in the interaction with TLR4/MD2. Two-dimensional nuclear Overhauser effect spectroscopy studies showed that upon glycol splitting, the heparin backbone gains extra adaptability that facilitates binding to proteins. However, unlike native heparin or glycol split non-anticoagulant heparin (NAH), hydrophobic derivatization of NAH forces sulfated iduronic acid residues to change configuration from a ²S₀ skew-boat to a ¹C₄ chair form. Whereas neither heparin nor NAH had any appreciable effect, NAHNP significantly inhibited lipopolysaccharide-induced activation of the NF-κB transcription factor. We showed that NAHNP binds to TLR4/MD2 with an affinity of 62.3 nM. In line with computational studies, biosensor-based structure-kinetic relationship studies demonstrated that 6-O-sulfo groups of d-glucosamine residue were essential in binding to arginines of both TLR4 and MD2 domains of the receptor complex. The desulfation of 6-O-sulfo groups decreases the association kinetics from 4.2 × 10⁴ M^-1 s^-1 to 3.8 × 10³ M^-1 s^-1, which results in a decreased affinity of 800 nM. Two aliphatic chains of NAHNP bound to the MD2 pocket similarly to lipopolysaccharide. A decrease in chain length resulted in a loss of inhibitory activity on NF-κB transcription and binding affinity to TLR4/MD2. In conclusion, the present study characterizes the immunosuppressive effect of aliphatic heparin derivatives and provides a promising strategy to develop selective immunosuppressants for acute and chronic inflammatory disorders.

Replacement of the L-iduronic acid unit of the anticoagulant pentasaccharide idraparinux by a 6-deoxy-L-talopyranose - Synthesis and conformational analysis

One critical part of the synthesis of heparinoid anticoagulants is the creation of the L-iduronic acid building block featured with unique conformational plasticity which is crucial for the anticoagulant activity. Herein, we studied whether a much more easily synthesizable sugar, the 6-deoxy-L-talose, built in a heparinoid oligosaccharide, could show a similar conformational plasticity, thereby can be a potential substituent of the L-idose. Three pentasaccharides related to the synthetic anticoagulant pentasaccharide idraparinux were prepared, in which the L-iduronate was replaced by a 6-deoxy-L-talopyranoside unit. The talo-configured building block was formed by C4 epimerisation of the commercially available L-rhamnose with high efficacy at both the monosaccharide and the disaccharide level. The detailed conformational analysis of these new derivatives, differing only in their methylation pattern, was performed and the conformationally relevant NMR parameters, such as proton-proton coupling constants and interproton distances were compared to the corresponding ones measured in idraparinux. The lack of anticoagulant activity of these novel heparin analogues could be explained by the biologically not favorable 1C4 chair conformation of their 6-deoxy-L-talopyranoside residues.

Predicting the Structures of Glycans, Glycoproteins, and Their Complexes

Complex carbohydrates are ubiquitous in nature, and together with proteins and nucleic acids they comprise the building blocks of life. But unlike proteins and nucleic acids, carbohydrates form nonlinear polymers, and they are not characterized by robust secondary or tertiary structures but rather by distributions of well-defined conformational states. Their molecular flexibility means that oligosaccharides are often refractory to crystallization, and nuclear magnetic resonance (NMR) spectroscopy augmented by molecular dynamics (MD) simulation is the leading method for their characterization in solution. The biological importance of carbohydrate-protein interactions, in organismal development as well as in disease, places urgency on the creation of innovative experimental and theoretical methods that can predict the specificity of such interactions and quantify their strengths. Additionally, the emerging realization that protein glycosylation impacts protein function and immunogenicity places the ability to define the mechanisms by which glycosylation impacts these features at the forefront of carbohydrate modeling. This review will discuss the relevant theoretical approaches to studying the three-dimensional structures of this fascinating class of molecules and interactions, with reference to the relevant experimental data and techniques that are key for validation of the theoretical predictions.

DSE promotes aggressive glioma cell phenotypes by enhancing HB-EGF/ErbB signaling

Remodeling of the extracellular matrix (ECM) in the tumor microenvironment promotes glioma progression. Chondroitin sulfate (CS) proteoglycans appear in the ECM and on the cell surface, and can be catalyzed by dermatan sulfate epimerase to form chondroitin sulfate/dermatan sulfate (CS/DS) hybrid chains. Dermatan sulfate epimerase 1 (DSE) is overexpressed in many types of cancer, and CS/DS chains mediate several growth factor signals. However, the role of DSE in gliomas has never been explored. In the present study, we determined the expression of DSE in gliomas by consulting a public database and conducting immunohistochemistry on a tissue array. Our investigation revealed that DSE was upregulated in gliomas compared with normal brain tissue. Furthermore, high DSE expression was associated with advanced tumor grade and poor survival. We found high DSE expression in several glioblastoma cell lines, and DSE expression directly mediated DS chain formation in glioblastoma cells. Knockdown of DSE suppressed the proliferation, migration, and invasion of glioblastoma cells. In contrast, overexpression of DSE in GL261 cells enhanced these malignant phenotypes and in vivo tumor growth. Interestingly, we found that DSE selectively regulated heparin-binding EGF-like growth factor (HB-EGF)-induced signaling in glioblastoma cells. Inhibiting epidermal growth factor receptor (EGFR) and ErbB2 with afatinib suppressed DSE-enhanced malignant phenotypes, establishing the critical role of the ErbB pathway in regulating the effects of DSE expression. This evidence indicates that upregulation of DSE in gliomas contributes to malignant behavior in cancer cells. We provide novel insight into the significance of DS chains in ErbB signaling and glioma pathogenesis.

Recognition and Conformational Properties of an Alternative Antithrombin Binding Sequence Obtained by Chemoenzymatic Synthesis

Heparin is a highly sulfated glycosaminoglycan (GAG) of natural origin used as an anticoagulant and antithrombotic drug. These properties are principally based on the binding and activation of antithrombin (AT) through the pentasaccharide sequence GlcNAc/NS,6S-GlcA-GlcNS,3,6S-IdoA2S-GlcNS,6S (AGA*IA). Literature data show that the population of the 2 S0 ring conformation of the 2-O-sulfo-α-l-iduronic acid (IdoA2S) motif correlates with the affinity and activation of AT. It was recently demonstrated that two synthetic AGA*IA-containing hexasaccharides (one G unit added at the reducing end), differing in the degree of sulfation of the IdoA unit, show comparable affinity and ability to activate AT, despite a different conformation of the IdoA residue. In this paper, the binding of these two glycans to AT was studied by isothermal titration microcalorimetry (ITC), transferred (tr-) NOESY, saturation transfer difference (STD) NMR spectroscopy and molecular dynamics (MD) simulations. Results indicated that both the IdoA2S and the IdoA units assume a 2 S0 conformation when bound with AT, and so present a common binding epitope for the two glycans, centred on the AGA*IA sequence.

A combinatorial approach towards the synthesis of non-hydrolysable triazole–iduronic acid hybrid inhibitors of human α-l-iduronidase: discovery of enzyme stabilizers for the potential treatment of MPSI

Synthesis and bioevaluation of substituent-diverse triazole–iduronic acid hybrid molecules are highlighted.

Computational Chemistry Tools in Glycobiology: Modelling of Carbohydrate–Protein Interactions

Molecular modelling provides a major impact in the field of glycosciences, helping in the characterisation of the molecular basis of the recognition between lectins from pathogens and human glycoconjugates, and in the design of glycocompounds with anti-infectious properties. The conformational properties of oligosaccharides are complex, and therefore, the simulation of these properties is a challenging task. Indeed, the development of suitable force fields is required for the proper simulation of important problems in glycobiology, such as the interatomic interactions responsible for oligosaccharide and glycoprotein dynamics, including O-linkages in oligo- and polysaccharides, and N- and O-linkages in glycoproteins. The computational description of representative examples is discussed, herein, related to biologically active oligosaccharides and their interaction with lectins and other proteins, and the new routes open for the design of glycocompounds with promising biological activities.

Insights into Hunter syndrome from the structure of iduronate-2-sulfatase

Hunter syndrome is a rare but devastating childhood disease caused by mutations in the IDS gene encoding iduronate-2-sulfatase, a crucial enzyme in the lysosomal degradation pathway of dermatan sulfate and heparan sulfate. These complex glycosaminoglycans have important roles in cell adhesion, growth, proliferation and repair, and their degradation and recycling in the lysosome is essential for cellular maintenance. A variety of disease-causing mutations have been identified throughout the IDS gene. However, understanding the molecular basis of the disease has been impaired by the lack of structural data. Here, we present the crystal structure of human IDS with a covalently bound sulfate ion in the active site. This structure provides essential insight into multiple mechanisms by which pathogenic mutations interfere with enzyme function, and a compelling explanation for severe Hunter syndrome phenotypes. Understanding the structural consequences of disease-associated mutations will facilitate the identification of patients that may benefit from specific tailored therapies.

Hunter syndrome is a lysosomal storage disease caused by mutations in the enzyme iduronate-2-sulfatase (IDS). Here, the authors present the IDS crystal structure and give mechanistic insights into mutations that cause Hunter syndrome.

Interpenetrating Conducting Hydrogel Materials for Neural Interfacing Electrodes

Conducting hydrogels (CHs) are an emerging technology in the field of medical electrodes and brain-machine interfaces. The greatest challenge to the fabrication of CH electrodes is the hybridization of dissimilar polymers (conductive polymer and hydrogel) to ensure the formation of interpenetrating polymer networks (IPN) required to achieve both soft and electroactive materials. A new hydrogel system is developed that enables tailored placement of covalently immobilized dopant groups within the hydrogel matrix. The role of immobilized dopant in the formation of CH is investigated through covalent linking of sulfonate doping groups to poly(vinyl alcohol) (PVA) macromers. These groups control the electrochemical growth of the conducting polymer poly(3,4-ethylenedioxythiophene) (PEDOT) and subsequent material properties. The effect of dopant density and interdopant spacing on the physical, electrochemical, and mechanical properties of the resultant CHs is examined. Cytocompatible PVA hydrogels with PEDOT penetration throughout the depth of the electrode are produced. Interdopant spacing is found to be the key factor in the formation of IPNs, with smaller interdopant spacing producing CH electrodes with greater charge storage capacity and lower impedance due to increased PEDOT growth throughout the network. This approach facilitates tailorable, high-performance CH electrodes for next generation, low impedance neuroprosthetic devices.

Pyranose ring conformations in mono- and oligosaccharides: a combined MD and DFT approach

A two-step computational protocol is proposed to efficiently study the conformational properties of hexopyranoses with a special emphasis on their ring-inversion-properties. By applying it, the errors resulting from overestimating the contribution of the hydrogen bond-rich, low-energy structures that are not abundant in aqueous solutions are avoided.

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

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

Extension and validation of the GLYCAM force field parameters for modeling glycosaminoglycans

Glycosaminoglycans (GAGs) are an important class of carbohydrates that serve critical roles in blood clotting, tissue repair, cell migration and adhesion, and lubrication. The variable sulfation pattern and iduronate ring conformations in GAGs influence their polymeric structure and nature of interaction. This study characterizes several heparin-like GAG disaccharides and tetrasaccharides using NMR and molecular dynamics simulations to assist in the development of parameters for GAGs within the GLYCAM06 force field. The force field additions include parameters and charges for a transferable sulfate group for O- and N-sulfation, neutral (COOH) forms of iduronic and glucuronic acid, and Δ4,5-unsaturated uronate (ΔUA) residues. ΔUA residues frequently arise from the enzymatic digestion of heparin and heparin sulfate. Simulations of disaccharides containing ΔUA reveal that the presence of sulfation on this residue alters the relative populations of 1H2 and 2H1 ring conformations. Simulations of heparin tetrasaccharides containing N-sulfation in place of N-acetylation on glucosamine residues influence the ring conformations of adjacent iduronate residues.

Pharmacology of Heparin and Related Drugs

Heparin has been recognized as a valuable anticoagulant and antithrombotic for several decades and is still widely used in clinical practice for a variety of indications. The anticoagulant activity of heparin is mainly attributable to the action of a specific pentasaccharide sequence that acts in concert with antithrombin, a plasma coagulation factor inhibitor. This observation has led to the development of synthetic heparin mimetics for clinical use. However, it is increasingly recognized that heparin has many other pharmacological properties, including but not limited to antiviral, anti-inflammatory, and antimetastatic actions. Many of these activities are independent of its anticoagulant activity, although the mechanisms of these other activities are currently less well defined. Nonetheless, heparin is being exploited for clinical uses beyond anticoagulation and developed for a wide range of clinical disorders. This article provides a "state of the art" review of our current understanding of the pharmacology of heparin and related drugs and an overview of the status of development of such drugs.

Uncovering the Relationship between Sulphation Patterns and Conformation of Iduronic Acid in Heparan Sulphate

The L-iduronic acid (IdoA) residue is a critically important structural component in heparan sulphate polysaccharide for the biological functions. The pyranose ring of IdoA is present in (1)C4-chair, (2)SO-skew boat, and less frequently, in (4)C1-chair conformations. Here, we analyzed the conformation of IdoA residue in eight hexasaccharides by NMR. The data demonstrate a correlation between the conformation of IdoA and sulphations in the surrounding saccharide residues. For the 2-O-sulpho IdoA residue, a high degree of sulphation on neighboring residues drives ring dynamics towards the (2)SO-skew boat conformer. In contrast, the nonsulphated IdoA residue is pushed towards the (1)C4-chair conformer when the neighboring residues are highly sulphated. Our data suggest that the conformation of IdoA is regulated by the sulphation pattern of nearby saccharides that is genetically controlled by the heparan sulphate biosynthetic pathway.

Xyloside-primed Chondroitin Sulfate/Dermatan Sulfate from Breast Carcinoma Cells with a Defined Disaccharide Composition Has Cytotoxic Effects in Vitro

We previously reported that the xyloside 2-(6-hydroxynaphthyl) β-d-xylopyranoside (XylNapOH), in contrast to 2-naphthyl β-d-xylopyranoside (XylNap), specifically reduces tumor growth both in vitro and in vivo Although there are indications that this could be mediated by the xyloside-primed glycosaminoglycans (GAGs) and that these differ in composition depending on xyloside and cell type, detailed knowledge regarding a structure-function relationship is lacking. In this study we isolated XylNapOH- and XylNap-primed GAGs from a breast carcinoma cell line, HCC70, and a breast fibroblast cell line, CCD-1095Sk, and demonstrated that both XylNapOH- and XylNap-primed chondroitin sulfate/dermatan sulfate GAGs derived from HCC70 cells had a cytotoxic effect on HCC70 cells and CCD-1095Sk cells. The cytotoxic effect appeared to be mediated by induction of apoptosis and was inhibited in a concentration-dependent manner by the XylNap-primed heparan sulfate GAGs. In contrast, neither the chondroitin sulfate/dermatan sulfate nor the heparan sulfate derived from CCD-1095Sk cells primed on XylNapOH or XylNap had any effect on the growth of HCC70 cells or CCD-105Sk cells. These observations were related to the disaccharide composition of the XylNapOH- and XylNap-primed GAGs, which differed between the two cell lines but was similar when the GAGs were derived from the same cell line. To our knowledge this is the first report on cytotoxic effects mediated by chondroitin sulfate/dermatan sulfate.

Revision of the GROMOS 56A6(CARBO) force field: Improving the description of ring-conformational equilibria in hexopyranose-based carbohydrates chains

This article describes a revised version 56A6(CARBO_R) of the GROMOS 56A6(CARBO) force field for hexopyranose-based carbohydrates. The simulated properties of unfunctionalized hexopyranoses are unaltered with respect to 56A6CARBO . In the context of both O1 -alkylated hexopyranoses and oligosaccharides, the revision stabilizes the regular (4) C1 chair for α-anomers, with the opposite effect for β-anomers. As a result, spurious ring inversions observed in α(1→4)-linked chains when using the original 56A6(CARBO) force field are alleviated. The (4) C1 chair is now the most stable conformation for all d-hexopyranose residues, irrespective of the linkage type and anomery, and of the position of the residue along the chain. The methylation of a d-hexopyranose leads to a systematic shift in the ring-inversion free energy ((4) C1 to (1) C4 ) by 7-8 kJ mol(-1), positive for the α-anomers and negative for the β-anomers, which is qualitatively compatible with the expected enhancement of the anomeric effect upon methylation at O1. The ring-inversion free energies for residues within chains are typically smaller in magnitude compared to those of the monomers, and correlate rather poorly with the latter. This suggests that the crowding of ring substituents upon chain formation alters the ring flexibility in a nonsystematic fashion. In general, the description of carbohydrate chains afforded by 56A6(CARBO_R) suggests a significant extent of ring flexibility, i.e., small but often non-negligible equilibrium populations of inverted chairs, and challenges the "textbook" picture of conformationally locked carbohydrate rings.

Solution Structure of Heparin Pentasaccharide: NMR and DFT Analysis

High-resolution NMR and density functional theory (DFT) calculations have been applied to analysis of heparin pentasaccharide 3D structure in aqueous solution. The fully optimized molecular geometry of two pentasaccharide conformations (differing from each other in the form, one (1)C4 and the other (2)S0, of the sulfated iduronic acid residue) were obtained using the B3LYP/6-311+G(d,p) level of theory in the presence of solvent, the latter included as explicit water molecules. The presented approach enabled insight into variations of the bond lengths, bond angles, and torsion angles, formations of intra- and intermolecular hydrogen bonds, and ionic interactions in the two pentasaccharide conformations. A rather complex hydrogen bond network is formed, including inter-residue and intraresidue bonds between the NH group in the GlcN,3,6S with oxygens linked to C-2 at the IdoA2S residue and the glycosidic O-1 and the neighboring OSO3(-) group linked to C-3 in the same residue. On the other hand, because the first hydration shell is strongly influenced by strong ion-ion and ion-dipole interactions between sodium ions, sulfates, carboxylates, and -OH groups, ionic interactions play an important role in the stabilization of the 3D structure. The DFT-computed three-bond proton-proton coupling constants also showed that best agreement with experiment was obtained with a weighted average of 15:85 ((1)C4/(2)S0) of the sulfated iduronic acid forms indicating that the ratio is even more shifted toward the (2)S0 form than previously supposed. The DFT-computed pentasaccharide conformation differs from the previously published data, with the main changes at the glycosidic linkages, namely, the ψ1 torsion angles and the ϕ3 angle. The comparison of the glycosidic linkage torsion angle values in solution with the antithrombin-pentasaccharide complex also indicates that the pentasaccharide conformation changes upon binding to antithrombin III. The data supports the assumption that the protein selects the more populated (2)S0 conformer of heparin pentasaccharide and, consequently, the binding process of heparin pentasaccharide with antithrombin III is energetically more favorable than formerly expected.

In Silico Analysis of the Metabolic Potential and Niche Specialization of Candidate Phylum "Latescibacteria" (WS3)

The “Latescibacteria” (formerly WS3), member of the Fibrobacteres–Chlorobi–Bacteroidetes (FCB) superphylum, represents a ubiquitous candidate phylum found in terrestrial, aquatic, and marine ecosystems. Recently, single-cell amplified genomes (SAGs) representing the “Latescibacteria” were obtained from the anoxic monimolimnion layers of Sakinaw Lake (British Columbia, Canada), and anoxic sediments of a coastal lagoon (Etoliko lagoon, Western Greece). Here, we present a detailed in-silico analysis of the four SAGs to gain some insights on their metabolic potential and apparent ecological roles. Metabolic reconstruction suggests an anaerobic fermentative mode of metabolism, as well as the capability to degrade multiple polysaccharides and glycoproteins that represent integral components of green (Charophyta and Chlorophyta) and brown (Phaeophycaea) algae cell walls (pectin, alginate, ulvan, fucan, hydroxyproline-rich glycoproteins), storage molecules (starch and trehalose), and extracellular polymeric substances (EPSs). The analyzed SAGs also encode dedicated transporters for the uptake of produced sugars and amino acids/oligopeptides, as well as an extensive machinery for the catabolism of all transported sugars, including the production of a bacterial microcompartment (BMC) to sequester propionaldehyde, a toxic intermediate produced during fucose and rhamnose metabolism. Finally, genes for the formation of gas vesicles, flagella, type IV pili, and oxidative stress response were found, features that could aid in cellular association with algal detritus. Collectively, these results indicate that the analyzed “Latescibacteria” mediate the turnover of multiple complex organic polymers of algal origin that reach deeper anoxic/microoxic habitats in lakes and lagoons. The implications of such process on our understanding of niche specialization in microbial communities mediating organic carbon turnover in stratified water bodies are discussed.

Proteoglycans and their heterogeneous glycosaminoglycans at the atomic scale

Proteoglycan spatiotemporal organization underpins extracellular matrix biology, but atomic scale glimpses of this microarchitecture are obscured by glycosaminoglycan size and complexity. To overcome this, multimicrosecond aqueous simulations of chondroitin and dermatan sulfates were abstracted into a prior coarse-grained model, which was extended to heterogeneous glycosaminoglycans and small leucine-rich proteoglycans. Exploration of relationships between sequence and shape led to hypotheses that proteoglycan size is dependent on glycosaminoglycan unit composition but independent of sequence permutation. Uronic acid conformational equilibria were modulated by adjacent hexosamine sulfonation and iduronic acid increased glycosaminoglycan chain volume and rigidity, while glucuronic acid imparted chain plasticity. Consequently, block copolymeric glycosaminoglycans contained microarchitectures capable of multivalent binding to growth factors and collagen, with potential for interactional synergy at greater chain number. The described atomic scale views of proteoglycans and heterogeneous glycosaminoglycans provide structural routes to understanding their fundamental signaling and mechanical biological roles and development of new biomaterials.