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.

LectomeXplore, an update of UniLectin for the discovery of carbohydrate-binding proteins based on a new lectin classification.

Lectins are non-covalent glycan-binding proteins mediating cellular interactions but their annotation in newly sequenced organisms is lacking. The limited size of functional domains and the low level of sequence similarity challenge usual bioinformatics tools. The identification of lectin domains in proteomes requires the manual curation of sequence alignments based on structural folds. A new lectin classification is proposed. It is built on three levels: (i) 35 lectin domain folds, (ii) 109 classes of lectins sharing at least 20% sequence similarity and (iii) 350 families of lectins sharing at least 70% sequence similarity. This information is compiled in the UniLectin platform that includes the previously described UniLectin3D database of curated lectin 3D structures. Since its first release, UniLectin3D has been updated with 485 additional 3D structures. The database is now complemented by two additional modules: PropLec containing predicted ß-propeller lectins and LectomeXplore including predicted lectins from sequences of the NBCI-nr and UniProt for every curated lectin class. UniLectin is accessible at https://www.unilectin.eu/.

Singular Interface Dynamics of the SARS-CoV-2 Delta Variant Explained with Contact Perturbation Analysis.

Emerging SARS-CoV-2 variants raise concerns about our ability to withstand the Covid-19 pandemic, and therefore, understanding mechanistic differences of those variants is crucial. In this study, we investigate disparities between the SARS-CoV-2 wild type and five variants that emerged in late 2020, focusing on the structure and dynamics of the spike protein interface with the human angiotensin-converting enzyme 2 (ACE2) receptor, by using crystallographic structures and extended analysis of microsecond molecular dynamics simulations. Dihedral angle principal component analysis (PCA) showed the strong similarities in the spike receptor binding domain (RBD) dynamics of the Alpha, Beta, Gamma, and Delta variants, in contrast with those of WT and Epsilon. Dynamical perturbation networks and contact PCA identified the peculiar interface dynamics of the Delta variant, which cannot be directly imputable to its specific L452R and T478K mutations since those residues are not in direct contact with the human ACE2 receptor. Our outcome shows that in the Delta variant the L452R and T478K mutations act synergistically on neighboring residues to provoke drastic changes in the spike/ACE2 interface; thus a singular mechanism of action eventually explains why it dominated over preceding variants.

UniLectin3D, a database of carbohydrate binding proteins with curated information on 3D structures and interacting ligands.

Lectins, and related receptors such as adhesins and toxins, are glycan-binding proteins from all origins that decipher the glycocode, i.e. the structural information encoded in the conformation of complex carbohydrates present on the surface of all cells. Lectins are still poorly classified and annotated, but since their functions are based on ligand recognition, their 3D-structures provide a solid foundation for characterization. UniLectin3D is a curated database that classifies lectins on origin and fold, with cross-links to literature, other databases in glycosciences and functional data such as known specificity. The database provides detailed information on lectins, their bound glycan ligands, and features their interactions using the Protein-Ligand Interaction Profiler (PLIP) server. Special care was devoted to the description of the bound glycan ligands with the use of simple graphical representation and numerical format for cross-linking to other databases in glycoscience. We conceived the design of the database architecture and the navigation tools to account for all organisms, as well as to search for oligosaccharide epitopes complexed within specified binding sites. UniLectin3D is accessible at https://www.unilectin.eu/unilectin3D.

MatrixDB: integration of new data with a focus on glycosaminoglycan interactions.

MatrixDB (http://matrixdb.univ-lyon1.fr/) is an interaction database focused on biomolecular interactions established by extracellular matrix (ECM) proteins and glycosaminoglycans (GAGs). It is an active member of the International Molecular Exchange (IMEx) consortium (https://www.imexconsortium.org/). It has adopted the HUPO Proteomics Standards Initiative standards for annotating and exchanging interaction data, either at the MIMIx (The Minimum Information about a Molecular Interaction eXperiment) or IMEx level. The following items related to GAGs have been added in the updated version of MatrixDB: (i) cross-references of GAG sequences to the GlyTouCan database, (ii) representation of GAG sequences in different formats (IUPAC and GlycoCT) and as SNFG (Symbol Nomenclature For Glycans) images and (iii) the GAG Builder online tool to build 3D models of GAG sequences from GlycoCT codes. The database schema has been improved to represent n-ary experiments. Gene expression data, imported from Expression Atlas (https://www.ebi.ac.uk/gxa/home), quantitative ECM proteomic datasets (http://matrisomeproject.mit.edu/ecm-atlas), and a new visualization tool of the 3D structures of biomolecules, based on the PDB Component Library and LiteMol, have also been added. A new advanced query interface now allows users to mine MatrixDB data using combinations of criteria, in order to build specific interaction networks related to diseases, biological processes, molecular functions or publications.

Mechanism of activation of plant monogalactosyldiacylglycerol synthase 1 (MGD1) by phosphatidylglycerol.

Mono- and digalactosyldiacylglycerol are essential galactolipids for the biogenesis of plastids and functioning of the photosynthetic machinery. In Arabidopsis, the first step of galactolipid synthesis is catalyzed by monogalactosyldiacylglycerol synthase 1 (MGD1), a monotopic protein located in the inner envelope membrane of chloroplasts, which transfers a galactose residue from UDP-galactose to diacylglycerol (DAG). MGD1 needs anionic lipids such as phosphatidylglycerol (PG) to be active, but the mechanism by which PG activates MGD1 is still unknown. Recent studies shed light on the catalytic mechanism of MGD1 and on the possible PG binding site. Particularly, Pro189 was identified as a potential residue implicated in PG binding and His155 as the putative catalytic residue. In the present study, using a multifaceted approach (Langmuir membrane models, atomic force microscopy, molecular dynamics; MD), we investigated the membrane binding properties of native MGD1 and mutants (P189A and H115A). We demonstrated that both residues are involved in PG binding, thus suggesting the existence of a PG-His catalytic dyad that should facilitate deprotonation of the nucleophile hydroxyl group of DAG acceptor. Interestingly, MD simulations showed that MGD1 induces a reorganization of lipids by attracting DAG molecules to create an optimal platform for binding.

Computational tools for drawing, building and displaying carbohydrates: a visual guide.

Drawing and visualisation of molecular structures are some of the most common tasks carried out in structural glycobiology, typically using various software. In this perspective article, we outline developments in the computational tools for the sketching, visualisation and modelling of glycans. The article also provides details on the standard representation of glycans, and glycoconjugates, which helps the communication of structure details within the scientific community. We highlight the comparative analysis of the available tools which could help researchers to perform various tasks related to structure representation and model building of glycans. These tools can be useful for glycobiologists or any researcher looking for a ready to use, simple program for the sketching or building of glycans.

A pipeline to translate glycosaminoglycan sequences into 3D models. Application to the exploration of glycosaminoglycan conformational space.

Mammalian glycosaminoglycans are linear complex polysaccharides comprising heparan sulfate, heparin, dermatan sulfate, chondroitin sulfate, keratan sulfate and hyaluronic acid. They bind to numerous proteins and these interactions mediate their biological activities. GAG-protein interaction data reported in the literature are curated mostly in MatrixDB database (http://matrixdb.univ-lyon1.fr/). However, a standard nomenclature and a machine-readable format of GAGs together with bioinformatics tools for mining these interaction data are lacking. We report here the building of an automated pipeline to (i) standardize the format of GAG sequences interacting with proteins manually curated from the literature, (ii) translate them into the machine-readable GlycoCT format and into SNFG (Symbol Nomenclature For Glycan) images and (iii) convert their sequences into a format processed by a builder generating three-dimensional structures of polysaccharides based on a repertoire of conformations experimentally validated by data extracted from crystallized GAG-protein complexes. We have developed for this purpose a converter (the CT23D converter) to automatically translate the GlycoCT code of a GAG sequence into the input file required to construct a three-dimensional model.

Virtual Screening Against Carbohydrate-Binding Proteins: Evaluation and Application to Bacterial Burkholderia ambifaria Lectin.

Bacterial adhesion to human epithelia via lectins constitutes a therapeutic opportunity to prevent infection. Specifically, BambL (the lectin from Burkholderia ambifaria) is implicated in cystic fibrosis, where lectin-mediated bacterial adhesion to fucosylated lung epithelia is suspected to play an important role. We employed structure-based virtual screening to identify inhibitors of BambL-saccharide interaction with potential therapeutic value. To enable such discovery, a virtual screening protocol was iteratively developed via 194 retrospective screening protocols against 4 bacterial lectins (BambL, BC2L-A, FimH, and LecA) with known ligands. Specific attention was given to the rigorous evaluation of retrospective screening, including calculation of analytical errors for enrichment metrics. The developed virtual screening workflow used crystallographic constraints, pharmacophore filters, and a final manual selection step. The protocol was applied to BambL, predicting 15 active compounds from virtual libraries of approximately 7 million compounds. Experimental validation using fluorescence polarization confirmed micromolar inhibitory activity for two compounds, which were further characterized by isothermal titration calorimetry and surface plasmon resonance. Subsequent testing against LecB from Pseudomonas aeruginosa demonstrated binding specificity of one of the hit compounds. This report demonstrates the utility of virtual screening protocols, integrating ligand-based pharmacophore filtering and structure-based constraints, in the search for bacterial lectin inhibitors.

Glycoscience@Synchrotron: Synchrotron radiation applied to structural glycoscience.

Synchrotron radiation is the most versatile way to explore biological materials in different states: monocrystalline, polycrystalline, solution, colloids and multiscale architectures. Steady improvements in instrumentation have made synchrotrons the most flexible intense X-ray source. The wide range of applications of synchrotron radiation is commensurate with the structural diversity and complexity of the molecules and macromolecules that form the collection of substrates investigated by glycoscience. The present review illustrates how synchrotron-based experiments have contributed to our understanding in the field of structural glycobiology. Structural characterization of protein-carbohydrate interactions of the families of most glycan-interacting proteins (including glycosyl transferases and hydrolases, lectins, antibodies and GAG-binding proteins) are presented. Examples concerned with glycolipids and colloids are also covered as well as some dealing with the structures and multiscale architectures of polysaccharides. Insights into the kinetics of catalytic events observed in the crystalline state are also presented as well as some aspects of structure determination of protein in solution.

Databases of Conformations and NMR Structures of Glycan Determinants.

The present study reports a comprehensive nuclear magnetic resonance (NMR) characterization and a systematic conformational sampling of the conformational preferences of 170 glycan moieties of glycosphingolipids as produced in large-scale quantities by bacterial fermentation. These glycans span across a variety of families including the blood group antigens (A, B and O), core structures (Types 1, 2 and 4), fucosylated oligosaccharides (core and lacto-series), sialylated oligosaccharides (Types 1 and 2), Lewis antigens, GPI-anchors and globosides. A complementary set of about 100 glycan determinants occurring in glycoproteins and glycosaminoglycans has also been structurally characterized using molecular mechanics-based computation. The experimental and computational data generated are organized in two relational databases that can be queried by the user through a user-friendly search engine. The NMR ((1)H and (13)C, COSY, TOCSY, HMQC, HMBC correlation) spectra and 3D structures are available for visualization and download in commonly used structure formats. Emphasis has been given to the use of a common nomenclature for the structural encoding of the carbohydrates and each glycan molecule is described by four different types of representations in order to cope with the different usages in chemistry and biology. These web-based databases were developed with non-proprietary software and are open access for the scientific community available at http://glyco3d.cermav.cnrs.fr.

Three-dimensional representations of complex carbohydrates and polysaccharides--SweetUnityMol: a video game-based computer graphic software.

A molecular visualization program tailored to deal with the range of 3D structures of complex carbohydrates and polysaccharides, either alone or in their interactions with other biomacromolecules, has been developed using advanced technologies elaborated by the video games industry. All the specific structural features displayed by the simplest to the most complex carbohydrate molecules have been considered and can be depicted. This concerns the monosaccharide identification and classification, conformations, location in single or multiple branched chains, depiction of secondary structural elements and the essential constituting elements in very complex structures. Particular attention was given to cope with the accepted nomenclature and pictorial representation used in glycoscience. This achievement provides a continuum between the most popular ways to depict the primary structures of complex carbohydrates to visualizing their 3D structures while giving the users many options to select the most appropriate modes of representations including new features such as those provided by the use of textures to depict some molecular properties. These developments are incorporated in a stand-alone viewer capable of displaying molecular structures, biomacromolecule surfaces and complex interactions of biomacromolecules, with powerful, artistic and illustrative rendering methods. They result in an open source software compatible with multiple platforms, i.e., Windows, MacOS and Linux operating systems, web pages, and producing publication-quality figures. The algorithms and visualization enhancements are demonstrated using a variety of carbohydrate molecules, from glycan determinants to glycoproteins and complex protein-carbohydrate interactions, as well as very complex mega-oligosaccharides and bacterial polysaccharides and multi-stranded polysaccharide architectures.

Negatively charged glyconanoparticles modulate and stabilize the secondary structures of a gp120 V3 loop peptide: toward fully synthetic HIV vaccine candidates.

The third variable region (V3 peptide) of the HIV-1 gp120 is a major immunogenic domain of HIV-1. Controlling the formation of the immunologically active conformation is a crucial step to the rational design of fully synthetic candidate vaccines. Herein, we present the modulation and stabilization of either the α-helix or ß-strand conformation of the V3 peptide by conjugation to negatively charged gold glyconanoparticles (GNPs). The formation of the secondary structure can be triggered by the variation of the buffer concentration and/or pH as indicated by circular dichoism. The peptide on the GNPs shows increased stability toward peptidase degradation as compared to the free peptide. Moreover, only the V3ß-GNPs bind to the anti-V3 human broadly neutralizing mAb 447-52D as demonstrated by surface plasmon resonance (SPR). The strong binding of V3ß-GNPs to the 447-52D mAb was the starting point to address its study as immunogen. V3ß-GNPs elicit antibodies in rabbits that recognize a recombinant gp120 and the serum displayed low but consistent neutralizing activity. These results open up the way for the design of new fully synthetic HIV vaccine candidates.

POLYS 2.0: An open source software package for building three-dimensional structures of polysaccharides.

This article describes an update of POLYS, the POLYSaccharide builder, for generating three-dimensional structures of polysaccharides and complex carbohydrates (Engelsen et al., Biopolymers 1996, 39, 417-433). POLYS is written in portable ANSI C and is now released under an open source license. Using this software, complex branched carbohydrate structures and polysaccharides can be constructed from their primary structure and the relevant monosaccharides stored in database containing information on optimized glycosidic linkage geometries. The constructed three-dimensional structures are described as Cartesian coordinate files which can be used as input to other molecular modeling software. The new version of POLYS includes a large database of monosaccharides and a helical generator to build and optimize regular single helix or double helix structures. To demonstrate the efficiency of POLYS to build carbohydrate structures, four examples of increasing complexity are presented in the manuscript, from simple alpha glucans over complex starch fragments and the double helical structure of amylopectin to the mega-oligosaccharide RhamnoGalacturonan II.

Influence of amylopectin structure and degree of phosphorylation on the molecular composition of potato starch lintners.

Morphology, molecular structure, and thermal properties of potato starch granules with low to high phosphate content were studied as an effect of mild acid hydrolysis (lintnerization) to 80% solubilization at two temperatures (25 and 45°C). Light microscopy showed that the lintners contained apparently intact granules, which disintegrated into fragments upon dehydration. Transmission electron microscopy of rehydrated lintners revealed lacy networks of smaller subunits. The molecular composition of the lintners suggested that they largely consisted of remnants of crystalline lamellae. When lintnerization was performed at 45°C, the lintners contained more of branched dextrins compared to 25°C in both low and intermediate phosphate-containing samples. High-phosphate-containing starch was, however, unaffected by temperature and this was probably due to an altered amylopectin structure rather than the phosphate content. After lintnerization, the melting endotherms were broad with decreased onset and increased peak melting temperatures. The relative crystallinity was lower in lintners prepared at 45°C. A hypothesis that combines the kinetics of lintnerization with the molecular and thermal characteristics of the lintners is presented.

Effect of distal sugars and interglycosidic linkage on the N-glycoprotein linkage region conformation: synthesis and X-ray crystallographic investigation of ß-1-N-alkanamide derivatives of cellobiose and maltose as disaccharide analogs of the conserved chitobiosylasparagine linkage.

The linkage region constituents, 2-deoxy-2-acetamido-ß-D-glucopyranose (GlcNAc) and L-asparagine (Asn) are conserved in the N-glycoproteins of all eukaryotes. Elucidation of the structure and conformation of the linkage region of glycoproteins is important to understand the presentation and dynamics of the carbohydrate chain at the protein/cell surface. Earlier crystallographic studies using monosaccharide models and analogs of N-glycoprotein linkage region have shown that the N-glycosidic torsion, ϕN, is more influenced by the structural variation in the sugar part than that of the aglycon moiety. To access the influence of distal sugar as well as interglycosidic linkage (α or ß) on the N-glycosidic torsion angles, cellobiosyl and maltosyl alkanamides have been synthesized and structural features of seven of these analogs have been characterized by X-ray crystallography. Comparative analysis of the seven disaccharide analogs with the reported monosaccharide analogs showed that the ϕN value of cellobiosyl analogs deviate ~9° with respect to GlcßNHAc. In the case of maltosyl analogs, deviation is more than 18°. These deviations indicate that the N-glycosidic torsion is influenced by addition of distal sugar as well as with respect to inter glycosidic linkage (α or ß); it is less influenced by changes occurring at the aglycon. The χ2 value of alkanamide derived from glucose, cellobiose and maltose exhibit a large range of variations (from 1.6° to -109.9°). This large span of χ2 value suggests the greater degree of rotational freedom around C1'-C2' bond which is restricted in GlcNAc alkanamides. The present finding explicitly proved the importance of molecular architecture in the N-glycoproteins linkage region to maintain the linearity, planarity and rigidity. These factors are necessary for N-glycan to serve role in inter- as well as intramolecular carbohydrate-protein interactions.

Glycosaminoglycan-mimetic infernan grafted with poly(N-isopropylacrylamide): Toward a thermosensitive polysaccharide.

Glycosaminoglycans (GAGs) are essential constituents of the cell surface and extracellular matrix, where they are involved in several cellular processes through their interactions with various proteins. For successful tissue regeneration, developing an appropriate matrix supporting biological activities of cells in a similar manner than GAGs remains still challenging. In this context, this study aims to design a thermosensitive polysaccharide that could further be used as hydrogel for tissue engineering applications. For this purpose, infernan, a marine bacterial exopolysaccharide (EPS) endowed with GAG-mimetic properties was grafted with a thermosensitive polymer, poly(N-isopropylacrylamide) (pNIPAM). Eight grafted polysaccharides were obtained by varying EPS/pNIPAM molar ratio and the molecular weight of pNIPAM. Their physicochemical characteristics and their thermosensitive properties were determined using a multi-technique, experimental approach. In parallel, molecular dynamics and Monte Carlo simulations were applied at two different scales to elucidate, respectively, the molecular conformation of grafted infernan chain and their ability to form an infinite network undergoing a sol-gel transition near the percolation, a necessary condition in hydrogel formation. It comes out from this study that thermosensitive infernan was successfully developed and its potential use in tissue regeneration as a hydrogel scaffold will further be assessed.

Glycosaminoglycans exhibit distinct interactions and signaling with BMP2 according to their nature and localization.

The role of glycosaminoglycans (GAGs) in modulating bone morphogenetic protein (BMP) signaling represents a recent and underexplored area. Conflicting reports suggest a dual effect: some indicate a positive influence, while others demonstrate a negative impact. This duality suggests that the localization of GAGs (either at the cell surface or within the extracellular matrix) or the specific type of GAG may dictate their signaling role. The precise sulfation patterns of heparan sulfate (HS) responsible for BMP2 binding remain elusive. BMP2 exhibits a preference for binding to HS over other GAGs. Using well-characterized biomaterials mimicking the extracellular matrix, our research reveals that HS promotes BMP2 signaling in the extracellular space, contrary to chondroitin sulfate (CS), which enhances BMP2 bioactivity at the cell surface. Further observations indicate that a central IdoA (2S)-GlcNS (6S) tri-sulfated motif within HS hexasaccharides enhances binding. Nevertheless, BMP2 exhibits a degree of adaptability to various HS sulfation types and sequences. Molecular dynamic simulations attribute this adaptability to the BMP2 N-terminal end flexibility. Our findings illustrate the complex interplay between GAGs and BMP signaling, highlighting the importance of localization and specific sulfation patterns. This understanding has implications for the development of biomaterials with tailored properties for therapeutic applications targeting BMP signaling pathways.

A biological guide to glycosaminoglycans: current perspectives and pending questions.

Mammalian glycosaminoglycans (GAGs), except hyaluronan (HA), are sulfated polysaccharides that are covalently attached to core proteins to form proteoglycans (PGs). This article summarizes key biological findings for the most widespread GAGs, namely HA, chondroitin sulfate/dermatan sulfate (CS/DS), keratan sulfate (KS), and heparan sulfate (HS). It focuses on the major processes that remain to be deciphered to get a comprehensive view of the mechanisms mediating GAG biological functions. They include the regulation of GAG biosynthesis and postsynthetic modifications in heparin (HP) and HS, the composition, heterogeneity, and function of the tetrasaccharide linkage region and its role in disease, the functional characterization of the new PGs recently identified by glycoproteomics, the selectivity of interactions mediated by GAG chains, the display of GAG chains and PGs at the cell surface and their impact on the availability and activity of soluble ligands, and on their move through the glycocalyx layer to reach their receptors, the human GAG profile in health and disease, the roles of GAGs and particular PGs (syndecans, decorin, and biglycan) involved in cancer, inflammation, and fibrosis, the possible use of GAGs and PGs as disease biomarkers, and the design of inhibitors targeting GAG biosynthetic enzymes and GAG-protein interactions to develop novel therapeutic approaches.

Influence of ligand presentation density on the molecular recognition of mannose-functionalised glyconanoparticles by bacterial lectin BC2L-A.

Polyvalent carbohydrate-protein interactions play a key role in bio- and pathological processes, including cell-cell communication and pathogen invasion. In order to study, control and manipulate these interactions gold nanoparticles have been employed as a 3D scaffold, presenting carbohydrate ligands in a multivalent fashion for use as high affinity binding partners and a model system for oligosaccharide presentation at biomacromolecular surfaces. In this study, the binding of a series of mannose-functionalised gold nanoparticles to the dimeric BC2L-A lectin from Burkholderia cenocepacia has been evaluated. BC2L-A is known to exhibit a high specificity for (oligo)mannosides. Due to the unique structure and binding nature of this lectin, it provides a useful tool to study (oligo)saccharides presented on multivalent scaffolds. Surface plasmon resonance and isothermal titration calorimetric assays were used to investigate the effect of ligand presentation density towards binding to the bacterial lectin. We show how a combination of structural complementarities between ligand presentation and lectin architecture and statistical re-binding effects are important for increasing the avidity of multivalent ligands for recognition by their protein receptors; further demonstrating the application of glyconanotechnology towards fundamental glycobiology research as well as a potential towards biomedical diagnostics and therapeutic treatments.