Predicting the Structures of Glycans, Glycoproteins, and Their Complexes

Nonconventional NMR Spin-Coupling Constants in Oligosaccharide Conformational Modeling: Structural Dependencies Determined from Density Functional Theory Calculations

Nonconventional NMR spin-coupling constants were investigated to determine their potential as conformational constraints in MA'AT modeling of the O-glycosidic linkages of oligosaccharides. Four (¹ J _C1',H1', ¹ J _C1',C2', ² J _C1',H2', and ² J _C2',H1') and eight (¹ J _C4,H4, ¹ J _C3,C4, ¹ J _C4,C5, ² J _C3,H4, ² J _C4,H3, ² J _C5,H4, ² J _C4,H5, and ² J _C3,C5) spin-couplings in methyl β-d-galactopyranosyl-(1→4)-β-d-glucopyranoside (methyl β-lactoside) were calculated using density functional theory (DFT) to determine their dependencies on O-glycosidic linkage C-O torsion angles, ϕ and ψ, respectively. Long-range ⁴ J _H1',H4 was also examined as a potential conformational constraint of either ϕ or ψ. Secondary effects of exocyclic (hydroxyl) C-O bond rotation within or proximal to these coupling pathways were investigated. Based on the findings of methyl β-lactoside, analogous J-couplings were studied in five additional two-bond O-glycosidic linkages [βGlcNAc-(1→4)-βMan, 2-deoxy-βGlc-(1→4)-βGlc, αMan-(1→3)-βMan, αMan-(1→2)-αMan, and βGlcNAc(1→2)-αMan] to determine whether the coupling behaviors observed in methyl β-lactoside were more broadly observed. Of the 13 nonconventional J-couplings studied, 7 exhibit properties that may be useful in future MA'AT modeling of O-glycosidic linkages, none of which involve coupling pathways that include the linkage C-O bonds. The findings also provide new insights into the general effects of exocyclic C-O bond conformation on the magnitude of experimental spin-couplings in saccharides and other hydroxyl-containing molecules.

Molecular simulations of complex carbohydrates and glycoconjugates

Complex carbohydrates (glycans) are the most abundant and versatile biopolymers in nature. The broad diversity of biochemical functions that carbohydrates cover is a direct consequence of the variety of 3D architectures they can adopt, displaying branched or linear arrangements, widely ranging in sizes, and with the highest diversity of building blocks of any other natural biopolymer. Despite this unparalleled complexity, a common denominator can be found in the glycans' inherent flexibility, which hinders experimental characterization, but that can be addressed by high-performance computing (HPC)-based molecular simulations. In this short review, I present and discuss the state-of-the-art of molecular simulations of complex carbohydrates and glycoconjugates, highlighting methodological strengths and weaknesses, important insights through emblematic case studies, and suggesting perspectives for future developments.

Glycoprotein molecular dynamics analysis: SARS-CoV-2 spike glycoprotein case study

Molecular Dynamics (MD) is a method used to calculate the movement of atoms and molecules broadly applied to several aspects of science. It involves computational simulation, which makes it, at first glance, not easily accessible. The rise of several automated tools to perform molecular simulations has allowed researchers to navigate through the various steps of MD. This enables to elucidate structural properties of proteins that could not be analyzed otherwise, such as the impact of glycosylation. Glycosylation dictates the physicochemical and biological properties of a protein modulating its solubility, stability, resistance to proteolysis, interaction partners, enzymatic activity, binding and recognition. Given the high conformational and compositional diversity of the glycan chains, assessing their influence on the protein structure is challenging using conventional analytical techniques. In this manuscript, we present a step-by-step workflow to build and perform MD analysis of glycoproteins focusing on the SPIKE glycoprotein of SARS-CoV-2 to appraise the impact of glycans in structure stabilization and antibody occlusion.

Experimental and computational characterization of dynamic biomolecular interaction systems involving glycolipid glycans

On cell surfaces, carbohydrate chains that modify proteins and lipids mediate various biological functions, which are exerted not only through carbohydrate-protein interactions but also through carbohydrate-carbohydrate interactions. These glycans exhibit considerable degrees of conformational variability and often form clusters providing multiple binding sites. The integration of nuclear magnetic resonance spectroscopy and molecular dynamics simulation has made it possible to delineate the dynamical structures of carbohydrate chains. This approach has facilitated the remodeling of oligosaccharide conformational space in the prebound state to improve protein-binding affinity and has been applied to visualize dynamic carbohydrate-carbohydrate interactions that control glycoprotein-glycoprotein complex formation. Functional glycoclusters have been characterized by experimental and computational approaches applied to various model membranes and artificial self-assembling systems. This line of investigation has provided dynamic views of molecular assembling on glycoclusters, giving mechanistic insights into physiological and pathological molecular events on cell surfaces as well as clues for the design and creation of molecular systems exerting improved glycofunctions. Further development and accumulation of such studies will allow detailed understanding and artificial control of the "glycosynapse" foreseen by Dr. Sen-itiroh Hakomori.

Carbohydrate anchored lipid nanoparticles

Nanotechnology has been a dynamic field for formulation scientists with multidisciplinary research being conducted worldwide. Advancements in development of functional nanosystems have led to evolution of breakthrough technologies. Lipidic nanosystems, in particular, are highly preferred owing to their non-immunogenic safety profiles along with a range of versatile intrinsic properties. Surface modification of lipid nanoparticles by anchoring carbohydrates to these systems is one such attractive drug delivery technology. Carbohydrates confer interesting properties to the nanosystems such as stealth, biostability, bioavailability, reduced toxicity due to decreased immunogenic response, targeting potential as well as ease of commercial availability. The carbohydrate anchored systems can be developed using methods such as adsorption, incorporation (nanoprecipitation or solvent displacement method), crosslinking and grafting. Current review provides a detailed overview of potential lipid based nanoparticulate systems with an emphasis on liposomes, solid lipid nanoparticles, nanostructures lipid carriers and micelles. Review further explores basics of surface modification, methods applied therein, advantages of carbohydrates as surface modifiers, their versatile applications, techniques for characterization of carbohydrate anchored systems and vital regulatory aspects concerned with these specialized systems.

Carbohydrate structure database (CSDB) oligosaccharide conformation tool

Population analysis in terms of glycosidic torsion angles is frequently used to reveal preferred conformers of glycans. However, due to high structural diversity and flexibility of carbohydrates, conformational characterization of complex glycans can be a challenging task. Herein we present a conformation module of oligosaccharide fragments occurring in natural glycan structures developed on the platform of the Carbohydrate Structure Database (CSDB). Currently, this module deposits free energy surface and conformer abundance maps plotted as a function of glycosidic torsions for 194 inter-residue bonds. Data are automatically and continuously derived from explicit-solvent molecular dynamics (MD) simulations. The module was also supplemented with high-temperature MD data of saccharides (2403 maps) provided by GlycoMapsDB (hosted by GLYCOSCIENCES.de project). Conformational data defined by up to four torsional degrees of freedom can be freely explored using a web interface of the module available at http://csdb.glycoscience.ru/database/core/search_conf.html.

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.

Comparative Molecular Modelling of Capsular Polysaccharide Conformations in Streptococcus suis Serotypes 1, 2, 1/2 and 14 Identifies Common Epitopes for Antibody Binding

Streptococcus suis is an encapsulated, commensal, potentially pathogenic bacterium that infects swine globally and causes sporadic life-threatening zoonotic septicemia and meningitis infections in humans. The capsular polysaccharide is a primary virulence factor for S. suis. As S. suis serotype 2 is the most prevalent serotype globally, the serotype 2 CPS is the primary target of current efforts to develop an effective glycoconjugate veterinary vaccine against S. suis. Possible cross-protection with related serotypes would broaden the coverage of a vaccine. The CPS in serotypes 2 and 1/2 differ at a single residue (Gal versus GalNAc), and both are similar to serotypes 1 and 14: all contain a terminal sialic acid on a side chain. However, despite this similarity, there is complex pattern of cross-protection for these serotypes, with varying estimations of the importance of sialic acid in a protective epitope. Further, a pentasaccharide without the terminal sialic acid has been identified as minimal epitope for serotype 2. Here we use molecular simulation to model the molecule conformations of the CPS in serotypes 2, 1/2, 1 and 14, as well as three vaccine candidate oligosaccharides. The common epitopes we identify assist in rationalizing the apparently contradictory immunological data and provide a basis for rational design of S. suis vaccines in the future.

Host-Dependent Producibility of Recombinant Cypridina noctiluca Luciferase With Glycosylation Defects

Cypridina noctiluca luciferase (CLuc) is a secreted luminescent protein that reacts with its substrate (Cypridina luciferin) to emit light. CLuc is known to be a thermostable protein and has been used for various research applications, including in vivo imaging and high-throughput reporter assays. Previously, we produced a large amount of recombinant CLuc for crystallographic analysis. However, this recombinant protein did not crystallize, probably due to heterogeneous N-glycan modifications. In this study, we produced recombinant CLuc without glycan modifications by introducing mutations at the N-glycan modification residues using mammalian Expi293F cells, silkworms, and tobacco Bright Yellow-2 cells. Interestingly, recombinant CLuc production depended heavily on the expression hosts. Among these selected hosts, we found that Expi293F cells efficiently produced the recombinant mutant CLuc without significant effects on its luciferase activity. We confirmed the lack of N-glycan modifications for this mutant protein by mass spectrometry analysis but found slight O-glycan modifications that we estimated were about 2% of the ion chromatogram peak area for the detected peptide fragments. Moreover, by using CLuc deletion mutants during the investigation of O-glycan modifications, we identified amino acid residues important to the luciferase activity of CLuc. Our results provide invaluable information related to CLuc function and pave the way for its crystallographic analysis.

Computational Prediction of N- and O-Linked Glycosylation Sites for Human and Mouse Proteins

Protein glycosylation is one of the most complex posttranslational modifications (PTM) that play a fundamental role in protein function. Identification and annotation of these sites using experimental approaches are challenging and time consuming. Hence, there is a demand to build fast and efficient computational methods to address this problem. Here, we present the SPRINT-Gly framework containing the largest dataset and a prediction model of glycosylation sites for a given protein sequence. In this framework, we construct a large dataset containing N- and O-linked glycosylation sites of human and mouse proteins, collected from different sources. We then introduce the SPRINT-Gly method to predict putative N- and O-linked sites. SPRINT-Gly is a machine learning-based approach consisting of a number of trained predictive models for glycosylation sites in both human and mouse proteins, separately. The method is built by incorporating sequence-based, predicted structural, and physicochemical information of the neighboring residues of each N- and O-linked glycosylation site and by training deep learning neural network and support vector machine as classifiers. SPRINT-Gly outperformed other existing methods by achieving 18% and 50% higher Matthew's correlation coefficient for N- and O-linked glycosylation site prediction, respectively. SPRINT-Gly is publicly available as an online and stand-alone predictor at https://sparks-lab.org/server/sprint-gly/ .

Pyranose Ring Puckering Thermodynamics for Glycan Monosaccharides Associated with Vertebrate Proteins

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

Construction of InstantPC derivatized glycan GU database: A foundation work for high-throughput and high-sensitivity glycomic analysis

Glycosylation is well-recognized as a critical quality attribute of biotherapeutics being routinely monitored to ensure desired product quality, safety, and efficacy. Additionally, as one of the most prominent and complex post-translational modifications, glycosylation plays a key role in disease manifestation. Changes in glycosylation may serve as a specific and sensitive biomarker for disease diagnostics and prognostics. However, the conventional 2-aminobenzamide based N-glycosylation analysis procedure is time-consuming and insensitive, with poor reproducibility. We have evaluated an innovative streamlined 96-well-plate-based platform utilizing InstantPC label for high-throughput, high-sensitivity glycan profiling, which is user-friendly, robust, and ready for automation. However, the limited availability of InstantPC labelled glycan standards has significantly hampered the applicability and transferability of this platform for expedited glycan structural profiling. To address this challenge, we have constructed a detailed InstantPC labelled glycan glucose unit database through analysis of human serum and a variety of other glycoproteins from various sources. Following preliminary hydrophilic interaction liquid chromatography with fluorescence detection separation and analysis, glycoproteins with complex glycan profiles were subjected to further fractionation by weak anion exchange hydrophilic interaction liquid chromatography and exoglycosidase sequential digestion for cross-validation of the glycan assignment. Hydrophilic interaction ultra-performance liquid chromatography coupled with electrospray ionization mass spectrometry was subsequently utilised for glycan fragmentation and accurate glycan mass confirmation. The constructed InstantPC glycan GU database is accurate and robust. It is believed that this database will enhance the application of the developed platform for high-throughput, high-sensitivity glycan profiling, and eventually advance glycan-based biopharmaceutical production and disease biomarker discovery.

Rationalizing the Stereoelectronic Influence of Interglycosidic Bond Conformations on the Reactivity of 1,4-O-Linked Disaccharide Donors

Disaccharide donors are key precursors in convergent glycan synthesis strategies. Unexpectedly, we observed that disaccharide thioglycosyl donors containing 1,4-O-linked α-glycosidic bonds are much more reactive than their β-analogues with the same protecting group pattern. Herein, we rationalized that such a difference in their reactivity is attributed to the conformation of the 1,4-O-interglycosidic bond which is controlled by anomeric and exo-anomeric effects. Moreover, the conformational preferences of these donors are dictated by the dihedral angles ϕ and ψ of their interglycosidic linkages and the torsional angle ω of their side chain along the C5-C6 bond. This fundamental research clarifies how the long-range stereoelectronic effects from the nonreducing end sugar can influence the reactivity of the leaving group at the reducing end and the behavior of disaccharide donors thereof.

Interplay of Trans- and Cis-Interactions of Glycolipids in Membrane Adhesion

Glycolipids mediate stable membrane adhesion of potential biological relevance. In this article, we investigate the trans- and cis-interactions of glycolipids in molecular dynamics simulations and relate these interactions to the glycolipid-induced average separations of membranes obtained from neutron scattering experiments. We find that the cis-interactions between glycolipids in the same membrane leaflet tend to strengthen the trans-interactions between glycolipids in apposing leaflets. The trans-interactions of the glycolipids in our simulations require local membrane separations that are significantly smaller than the average membrane separations in the neutron scattering experiments, which indicates an important role of membrane shape fluctuations in glycolipid trans-binding. Simulations at the experimentally measured average membrane separations provide a molecular picture of the interplay between glycolipid attraction and steric repulsion of the fluctuating membranes probed in the experiments.

The Flexibility of Oligosaccharides Unveiled Through Residual Dipolar Coupling Analysis

The intrinsic flexibility of glycans complicates the study of their structures and dynamics, which are often important for their biological function. NMR has provided insights into the conformational, dynamic and recognition features of glycans, but suffers from severe chemical shift degeneracy. We employed labelled glycans to explore the conformational behaviour of a β(1-6)-Glc hexasaccharide model through residual dipolar couplings (RDCs). RDC delivered information on the relative orientation of specific residues along the glycan chain and provided experimental clues for the existence of certain geometries. The use of two different aligning media demonstrated the adaptability of flexible oligosaccharide structures to different environments.

Structural dynamics and vibrational feature of N-Acetyl-d-glucosamine in aqueous solution

Molecular dynamics simulations and DFT calculations were performed for the demonstration of the structural dynamics and vibrational feature of N-Acetyl-d-glucosamine (NAG) in solution phase. The interactions between NAG and solvent molecules were evaluated through spatial distribution function and radial distribution function, and the preferred conformations of NAG in aqueous solution were revealed by cluster analysis. Results from normal mode analysis show that the solvent induced structural fluctuation of NAG could be reflected in the vibrational feature of specific chromophores, thus we can evaluate the molecular structure with the help of its vibrational signature based on the built correlation between molecular structure and vibrational frequencies of specific groups.

Targeted Chemical Modifications Identify Key Features of Carbohydrate Assemblies and Generate Tailored Carbohydrate Materials

The molecular level description of carbohydrate assemblies is hampered by their structural complexity and the lack of suitable analytical methods. Here, we employed systematic chemical modifications to identify key non-covalent interactions that triggered the supramolecular assembly of a disaccharide model. While some modifications disrupted the supramolecular organization, others were tolerated, delivering important information on the aggregation process. The screening identified new geometries, including nanotubes, and twisted ribbons that were characterized with electron tomography and electron diffraction (ED) methods. This work demonstrates that the combination of synthetic chemistry and ED methods is a powerful tool to draw correlations between the molecular structure and the nanoscale architecture of carbohydrate assemblies.

N-Glycosylation Enhances Conformational Flexibility of Protein Disulfide Isomerase Revealed by Microsecond Molecular Dynamics and Markov State Modeling

Secreted proteins of eukaryotes are decorated with branched carbohydrate oligomers called glycans. This fact is only starting to be considered for in silico investigations of protein dynamics. Using all-atom molecular dynamics (MD) simulations and Markov state modeling (MSM), we unveil the influence of glycans on the conformational flexibility of the multidomain protein disulfide isomerase (PDI), which is a ubiquitous chaperone in the endoplasmic reticulum (ER). Yeast PDI (yPDI) from Saccharomyces cerevisiae is glycosylated at asparagine side chains and the knowledge of its five modified sites enables a realistic computational modeling. We compare simulations of glycosylated and unglycosylated yPDI and find that the presence of glycan-glycan and glycan-protein interactions influences the flexibility of PDI in different ways. For example, glycosylation reduces interdomain interactions, shifting the conformational ensemble toward more open, extended structures. In addition, we compare our results on yPDI with structural information of homologous proteins such as human PDI (hPDI), which is natively unglycosylated. Interestingly, hPDI lacks a surface recess that is present in yPDI. We find that glycosylation of yPDI facilitates its catalytic site to reach close to this surface recess. Hence, this might point to a possible functional relevance of glycosylation in yeast to act on substrates, while glycosylation seems redundant for the human homologous protein. We conclude that glycosylation is fundamental for protein dynamics, making it a necessity for a truthful representation of the flexibility and function in in silico studies of glycoproteins.

A Computational Perspective on Molecular Recognition by Galectins

This article presents an overview of recent computational studies dedicated to the analysis of binding between galectins and small-molecule ligands. We first present a summary of the most popular simulation techniques adopted for calculating binding poses and binding energies, and then discuss relevant examples reported in the literature for the three main classes of galectins (dimeric, tandem and chimera). We show that simulation of galectin-ligand interactions is a mature field which has proven invaluable for completing and unraveling experimental observations. Future perspectives to further improve the accuracy and cost-effectiveness of existing computational approaches will involve the development of new schemes to account for solvation and entropy effects, which represent the main current limitations to the accuracy of computational results.

Development and Evaluation of GlycanDock: A Protein-Glycoligand Docking Refinement Algorithm in Rosetta

Carbohydrate chains are ubiquitous in the complex molecular processes of life. These highly diverse chains are recognized by a variety of protein receptors, enabling glycans to regulate many biological functions. High-resolution structures of protein-glycoligand complexes reveal the atomic details necessary to understand this level of molecular recognition and inform application-focused scientific and engineering pursuits. When experimental challenges hinder high-throughput determination of quality structures, computational tools can, in principle, fill the gap. In this work, we introduce GlycanDock, a residue-centric protein-glycoligand docking refinement algorithm developed within the Rosetta macromolecular modeling and design software suite. We performed a benchmark docking assessment using a set of 109 experimentally determined protein-glycoligand complexes as well as 62 unbound protein structures. The GlycanDock algorithm can sample and discriminate among protein-glycoligand models of native-like structural accuracy with statistical reliability from starting structures of up to 7 Å root-mean-square deviation in the glycoligand ring atoms. We show that GlycanDock-refined models qualitatively replicated the known binding specificity of a bacterial carbohydrate-binding module. Finally, we present a protein-glycoligand docking pipeline for generating putative protein-glycoligand complexes when only the glycoligand sequence and unbound protein structure are known. In combination with other carbohydrate modeling tools, the GlycanDock docking refinement algorithm will accelerate research in the glycosciences.

Molecular dynamic simulation: Structural insights of multi-stranded curdlan in aqueous solution

In this work, by using molecular dynamic simulation we provide microscale structure information which helps to reveal the molecular mechanisms concerning the multi-chain conformational behavior of short curdlan. Through simulations starting with different conformations of curldan dodecasaccharides, it is found that the right-handed triple helix is thermodynamically the most stable conformation in aqueous solutions, which is well maintained and stabilized by an inter-strand hydrogen bonding network of the C₂ hydroxyls. Unlike any predicted forms, the inter-strand hydrogen bonds exhibit a left-handed double helix pattern with preferred global orientations. Temperature REMD results suggest that the formation of triple helix is temperature sensitive, but the already formed triple helix is not. Investigation of curdlan with numbers of repeating units from 3 to 12 captures a critical value of 6, which in a way elucidates the relationship between the formation of triple helix and the chain length.

A Bittersweet Computational Journey among Glycosaminoglycans

Glycosaminoglycans (GAGs) are linear polysaccharides. In proteoglycans (PGs), they are attached to a core protein. GAGs and PGs can be found as free molecules, associated with the extracellular matrix or expressed on the cell membrane. They play a role in the regulation of a wide array of physiological and pathological processes by binding to different proteins, thus modulating their structure and function, and their concentration and availability in the microenvironment. Unfortunately, the enormous structural diversity of GAGs/PGs has hampered the development of dedicated analytical technologies and experimental models. Similarly, computational approaches (in particular, molecular modeling, docking and dynamics simulations) have not been fully exploited in glycobiology, despite their potential to demystify the complexity of GAGs/PGs at a structural and functional level. Here, we review the state-of-the art of computational approaches to studying GAGs/PGs with the aim of pointing out the “bitter” and “sweet” aspects of this field of research. Furthermore, we attempt to bridge the gap between bioinformatics and glycobiology, which have so far been kept apart by conceptual and technical differences. For this purpose, we provide computational scientists and glycobiologists with the fundamentals of these two fields of research, with the aim of creating opportunities for their combined exploitation, and thereby contributing to a substantial improvement in scientific knowledge.

Stereoelectronic effects in stabilizing protein-N-glycan interactions revealed by experiment and machine learning

The energetics of protein-carbohydrate interactions, central to many life processes, cannot yet be manipulated predictably. This is mostly due to an incomplete quantitative understanding of the enthalpic and entropic basis of these interactions in aqueous solution. Here, we show that stereoelectronic effects contribute to stabilizing protein-N-glycan interactions in the context of a cooperatively folding protein. Double-mutant cycle analyses of the folding data from 52 electronically varied N-glycoproteins demonstrate an enthalpy-entropy compensation depending on the electronics of the interacting side chains. Linear and nonlinear models obtained using quantum mechanical calculations and machine learning explain up to 79% and 97% of the experimental interaction energy variability, as inferred from the R² value of the respective models. Notably, the protein-carbohydrate interaction energies strongly correlate with the molecular orbital energy gaps of the interacting substructures. This suggests that stereoelectronic effects must be given a greater weight than previously thought for accurately modelling the short-range dispersive van der Waals interactions between the N-glycan and the protein.

Stereoelectronic effects in stabilizing protein-N-glycan interactions revealed by experiment and machine learning

The energetics of protein-carbohydrate interactions, central to many life processes, cannot yet be manipulated predictably. This is mostly due to an incomplete quantitative understanding of the enthalpic and entropic basis of these interactions in aqueous solution. Here, we show that stereoelectronic effects contribute to stabilizing protein-N-glycan interactions in the context of a cooperatively folding protein. Double-mutant cycle analyses of the folding data from 52 electronically varied N-glycoproteins demonstrate an enthalpy-entropy compensation depending on the electronics of the interacting side chains. Linear and nonlinear models obtained using quantum mechanical calculations and machine learning explain up to 79% and 97% of the experimental interaction energy variability, as inferred from the R² value of the respective models. Notably, the protein-carbohydrate interaction energies strongly correlate with the molecular orbital energy gaps of the interacting substructures. This suggests that stereoelectronic effects must be given a greater weight than previously thought for accurately modelling the short-range dispersive van der Waals interactions between the N-glycan and the protein.

Glycan-protein interactions determine kinetics of N-glycan remodeling

A hallmark of N-linked glycosylation in the secretory compartments of eukaryotic cells is the sequential remodeling of an initially uniform oligosaccharide to a site-specific, heterogeneous ensemble of glycostructures on mature proteins. To understand site-specific processing, we used protein disulfide isomerase (PDI), a model protein with five glycosylation sites, for molecular dynamics (MD) simulations and compared the result to a biochemical in vitro analysis with four different glycan processing enzymes. As predicted by an analysis of the accessibility of the N-glycans for their processing enzymes derived from the MD simulations, N-glycans at different glycosylation sites showed different kinetic properties for the processing enzymes. In addition, altering the tertiary structure of the glycoprotein PDI affected its N-glycan remodeling in a site-specific way. We propose that the observed differential N-glycan reactivities depend on the surrounding protein tertiary structure and lead to different glycan structures in the same protein through kinetically controlled processing pathways.

Atomistic glycoprotein simulations reveal a site-specific availability of glycan substrates in time-resolved mass spectrometry of maturating enzyme kinetics.

A perspective on the PDB's impact on the field of glycobiology

Structures deposited in the Protein Data Bank (PDB) facilitate our understanding of many biological processes including those that fall under the general category of glycobiology. However, structure-based studies of how glycans affect protein structure, how they are synthesized, and how they regulate other biological processes remain challenging. Despite the abundant presence of glycans on proteins and the dense layers of glycans that surround most of our cells, structures containing glycans are underrepresented in the PDB. There are sound reasons for this, including difficulties in producing proteins with well-defined glycosylation and the tendency of mobile and heterogeneous glycans to inhibit crystallization. Nevertheless, the structures we do find in the PDB, even some of the earliest deposited structures, have had an impact on our understanding of function. I highlight a few examples in this review and point to some promises for the future. Promises include new structures from methodologies, such as cryo-EM, that are less affected by the presence of glycans and experiment-aided computational methods that build on existing structures to provide insight into the many ways glycans affect biological function.

Laboratory Observation of, Astrochemical Search for, and Structure of Elusive Erythrulose in the Interstellar Medium

Rotational spectroscopy provides the most powerful means of identifying molecules of biological interest in the interstellar medium (ISM), but despite their importance, the detection of carbohydrates has remained rather elusive. Here, we present a comprehensive Fourier transform rotational spectroscopic study of elusive erythrulose, a sugar building block likely to be present in the ISM, employing a novel method of transferring the hygroscopic oily carbohydrate into the gas phase. The high sensitivity of the experiment allowed the rotational spectra of all monosubstituted isotopologue species of ¹³C-¹²C₃H₈O₄ to be recorded, which, together with quantum chemical calculations, enabled us to determine their equilibrium geometries (r_e^SE) with great precision. Searches employing the new experimental data for erythrulose have been undertaken in different ISM regions, so far including the cold areas Barnard 1, the pre-stellar core TMC-1, Sagittarius B2. Although no lines of erythrulose were found, this data will serve to enable future searches and possible detections in other ISM regions.

Systematic Structural Characterization of Chitooligosaccharides Enabled by Automated Glycan Assembly

Chitin, a polymer composed of β(1-4)-linked N-acetyl-glucosamine monomers, and its partially deacetylated analogue chitosan, are abundant biopolymers with outstanding mechanical as well as elastic properties. Their degradation products, chitooligosaccharides (COS), can trigger the innate immune response in humans and plants. Both material and biological properties are dependent on polymer length, acetylation, as well as the pH. Without well-defined samples, a complete molecular description of these factors is still missing. Automated glycan assembly (AGA) enabled rapid access to synthetic well-defined COS. Chitin-cellulose hybrid oligomers were prepared as important tools for a systematic structural analysis. Intramolecular interactions, identified by molecular dynamics simulations and NMR analysis, underscore the importance of the chitosan amino group for the stabilization of specific geometries.

Carbohydrate-Protein Interactions: Advances and Challenges

A carbohydrate, also called saccharide in biochemistry, is a biomolecule consisting of carbon (C), hydrogen (H) and oxygen (O) atoms. For example, sugars are low molecular-weight carbohydrates, and starches are high molecular-weight carbohydrates. Carbohydrates are the most abundant organic substances in nature and essential constituents of all living things. Protein-carbohydrate interactions play important roles in many biological processes, such as cell growth, differentiation, and aggregation. They also have broad applications in pharmaceutical drug design. In this review, we will summarize the characteristic features of protein-carbohydrate interactions and review the computational methods for structure prediction, energy calculations, and kinetic studies of protein-carbohydrate complexes. Finally, we will discuss the challenges in this field.

Exploring the Molecular Conformation Space by Soft Molecule-Surface Collision

Biomolecules function by adopting multiple conformations. Such dynamics are governed by the conformation landscape whose study requires characterization of the ground and excited conformation states. Here, the conformational landscape of a molecule is sampled by exciting an initial gas-phase molecular conformer into diverse conformation states, using soft molecule-surface collision (0.5-5.0 eV). The resulting ground and excited molecular conformations, adsorbed on the surface, are imaged at the single-molecule level. This technique permits the exploration of oligosaccharide conformations, until now, limited by the high flexibility of oligosaccharides and ensemble-averaged analytical methods. As a model for cellulose, cellohexaose chains are observed in two conformational extremes, the typical "extended" chain and the atypical "coiled" chain-the latter identified as the gas-phase conformer preserved on the surface. Observing conformations between these two extremes reveals the physical properties of cellohexaose, behaving as a rigid ribbon that becomes flexible when twisted. The conformation space of any molecule that can be electrosprayed can now be explored.

Quantification of the Resilience and Vulnerability of HIV-1 Native Glycan Shield at Atomistic Detail

Dense surface glycosylation on the HIV-1 envelope (Env) protein acts as a shield from the adaptive immune system. However, the molecular complexity and flexibility of glycans make experimental studies a challenge. Here we have integrated high-throughput atomistic modeling of fully glycosylated HIV-1 Env with graph theory to capture immunologically important features of the shield topology. This is the first complete all-atom model of HIV-1 Env SOSIP glycan shield that includes both oligomannose and complex glycans, providing physiologically relevant insights of the glycan shield. This integrated approach including quantitative comparison with cryo-electron microscopy data provides hitherto unexplored details of the native shield architecture and its difference from the high-mannose glycoform. We have also derived a measure to quantify the shielding effect over the antigenic protein surface that defines regions of relative vulnerability and resilience of the shield and can be harnessed for rational immunogen design.

Supramolecular Assembly and Chirality of Synthetic Carbohydrate Materials

Hierarchical carbohydrate architectures serve multiple roles in nature. Hardly any correlations between the carbohydrate chemical structures and the material properties are available due to the lack of standards and suitable analytic techniques. Therefore, designer carbohydrate materials remain highly unexplored, as compared to peptides and nucleic acids. A synthetic D-glucose disaccharide, DD, was chosen as a model to explore carbohydrate materials. Microcrystal electron diffraction (MicroED), optimized for oligosaccharides, revealed that DD assembled into highly crystalline left-handed helical fibers. The supramolecular architecture was correlated to the local crystal organization, allowing for the design of the enantiomeric right-handed fibers, based on the L-glucose disaccharide, LL, or flat lamellae, based on the racemic mixture. Tunable morphologies and mechanical properties suggest the potential of carbohydrate materials for nanotechnology applications.

GlycoTorch Vina: Docking Designed and Tested for Glycosaminoglycans

Glycosaminoglycans (GAGs) are a family of anionic carbohydrates that play an essential role in the physiology and pathology of all eukaryotic life forms. Experimental determination of GAG-protein complexes is challenging due to their difficult isolation from biological sources, natural heterogeneity, and conformational flexibility-including possible ring puckering of sulfated iduronic acid from ¹C₄ to ²S_O conformation. To overcome these challenges, we present GlycoTorch Vina (GTV), a molecular docking tool based on the carbohydrate docking program VinaCarb (VC). Our program is unique in that it contains parameters to model ²S_O sugars while also supporting glycosidic linkages specific to GAGs. We discuss how crystallographic models of carbohydrates can be biased by the choice of refinement software and structural dictionaries. To overcome these variations, we carefully curated 12 of the best available GAG and GAG-like crystal structures (ranging from tetra- to octasaccharides or longer) obtained from the PDB-REDO server and refined using the same protocol. Both GTV and VC produced pose predictions with a mean root-mean-square deviation (RMSD) of 3.1 Å from the native crystal structure-a statistically significant improvement when compared to AutoDock Vina (4.5 Å) and the commercial software Glide (5.9 Å). Examples of how real-space correlation coefficients can be used to better assess the accuracy of docking pose predictions are given. Comparisons between statistical distributions of empirical "salt bridge" interactions, relevant to GAGs, were compared to density functional theory (DFT) studies of model salt bridges, and water-mediated salt bridges; however, there was generally a poor agreement between these data. Water bridges appear to play an important, yet poorly understood, role in the structures of GAG-protein complexes. To aid in the rapid prototyping of future pose scoring functions, we include a module that allows users to include their own torsional and nonbonded parameters.

Glycan-mediated functional assembly of IL-1RI: structural insights into completion of the current description for immune response

Interleukin 1 Receptor type I (IL-1RI) is a multi-domain transmembrane receptor that triggers the inflammatory response. Understanding its detailed mechanism of action is crucial for treating immune disorders. IL-1RI is activated upon formation of its functional assembly that occurs by binding of the IL-1 cytokine and the accessory protein (Il-1RAcP) to it. X-ray crystallography, small-Angle X-ray Scattering and molecular dynamics simulation studies showed that IL-1RI adopts two types of 'compact' and 'extended' conformational states in its dynamical pattern. Furthermore, glycosylation has shown to play a critical role in its activation process. Here, classical and accelerated atomistic molecular dynamics were carried out to examine the role of full glycosylation of IL-1RI and IL-1RAcP in arrangement of the functional assembly. Simulations showed that the 'compact' and 'extended' IL-1RI form two types of 'cytokine-inaccessible-non-signaling' and 'cytokine-accessible-signaling' assemblies with the IL-1RacP, respectively that are both abiding in the presence of glycans. Suggesting that the cytokine binding to IL-1RI is not required for the formation of IL-1RI-IL-1RAcP complex and the 'compact' complex could act as a down-regulatory mechanism. The 'extended' complex is maintained by formation of several persistent hydrogen bonds between the IL-1RI-IL-1RAcP inter-connected glycans. Taken together, it was shown that full glycosylation regulates formation of the IL-1RI functional assembly and play critical role in cytokine biding and triggering the IL-1RI involved downstream pathways in the cell.Communicated by Ramaswamy H. Sarma.

NMR of glycoproteins: profiling, structure, conformation and interactions

In glycoproteins, carbohydrates are responsible for the selective interaction and tight regulation of cellular processes, constituting the main information transducer interface in protein-glycoprotein interactions. Increasing experimental and computational evidence suggest that such interactions often induce allosteric changes in the host protein, underlining the importance of studying intact glycoproteins. Technical issues have precluded such studies for years but, nowadays, a promising era is emerging where NMR spectroscopy, among other techniques, allows the characterization of the composition, structure and segmental dynamics of glycoproteins. In this review, we discuss such advances and highlight some selected examples. This novel technology unravels multiple new functional mechanisms, subtly hidden within the sugar code.

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.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Comparison of Methods for Bulk Automated Simulation of Glycosidic Bond Conformations

Six empirical force fields were tested for applicability to calculations for automated carbohydrate database filling. They were probed on eleven disaccharide molecules containing representative structural features from widespread classes of carbohydrates. The accuracy of each method was queried by predictions of nuclear Overhauser effects (NOEs) from conformational ensembles obtained from 50 to 100 ns molecular dynamics (MD) trajectories and their comparison to the published experimental data. Using various ranking schemes, it was concluded that explicit solvent MM3 MD yielded non-inferior NOE accuracy with newer GLYCAM-06, and ultimately PBE0-D3/def2-TZVP (Triple-Zeta Valence Polarized) Density Functional Theory (DFT) simulations. For seven of eleven molecules, at least one empirical force field with explicit solvent outperformed DFT in NOE prediction. The aggregate of characteristics (accuracy, speed, and compatibility) made MM3 dynamics with explicit solvent at 300 K the most favorable method for bulk generation of disaccharide conformation maps for massive database filling.

Analysis of the SARS-CoV-2 spike protein glycan shield reveals implications for immune recognition

Here we have generated 3D structures of glycoforms of the spike (S) glycoprotein from SARS-CoV-2, based on reported 3D structures and glycomics data for the protein produced in HEK293 cells. We also analyze structures for glycoforms representing those present in the nascent glycoproteins (prior to enzymatic modifications in the Golgi), as well as those that are commonly observed on antigens present in other viruses. These models were subjected to molecular dynamics (MD) simulation to determine the extent to which glycan microheterogeneity impacts the antigenicity of the S glycoprotein. Lastly, we have identified peptides in the S glycoprotein that are likely to be presented in human leukocyte antigen (HLA) complexes, and discuss the role of S protein glycosylation in potentially modulating the innate and adaptive immune response to the SARS-CoV-2 virus or to a related vaccine. The 3D structures show that the protein surface is extensively shielded from antibody recognition by glycans, with the notable exception of the ACE2 receptor binding domain, and also that the degree of shielding is largely insensitive to the specific glycoform. Despite the relatively modest contribution of the glycans to the total molecular weight of the S trimer (17% for the HEK293 glycoform) they shield approximately 40% of the protein surface.

Dissecting the conformation of glycans and their interactions with proteins

The use of in silico strategies to develop the structural basis for a rational optimization of glycan-protein interactions remains a great challenge. This problem derives, in part, from the lack of technologies to quantitatively and qualitatively assess the complex assembling between a glycan and the targeted protein molecule. Since there is an unmet need for developing new sugar-targeted therapeutics, many investigators are searching for technology platforms to elucidate various types of molecular interactions within glycan-protein complexes and aid in the development of glycan-targeted therapies. Here we discuss three important technology platforms commonly used in the assessment of the complex assembly of glycosylated biomolecules, such as glycoproteins or glycosphingolipids: Biacore analysis, molecular docking, and molecular dynamics simulations. We will also discuss the structural investigation of glycosylated biomolecules, including conformational changes of glycans and their impact on molecular interactions within the glycan-protein complex. For glycoproteins, secreted protein acidic and rich in cysteine (SPARC), which is associated with various lung disorders, such as chronic obstructive pulmonary disease (COPD) and lung cancer, will be taken as an example showing that the core fucosylation of N-glycan in SPARC regulates protein-binding affinity with extracellular matrix collagen. For glycosphingolipids (GSLs), Globo H ceramide, an important tumor-associated GSL which is being actively investigated as a target for new cancer immunotherapies, will be used to demonstrate how glycan structure plays a significant role in enhancing angiogenesis in tumor microenvironments.

Weak carbohydrate-carbohydrate interactions in membrane adhesion are fuzzy and generic

Carbohydrates such as the trisaccharide motif LeX are key constituents of cell surfaces. Despite intense research, the interactions between carbohydrates of apposing cells or membranes are not well understood. In this article, we investigate carbohydrate-carbohydrate interactions in membrane adhesion as well as in solution with extensive atomistic molecular dynamics simulations that exceed the simulation times of previous studies by orders of magnitude. For LeX, we obtain association constants of soluble carbohydrates, adhesion energies of lipid-anchored carbohydrates, and maximally sustained forces of carbohydrate complexes in membrane adhesion that are in good agreement with experimental results in the literature. Our simulations thus appear to provide a realistic, detailed picture of LeX-LeX interactions in solution and during membrane adhesion. In this picture, the LeX-LeX interactions are fuzzy, i.e. LeX pairs interact in a large variety of short-lived, bound conformations. For the synthetic tetrasaccharide Lac 2, which is composed of two lactose units, we observe similarly fuzzy interactions and obtain association constants of both soluble and lipid-anchored variants that are comparable to the corresponding association constants of LeX. The fuzzy, weak carbohydrate-carbohydrate interactions quantified in our simulations thus appear to be a generic feature of small, neutral carbohydrates such as LeX and Lac 2.

Molecular dynamic simulation: Conformational properties of single-stranded curdlan in aqueous solution

Recently, molecular dynamic simulation technique has been proved to be a powerful tool providing structural insights for better understanding the functionality of carbohydrates. Here, by using molecular dynamic simulation method we investigated the detailed conformational properties of the single-stranded curdlan with 12 glucose units. The results showed that the right-handed 6/1 helix structure was thermodynamically the most stable conformation in solution. The formation of the helix conformation was determined by many factors such as the glycosidic linkage, explicit water solvation and hydrogen bonds. When temperature was increased, the representative helix conformation was found becoming unstable giving rise to metastable conformations because when water mobility was accelerated with temperature, the hydrogen bonding strength between the curdlan chain and water went down, breaking the continuity of the hydrogen bonding network of water and hydroxyl groups. When the number of repeating glucose units varied from 6 to 24, the major helix conformation remained, but the conformational properties of longer chains were more apparently aff ;ected by chain flexibility.

Conformational Studies of Oligosaccharides

The conformation of a molecule strongly affects its function, as demonstrated for peptides and nucleic acids. This correlation is much less established for carbohydrates, the most abundant organic materials in nature. Recent advances in synthetic and analytical techniques have enabled the study of carbohydrates at the molecular level. Recurrent structural features were identified as responsible for particular biological activities or material properties. In this Minireview, recent achievements in the structural characterization of carbohydrates, enabled by systematic studies of chemically defined oligosaccharides, are discussed. These findings can guide the development of more potent glycomimetics. Synthetic carbohydrate materials by design can be envisioned.

N-glycosylation of High Mobility Group Box 1 protein (HMGB1) modulates the interaction with glycyrrhizin: A molecular modeling study

High Mobility Group Box 1 protein (HMGB1) is an abundant protein with multiple functions in cells, acting as a DNA chaperone and damage-associated molecular pattern molecule. It represents an attractive target for the treatment of inflammatory diseases and cancers. The plant natural product glycyrrhizin (GLR) is a well-characterized ligand of HMGB1 and a drug used to treat diverse liver and skin diseases. The drug is known to bind to each of the two adjacent HMG boxes of the non-glycosylated protein. In cells, HMGB1 is N-glycosylated at three asparagine residues located in boxes A and B, and these N-glycans are essential for the nucleocytoplasmic transport of the protein. But the impact of the N-glycans on drug binding is unknown. Here we have investigated the effect of the N-glycosylation of HMGB1 on its interaction with GLR using molecular modelling, after incorporation of three N-glycans on a Human HMGB1 structure (PDB code 2YRQ). Sialylated bi-antennary N-glycans were introduced on the protein and exposed in a folded or an extended conformation for the drug binding study. The docking of the drug was performed using both 18α- and 18β-epimers of GLR and the conformations and potential energy of interaction (ΔE) of the different drug-protein complexes were compared. The N-glycans do not shield the drug binding sites on boxes A and B but can modulate the drug-protein interaction, via both direct and indirect effects. The calculations indicate that binding of 18α/β-GLR to the HMG box is generally reduced when the protein is N-glycosylated vs. the non-glycosylated protein. In particular, the N-glycans in an extended configuration significantly weaken the binding of GLR to box-B. The effects of the N-glycans are mostly indirect, but in one case a direct contact with the drug, via a carbohydrate-carbohydrate interaction, was observed with 18β-GLR bound to Box-B of glycosylated HMGB1. For the first time, it is shown (at least in silico) that N-glycosylation, one of the many post-translational modifications of HMGB1, can affect drug binding.

Analysis of the SARS-CoV-2 spike protein glycan shield: implications for immune recognition

Here we have generated 3D structures of glycoforms of the spike (S) glycoprotein from SARS-CoV-2, based on reported 3D structures and glycomics data for the protein produced in HEK293 cells. We also analyze structures for glycoforms representing those present in the nascent glycoproteins (prior to enzymatic modifications in the Golgi), as well as those that are commonly observed on antigens present in other viruses. These models were subjected to molecular dynamics (MD) simulation to determine the extent to which glycan microheterogeneity impacts the antigenicity of the S glycoprotein. Lastly, we have identified peptides in the S glycoprotein that are likely to be presented in human leukocyte antigen (HLA) complexes, and discuss the role of S protein glycosylation in potentially modulating the adaptive immune response to the SARS-CoV-2 virus or to a related vaccine. The 3D structures show that the protein surface is extensively shielded from antibody recognition by glycans, with the exception of the ACE2 receptor binding domain, and also that the degree of shielding is largely insensitive to the specific glycoform. Despite the relatively modest contribution of the glycans to the total molecular weight (17% for the HEK293 glycoform) the level of surface shielding is disproportionately high at 42%.