Multifaceted Computational Modeling in Glycoscience

Glycoconjugates: Advances in modern medicines and human health

Glycans and their glycoconjugates are complex biomolecules that are crucial for various biological processes. Glycoconjugates are found in all domains of life. They are covalently linked to key biomolecules such as proteins and lipids to play a pivotal role in cell signaling, adhesion, and recognition. The diversity of glycan structures and the associated complexity of glycoconjugates is the reason for their role in intricate biosynthetic pathways. Glycoconjugates play an important role in various diseases where they are actively involved in the immune response as well as in the pathogenicity of infectious diseases. In addition, various autoimmune diseases have been linked to glycosylation defects of different biomolecules, making them an important molecule in the field of medicine. The glycoconjugates have been explored for the development of therapeutics and vaccines, representing a breakthrough in medical science. They also hold significance in research studies to understand the mechanisms behind various biological processes. Finally, glycoconjugates have found an emerging role in various industrial and environmental applications which have been discussed here.

Unravelling carbohydrate binding module 21 (CBM21) dynamics of interaction with amylose

The carbohydrate binding module 21 (CBM21) from Rhizopus oryzae is a dual-site CBM proposed to disrupt polysaccharide structures. Additionally, it serves as a purification tag in industry. CBM21 crystal structure features a Glc residue in an unusual 1S3 conformation, whose relevance for the CBM mechanism of action is unclear. In this context, we seek to contribute for the understanding of CBM21 mechanism of action by: i) investigating the role of the 1S3 conformation on carbohydrate recognition, and ii) characterize the protein-carbohydrate binding dynamics using molecular dynamics and metadynamics simulations at MM and QM/MM levels. Results indicate the 1S3 Glc conformation is unlikely to occur under biological conditions, being originated from the crystallographic environment. CBM21 binding to small ligands appears transient and unstable, while protein dimerization and polysaccharide chain size influence complex stability. In interactions with amylose, CBM21 exhibits a repeated unbinding followed by re-binding, while simultaneously alternating between binding sites I and II. These results suggest that CBM21 acts through transient interactions, directing carbohydrates to the catalytic center rather than forming strong and long-lasting bonds with carbohydrates. Accordingly, we expect such atomistic depiction of CBM21 mechanism could aid in CBM design targeting biotechnological applications.

Structural characterization, molecular dynamic simulation, and conformational visualization of a water-soluble glucan with high molecular weight from Gastrodia elata Blume

Polysaccharides have been widely used in the development of natural drugs and health food. However, polysaccharide characterization lags due to inherently complicated features and the limitations of existing detection approaches. We aimed to provide new insight into the fine structure and conformational visualization of polysaccharides from Gastrodia elata Blume, a medicinal and edible plant. A water-soluble polysaccharide (GEP2-6) with the high molecular weight of 2.7 × 106 Da was first obtained, and its purity reached 99.2 %. Chemical and spectroscopic analyses jointly revealed that GEP2-6 was a glucan linked by α-(1 → 4) and α-(1 → 6) glycosidic bonds. After enzymolysis, the local structure of GEP2-6 included α-1,4-Glcp, α-1,6-Glcp, α-1,4,6-Glcp, and α-1-Glcp at a molar ratio of 31.27∶1.32∶1.08∶0.93. The glycosidic linkage pattern of repeating units was further simulated by a glycan database and spatial examination software. The good dissolution performance was interpreted by dynamics simulation and practical molecular characteristics. Spherical flexible chains and the porous stable conformation were corroborated using atomic force microscopy. In addition, GEP2-6 could effectively scavenge DPPH and hydroxyl radicals as a promising natural antioxidant. These efforts will contribute to the expansion of clinical applications of this G. elata polysaccharide and the structural elucidation for macromolecular polysaccharides combined with traditional and modern analysis techniques.

Glycans modulate lipid binding in Lili-Mip lipocalin protein: insights from molecular simulations and protein network analyses

The unique viviparous Pacific Beetle cockroaches provide nutrition to their embryo by secreting milk proteins Lili-Mip, a lipid-binding glycoprotein that crystallises in-vivo. The resolved in-vivo crystal structure of variably glycosylated Lili-Mip shows a classical Lipocalin fold with an eight-stranded antiparallel beta-barrel enclosing a fatty acid. The availability of physiologically unaltered glycoprotein structure makes Lili-Mip a very attractive model system to investigate the role of glycans on protein structure, dynamics, and function. Towards that end, we have employed all-atom molecular dynamics simulations on various glycosylated stages of a bound and free Lili-Mip protein and characterised the impact of glycans and the bound lipid on the dynamics of this glycoconjugate. Our work provides important molecular-level mechanistic insights into the role of glycans in the nutrient storage function of the Lili-Mip protein. Our analyses show that the glycans stabilise spatially proximal residues and regulate the low amplitude opening motions of the residues at the entrance of the binding pocket. Glycans also preserve the native orientation and conformational flexibility of the ligand. However, we find that either deglycosylation or glycosylation with high-mannose and paucimannose on the core glycans, which better mimic the natural insect glycosylation state, significantly affects the conformation and dynamics. A simple but effective distance- and correlation-based network analysis of the protein also reveals the key residues regulating the barrel's architecture and ligand binding characteristics in response to glycosylation.

Many Roles of Carbohydrates: A Computational Spotlight on the Coronavirus S Protein Binding

Glycosylation is one of the post-translational modifications with more than 50% of human proteins being glycosylated. The exact nature and chemical composition of glycans are inaccessible to X-ray or cryo-electron microscopy imaging techniques. Therefore, computational modeling studies and molecular dynamics must be used as a "computational microscope". The spike (S) protein of SARS-CoV-2 is heavily glycosylated, and a few glycans play a more functional role "beyond shielding". In this mini-review, we discuss computational investigations of the roles of specific S-protein and ACE2 glycans in the overall ACE2-S protein binding. We highlight different functions of specific glycans demonstrated in myriad computational models and simulations in the context of the SARS-CoV-2 virus binding to the receptor. We also discuss interactions between glycocalyx and the S protein, which may be utilized to design prophylactic polysaccharide-based therapeutics targeting the S protein. In addition, we underline the recent emergence of coronavirus variants and their impact on the S protein and its glycans.

Molecular Gas-Phase Conformational Ensembles

Accurately determining the global minima of a molecular structure is important in diverse scientific fields, including drug design, materials science, and chemical synthesis. Conformational search engines serve as valuable tools for exploring the extensive conformational space of molecules and for identifying energetically favorable conformations. In this study, we present a comparison of Auto3D, CREST, Balloon, and ETKDG (from RDKit), which are freely available conformational search engines, to evaluate their effectiveness in locating global minima. These engines employ distinct methodologies, including machine learning (ML) potential-based, semiempirical, and force field-based approaches. To validate these methods, we propose the use of collisional cross-section (CCS) values obtained from ion mobility-mass spectrometry studies. We hypothesize that experimental gas-phase CCS values can provide experimental evidence that we likely have the global minimum for a given molecule. To facilitate this effort, we used our gas-phase conformation library (GPCL) which currently consists of the full ensembles of 20 small molecules and can be used by the community to validate any conformational search engine. Further members of the GPCL can be readily created for any molecule of interest using our standard workflow used to compute CCS values, expanding the ability of the GPCL in validation exercises. These innovative validation techniques enhance our understanding of the conformational landscape and provide valuable insights into the performance of conformational generation engines. Our findings shed light on the strengths and limitations of each search engine, enabling informed decisions for their utilization in various scientific fields, where accurate molecular structure determination is crucial for understanding biological activity and designing targeted interventions. By facilitating the identification of reliable conformations, this study significantly contributes to enhancing the efficiency and accuracy of molecular structure determination, with particular focus on metabolite structure elucidation. The findings of this research also provide valuable insights for developing effective workflows for predicting the structures of unknown compounds with high precision.

Exploring Ring Conformation in Uronate Monosaccharides: Insights from Ab Initio Calculations and Classical Molecular Dynamics Simulations

The study focuses on the conformational properties of biologically relevant monosaccharides belonging to the group of uronates: α-l-iduronate, O2-sulfated-α-l-iduronate, and O2-sulfated-α-l-guluronate, either unfunctionalized or O1-methylated. We applied the previously proposed two-step methodology, combining classical MD simulations and subsequent ab initio (QM) calculations, performed on a rationally subsampled set of molecular configurations. We found that, regardless of the number of molecular configurations considered, the level of theory, and the weighting scheme applied, none of the QM approaches is capable of predicting the correct conformational equilibrium of sulfated iduronates as long as the tight counterion binding is not considered. Multicenter, ring-shape-specific binding of either Na+ or Ca2+ ions drastically shifts the conformational equilibrium of the pyranose ring in sulfated iduronates toward 1C4 but does not significantly affect the conformation of non-sulfated compounds. A similar shift is observed upon the protonation of carboxyl groups in all iduronates. In addition, we report a set of average J-coupling constant values related to vicinal protons bound to the pyranose ring of iduronates and corresponding to each of the three main groups of ring conformers, i.e., 4C1, B/S (boat/skew boat), and 1C4. In combination with the conformational energies or with the experimental data, these values allowed the relative proportions of the ring conformers to be estimated and the Karplus-type equations linking the 3JHH-coupling constants to the torsion angles within the pyranose ring to be refined.

Gas Phase Conformation of Trisaccharides and Core Pentasaccharide: A Three-Step Tree-Based Sampling and Quantum Mechanical Computational Approach

As an important component of N-linked glycoproteins, the core pentasaccharide is highly crucial to the potential application prospect of glycoprotein. However, the gas phase conformation study is a challenging one due to the size and complexity of the molecule, together with the necessity to rely on quantum chemistry modeling for relevant energetics and structures. In this paper, the structures of the trisaccharides and core pentasaccharides in N-linked glycans in the gas phase were constructed by a three-step tree-based (TSTB) sampling. Since single point energies of all the conformers are calculated at the temperature of zero, it is necessary to evaluate the stability at a high temperature. We calculate the Gibbs free energies using the standard thermochemistry model (T = 298.15 K). For trimannose, the energetic ordering at 298.15 K can be strongly changed compared to 0 K. Moreover, two structures of trimannose with high energies at 0 K are considered to provide a much better match of IR vibration signatures with the low Gibbs free energies. On this basis, the core pentasaccharide was constructed in three ways. The building configurations of core pentasaccharide were optimized to obtain reasonable low-energy stable conformers. Fortunately, the lowest-energy structure of core pentasaccharide is eventually the minimum at 0 K and 298.15 K. Furthermore, spectrum analysis of core pentasaccharide was carried out. Although poorly resolved, its contour from the experiment was in qualitative correspondence with the computed IR spectrum associated with its minimum free energy structure. A large number of strongly and weakly hydrogen-bonded hydroxyl and acetylamino groups contribute to a highly congested set of overlapping bands. Compared with traditional conformation generators, the TSTB sampling is employed to efficiently and comprehensively obtain preferred conformers of larger saccharides with lower energy.

Mechanism of multivalent glycoconjugate-lectin interaction: An update

Lectins are predominantly oligomeric proteins with several binding sites per molecule. Glycoconjugates are their natural ligands, which often possess multiple binding epitopes. Thus, lectin-glycoconjugate interactions are mostly multivalent in nature. The mechanism of multivalent binding is fundamentally different from those described for monovalent interactions in textbooks and research papers. Over the years, binding studies that make use of different lectins and a variety of multivalent glycoconjugate ligands were conducted in order to understand the underlying principles of multivalency. Starting with seemingly simple synthetic multivalent analogs, systematic studies were carried out using natural glycoconjugate ligands with increasing valency and complexity. Those ligands included multivalent glycoproteins, polyvalent polysaccharides, including glycosaminoglycans, as well as supra-valent mucins and proteoglycans. Models and mechanisms of multivalent binding derived from quantitative data are summarized in the present updated review.

Computational modeling of protein-carbohydrate interactions: Current trends and future challenges

The article leads the reader through an up-to-date presentation of the concepts, developments, and main applications of computational modeling to study protein-carbohydrate interactions. It follows with the presentation of some current issues and perspectives arising from the expected evolution of generic methodological developments in deep learning, immersive analytics, and virtual reality for molecular visualization and data management. Such methodological developments for macromolecular interactions would greatly benefit a wide range of scientific endeavors in the field of carbohydrate chemistry and biochemistry, including the following interrelated efforts dealing with highly crowded media, with examples concerning glycoside transferases, the extracellular matrix, and the exploration of interactions between complex carbohydrates and intrinsically disordered proteins.

Chitosan Polysaccharides from a Polarizable Multiscale Approach

We report simulations of chitosan polysaccharides in the aqueous phase, at infinite dilute conditions and zero ionic strength. Those simulations are performed by means of a polarizable multiscale modeling scheme that relies on a polarizable all atom force field to model solutes and on a polarizable solvent coarse grained approach. Force field parameters are assigned only from quantum chemistry ab initio data. We simulate chitosan monomer units, dimers and 50-long chains. Regarding the 50-long chains we simulate three sets of ten randomly built chain replica at three different pH conditions (corresponding to different chain protonation states, the chain degree of deacetylation is 85%). Our simulations show the persistence length of 50-long chitosan chains at strong acidic conditions (pH <5) to be 24 ± 2 nm (at weak/negligible ionic strength conditions), and to be 1 order of magnitude shorter at usual pH conditions. Our simulation data support the most recent simulation and experimental studies devoted to chitosan polysaccharides in the aqueous phase.

Mechanisms of interactions in pattern-recognition of common glycostructures across pectin-derived heteropolysaccharides by Toll-like receptor 4

Complex pectin, originating from terrestrial plant cell walls has been attracting research attention as a promising source of a new innate immune modulator. Numerous bioactive polysaccharides associated with pectin are newly reported every year, but the general mechanism of their immunological action remains unclear owing to the complexity and heterogeneity of pectin. Herein, we systematically investigated the interactions in pattern-recognition for common glycostructures of pectic heteropolysaccharides (HPSs) by Toll-like receptors (TLRs). The compositional similarity of glycosyl residues derived from pectic HPS was confirmed by conducting systematic reviews, leading to molecular modeling of representative pectic segments. Via structural investigation, the inner concavity of leucine-rich repeats of TLR4 was predicted to act as a binding motif for carbohydrate recognition, and subsequent simulations predicted the binding modes and conformations. We experimentally demonstrated that pectic HPS exhibits the non-canonical and multivalent binding aspects for TLR4 resulting in receptor activation. Furthermore, we showed that pectic HPSs were selectively clustered with TLR4 during endocytosis, inducing downstream signals to cause phenotypic activation of macrophages. Overall, we have presented a better explanation for the pattern recognition of pectic HPS and further proposed an approach to understand the interaction between complex carbohydrates and proteins.

Multiscale QM/MM Simulations Identify the Roles of Asp239 and 1-OH···Nucleophile in Transition State Stabilization in Arabidopsis thaliana Cell-Wall Invertase 1

Arabidopsis thaliana cell-wall invertase 1 (AtCWIN1), a key enzyme in sucrose metabolism in plants, catalyzes the hydrolysis of sucrose into fructose and glucose. AtCWIN1 belongs to the glycoside hydrolase GH-J clan, where two carboxylate residues (Asp23 and Glu203 in AtCWIN1) are well documented as a nucleophile and an acid/base catalyst. However, details at the atomic level about the role of neighboring residues and enzyme-substrate interactions during catalysis are not fully understood. Here, quantum mechanical/molecular mechanical (QM/MM) free-energy simulations were carried out to clarify the origin of the observed decreased rates in Asp239Ala, Asp239Asn, and Asp239Phe in AtCWIN1 compared to the wild type and delineate the role of Asp239 in catalysis. The glycosylation and deglycosylation steps were considered in both wild type and mutants. Deglycosylation is predicted to be the rate-determining step in the reaction, with a calculated overall free-energy barrier of 15.9 kcal/mol, consistent with the experimental barrier (15.3 kcal/mol). During the reaction, the -1 furanosyl ring underwent a conformational change corresponding to 3E ↔ [E2]⧧ ↔ 1E according to the nomenclature of saccharide structures along the full catalytic reaction. Asp239 was found to stabilize not only the transition state but also the fructosyl-enzyme intermediate, which explains findings from previous structural and mutagenesis experiments. The 1-OH···nucleophile interaction has been found to provide an important contribution to the transition state stabilization, with a contribution of ∼7 kcal/mol, and affected glycosylation more significantly than deglycosylation. This study provides molecular insights that improve the current understanding of sucrose binding and hydrolysis in members of clan GH-J, which may benefit protein engineering research. Finally, a rationale on the sucrose inhibitor configuration in chicory 1-FEH IIa, proposed a long time ago in the literature, is also provided based on the QM/MM calculations.

Application of bioanalytical and computational methods in decoding the roles of glycans in host-pathogen interactions

Host-pathogen interactions (HPIs) are complex processes that require tight regulation. A common regulatory mechanism of HPIs is through glycans of either host cells or pathogens. Due to their diverse sequences, complex structures, and conformations, studies of glycans require highly sensitive and powerful tools. Recent improvements in technology have enabled the application of many bioanalytical techniques and modeling methods to investigate glycans and their mechanisms in HPIs. This mini-review highlights how these advances have been used to understand the role glycans play in HPIs in the past 2 years.

Diversity of Bioinspired Hydrogels: From Structure to Applications

Hydrogels are three-dimensional networks with a variety of structures and functions that have a remarkable ability to absorb huge amounts of water or biological fluids. They can incorporate active compounds and release them in a controlled manner. Hydrogels can also be designed to be sensitive to external stimuli: temperature, pH, ionic strength, electrical or magnetic stimuli, specific molecules, etc. Alternative methods for the development of various hydrogels have been outlined in the literature over time. Some hydrogels are toxic and therefore are avoided when obtaining biomaterials, pharmaceuticals, or therapeutic products. Nature is a permanent source of inspiration for new structures and new functionalities of more and more competitive materials. Natural compounds present a series of physico-chemical and biological characteristics suitable for biomaterials, such as biocompatibility, antimicrobial properties, biodegradability, and nontoxicity. Thus, they can generate microenvironments comparable to the intracellular or extracellular matrices in the human body. This paper discusses the main advantages of the presence of biomolecules (polysaccharides, proteins, and polypeptides) in hydrogels. Structural aspects induced by natural compounds and their specific properties are emphasized. The most suitable applications will be highlighted, including drug delivery, self-healing materials for regenerative medicine, cell culture, wound dressings, 3D bioprinting, foods, etc.

Glycosaminoglycans: What Remains To Be Deciphered?

Glycosaminoglycans (GAGs) are complex polysaccharides exhibiting a vast structural diversity and fulfilling various functions mediated by thousands of interactions in the extracellular matrix, at the cell surface, and within the cells where they have been detected in the nucleus. It is known that the chemical groups attached to GAGs and GAG conformations comprise "glycocodes" that are not yet fully deciphered. The molecular context also matters for GAG structures and functions, and the influence of the structure and functions of the proteoglycan core proteins on sulfated GAGs and vice versa warrants further investigation. The lack of dedicated bioinformatic tools for mining GAG data sets contributes to a partial characterization of the structural and functional landscape and interactions of GAGs. These pending issues will benefit from the development of new approaches reviewed here, namely (i) the synthesis of GAG oligosaccharides to build large and diverse GAG libraries, (ii) GAG analysis and sequencing by mass spectrometry (e.g., ion mobility-mass spectrometry), gas-phase infrared spectroscopy, recognition tunnelling nanopores, and molecular modeling to identify bioactive GAG sequences, biophysical methods to investigate binding interfaces, and to expand our knowledge and understanding of glycocodes governing GAG molecular recognition, and (iii) artificial intelligence for in-depth investigation of GAGomic data sets and their integration with proteomics.

Enzymatic Glyco-Modification of Synthetic Membrane Systems

The present report assesses the capability of a soluble glycosyltransferase to modify glycolipids organized in two synthetic membrane systems that are attractive models to mimic cell membranes: giant unilamellar vesicles (GUVs) and supported lipid bilayers (SLBs). The objective was to synthesize the Gb3 antigen (Galα1,4Galβ1,4Glcβ-Cer), a cancer biomarker, at the surface of these membrane models. A soluble form of LgtC that adds a galactose residue from UDP-Gal to lactose-containing acceptors was selected. Although less efficient than with lactose, the ability of LgtC to utilize lactosyl-ceramide as an acceptor was demonstrated on GUVs and SLBs. The reaction was monitored using the B-subunit of Shiga toxin as Gb3-binding lectin. Quartz crystal microbalance with dissipation analysis showed that transient binding of LgtC at the membrane surface was sufficient for a productive conversion of LacCer to Gb3. Molecular dynamics simulations provided structural elements to help rationalize experimental data.

Molecular Dynamics-Based Comparative Analysis of Chondroitin and Dermatan Sulfates

Glycosaminoglycans (GAGs) are a class of linear anionic periodic polysaccharides containing disaccharide repetitive units. These molecules interact with a variety of proteins in the extracellular matrix and so participate in biochemically crucial processes such as cell signalling affecting tissue regeneration as well as the onset of cancer, Alzheimer's or Parkinson's diseases. Due to their flexibility, periodicity and chemical heterogeneity, often termed "sulfation code", GAGs are challenging molecules both for experiments and computation. One of the key questions in the GAG research is the specificity of their intermolecular interactions. In this study, we make a step forward to deciphering the "sulfation code" of chondroitin sulfates-4,6 (CS4, CS6, where the numbers correspond to the position of sulfation in NAcGal residue) and dermatan sulfate (DS), which is different from CSs by the presence of IdoA acid instead of GlcA. We rigorously investigate two sets of these GAGs in dimeric, tetrameric and hexameric forms with molecular dynamics-based descriptors. Our data clearly suggest that CS4, CS6 and DS are substantially different in terms of their structural, conformational and dynamic properties, which contributes to the understanding of how these molecules can be different when they bind proteins, which could have practical implications for the GAG-based drug design strategies in the regenerative medicine.

The In Silico Characterization of Monocotyledonous α-l-Arabinofuranosidases on the Example of Maize

Plant α-l-arabinofuranosidases remove terminal arabinose from arabinose-containing substrates such as plant cell wall polysaccharides, including arabinoxylans, arabinogalactans, and arabinans. In plants, de-arabinosylation of cell wall polysaccharides accompanies different physiological processes such as fruit ripening and elongation growth. In this report, we address the diversity of plant α-l-arabinofuranosidases of the glycoside hydrolase (GH) family 51 through their phylogenetic analysis as well as their structural features. The CBM4-like domain at N-terminus was found to exist only in GH51 family proteins and was detected in almost 90% of plant sequences. This domain is similar to bacterial CBM4, but due to substitutions of key amino acid residues, it does not appear to be able to bind carbohydrates. Despite isoenzymes of GH51 being abundant, in particular in cereals, almost half of the GH51 proteins in Poales have a mutation of the acid/base residue in the catalytic site, making them potentially inactive. Open-source data on the transcription and translation of GH51 isoforms in maize were analyzed to discuss possible functions of individual isoenzymes. The results of homology modeling and molecular docking showed that the substrate binding site can accurately accommodate terminal arabinofuranose and that arabinoxylan is a more favorable ligand for all maize GH51 enzymes than arabinan.

Primary Structure of Glycans by NMR Spectroscopy

Glycans, carbohydrate molecules in the realm of biology, are present as biomedically important glycoconjugates and a characteristic aspect is that their structures in many instances are branched. In determining the primary structure of a glycan, the sugar components including the absolute configuration and ring form, anomeric configuration, linkage(s), sequence, and substituents should be elucidated. Solution state NMR spectroscopy offers a unique opportunity to resolve all these aspects at atomic resolution. During the last two decades, advancement of both NMR experiments and spectrometer hardware have made it possible to unravel carbohydrate structure more efficiently. These developments applicable to glycans include, inter alia, NMR experiments that reduce spectral overlap, use selective excitations, record tilted projections of multidimensional spectra, acquire spectra by multiple receivers, utilize polarization by fast-pulsing techniques, concatenate pulse-sequence modules to acquire several spectra in a single measurement, acquire pure shift correlated spectra devoid of scalar couplings, employ stable isotope labeling to efficiently obtain homo- and/or heteronuclear correlations, as well as those that rely on dipolar cross-correlated interactions for sequential information. Refined computer programs for NMR spin simulation and chemical shift prediction aid the structural elucidation of glycans, which are notorious for their limited spectral dispersion. Hardware developments include cryogenically cold probes and dynamic nuclear polarization techniques, both resulting in enhanced sensitivity as well as ultrahigh field NMR spectrometers with a ¹H NMR resonance frequency higher than 1 GHz, thus improving resolution of resonances. Taken together, the developments have made and will in the future make it possible to elucidate carbohydrate structure in great detail, thereby forming the basis for understanding of how glycans interact with other molecules.

Glycan Conformation in the Heavily Glycosylated Protein, CEACAM1

Glycans attached to glycoproteins can contribute to stability, mediate interactions with other proteins, and initiate signal transduction. Glycan conformation, which is critical to these processes, is highly variable and often depicted as sampling a multitude of conformers. These conformers can be generated by molecular dynamics simulations, and more inclusively by accelerated molecular dynamics, as well as other extended sampling methods. However, experimental assessments of the contribution that various conformers make to a native ensemble are rare. Here, we use long-range pseudo-contact shifts (PCSs) of NMR resonances from an isotopically labeled glycoprotein to identify preferred conformations of its glycans. The N-terminal domain from human Carcinoembryonic Antigen Cell Adhesion Molecule 1, hCEACAM1-Ig1, was used as the model glycoprotein in this study. It has been engineered to include a lanthanide-ion-binding loop that generates PCSs, as well as a homogeneous set of three ¹³C-labeled N-glycans. Analysis of the PCSs indicates that preferred glycan conformers have extensive contacts with the protein surface. Factors leading to this preference appear to include interactions between N-acetyl methyls of GlcNAc residues and hydrophobic surface pockets on the protein surface.

Glycoinformatics in the Artificial Intelligence Era

Artificial intelligence (AI) methods have been and are now being increasingly integrated in prediction software implemented in bioinformatics and its glycoscience branch known as glycoinformatics. AI techniques have evolved in the past decades, and their applications in glycoscience are not yet widespread. This limited use is partly explained by the peculiarities of glyco-data that are notoriously hard to produce and analyze. Nonetheless, as time goes, the accumulation of glycomics, glycoproteomics, and glycan-binding data has reached a point where even the most recent deep learning methods can provide predictors with good performance. We discuss the historical development of the application of various AI methods in the broader field of glycoinformatics. A particular focus is placed on shining a light on challenges in glyco-data handling, contextualized by lessons learnt from related disciplines. Ending on the discussion of state-of-the-art deep learning approaches in glycoinformatics, we also envision the future of glycoinformatics, including development that need to occur in order to truly unleash the capabilities of glycoscience in the systems biology era.

Multivalent glycocyclopeptides: conjugation methods and biological applications

Click chemistry was extensively used to decorate synthetic multivalent scaffolds with glycans to mimic the cell surface glycocalyx and to develop applications in glycosciences. Conjugation methods such as oxime ligation, copper(I)-catalyzed alkyne-azide cycloaddition, thiol-ene coupling, squaramide coupling or Lansbury aspartylation proved particularly suitable to achieve this purpose. This review summarizes the synthetic strategies that can be used either in a stepwise manner or in an orthogonal one-pot approach, to conjugate multiple copies of identical or different glycans to cyclopeptide scaffolds (namely multivalent glycocyclopeptides) having different size, valency, geometry and molecular composition. The second part of this review will describe the potential of these structures to interact with various carbohydrate binding proteins or to stimulate immunity against tumor cells.

Molecular Insight into the Self-Assembly Process of Cellulose Iβ Microfibril

The self-assembly process of β-D-glucose oligomers on the surface of cellulose Iβ microfibril involves crystallization, and this process is analyzed herein, in terms of the length and flexibility of the oligomer chain, by means of molecular dynamics (MD) simulations. The characterization of this process involves the structural relaxation of the oligomer, the recognition of the cellulose I microfibril, and the formation of several hydrogen bonds (HBs). This process is monitored on the basis of the changes in non-bonded energies and the interaction with hydrophilic and hydrophobic crystal faces. The oligomer length is considered a parameter for capturing insight into the energy landscape and its stability in the bound form with the cellulose I microfibril. We notice that the oligomer–microfibril complexes are more stable by increasing the number of hydrogen bond interactions, which is consistent with a gain in electrostatic energy. Our studies highlight the interaction with hydrophilic crystal planes on the microfibril and the acceptor role of the flexible oligomers in HB formation. In addition, we study by MD simulation the interaction between a protofibril and the cellulose I microfibril in solution. In this case, the main interaction consists of the formation of hydrogen bonds between hydrophilic faces, and those HBs involve donor groups in the protofibril.