Strategies for carbohydrate model building, refinement and validation

Monitoring carbohydrate 3D structure quality with the Privateer database

The remediation of the carbohydrate data of the Protein Data Bank (PDB) has brought numerous enhancements to the findability and interpretability of deposited glycan structures, yet crucial quality indicators are either missing or hard to find on the PDB pages. Without a way to access wider glycochemical context, problematic structures may be taken as fact by keen but inexperienced scientists. The Privateer software is a validation and analysis tool that provides access to a number of metrics and links to external experimental resources, allowing users to evaluate structures using carbohydrate-specific methods. Here, we present the Privateer database, a free resource that aims to complement the growing glycan content of the PDB.

Online carbohydrate 3D structure validation with the Privateer web app

Owing to the difficulties associated with working with carbohydrates, validating glycan 3D structures prior to deposition into the Protein Data Bank has become a staple of the structure-solution pipeline. The Privateer software provides integrative methods for the validation, analysis, refinement and graphical representation of 3D atomic structures of glycans, both as ligands and as protein modifiers. While Privateer is free software, it requires users to install any of the structural biology software suites that support it or to build it from source code. Here, the Privateer web app is presented, which is always up to date and available to be used online (https://privateer.york.ac.uk) without installation. This self-updating tool, which runs locally on the user's machine, will allow structural biologists to simply and quickly analyse carbohydrate ligands and protein glycosylation from a web browser whilst retaining all confidential information on their devices.

Highly functionalized diaminocyclopentanes: A new route to potent and selective inhibitors of human O-GlcNAcase

A new class of compounds inhibiting de-O-glycosylation of proteins has been identified. Highly substituted diaminocyclopentanes are impressively selective reversible non-transition state O-β-N-acetyl-d-glucosaminidase (O-GlcNAcase) inhibitors. The ease of preparative access and remarkable biological activities provide highly viable leads for the development of anti-tau-phosphorylation agents with a view to eventually ameliorating Alzheimer's disease.

Structural basis for Lewis antigen synthesis by the α1,3-fucosyltransferase FUT9

Mammalian cell surface and secreted glycoproteins exhibit remarkable glycan structural diversity that contributes to numerous physiological and pathogenic interactions. Terminal glycan structures include Lewis antigens synthesized by a collection of α1,3/4-fucosyltransferases (CAZy GT10 family). At present, the only available crystallographic structure of a GT10 member is that of the Helicobacter pylori α1,3-fucosyltransferase, but mammalian GT10 fucosyltransferases are distinct in sequence and substrate specificity compared with the bacterial enzyme. Here, we determined crystal structures of human FUT9, an α1,3-fucosyltransferase that generates Lewisx and Lewisy antigens, in complex with GDP, acceptor glycans, and as a FUT9-donor analog-acceptor Michaelis complex. The structures reveal substrate specificity determinants and allow prediction of a catalytic model supported by kinetic analyses of numerous active site mutants. Comparisons with other GT10 fucosyltransferases and GT-B fold glycosyltransferases provide evidence for modular evolution of donor- and acceptor-binding sites and specificity for Lewis antigen synthesis among mammalian GT10 fucosyltransferases.

Germline-targeting HIV-1 Env vaccination induces VRC01-class antibodies with rare insertions

Targeting germline (gl-) precursors of broadly neutralizing antibodies (bNAbs) is acknowledged as an important strategy for HIV-1 vaccines. The VRC01-class of bNAbs is attractive because of its distinct genetic signature. However, VRC01-class bNAbs often require extensive somatic hypermutation, including rare insertions and deletions. We describe a BG505 SOSIP trimer, termed GT1.2, to optimize binding to gl-CH31, the unmutated common precursor of the CH30-34 bNAb lineage that acquired a large CDRH1 insertion. The GT1.2 trimer activates gl-CH31 naive B cells in knock-in mice, and B cell responses could be matured by selected boosting immunogens to generate cross-reactive Ab responses. Next-generation B cell sequencing reveals selection for VRC01-class mutations, including insertions in CDRH1 and FWR3 at positions identical to VRC01-class bNAbs, as well as CDRL1 deletions and/or glycine substitutions to accommodate the N276 glycan. These results provide proof of concept for vaccine-induced affinity maturation of B cell lineages that require rare insertions and deletions.

Exploring multivalent carbohydrate-protein interactions by NMR

Nuclear Magnetic Resonance (NMR) has been widely employed to assess diverse features of glycan-protein molecular recognition events. Different types of qualitative and quantitative information at different degrees of resolution and complexity can be extracted from the proper application of the available NMR-techniques. In fact, affinity, structural, kinetic, conformational, and dynamic characteristics of the binding process are available. Nevertheless, except in particular cases, the affinity of lectin-sugar interactions is weak, mostly at the low mM range. This feature is overcome in biological processes by using multivalency, thus augmenting the strength of the binding. However, the application of NMR methods to monitor multivalent lectin-glycan interactions is intrinsically challenging. It is well known that when large macromolecular complexes are formed, the NMR signals disappear from the NMR spectrum, due to the existence of fast transverse relaxation, related to the large size and exchange features. Indeed, at the heart of the molecular recognition event, the associated free-bound chemical exchange process for both partners takes place in a particular timescale. Thus, these factors have to be considered and overcome. In this review article, we have distinguished, in a subjective manner, the existence of multivalent presentations in the glycan or in the lectin. From the glycan perspective, we have also considered whether multiple epitopes of a given ligand are presented in the same linear chain of a saccharide (i.e., poly-LacNAc oligosaccharides) or decorating different arms of a multiantennae scaffold, either natural (as in multiantennae N-glycans) or synthetic (of dendrimer or polymer nature). From the lectin perspective, the presence of an individual binding site at every monomer of a multimeric lectin may also have key consequences for the binding event at different levels of complexity.

The retaining β-Kdo glycosyltransferase WbbB uses a double-displacement mechanism with an intermediate adduct rearrangement step

WbbB, a lipopolysaccharide O-antigen synthesis enzyme from Raoultella terrigena, contains an N-terminal glycosyltransferase domain with a highly modified architecture that adds a terminal β-Kdo (3-deoxy-D-manno-oct-2-ulosonic acid) residue to the O-antigen saccharide, with retention of stereochemistry. We show, using mass spectrometry, that WbbB forms a covalent adduct between the catalytic nucleophile, Asp232, and Kdo. We also determine X-ray structures for the CMP-β-Kdo donor complex, for Kdo-adducts with D232N and D232C WbbB variants, for a synthetic disaccharide acceptor complex, and for a ternary complex with both a Kdo-adduct and the acceptor. Together, these structures show that the enzyme-linked Asp232-Kdo adduct rotates to reposition the Kdo into a second sub-site, which then transfers Kdo to the acceptor. Retaining glycosyltransferases were thought to use only the front-side S_Ni substitution mechanism; here we show that retaining glycosyltransferases can also potentially use double-displacement mechanisms, but incorporating an additional catalytic subsite requires rearrangement of the protein's architecture.

Multifaceted Computational Modeling in Glycoscience

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

A database of calculated solution parameters for the AlphaFold predicted protein structures

Recent spectacular advances by AI programs in 3D structure predictions from protein sequences have revolutionized the field in terms of accuracy and speed. The resulting “folding frenzy” has already produced predicted protein structure databases for the entire human and other organisms’ proteomes. However, rapidly ascertaining a predicted structure’s reliability based on measured properties in solution should be considered. Shape-sensitive hydrodynamic parameters such as the diffusion and sedimentation coefficients (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D_{t(20,w)}^{0}}$$\end{document}Dt(20,w)0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${s_{{\left( {{20},w} \right)}}^{{0}} }$$\end{document}s20,w0) and the intrinsic viscosity ([η]) can provide a rapid assessment of the overall structure likeliness, and SAXS would yield the structure-related pair-wise distance distribution function p(r) vs. r. Using the extensively validated UltraScan SOlution MOdeler (US-SOMO) suite, a database was implemented calculating from AlphaFold structures the corresponding \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${D_{t(20,w)}^{0}}$$\end{document}Dt(20,w)0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${s_{{\left( {{20},w} \right)}}^{{0}} }$$\end{document}s20,w0, [η], p(r) vs. r, and other parameters. Circular dichroism spectra were computed using the SESCA program. Some of AlphaFold’s drawbacks were mitigated, such as generating whenever possible a protein’s mature form. Others, like the AlphaFold direct applicability to single-chain structures only, the absence of prosthetic groups, or flexibility issues, are discussed. Overall, this implementation of the US-SOMO-AF database should already aid in rapidly evaluating the consistency in solution of a relevant portion of AlphaFold predicted protein structures.

Updated restraint dictionaries for carbohydrates in the pyranose form

Restraint dictionaries are used during macromolecular structure refinement to encapsulate intramolecular connectivity and geometric information. These dictionaries allow previously determined `ideal' values of features such as bond lengths, angles and torsions to be used as restraint targets. During refinement, restraints influence the model to adopt a conformation that agrees with prior observation. This is especially important when refining crystal structures of glycosylated proteins, as their resolutions tend to be worse than those of nonglycosylated proteins. Pyranosides, the overwhelming majority component in all forms of protein glycosylation, often display conformational errors in crystal structures. Whilst many of these flaws usually relate to model building, refinement issues may also have their root in suboptimal restraint dictionaries. In order to avoid subsequent misinterpretation and to improve the quality of all pyranose monosaccharide entries in the CCP4 Monomer Library, new dictionaries with improved ring torsion restraints, coordinates reflecting the lowest-energy ring pucker and updated geometry have been produced and evaluated. These new dictionaries are now part of the CCP4 Monomer Library and will be released with CCP4 version 8.0.

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.

C-Glycosylcrotylboronates for the Synthesis of Glycomimetics

The stereoselective synthesis of E- and Z- isomers of a C- mannosyl crotylpinacolboronate via Ni-promoted reactions on an allylic acetate and a diene precursor, respectively, is described. The E- and Z- isomers reacted with 1,2-O-isopropylidene glyceraldehyde in the presence or absence of (R)- and (S)- TRIP catalysts, to give predominantly 3,4-anti and 3,4-syn crotylation products, respectively, with moderate to high facial selectivity. These products were transformed to biologically relevant C-manno-disaccharides.

Molecular Recognition Insights of Sialic Acid Glycans by Distinct Receptors Unveiled by NMR and Molecular Modeling

All cells are decorated with a highly dense and complex structure of glycan chains, which are mostly attached to proteins and lipids. In this context, sialic acids are a family of nine-carbon acidic monosaccharides typically found at the terminal position of glycan chains, modulating several physiological and pathological processes. Sialic acids have many structural and modulatory roles due to their negative charge and hydrophilicity. In addition, the recognition of sialic acid glycans by mammalian cell lectins, such as siglecs, has been described as an important immunological checkpoint. Furthermore, sialic acid glycans also play a pivotal role in host-pathogen interactions. Various pathogen receptors exposed on the surface of viruses and bacteria are responsible for the binding to sialic acid sugars located on the surface of host cells, becoming a critical point of contact in the infection process. Understanding the molecular mechanism of sialic acid glycans recognition by sialic acid-binding proteins, present on the surface of pathogens or human cells, is essential to realize the biological mechanism of these events and paves the way for the rational development of strategies to modulate sialic acid-protein interactions in diseases. In this perspective, nuclear magnetic resonance (NMR) spectroscopy, assisted with molecular modeling protocols, is a versatile and powerful technique to investigate the structural and dynamic aspects of glycoconjugates and their interactions in solution at the atomic level. NMR provides the corresponding ligand and protein epitopes, essential for designing and developing potential glycan-based therapies. In this review, we critically discuss the current state of knowledge about the structural features behind the molecular recognition of sialic acid glycans by different receptors, naturally present on human cells or pathogens, disclosed by NMR spectroscopy and molecular modeling protocols.

Modernized uniform representation of carbohydrate molecules in the Protein Data Bank

Since 1971, the Protein Data Bank (PDB) has served as the single global archive for experimentally determined 3D structures of biological macromolecules made freely available to the global community according to the FAIR principles of Findability-Accessibility-Interoperability-Reusability. During the first 50 years of continuous PDB operations, standards for data representation have evolved to better represent rich and complex biological phenomena. Carbohydrate molecules present in more than 14,000 PDB structures have recently been reviewed and remediated to conform to a new standardized format. This machine-readable data representation for carbohydrates occurring in the PDB structures and the corresponding reference data improves the findability, accessibility, interoperability and reusability of structural information pertaining to these molecules. The PDB Exchange MacroMolecular Crystallographic Information File data dictionary now supports (i) standardized atom nomenclature that conforms to International Union of Pure and Applied Chemistry-International Union of Biochemistry and Molecular Biology (IUPAC-IUBMB) recommendations for carbohydrates, (ii) uniform representation of branched entities for oligosaccharides, (iii) commonly used linear descriptors of carbohydrates developed by the glycoscience community and (iv) annotation of glycosylation sites in proteins. For the first time, carbohydrates in PDB structures are consistently represented as collections of standardized monosaccharides, which precisely describe oligosaccharide structures and enable improved carbohydrate visualization, structure validation, robust quantitative and qualitative analyses, search for dendritic structures and classification. The uniform representation of carbohydrate molecules in the PDB described herein will facilitate broader usage of the resource by the glycoscience community and researchers studying glycoproteins.

Towards consistency in geometry restraints for carbohydrates in the pyranose form: modern dictionary generators reviewed

Macromolecular restrained refinement is nowadays the most used method for improving the agreement between an atomic structural model and experimental data. Restraint dictionaries, a key tool behind the success of the method, allow fine-tuning geometric properties such as distances and angles between atoms beyond simplistic expectations. Dictionary generators can provide restraint target estimates derived from different sources, from fully theoretical to experimental and any combination in between. Carbohydrates are stereochemically complex biomolecules and, in their pyranose form, have clear conformational preferences. As such, they pose unique problems to dictionary generators and in the course of this study, require special attention from software developers. Functional differences between restraint generators will be discussed, as well as the process of achieving consistent results with different software designs. The study will conclude a set of practical considerations, as well as recommendations for the generation of new restraint dictionaries, using the improved software alternatives discussed.

Neutralizing Antibodies Induced by First-Generation gp41-Stabilized HIV-1 Envelope Trimers and Nanoparticles

The immunogenicity of gp41-stabilized HIV-1 BG505 envelope (Env) trimers and nanoparticles (NPs) was recently assessed in mice and rabbits. Here, we combined Env-specific B-cell sorting and repertoire sequencing to identify neutralizing antibodies (NAbs) from immunized animals. A panel of mouse NAbs was isolated from mice immunized with a 60-meric I3-01 NP presenting 20 stabilized trimers. Three mouse NAbs potently neutralized BG505.T332N by recognizing a glycan epitope centered in the C3/V4 region on BG505 Env, as revealed by electron microscopy (EM), X-ray crystallography, and epitope mapping. A set of rabbit NAbs was isolated from rabbits immunized with a soluble trimer and a 24-meric ferritin NP presenting 8 trimers. Neutralization assays against BG505.T332N variants confirmed that potent rabbit NAbs targeted previously described glycan holes on BG505 Env and accounted for a significant portion of the autologous NAb response in both the trimer and ferritin NP groups. Last, we examined NAb responses that were induced by non-BG505 Env immunogens. We determined a 3.4-Å-resolution crystal structure for the clade C transmitted/founder (T/F) Du172.17 Env with a redesigned heptad repeat 1 (HR1) bend in gp41. This clade C Env, in a soluble trimer form and in a multivalent form with 8 trimers attached to ferritin NP, and the gp41-stabilized clade A Q482-d12 Env trimer elicited distinct NAb responses in rabbits, with notable differences in neutralization breadth. Although eliciting a broad NAb response remains a major challenge, our study provides valuable information on an HIV-1 vaccine design strategy that combines gp41 stabilization and NP display. IMPORTANCE Self-assembling protein nanoparticles (NPs) presenting BG505 envelope (Env) trimers can elicit tier 2 HIV-1-neutralizing antibody (NAb) responses more effectively than soluble trimers. In the present study, monoclonal NAbs were isolated from previously immunized mice and rabbits for structural and functional analyses, which revealed that potent mouse NAbs recognize the C3/V4 region and small NP-elicited rabbit NAbs primarily target known glycan holes on BG505 Env. This study validates the gp41 stabilization strategy for HIV-1 Env vaccine design and highlights the challenge in eliciting a broad NAb response.

Active site dynamics and catalytic mechanism in arabinan hydrolysis catalyzed by GH43 endo-arabinanase from QM/MM molecular dynamics simulation and potential energy surface

The endo-1,5-α-L-arabinanases, belonging to glycoside hydrolase family 43 (GH43), catalyse the hydrolysis of α-1,5-arabinofuranosidic bonds in arabinose-containing polysaccharides. These enzymes are proposed targets for industrial and medical applications. Here, molecular dynamics (MD), potential energy surface and free energy (potential of mean force) simulations are undertaken using hybrid quantum mechanical/molecular mechanical (QM/MM) potentials to understand the active site dynamics, catalytic mechanism and the electrostatic influence of active site residues of the GH43 endo-arabinanase from G. stearothermophilus. The calculated results give support to the single-displacement mechanism proposed for the inverting GH43 enzymes: first a proton is transferred from the general acid E201 to the substrate, followed by a nucleophilic attack by water, activated by the general base D27, on the anomer carbon. A conformational change (²E ↔E₃ ↔ ⁴E) in the -1 sugar ring is observed involving a transition state featuring an oxocarbenium ion character. Residues D87, K106, H271 are highlighted as potential targets for future mutation experiments in order to increase the efficiency of the reaction. To our knowledge, this is the first QM/MM study providing molecular insights into the glycosidic bond hydrolysis of a furanoside substrate by an inverting GH in solution.Communicated by Ramaswamy H. Sarma.

Relation between glycosidic linkage, structure and dynamics of α- and β-glucans in water

In a molecular dynamics simulation study of several oligosaccharides comprising of the very basic building block of carbohydrate, the α- or β-d glucopyranose units, linked by any one of the 1-3/1-4 or 1-6 glycosidic linkages, we compare and contrast their structural and dynamical properties. Results indicate that the litheness of the oligosaccharide chain is noticeably controlled by the composition, anomeric nature and glycosidic linkage type of the units. In mixed β 1-4/1-3 d-glucopyranosides, as those found in oats and barley, the ratio of the β 1-4 and β 1-3 linked residues is crucial in determining the structural and dynamical attributes. Principal component analysis (PCA) using the internal coordinates of torsion angles subtended by glycosidic oxygen atoms and subsequent K-means clustering of the dynamical space spanned by PC1 to PC2 point to the dynamical and structural disparity in the various types of oligosaccharides studied. The properties simulated in this work are meant to provide a systematic yet comparative understanding of the importance of linkage and anomericity on the oligosaccharide chain properties and are in line with some experimental structural attributes.

Structural basis for the design of bisubstrate inhibitors of protein kinase CK2 provided by complex structures with the substrate-competitive inhibitor heparin

The Ser/Thr kinase CK2, a member of the superfamily of eukaryotic protein kinases, has an acidophilic substrate profile with the substrate recognition sequence S/T-D/E-X-D/E, and it is inhibited by polyanionic substances like heparin. The latter, a highly sulphated glucosamino glycan composed mainly of repeating 2-O-sulpho-α-l-idopyranuronic acid/N,O6-disulpho-α-d-glucosamine disaccharide units, is the longest known substrate-competitive CK2 inhibitor. The structural basis of CK2's preference for anionic substrates and substrate-competitive inhibitors is only vaguely known which limits the value of the substrate-binding region for the structure-based development of CK2 bisubstrate inhibitors. Here, a tetragonal and a monoclinic co-crystal structure of CK2α, the catalytic subunit of CK2, with a decameric heparin fragment are described. In the tetragonal structure, the heparin molecule binds to the polybasic stretch at the beginning of CK2α's helix αC, whereas in the monoclinic structure it occupies the central substrate-recognition region around the P+1 loop. Together, the structures rationalize the inhibitory efficacy of heparin fragments as a function of chain length. The monoclinic CK2α/heparin structure, in which the heparin fragment is particularly well defined, is the first CK2 structure with an anionic inhibitor of considerable size at the central part of the substrate-recognition site. The bound heparin fragment is so close to the binding site of ATP-competitive inhibitors that it can guide the design of linkers and pave the way to efficient CK2 bisubstrate inhibitors in the future.

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.

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.

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.

Leveraging glycomics data in glycoprotein 3D structure validation with Privateer

The heterogeneity, mobility and complexity of glycans in glycoproteins have been, and currently remain, significant challenges in structural biology. These aspects present unique problems to the two most prolific techniques: X-ray crystallography and cryo-electron microscopy. At the same time, advances in mass spectrometry have made it possible to get deeper insights on precisely the information that is most difficult to recover by structure solution methods: the full-length glycan composition, including linkage details for the glycosidic bonds. The developments have given rise to glycomics. Thankfully, several large scale glycomics initiatives have stored results in publicly available databases, some of which can be accessed through API interfaces. In the present work, we will describe how the Privateer carbohydrate structure validation software has been extended to harness results from glycomics projects, and its use to greatly improve the validation of 3D glycoprotein structures.

A VH1-69 antibody lineage from an infected Chinese donor potently neutralizes HIV-1 by targeting the V3 glycan supersite

An oligomannose patch around the V3 base of HIV-1 envelope glycoprotein (Env) is recognized by multiple classes of broadly neutralizing antibodies (bNAbs). Here, we investigated the bNAb response to the V3 glycan supersite in an HIV-1-infected Chinese donor by Env-specific single B cell sorting, structural and functional studies, and longitudinal analysis of antibody and virus repertoires. Monoclonal antibodies 438-B11 and 438-D5 were isolated that potently neutralize HIV-1 with moderate breadth, are encoded by the VH1-69 germline gene, and have a disulfide-linked long HCDR3 loop. Crystal structures of Env-bound and unbound antibodies revealed heavy chain-mediated recognition of the glycan supersite with a unique angle of approach and a critical role of the intra-HCDR3 disulfide. The mechanism of viral escape was examined via single-genome amplification/sequencing and glycan mutations around the N332 supersite. Our findings further emphasize the V3 glycan supersite as a prominent target for Env-based vaccine design.

Current developments in Coot for macromolecular model building of Electron Cryo-microscopy and Crystallographic Data

Coot is a tool widely used for model building, refinement, and validation of macromolecular structures. It has been extensively used for crystallography and, more recently, improvements have been introduced to aid in cryo‐EM model building and refinement, as cryo‐EM structures with resolution ranging 2.5–4 A are now routinely available. Model building into these maps can be time‐consuming and requires experience in both biochemistry and building into low‐resolution maps. To simplify and expedite the model building task, and minimize the needed expertise, new tools are being added in Coot. Some examples include morphing, Geman‐McClure restraints, full‐chain refinement, and Fourier‐model based residue‐type‐specific Ramachandran restraints. Here, we present the current state‐of‐the‐art in Coot usage.

Crystal Structure of a GH3 β-Glucosidase from the Thermophilic Fungus Chaetomium thermophilum

Beta-glucosidases (β-glucosidases) have attracted considerable attention in recent years for use in various biotechnological applications. They are also essential enzymes for lignocellulose degradation in biofuel production. However, cost-effective biomass conversion requires the use of highly efficient enzymes. Thus, the search for new enzymes as better alternatives of the currently available enzyme preparations is highly important. Thermophilic fungi are nowadays considered as a promising source of enzymes with improved stability. Here, the crystal structure of a family GH3 β-glucosidase from the thermophilic fungus Chaetomium thermophilum (CtBGL) was determined at a resolution of 2.99 Å. The structure showed the three-domain architecture found in other β-glucosidases with variations in loops and linker regions. The active site catalytic residues in CtBGL were identified as Asp287 (nucleophile) and Glu517 (acid/base). Structural comparison of CtBGL with Protein Data Bank (PDB)-deposited structures revealed variations among glycosylated Asn residues. The enzyme displayed moderate glycosylation compared to other GH3 family β-glucosidases with similar structure. A new glycosylation site at position Asn504 was identified in CtBGL. Moreover, comparison with respect to several thermostability parameters suggested that glycosylation and charged residues involved in electrostatic interactions may contribute to the stability of the enzyme at elevated temperatures. The reported CtBGL structure provides additional insights into the family GH3 enzymes and could offer new ideas for further improvements in β-glucosidases for more efficient use in biotechnological applications regarding cellulose degradation.

The molecular basis of endolytic activity of a multidomain alginate lyase from Defluviitalea phaphyphila, a representative of a new lyase family, PL39

Alginate is a polymer containing two uronic acid epimers, β-d-mannuronate (M) and α-l-guluronate (G), and is a major component of brown seaweed that is depolymerized by alginate lyases. These enzymes have diverse specificity, cleaving the chain with endo- or exotype activity and with differential selectivity for the sequence of M or G at the cleavage site. Dp0100 is a 201-kDa multimodular, broad-specificity endotype alginate lyase from the marine thermophile Defluviitalea phaphyphila, which uses brown algae as a carbon source, converting it to ethanol, and bioinformatics analysis suggested that its catalytic domain represents a new polysaccharide lyase family, PL39. The structure of the Dp0100 catalytic domain, determined at 2.07 Å resolution, revealed that it comprises three regions strongly resembling those of the exotype lyase families PL15 and PL17. The conservation of key catalytic histidine and tyrosine residues belonging to the latter suggests these enzymes share mechanistic similarities. A complex of Dp0100 with a pentasaccharide, M₅, showed that the oligosaccharide is located in subsites -2, -1, +1, +2, and +3 in a long, deep canyon open at both ends, explaining the endotype activity of this lyase. This contrasted with the hindered binding sites of the exotype enzymes, which are blocked such that only one sugar moiety can be accommodated at the -1 position in the catalytic site. The biochemical and structural analyses of Dp0100, the first for this new class of endotype alginate lyases, have furthered our understanding of the structure-function and evolutionary relationships within this important class of enzymes.

Novel NMR Avenues to Explore the Conformation and Interactions of Glycans

This perspective article is focused on the presentation of the latest advances in NMR methods and applications that are behind the exciting achievements in the understanding of glycan receptors in molecular recognition events. Different NMR-based methodologies are discussed along with their applications to scrutinize the conformation and dynamics of glycans as well as their interactions with protein receptors.

A Million Crystal Structures: The Whole Is Greater than the Sum of Its Parts

The founding in 1965 of what is now called the Cambridge Structural Database (CSD) has reaped dividends in numerous and diverse areas of chemical research. Each of the million or so crystal structures in the database was solved for its own particular reason, but collected together, the structures can be reused to address a multitude of new problems. In this Review, which is focused mainly on the last 10 years, we chronicle the contribution of the CSD to research into molecular geometries, molecular interactions, and molecular assemblies and demonstrate its value in the design of biologically active molecules and the solid forms in which they are delivered. Its potential in other commercially relevant areas is described, including gas storage and delivery, thin films, and (opto)electronics. The CSD also aids the solution of new crystal structures. Because no scientific instrument is without shortcomings, the limitations of CSD research are assessed. We emphasize the importance of maintaining database quality: notwithstanding the arrival of big data and machine learning, it remains perilous to ignore the principle of garbage in, garbage out. Finally, we explain why the CSD must evolve with the world around it to ensure it remains fit for purpose in the years ahead.

Building and rebuilding N-glycans in protein structure models

N-Glycosylation is one of the most common post-translational modifications and is implicated in, for example, protein folding and interaction with ligands and receptors. N-Glycosylation trees are complex structures of linked carbohydrate residues attached to asparagine residues. While carbohydrates are typically modeled in protein structures, they are often incomplete or have the wrong chemistry. Here, new tools are presented to automatically rebuild existing glycosylation trees, to extend them where possible, and to add new glycosylation trees if they are missing from the model. The method has been incorporated in the PDB-REDO pipeline and has been applied to build or rebuild 16 452 carbohydrate residues in 11 651 glycosylation trees in 4498 structure models, and is also available from the PDB-REDO web server. With better modeling of N-glycosylation, the biological function of this important modification can be better and more easily understood.

Capturing the inherent structural dynamics of the HIV-1 envelope glycoprotein fusion peptide

The N-terminal fusion peptide (FP) of the human immunodeficiency virus (HIV)-1 envelope glycoprotein (Env) gp41 subunit plays a critical role in cell entry. However, capturing the structural flexibility in the unbound FP is challenging in the native Env trimer. Here, FP conformational isomerism is observed in two crystal structures of a soluble clade B transmitted/founder virus B41 SOSIP.664 Env with broadly neutralizing antibodies (bNAbs) PGT124 and 35O22 to aid in crystallization and that are not specific for binding to the FP. Large rearrangements in the FP and fusion peptide proximal region occur around M530, which remains anchored in the tryptophan clasp (gp41 W623, W628, W631) in the B41 Env prefusion state. Further, we redesigned the FP at position 518 to reinstate the bNAb VRC34.01 epitope. These findings provide further structural evidence for the dynamic nature of the FP and how a bNAb epitope can be restored during vaccine design.

The fusion peptide (FP) of HIV envelope (Env) is critical in the cell entry process. Here, Kumar et al. present crystal structures of B41 SOSIP.664 Env trimer and show the dynamic nature of the FP and proximal region, which likely relates to conformational rearrangements required for membrane fusion.

Automatically Fixing Errors in Glycoprotein Structures with Rosetta

Recent advances in single-particle cryo-electron microscopy (cryoEM) have resulted in determination of an increasing number of protein structures with resolved glycans. However, existing protocols for the refinement of glycoproteins at low resolution have failed to keep up with these advances. As a result, numerous deposited structures contain glycan stereochemical errors. Here, we describe a Rosetta-based approach for both cryoEM and X-ray crystallography refinement of glycoproteins that is capable of correcting conformational and configurational errors in carbohydrates. Building upon a previous Rosetta framework, we introduced additional features and score terms enabling automatic detection, setup, and refinement of glycan-containing structures. We benchmarked this approach using 12 crystal structures and showed that glycan geometries can be automatically improved while maintaining good fit to the crystallographic data. Finally, we used this method to refine carbohydrates of the human coronavirus NL63 spike glycoprotein and of an HIV envelope glycoprotein, demonstrating its usefulness for cryoEM refinement.

•

New method for refinement of carbohydrates with low-resolution electron density

•

Improved physical geometry of glycans in protein structures

•

Compatible with cryoEM and X-ray crystallography data

The recognition of glycans by protein receptors. Insights from NMR spectroscopy

Carbohydrates (glycans, saccharides, sugars) are everywhere. In fact, glycan-protein interactions are involved in many essential processes of life and disease. The understanding of the key structural details at the atomic and molecular level is of paramount importance to effectively design molecules for therapeutic purposes. Different approximations may be employed to decipher these molecular recognition processes with high resolution. Advances in cryo-electron microscopy are providing exquisite details on different biological mechanisms involving sugars, while better and better protocols for structural refinement in the application of X-ray methods for protein-sugar complexes and glycoproteins are also permitting fantastic advances in the glycoscience arena. Alternatively, NMR spectroscopy remains as one of the most rewarding techniques to explore protein-carbohydrate interactions. In fact, given the intrinsic dynamic nature of saccharides, NMR can afford exquisite structural information at the atomic detail, not accessible by other techniques. However, the access to this information is sometimes intricate, and requires careful analysis and well-defined strategies. In this review, we have highlighted these issues and presented an overview of different modern NMR approaches with a focus on the latest developments and challenges.

Crystal structure of highly glycosylated human leukocyte elastase in complex with an S2' site binding inhibitor

Glycosylated human leukocyte elastase (HLE) was crystallized and structurally analysed in complex with a 1,3-thiazolidine-2,4-dione derivative that had been identified as an HLE inhibitor in preliminary studies. In contrast to previously described HLE structures with small-molecule inhibitors, in this structure the inhibitor does not bind to the S1 and S2 substrate-recognition sites; rather, this is the first HLE structure with a synthetic inhibitor in which the S2' site is blocked that normally binds the second side chain at the C-terminal side of the scissile peptide bond in a substrate protein. The inhibitor also induces the formation of crystalline HLE dimers that block access to the active sites and that are also predicted to be stable in solution. Neither such HLE dimers nor the corresponding crystal packing have been observed in previous HLE crystal structures. This novel crystalline environment contributes to the observation that comparatively large parts of the N-glycan chains of HLE are defined by electron density. The final HLE structure contains the largest structurally defined carbohydrate trees among currently available HLE structures.

Carbohydrate structure hits the groove

An introduction to the Acta Cryst. F special issue on glycoproteins and protein–carbohydrate complexes in which the contents, the current state of the field and the future of glycan structural biology are briefly discussed.

Structural analysis of glycoproteins: building N-linked glycans with Coot

Coot is a graphics application that is used to build or manipulate macromolecular models; its particular forte is manipulation of the model at the residue level. The model-building tools of Coot have been combined and extended to assist or automate the building of N-linked glycans. The model is built by the addition of monosaccharides, placed by variation of internal coordinates. The subsequent model is refined by real-space refinement, which is stabilized with modified and additional restraints. It is hoped that these enhanced building tools will help to reduce building errors of N-linked glycans and improve our knowledge of the structures of glycoproteins.

X-ray and EM structures of a natively glycosylated HIV-1 envelope trimer

The structural and biochemical characterization of broadly neutralizing anti-HIV-1 antibodies (bNAbs) has been essential in guiding the design of potential vaccines to prevent infection by HIV-1. While these studies have revealed critical mechanisms by which bNAbs recognize and/or accommodate N-glycans on the trimeric envelope glycoprotein (Env), they have been limited to the visualization of high-mannose glycan forms only, since heterogeneity introduced from the presence of complex glycans makes it difficult to obtain high-resolution structures. 3.5 and 3.9 Å resolution crystal structures of the HIV-1 Env trimer with fully processed and native glycosylation were solved, revealing a glycan shield of high-mannose and complex-type N-glycans that were used to define the complete epitopes of two bNAbs. Here, the refinement of the N-glycans in the crystal structures is discussed and comparisons are made with glycan densities in glycosylated Env structures derived by single-particle cryo-electron microscopy.

Keep it together: restraints in crystallographic refinement of macromolecule-ligand complexes

A short introduction is provided to the concept of restraints in macromolecular crystallographic refinement. A typical ligand restraint-generation process is then described, covering types of input, the methodology and the mechanics behind the software in general terms, how this has evolved over recent years and what to look for in the output. Finally, the currently available restraint-generation software is compared, concluding with some thoughts for the future.

Ligand fitting with CCP4

Crystal structures of protein-ligand complexes are often used to infer biology and inform structure-based drug discovery. Hence, it is important to build accurate, reliable models of ligands that give confidence in the interpretation of the respective protein-ligand complex. This paper discusses key stages in the ligand-fitting process, including ligand binding-site identification, ligand description and conformer generation, ligand fitting, refinement and subsequent validation. The CCP4 suite contains a number of software tools that facilitate this task: AceDRG for the creation of ligand descriptions and conformers, Lidia and JLigand for two-dimensional and three-dimensional ligand editing and visual analysis, Coot for density interpretation, ligand fitting, analysis and validation, and REFMAC5 for macromolecular refinement. In addition to recent advancements in automatic carbohydrate building in Coot (LO/Carb) and ligand-validation tools (FLEV), the release of the CCP4i2 GUI provides an integrated solution that streamlines the ligand-fitting workflow, seamlessly passing results from one program to the next. The ligand-fitting process is illustrated using instructive practical examples, including problematic cases such as post-translational modifications, highlighting the need for careful analysis and rigorous validation.