Impact of sialic acids on the molecular dynamic of bi-antennary and tri-antennary glycans

Sialic acids cleavage induced by elastin-derived peptides impairs the interaction between insulin and its receptor in adipocytes 3T3-L1

The insulin receptor (IR) plays an important role in insulin signal transduction, the defect of which is believed to be the root cause of type 2 diabetes. In 3T3-L1 adipocytes as in other cell types, the mature IR is a heterotetrameric cell surface glycoprotein composed of two α subunits and two β subunits. Our objective in our study, is to understand how the desialylation of N-glycan chains, induced by elastin-derived peptides, plays a major role in the function of the IR. Using the 3T3-L1 adipocyte line, we show that removal of the sialic acid from N-glycan chains (N893 and N908), induced by the elastin receptor complex (ERC) and elastin derived-peptides (EDPs), leads to a decrease in the autophosphorylation activity of the insulin receptor. We demonstrate by molecular dynamics approaches that the absence of sialic acids on one of these two sites is sufficient to generate local and general modifications of the structure of the IR. Biochemical approaches highlight a decrease in the interaction between insulin and its receptor when ERC sialidase activity is induced by EDPs. Therefore, desialylation by EDPs is synonymous with a decrease of IR sensitivity in adipocytes and could thus be a potential source of insulin resistance associated with diabetic conditions.

Differential MMP-14 targeting by biglycan, decorin, fibromodulin, and lumican unraveled by in silico approach

Small leucine-rich proteoglycans (SLRPs) are major regulators of extracellular matrix assembly and cell signaling. Lumican, a member of the SLRPs family, and its derived peptides were shown to possess antitumor activity by interacting directly with the catalytic domain of MMP-14 leading to the inhibition of its activity. The aim of the present report was to characterize by in silico three-dimensional (3D) modeling the structure and the dynamics of four SLRPs including their core protein and their specific polysaccharide chains to assess their capacity to bind to MMP-14 and to regulate its activity. Molecular docking experiments were performed to identify the specific amino acids of MMP-14 interacting with each of the four SLRPs. The inhibition of each SLRP (100 nM) on MMP-14 activity was measured and the constants of inhibition (K_i) were evaluated. The impact of the number of glycan chains, structures, and dynamics of lumican on the interaction with MMP-14 was assessed by molecular dynamics simulations. Molecular docking analysis showed that all SLRPs bind to MMP-14 through their concave face, but in different regions of the catalytic domain of MMP-14. Each SLRPs inhibited significantly the MMP-14 activity. Finally, molecular dynamics showed the role of glycan chains in interaction with MMP-14 and shielding effect of SLRPs. Altogether, the results demonstrated that each SLRP exhibited inhibition of MMP-14 activity. However, the differential targeting of MMP-14 by the SLRPs was shown to be related not only to the core protein conformation but also to the glycan chain structures and dynamics.

A Multiscale Molecular Dynamic Analysis Reveals the Effect of Sialylation on EGFR Clustering in a CRISPR/Cas9-Derived Model

Bacterial and viral pathogens can modulate the glycosylation of key host proteins to facilitate pathogenesis by using various glycosidases, particularly sialidases. Epidermal growth factor receptor (EGFR) signaling is activated by ligand-induced receptor dimerization and oligomerization. Ligand binding induces conformational changes in EGFR, leading to clusters and aggregation. However, information on the relevance of EGFR clustering in the pattern of glycosylation during bacterial and viral invasion remains unclear. In this study, (1) we established CRISPR/Cas9-mediated GFP knock-in (EGFP-KI) HeLa cells expressing fluorescently tagged EGFR at close to endogenous levels to study EGF-induced EGFR clustering and molecular dynamics; (2) We studied the effect of sialylation on EGF-induced EGFR clustering and localization in live cells using a high content analysis platform and raster image correlation spectroscopy (RICS) coupled with a number and brightness (N&B) analysis; (3) Our data reveal that the removal of cell surface sialic acids by sialidase treatment significantly decreases EGF receptor clustering with reduced fluorescence intensity, number, and area of EGFR-GFP clusters per cell upon EGF stimulation. Sialylation appears to mediate EGF-induced EGFR clustering as demonstrated by the change of EGFR-GFP clusters in the diffusion coefficient and molecular brightness, providing new insights into the role of sialylation in EGF-induced EGFR activation; and (4) We envision that the combination of CRISPR/Cas9-mediated fluorescent tagging of endogenous proteins and fluorescence imaging techniques can be the method of choice for studying the molecular dynamics and interactions of proteins in live cells.

Haptoglobin polymorphism affects its N-glycosylation pattern in serum

•

Polymorphism affects glycosylation pattern of haptoglobin in healthy population.

•

Sample phenotype classification was done based on the number and type of α-chains.

•

Glycoproteomic analyses of haptoglobin were done using enzyme-assisted LC-MS/MS.

•

Significant differences were obseerved in branching, sialylation and fucosylation.

Introduction

Haptoglobin (Hp) is an abundant acute-phase protein secreted mainly by the liver into the bloodstream. There are three Hp protein phenotypes (Hp type 1–1, 2–1, and 2–2), which differ in the number of α- and β-chains, type of α-chain (the β-chain type remains the same in all the Hp phenotypes), and the polymers that they form via disulfide bonds. Hp has four N-glycosylation sites on the β-chain. Glycosylation of Hp has been reported frequently as a potential glycobiomarker for many diseases; however, whether Hp polymorphism affects its glycosylation has not yet been addressed extensively or in depth.

Objectives

This study investigated the differences between the glycosylation patterns of Hp phenotypes using serum from 12 healthy individuals (four for each Hp phenotype).

Method

An efficient method for isolating Hp from serum was established and subsequently the Hp phenotype of each sample was characterized by immunoblotting. Then, LC-MS/MS analysis of isolated Hp after treatment with three exoglycosidases (sialidase, α2-3 neuraminidase, Endo F3) was performed to characterize the glycosylation pattern of Hp for each individual sample.

Results

The data reveal significant differences among the branching, sialylation, and fucosylation of Hp types, documenting the effect of Hp polymorphism on its glycosylation.

Conclusion

Overall, the study suggests that Hp phenotype characterization should be considered during the investigation of Hp glycosylation.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

How molecular modelling can better broaden the understanding of glycosylations

Glycosylations are among the most ubiquitous post-translational modifications (PTMs) in proteins, and the effects of their perturbations are seen in various diseases such as cancers, diabetes and arthritis to name a few. Yet they remain one of the most enigmatic aspects of protein structure and function. On the other hand, molecular modelling techniques have been rapidly bridging this knowledge gap since the last decade. In this review, we discuss how these techniques have proven to be indispensable for a better understanding of the role of glycosylations in glycoprotein structure and function.

Structure-guided mutagenesis of a mucin-selective metalloprotease from Akkermansia muciniphila alters substrate preferences

Akkermansia muciniphila, a mucin-degrading microbe found in the human gut, is often associated with positive health outcomes. The abundance of A. muciniphila is modulated by the presence and accessibility of nutrients, which can be derived from diet or host glycoproteins. In particular, the ability to degrade host mucins, a class of proteins carrying densely O-glycosylated domains, provides a competitive advantage in the sustained colonization of niche mucosal environments. Although A. muciniphila is known to rely on mucins as a carbon and nitrogen source, the enzymatic machinery used by this microbe to process mucins in the gut is not yet fully characterized. Here, we focus on the mucin-selective metalloprotease, Amuc_0627 (AM0627), which is known to cleave between adjacent residues carrying truncated core 1 O-glycans. We showed that this enzyme is capable of degrading purified mucin 2 (MUC2), the major protein component of mucus in the gut. An X-ray crystal structure of AM0627 (1.9 Å resolution) revealed O-glycan–binding residues that are conserved between structurally characterized enzymes from the same family. We further rationalized the substrate cleavage motif using molecular modeling to identify nonconserved glycan-interacting residues. We conclude that mutagenesis of these residues resulted in altered substrate preferences down to the glycan level, providing insight into the structural determinants of O-glycan recognition.

Conformational preferences of triantennary and tetraantennary hybrid N-glycans in aqueous solution: Insights from 20 μs long atomistic molecular dynamic simulations

In the current study, we have investigated the conformational dynamics of a triantennary (N-glycan1) and tetraantennary (N-glycan2) hybrid N-glycans found on the surface of the HIV glycoprotein using 20 μs long all-atom molecular dynamics (MD) simulations. The main objective of the present study is to elucidate the influence of adding a complex branch on the overall glycan structural dynamics. Our investigation suggests that the average RMSD value increases when a complex branch is added to N-glycan1. However, the RMSD distribution is relatively wider in the case of N-glycan1 compared to N-glycan2, which indicates that multiple complex branches restrict the conformational variability of glycans. A similar observation is obtained from the principal component analysis of both glycans. All the puckering states (⁴C₁ to ¹C₄) of each monosaccharide except mannose are sampled in our simulations, although the ⁴C₁ chair form is energetically more favorable than ¹C₄. In N-glycan1, the 1-6 linkage in the mannose branch [Man(9)-α(1-6)-Man(5)] stays in the gauche-gauche cluster, whereas it moves towards trans-gauche in N-glycan2. For both glycans, mannose branches are more flexible than the complex branches, and adding a complex branch does not influence the dynamics of the mannose branches. We have noticed that the end-to-end distance of the complex branch shortens by ∼ 10 Å in the presence of another complex branch. This suggests that in the presence of an additional complex branch, the other complex branch adopts a close folded structure. All these conformational changes involve the selective formation of inter-residue and water-mediated hydrogen-bond networks.

Native Ion Mobility-Mass Spectrometry-Enabled Fast Structural Interrogation of Labile Protein Surface Modifications at the Intact Protein Level

Protein sialylation has been closely linked to many diseases including Alzheimer's disease (AD). It is also broadly implicated in therapeutics operating in a pattern-dependent (e.g., Neu5Ac vs Neu5Gc) manner. However, how the sialylation pattern affects the AD-associated, transferrin-assisted iron/Aβ cellular uptake process remains largely ill-defined. Herein, we report the use of native ion mobility-mass spectrometry (IM-MS)-based fast structural probing methodology, enabling well-controlled, synergistic, and in situ manipulation of mature glycoproteins and attached sialic acids. IM-MS-centered experiments enable the combinatorial interrogation of sialylation effects on Aβ cytotoxicity and the chemical, conformational, and topological stabilities of transferrin. Cell viability experiments suggest that Neu5Gc replacement enhances the transferrin-assisted, iron loading-associated Aβ cytotoxicity. Native gel electrophoresis and IM-MS reveal that sialylation stabilizes transferrin conformation but inhibits its dimerization. Collectively, IM-MS is adapted to capture key sialylation intermediates involved in fine-tuning AD-associated glycoprotein structural microheterogeneity. Our results provide the molecular basis for the importance of sustaining moderate TF sialylation levels, especially Neu5Ac, in promoting iron cellular transportation and rescuing iron-enhanced Aβ cytotoxicity.

Effects of changes in glycan composition on glycoprotein dynamics: example of N-glycans on insulin receptor

Glycosylation is among the most common post-translational modifications in proteins, although it is observed in only about 10% of all the protein structures in protein data bank (PDB). Modifications of sugar composition in glycoproteins profoundly impact the overall physiology of the organism. One such example is the development of insulin resistance, which has been attributed to the removal of sialic acid residues from N-glycans of insulin receptor (IR) from various experimental studies. How such modifications affect the glycan-glycoprotein dynamics, and ultimately their function is not clearly understood to date. In this study, we performed molecular dynamics simulations of glycans in different environments. We studied the effects of removal of sialic acid on the glycan, as well as on the dynamics of leucine-rich repeat L1 domain of the IR ectodomain. We observed perturbations in L1 domain dynamics as a result of the removal of sialic acid. The perturbations include an increase in the flexibility of insulin-binding residues, which may affect insulin binding with IR. These changes are accompanied by perturbations in glycan-protein interactions and perturbation of long-range allosteric dynamics. Our observations will further aid in understanding the role of sugars in maintaining homeostasis and how changes in glycan composition may lead to perturbations in homeostasis, ultimately leading to conditions such as insulin resistance.

Elastic fibers: formation, function, and fate during aging and disease

Elastic fibers are extracellular components of higher vertebrates and confer elasticity and resilience to numerous tissues and organs such as large blood vessels, lungs, and skin. Their formation and maturation take place in a complex multistage process called elastogenesis. It requires interactions between very different proteins but also other molecules and leads to the deposition and crosslinking of elastin's precursor on a scaffold of fibrillin-rich microfibrils. Mature fibers are exceptionally resistant to most influences and, under healthy conditions, retain their biomechanical function over the life of the organism. However, due to their longevity, they accumulate damages during aging. These are caused by proteolytic degradation, formation of advanced glycation end products, calcification, oxidative damage, aspartic acid racemization, lipid accumulation, carbamylation, and mechanical fatigue. The resulting changes can lead to diminution or complete loss of elastic fiber function and ultimately affect morbidity and mortality. Particularly, the production of elastokines has been clearly shown to influence several life-threatening diseases. Moreover, the structure, distribution, and abundance of elastic fibers are directly or indirectly influenced by a variety of inherited pathological conditions, which mainly affect organs and tissues such as skin, lungs, or the cardiovascular system. A distinction can be made between microfibril-related inherited diseases that are the result of mutations in diverse microfibril genes and indirectly affect elastogenesis, and elastinopathies that are linked to changes in the elastin gene. This review gives an overview on the formation, structure, and function of elastic fibers and their fate over the human lifespan in health and disease.

Identification of Isomeric N-Glycans by Conformer Distribution Fingerprinting using Ion Mobility Mass Spectrometry

Glycans possess unparalleled structural complexity arising from chemically similar monosaccharide building blocks, configurations of anomeric linkages and different branching patterns, potentially giving rise to many isomers. This level of complexity is one of the main reasons that identification of exact glycan structures in biological samples still lags behind that of other biomolecules. Here, we introduce a methodology to identify isomeric N-glycans by determining gas phase conformer distributions (CDs) by measuring arrival time distributions (ATDs) using drift-tube ion mobility spectrometry-mass spectrometry. Key to the approach is the use of a range of well-defined synthetic glycans that made it possible to investigate conformer distributions in the gas phase of isomeric glycans in a systematic manner. In addition, we have computed CD fingerprints by molecular dynamics (MD) simulation, which compared well with experimentally determined CDs. It supports that ATDs resemble conformational populations in the gas phase and offer the prospect that such an approach can contribute to generating a library of CCS distributions (CCSDs) for structure identification.

Protein glycosylation in head and neck cancers: From diagnosis to treatment

Glycosylation is the most common post-translational modification (PTM) of proteins. Malignant tumour cells frequently undergo an alteration in surface protein glycosylation. This phenomenon is also common in cancers of the head and neck, most of which are squamous cell carcinomas (HNSCC). It affects cell functions, including proliferation, motility and invasiveness, thus increasing the propensity to metastasise. HNSCC represents the sixth most frequent malignancy worldwide. These neoplasms, which arise from the mucous membranes of the various anatomical subsites of the upper aero-digestive tract, are heterogeneous in terms of aetiology and clinico-pathologic features. With current treatments, only about 50% of HNSCC patients survive beyond 5-years. Therefore, there is the pressing need to dissect NHSCC heterogeneity to inform treatment choices. In particular, reliable biomarkers of predictive and prognostic value are eagerly needed. This review describes the current state of the art and bio-pathological meaning of glycosylation signatures associated with HNSCC and explores the possible role of tumour specific glycoproteins as potential biomarkers and attractive therapeutic targets. We have also compiled data relating to altered glycosylation and the nature of glycoproteins as tools for the identification of circulating tumour cells (CTCs) in the new era of liquid biopsy.

Umbrella Visualization: A method of analysis dedicated to glycan flexibility with UnityMol

N-glycosylation is a post-translational modification heavily impacting protein functions. Some alterations of glycosylation, such as sialic acid hydrolysis, are related to protein dysfunction. Because of their high flexibility and the many reactive groups of the glycan chains, studying glycans with in vitro methods is a challenging task. Molecular dynamics is a useful tool and probably the only one in biology able to overcome this problem and gives access to conformational information through exhaustive sampling. To better decipher the impact of N-glycans, the analysis and visualization of their influence over time on protein structure is a prerequisite. We developed the Umbrella Visualization, a graphical method that assigns the glycan intrinsic flexibility during a molecular dynamics trajectory. The density plot generated by this method brought relevant informations regarding glycans dynamics and flexibility, but needs further development in order to integrate an accurate description of the protein topology and its interactions. We propose here to transform this analysis method into a visualization mode in UnityMol. UnityMol is a molecular editor, viewer and prototyping platform, coded in C#. The new representation of glycan chains presented in this study takes into account both the main positions adopted by each antenna of a glycan and their statistical relevance. By displaying the collected data on the protein surface, one is then able to investigate the protein/glycan interactions.

Quantification of the α2-6 Sialic Acid Linkage in Branched N-Glycan Structures with Capillary Nanogel Electrophoresis

Sialylation and sialic acid linkage in N-glycans are markers of disease but are analytically challenging to quantify. A capillary electrophoresis method is reported that integrates a unique combination of enzymes and lectins to modify sialylated N-glycans in real time in the capillary so that N-glycan structures containing α2-6-linked sialic acid are easily separated, detected, and quantified. In this study, N-glycans were sequentially cleaved by enzymes at the head of the separation capillary so that the presence of α2-6-linked sialic acids corresponded to a shift in the analyte migration time in a manner that enabled interpretation of the N-glycan structure. Following injection, only afucosylated N-glycan structures were passed through enzyme zones that contained α2-3 sialidase, followed by β1-3,4 galactosidase, which cleaved any terminal α2-3-linked sialic acid and underlying galactose yielding a terminal N-acetyl glucosamine. With this treatment complete, a third zone of α2-3,6,8 sialidase converted the remaining α2-6-linked sialic acid to terminal galactose. With these enzyme processing steps the α2-6-linked sialic acid residues on an N-glycan correlated directly to the number of terminal galactose residues that remained. The number of terminal galactose residues could be interpreted as a stepwise decrease in the migration time. Complex N-glycans from α-1-acid glycoprotein were analyzed using this approach, revealing that a limited number of α2-6-linked sialic acids were present with biantennary, triantennary, and tetraantennary N-glycans of α-1-acid glycoprotein generally containing 0 or 1 α2-6-linked sialic acid.

Elastic fibers and elastin receptor complex: Neuraminidase-1 takes the center stage

Extracellular matrix (ECM) has for a long time being considered as a simple architectural support for cells. It is now clear that ECM presents a fundamental influence on cells driving their phenotype and fate. This complex network is highly specialized and the different classes of macromolecules that comprise the ECM determine its biological functions. For instance, collagens are responsible for the tensile strength of tissues, proteoglycans and glycosaminoglycans are essential for hydration and resistance to compression, and glycoproteins such as laminins facilitate cell attachment. The largest structures of the ECM are the elastic fibers found in abundance in tissues suffering high mechanical constraints such as skin, lungs or arteries. These structures present a very complex composition whose core is composed of elastin surrounded by a microfibrils mantle. Elastogenesis is a tightly regulated process involving the sialidase activity of the Neuraminidase-1 (Neu-1) sub-unit of the Elastin Receptor Complex. Interestingly, Neu-1 subunit also serves as a sensor of elastin degradation via its ability to transmit elastin-derived peptides signaling. Finally, reports showing that neuraminidase activity is able to regulate TGF-β activation raises questions about a possible role for Neu-1 in elastic fibers remodeling. In this mini review, we develop the concept of the regulation of the whole life of elastic fibers through an original scope, the key role of Neu-1 sialidase enzymatic activity.

Sialic acid as a target for the development of novel antiangiogenic strategies

Sialic acid is associated with glycoproteins and gangliosides of eukaryotic cells. It regulates various molecular interactions, being implicated in inflammation and cancer, where its expression is regulated by sialyltransferases and sialidases. Angiogenesis, the formation of new capillaries, takes place during inflammation and cancer, and represents the outcome of several interactions occurring at the endothelial surface among angiogenic growth factors, inhibitors, receptors, gangliosides and cell-adhesion molecules. Here, we elaborate on the evidences that many structures involved in angiogenesis are sialylated and that their interactions depend on sialic acid with implications in angiogenesis itself, inflammation and cancer. We also discuss the possibility to exploit sialic acid as a target for the development of novel antiangiogenic drugs.

Selective Inhibitors of Human Neuraminidase 1 (NEU1)

Inhibitors of human neuraminidase enzymes (NEU) are recognized as important tools for the study of the biological functions of NEU and will be potent tools for elucidating the role of these enzymes in regulating the repertoire of cellular glycans. Here we report the discovery of selective inhibitors of the human neuraminidase 1 (NEU1) and neuraminidase 2 (NEU2) enzymes with exceptional potency. A library of modified 2-deoxy-2,3-didehydro- N-acetylneuraminic acid (DANA) analogues, with variability in the C5- or C9-position, were synthesized and evaluated against four human neuraminidase isoenyzmes (NEU1-4). Hydrophobic groups with an amide linker at the C5 and C9 positions were well accommodated by NEU1, and a hexanamido group was found to give the best potency at both positions. While the C5-hexanamido-C9-hexanamido-DANA analogue did not show synergistic improvements for combined modification, an extended alkylamide at an individual position combined with a smaller group at the second gave increased potency. The best NEU1 inhibitor identified was a C5-hexanamido-C9-acetamido-DANA that had a K_i of 53 ± 5 nM and 340-fold selectivity over other isoenzymes. Additionally, we demonstrated that C5-modifications combined with a C4-guandino group provided the most potent NEU2 inhibitor reported, with a K_i of 1.3 ± 0.2 μM and 7-fold selectivity over other NEU isoenzymes.

Production of Elastin-Derived Peptides Contributes to the Development of Nonalcoholic Steatohepatitis

Affecting more than 30% of the Western population, nonalcoholic fatty liver disease (NAFLD) is the most common liver disease and can lead to multiple complications, including nonalcoholic steatohepatitis (NASH), cancer, hypertension, and atherosclerosis. Insulin resistance and obesity are described as potential causes of NAFLD. However, we surmised that factors such as extracellular matrix remodeling of large blood vessels, skin, or lungs may also participate in the progression of liver diseases. We studied the effects of elastin-derived peptides (EDPs), biomarkers of aging, on NAFLD progression. We evaluated the consequences of EDP accumulation in mice and of elastin receptor complex (ERC) activation on lipid storage in hepatocytes, inflammation, and fibrosis development. The accumulation of EDPs induces hepatic lipogenesis (i.e., SREBP1c and ACC), inflammation (i.e., Kupffer cells, IL-1β, and TGF-β), and fibrosis (collagen and elastin expression). These effects are induced by inhibition of the LKB1-AMPK pathway by ERC activation. In addition, pharmacological inhibitors of EDPs demonstrate that this EDP-driven lipogenesis and fibrosis relies on engagement of the ERC. Our data reveal a major role of EDPs in the development of NASH, and they provide new clues for understanding the relationship between NAFLD and vascular aging.