Computational glycoscience: characterizing the spatial and temporal properties of glycans and glycan-protein complexes

Conformational dynamics of α-1 acid glycoprotein (AGP) in cancer: A comparative study of glycosylated and unglycosylated AGP

α-1 acid glycoprotein (AGP) is one of the most abundant plasma proteins. It fulfills two important functions: immunomodulation, and binding to various drugs and receptors. These different functions are closely associated and modulated via changes in glycosylation and cancer missense mutations. From a structural point of view, glycans alter the local biophysical properties of the protein leading to a diverse ligand-binding spectrum. However, glycans can typically not be observed in the resolved X-ray crystallography structure of AGP due to their high flexibility and microheterogeneity, so limiting our understanding of AGP's conformational dynamics 70 years after its discovery. We here investigate how mutations and glycosylation interfere with AGP's conformational dynamics changing its biophysical behavior, by using molecular dynamics (MD) simulations and sequence-based dynamics predictions. The MD trajectories show that glycosylation decreases the local backbone flexibility of AGP and increases the flexibility of distant regions through allosteric effects. We observe that mutations near the glycosylation site affect glycan's conformational preferences. Thus, we conclude that mutations control glycan dynamics which modulates the protein's backbone flexibility directly affecting its accessibility. These findings may assist in the drug design targeting AGP's glycosylation and mutations in cancer.

Structural Analysis of Oligosaccharides and Glycoconjugates Using NMR

Carbohydrate chains play critical roles in cellular recognition and subsequent signal transduction in the nervous system. Furthermore, gangliosides are targets for various amyloidogenic proteins associated with neurodegenerative disorders. To better understand the molecular mechanisms underlying these biological phenomena, atomic views are essential to delineate dynamic biomolecular interactions. Nuclear magnetic resonance (NMR) spectroscopy provides powerful tools for studying structures, dynamics, and interactions of biomolecules at the atomic level. This chapter describes the basics of solution NMR techniques and their applications to the analysis of 3D structures and interactions of glycoconjugates in the nervous system.

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.

Glycans in Virus-Host Interactions: A Structural Perspective

Many interactions between microbes and their hosts are driven or influenced by glycans, whose heterogeneous and difficult to characterize structures have led to an underappreciation of their role in these interactions compared to protein-based interactions. Glycans decorate microbe glycoproteins to enhance attachment and fusion to host cells, provide stability, and evade the host immune system. Yet, the host immune system may also target these glycans as glycoepitopes. In this review, we provide a structural perspective on the role of glycans in host-microbe interactions, focusing primarily on viral glycoproteins and their interactions with host adaptive immunity. In particular, we discuss a class of topological glycoepitopes and their interactions with topological mAbs, using the anti-HIV mAb 2G12 as the archetypical example. We further offer our view that structure-based glycan targeting strategies are ready for application to viruses beyond HIV, and present our perspective on future development in this area.

Computational prediction method to decipher receptor-glycoligand interactions in plant immunity

Microbial and plant cell walls have been selected by the plant immune system as a source of microbe- and plant damage-associated molecular patterns (MAMPs/DAMPs) that are perceived by extracellular ectodomains (ECDs) of plant pattern recognition receptors (PRRs) triggering immune responses. From the vast number of ligands that PRRs can bind, those composed of carbohydrate moieties are poorly studied, and only a handful of PRR/glycan pairs have been determined. Here we present a computational screening method, based on the first step of molecular dynamics simulation, that is able to predict putative ECD-PRR/glycan interactions. This method has been developed and optimized with Arabidopsis LysM-PRR members CERK1 and LYK4, which are involved in the perception of fungal MAMPs, chitohexaose (1,4-β-d-(GlcNAc)6 ) and laminarihexaose (1,3-β-d-(Glc)6 ). Our in silico results predicted CERK1 interactions with 1,4-β-d-(GlcNAc)6 whilst discarding its direct binding by LYK4. In contrast, no direct interaction between CERK1/laminarihexaose was predicted by the model despite CERK1 being required for laminarihexaose immune activation, suggesting that CERK1 may act as a co-receptor for its recognition. These in silico results were validated by isothermal titration calorimetry binding assays between these MAMPs and recombinant ECDs-LysM-PRRs. The robustness of the developed computational screening method was further validated by predicting that CERK1 does not bind the DAMP 1,4-β-d-(Glc)6 (cellohexaose), and then probing that immune responses triggered by this DAMP were not impaired in the Arabidopsis cerk1 mutant. The computational predictive glycan/PRR binding method developed here might accelerate the discovery of protein-glycan interactions and provide information on immune responses activated by glycoligands.

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.

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.

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.

In-depth structural analysis of glycans in the gas phase

Although there have been substantial improvements in glycan analysis over the past decade, the lack of both high-resolution and high-throughput methods hampers progress in glycomics. This perspective article highlights the current developments of liquid chromatography, mass spectrometry, ion-mobility spectrometry and cryogenic IR spectroscopy for glycan analysis and gives a critical insight to their individual strengths and limitations. Moreover, we discuss a novel concept in which ion mobility-mass spectrometry and cryogenic IR spectroscopy is combined in a single instrument such that datasets consisting of m/z, collision cross sections and IR fingerprints can be obtained. This multidimensional data will then be compared to a comprehensive reference library of intact glycans and their fragments to accurately identify unknown glycans on a high-throughput scale with minimal sample requirements. Due to the complementarity of the obtained information, this novel approach is highly diagnostic and also suitable for the identification of larger glycans; however, the workflow and instrumentation is straightforward enough to be implemented into a user-friendly setup.

Molecular Dynamics of Gangliosides

Computational methodologies have immense potential to delineate the dynamic conformations of glycoconjugates including gangliosides, thereby characterizing the conformational adaptability of their glycans upon interacting with various target proteins. Replica-exchange molecular dynamics simulations have been employed to effectively explore the vast conformational spaces of large, branched carbohydrate moieties. When experimentally validated using NMR, molecular simulations can provide dynamical views of molecular recognition events involving the ganglioside glycans.

Identification of Rare Lewis Oligosaccharide Conformers in Aqueous Solution Using Enhanced Sampling Molecular Dynamics

Determining the conformations accessible to carbohydrate ligands in aqueous solution is important for understanding their biological action. In this work, we evaluate the conformational free-energy surfaces of Lewis oligosaccharides in explicit aqueous solvent using a multidimensional variant of the swarm-enhanced sampling molecular dynamics (msesMD) method; we compare with multi-microsecond unbiased MD simulations, umbrella sampling, and accelerated MD approaches. For the sialyl Lewis A tetrasaccharide, msesMD simulations in aqueous solution predict conformer landscapes in general agreement with the other biased methods and with triplicate unbiased 10 μs trajectories; these simulations find a predominance of closed conformer and a range of low-occupancy open forms. The msesMD simulations also suggest closed-to-open transitions in the tetrasaccharide are facilitated by changes in ring puckering of its GlcNAc residue away from the ⁴C₁ form, in line with previous work. For sialyl Lewis X tetrasaccharide, msesMD simulations predict a minor population of an open form in solution corresponding to a rare lectin-bound pose observed crystallographically. Overall, from comparison with biased MD calculations, we find that triplicate 10 μs unbiased MD simulations may not be enough to fully sample glycan conformations in aqueous solution. However, the computational efficiency and intuitive approach of the msesMD method suggest potential for its application in glycomics as a tool for analysis of oligosaccharide conformation.

A molecular dynamics-based algorithm for evaluating the glycosaminoglycan mimicking potential of synthetic, homogenous, sulfated small molecules

Glycosaminoglycans (GAGs) are key natural biopolymers that exhibit a range of biological functions including growth and differentiation. Despite this multiplicity of function, natural GAG sequences have not yielded drugs because of problems of heterogeneity and synthesis. Recently, several homogenous non-saccharide glycosaminoglycan mimetics (NSGMs) have been reported as agents displaying major therapeutic promise. Yet, it remains unclear whether sulfated NSGMs structurally mimic sulfated GAGs. To address this, we developed a three-step molecular dynamics (MD)-based algorithm to compare sulfated NSGMs with GAGs. In the first step of this algorithm, parameters related to the range of conformations sampled by the two highly sulfated molecules as free entities in water were compared. The second step compared identity of binding site geometries and the final step evaluated comparable dynamics and interactions in the protein-bound state. Using a test case of interactions with fibroblast growth factor-related proteins, we show that this three-step algorithm effectively predicts the GAG structure mimicking property of NSGMs. Specifically, we show that two unique dimeric NSGMs mimic hexameric GAG sequences in the protein-bound state. In contrast, closely related monomeric and trimeric NSGMs do not mimic GAG in either the free or bound states. These results correspond well with the functional properties of NSGMs. The results show for the first time that appropriately designed sulfated NSGMs can be good structural mimetics of GAGs and the incorporation of a MD-based strategy at the NSGM library screening stage can identify promising mimetics of targeted GAG sequences.

Structural Analysis of the Glycosylated Intact HIV-1 gp120-b12 Antibody Complex Using Hydroxyl Radical Protein Footprinting

Glycoprotein gp120 is a surface antigen and virulence factor of human immunodeficiency virus 1. Broadly neutralizing antibodies (bNAbs) that react to gp120 from a variety of HIV isolates offer hope for the development of broadly effective immunogens for vaccination purposes, if the interactions between gp120 and bNAbs can be understood. From a structural perspective, gp120 is a particularly difficult system because of its size, the presence of multiple flexible regions, and the large amount of glycosylation, all of which are important in gp120–bNAb interactions. Here, the interaction of full-length, glycosylated gp120 with bNAb b12 is probed using high-resolution hydroxyl radical protein footprinting (HR-HRPF) by fast photochemical oxidation of proteins. HR-HRPF allows for the measurement of changes in the average solvent accessible surface area of multiple amino acids without the need for measures that might alter the protein conformation, such as mutagenesis. HR-HRPF of the gp120–b12 complex coupled with computational modeling shows a novel extensive interaction of the V1/V2 domain, probably with the light chain of b12. Our data also reveal HR-HRPF protection in the C3 domain caused by interaction of the N330 glycan with the b12 light chain. In addition to providing information about the interactions of full-length, glycosylated gp120 with b12, this work serves as a template for the structural interrogation of full-length glycosylated gp120 with other bNAbs to better characterize the interactions that drive the broad specificity of the bNAb.

Binding interaction between rice glutelin and amylose: Hydrophobic interaction and conformational changes

The interaction of rice glutelin (RG) with amylose was characterized by spectroscopic and molecular docking studies. The intrinsic fluorescence of RG increased upon the addition of amylose. The binding sites, binding constant and thermodynamic features indicated that binding process was spontaneous and the main driving force of the interaction was hydrophobic interaction. The surface hydrophobicity of RG decreased with increasing amount of amylose. Furthermore, synchronous fluorescence and circular dichroism (CD) spectra provided data concerning conformational and micro-environmental changes of RG. With the concentration of amylose increasing, the polarity around the tyrosine residues increased while the hydrophobicity decreased. Alteration of protein conformation was observed with increasing of α-helix and reducing of β-sheet. Finally, a visual representation of two binding sites located in the amorphous area of RG was presented by molecular modeling studies.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Revealing the Mechanisms of Protein Disorder and N-Glycosylation in CD44-Hyaluronan Binding Using Molecular Simulation

The extracellular N-terminal hyaluronan binding domain (HABD) of CD44 is a small globular domain that confers hyaluronan (HA) binding functionality to this large transmembrane glycoprotein. When recombinantly expressed by itself, HABD exists as a globular water-soluble protein that retains the capacity to bind HA. This has enabled atomic-resolution structural biology experiments that have revealed the structure of HABD and its binding mode with oligomeric HA. Such experiments have also pointed to an order-to-disorder transition in HABD that is associated with HA binding. However, it had remained unclear how this structural transition was involved in binding since it occurs in a region of HABD distant from the HA-binding site. Furthermore, HABD is known to be N-glycosylated, and such glycosylation can diminish HA binding when the associated N-glycans are capped with sialic acid residues. The intrinsic flexibility of disordered proteins and of N-glycans makes it difficult to apply experimental structural biology approaches to probe the molecular mechanisms of how the order-to-disorder transition and N-glycosylation can modulate HA binding by HABD. We review recent results from molecular dynamics simulations that provide atomic-resolution mechanistic understanding of such modulation to help bridge gaps between existing experimental binding and structural biology data. Findings from these simulations include: Tyr42 may function as a molecular switch that converts the HA-binding site from a low affinity to a high affinity state; in the partially disordered form of HABD, basic amino acids in the C-terminal region can gain sufficient mobility to form direct contacts with bound HA to further stabilize binding; and terminal sialic acids on covalently attached N-glycans can form charge-paired hydrogen bonding interactions with basic amino acids that could otherwise bind to HA, thereby blocking HA binding to glycosylated CD44 HABD.

All-Atom Internal Coordinate Mechanics (ICM) Force Field for Hexopyranoses and Glycoproteins

We present an extension of the all-atom internal-coordinate force field, ICMFF, that allows for simulation of heterogeneous systems including hexopyranose saccharides and glycan chains in addition to proteins. A library of standard glycan geometries containing α- and β-anomers of the most common hexapyranoses, i.e., d-galactose, d-glucose, d-mannose, d-xylose, l-fucose, N-acetylglucosamine, N-acetylgalactosamine, sialic, and glucuronic acids, is created based on the analysis of the saccharide structures reported in the Cambridge Structural Database. The new force field parameters include molecular electrostatic potential-derived partial atomic charges and the torsional parameters derived from quantum mechanical data for a collection of minimal molecular fragments and related molecules. The ϕ/ψ torsional parameters for different types of glycosidic linkages are developed using model compounds containing the key atoms in the full carbohydrates, i.e., glycosidic-linked tetrahydropyran-cyclohexane dimers. Target data for parameter optimization include two-dimensional energy surfaces corresponding to the ϕ/ψ glycosidic dihedral angles in the disaccharide analogues, as determined by quantum mechanical MP2/6-31G** single-point energies on HF/6-31G** optimized structures. To achieve better agreement with the observed geometries of glycosidic linkages, the bond angles at the O-linkage atoms are added to the internal variable set and the corresponding bond bending energy term is parametrized using quantum mechanical data. The resulting force field is validated on glycan chains of 1-12 residues from a set of high-resolution X-ray glycoprotein structures based on heavy atom root-mean-square deviations of the lowest-energy glycan conformations generated by the biased probability Monte Carlo (BPMC) molecular mechanics simulations from the native structures. The appropriate BPMC distributions for monosaccharide-monosaccharide and protein-glycan linkages are derived from the extensive analysis of conformational properties of glycoprotein structures reported in the Protein Data Bank. Use of the BPMC search leads to significant improvements in sampling efficiency for glycan simulations. Moreover, good agreement with the X-ray glycoprotein structures is achieved for all glycan chain lengths. Thus, average/median RMSDs are 0.81/0.68 Å for one-residue glycans and 1.32/1.47 Å for three-residue glycans. RMSD from the native structure for the lowest-energy conformation of the 12-residue glycan chain (PDB ID 3og2) is 1.53 Å. Additionally, results obtained for free short oligosaccharides using the new force field are in line with the available experimental data, i.e., the most populated conformations in solution are predicted to be the lowest energy ones. The newly developed parameters allow for the accurate modeling of linear and branched hexopyranose glycosides in heterogeneous systems.

Computational tools for epitope vaccine design and evaluation

Rational approaches will be required to develop universal vaccines for viral pathogens such as human immunodeficiency virus, hepatitis C virus, and influenza, for which empirical approaches have failed. The main objective of a rational vaccine strategy is to design novel immunogens that are capable of inducing long-term protective immunity. In practice, this requires structure-based engineering of the target neutralizing epitopes and a quantitative readout of vaccine-induced immune responses. Therefore, computational tools that can facilitate these two areas have played increasingly important roles in rational vaccine design in recent years. Here we review the computational techniques developed for protein structure prediction and antibody repertoire analysis, and demonstrate how they can be applied to the design and evaluation of epitope vaccines.

Crystal structure of a fully glycosylated HIV-1 gp120 core reveals a stabilizing role for the glycan at Asn262

The crystal structure of a fully glycosylated HIV-1 gp120 core in complex with CD4 receptor and Fab 17b at 4.5-Å resolution reveals 9 of the 15 N-linked glycans of core gp120 to be partially ordered. The glycan at position Asn262 had the most extensive and well-ordered electron density, and a GlcNAc(2)Man(7) was modeled. The GlcNAc stem of this glycan is largely buried in a cleft in gp120, suggesting a role in gp120 folding and stability. Its arms interact with the stems of neighboring glycans from the oligomannose patch, which is a major target for broadly neutralizing antibodies.

Calculating binding free energies for protein-carbohydrate complexes

A variety of computational techniques may be applied to compute theoretical binding free energies for protein-carbohydrate complexes. Elucidation of the intermolecular interactions, as well as the thermodynamic effects, that contribute to the relative strength of receptor binding can shed light on biomolecular recognition, and the resulting initiation or inhibition of a biological process. Three types of free energy methods are discussed here, including MM-PB/GBSA, thermodynamic integration, and a non-equilibrium alternative utilizing SMD. Throughout this chapter, the well-known concanavalin A lectin is employed as a model system to demonstrate the application of these methods to the special case of carbohydrate binding.

Conformational analysis of oligosaccharides and polysaccharides using molecular dynamics simulations

Complex carbohydrates usually have a large number of rotatable bonds and consequently a large number of theoretically possible conformations can be generated (combinatorial explosion). The application of systematic search methods for conformational analysis of carbohydrates is therefore limited to disaccharides and trisaccharides in a routine analysis. An alternative approach is to use Monte-Carlo methods or (high-temperature) molecular dynamics (MD) simulations to explore the conformational space of complex carbohydrates. This chapter describes how to use MD simulation data to perform a conformational analysis (conformational maps, hydrogen bonds) of oligosaccharides and how to build realistic 3D structures of large polysaccharides using Conformational Analysis Tools (CAT).

Carbohydrate Microarrays

Carbohydrates, like nucleic acids and proteins, are essential biological molecules. Owing to their intrinsic physicochemical properties, carbohydrates are capable of generating structural diversity in a multitude of ways and are prominently displayed on the surfaces of cell membranes or on the exposed regions of macromolecules. Recent studies highlight that carbohydrate moieties are critical for molecular recognition, cell-cell interactions, and cell signaling in many physiological and pathological processes, and for biocommunication between microbes and host species. Modern carbohydrate microarrays emerged in 2002 and brought in new high-throughput tools for “glyco code” exploration. In this section, some basic concepts of sugar chain diversity, glyco-epitope recognition, and the evolving area of glyco-epitomics and biomarker discovery are discussed. Two complementary technologies, carbohydrate antigen arrays and photogenerated glyco-chips, serve as models to illustrate how to apply carbohydrate microarrays to address biomedical questions.

Glyco3D: a portal for structural glycosciences

The present work describes, in a detailed way, a family of databases covering the three-dimensional features of monosaccharides, disaccharides, oligosaccharides, polysaccharides, glycosyltransferases, lectins, monoclonal antibodies against carbohydrates, and glycosaminoglycan-binding proteins. These databases have been developed with non-proprietary software, and they are open freely to the scientific community. They are accessible through the common portal called "Glyco3D" http://www.glyco3d.cermav.cnrs.fr. The databases are accompanied by a user-friendly graphical user interface (GUI) which offers several search options. All three-dimensional structures are available for visual consultations (with basic measurements possibilities) and can be downloaded in commonly used formats for further uses.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Crossroads between Bacterial and Mammalian Glycosyltransferases

Bacterial glycosyltransferases (GT) often synthesize the same glycan linkages as mammalian GT; yet, they usually have very little sequence identity. Nevertheless, enzymatic properties, folding, substrate specificities, and catalytic mechanisms of these enzyme proteins may have significant similarity. Thus, bacterial GT can be utilized for the enzymatic synthesis of both bacterial and mammalian types of complex glycan structures. A comparison is made here between mammalian and bacterial enzymes that synthesize epitopes found in mammalian glycoproteins, and those found in the O antigens of Gram-negative bacteria. These epitopes include Thomsen–Friedenreich (TF or T) antigen, blood group O, A, and B, type 1 and 2 chains, Lewis antigens, sialylated and fucosylated structures, and polysialic acids. Many different approaches can be taken to investigate the substrate binding and catalytic mechanisms of GT, including crystal structure analyses, mutations, comparison of amino acid sequences, NMR, and mass spectrometry. Knowledge of the protein structures and functions helps to design GT for specific glycan synthesis and to develop inhibitors. The goals are to develop new strategies to reduce bacterial virulence and to synthesize vaccines and other biologically active glycan structures.

Predicting the Origins of Anti-Blood Group Antibody Specificity: A Case Study of the ABO A- and B-Antigens

The ABO blood group system is the most important blood type system in human transfusion medicine. Here, we explore the specificity of antibody recognition toward ABO blood group antigens using computational modeling and biolayer interferometry. Automated docking and molecular dynamics simulations were used to explore the origin of the specificity of an anti-blood group A antibody variable fragment (Fv AC1001). The analysis predicts a number of Fv-antigen interactions that contribute to affinity, including a hydrogen bond between a His^L49 and the carbonyl moiety of the GalNAc in antigen A. This interaction was consistent with the dependence of affinity on pH, as measured experimentally; at lower pH there is an increase in binding affinity. Binding energy calculations provide unique insight into the origin of interaction energies at a per-residue level in both the scFv and the trisaccharide antigen. The calculations indicate that while the antibody can accommodate both blood group A and B antigens in its combining site, the A antigen is preferred by 4 kcal/mol, consistent with the lack of binding observed for the B antigen.

A guide into glycosciences: How chemistry, biochemistry and biology cooperate to crack the sugar code

BACKGROUND

The most demanding challenge in research on molecular aspects within the flow of biological information is posed by the complex carbohydrates (glycan part of cellular glycoconjugates). How the 'message' encoded in carbohydrate 'letters' is 'read' and 'translated' can only be unraveled by interdisciplinary efforts.

SCOPE OF REVIEW

This review provides a didactic step-by-step survey of the concept of the sugar code and the way strategic combination of experimental approaches characterizes structure-function relationships, with resources for teaching.

MAJOR CONCLUSIONS

The unsurpassed coding capacity of glycans is an ideal platform for generating a broad range of molecular 'messages'. Structural and functional analyses of complex carbohydrates have been made possible by advances in chemical synthesis, rendering production of oligosaccharides, glycoclusters and neoglycoconjugates possible. This availability facilitates to test the glycans as ligands for natural sugar receptors (lectins). Their interaction is a means to turn sugar-encoded information into cellular effects. Glycan/lectin structures and their spatial modes of presentation underlie the exquisite specificity of the endogenous lectins in counterreceptor selection, that is, to home in on certain cellular glycoproteins or glycolipids.

GENERAL SIGNIFICANCE

Understanding how sugar-encoded 'messages' are 'read' and 'translated' by lectins provides insights into fundamental mechanisms of life, with potential for medical applications.

Importance of ligand conformational energies in carbohydrate docking: Sorting the wheat from the chaff

Docking algorithms that aim to be applicable to a broad range of ligands suffer reduced accuracy because they are unable to incorporate ligand-specific conformational energies. Here, we develop a set of Carbohydrate Intrinsic (CHI) energy functions that quantify the conformational properties of oligosaccharides, based on the values of their glycosidic torsion angles. The relative energies predicted by the CHI energy functions mirror the conformational distributions of glycosidic linkages determined from a survey of oligosaccharide-protein complexes in the protein data bank. Addition of CHI energies to the standard docking scores in Autodock 3, 4.2, and Vina consistently improves pose ranking of oligosaccharides docked to a set of anticarbohydrate antibodies. The CHI energy functions are also independent of docking algorithm, and with minor modifications, may be incorporated into both theoretical modeling methods, and experimental NMR or X-ray structure refinement programs.

NMR-based dynamics of free glycosaminoglycans in solution

Dynamical behaviors of glycosaminoglycans, as here illustrated with a hyaluronan oligosaccharide, are key regulators of biological functions.

The influence of N-linked glycans on the molecular dynamics of the HIV-1 gp120 V3 loop

N-linked glycans attached to specific amino acids of the gp120 envelope trimer of a HIV virion can modulate the binding affinity of gp120 to CD4, influence coreceptor tropism, and play an important role in neutralising antibody responses. Because of the challenges associated with crystallising fully glycosylated proteins, most structural investigations have focused on describing the features of a non-glycosylated HIV-1 gp120 protein. Here, we use a computational approach to determine the influence of N-linked glycans on the dynamics of the HIV-1 gp120 protein and, in particular, the V3 loop. We compare the conformational dynamics of a non-glycosylated gp120 structure to that of two glycosylated gp120 structures, one with a single, and a second with five, covalently linked high-mannose glycans. Our findings provide a clear illustration of the significant effect that N-linked glycosylation has on the temporal and spatial properties of the underlying protein structure. We find that glycans surrounding the V3 loop modulate its dynamics, conferring to the loop a marked propensity towards a more narrow conformation relative to its non-glycosylated counterpart. The conformational effect on the V3 loop provides further support for the suggestion that N-linked glycosylation plays a role in determining HIV-1 coreceptor tropism.

Conference report: "Functional Glycomics in HIV Type 1 Vaccine Design" workshop report, Bethesda, Maryland, April 30-May 1, 2012

A vital part of the renewed hope for a vaccine against the human immunodeficiency virus (HIV-1) is based on recent studies that have highlighted major sites of HIV-1 vulnerability that could be effectively targeted by a preventive vaccine. One of these potential vulnerabilities includes the dense cluster of carbohydrates surrounding HIV-1's envelope glycoproteins gp120 and gp41, typically referred to as the "glycan shield." Recent data from several laboratories have shown that glycans on the HIV-1 envelope form key epitopes for broadly neutralizing antibodies (bNAb). Moreover, HIV-1 envelope glycans play an important role in viral transmission, antigenicity, and immunogenicity. The recent availability of novel tools and technologies has now allowed investigators to leverage glycomic structure-function relationships in the design of candidate HIV-1 vaccines. Additionally, glycans modulate the immune response, playing an essential role in Fc receptor and complement activity. To promote cross-disciplinary collaboration and promote synergistic HIV-1- glycomics research, the National Institutes of Health (NIH) cosponsored and convened a 1.5-day workshop entitled "Functional Glycomics in HIV-1 Vaccine Design." The meeting focused on the role of glycan interactions with neutralizing antibodies, the influence of immunoglobulin G (IgG) Fc receptor glycosylation, newly available glycomics technologies, and how new information on the role of glycans could be applied in HIV-1 immunogen design strategies. This report summarizes the discussions of this workshop.

Protein-carbohydrate interactions studied by NMR: from molecular recognition to drug design

Diseases that result from infection are, in general, a consequence of specific interactions between a pathogenic organism and the cells. The study of host-pathogen interactions has provided insights for the design of drugs with therapeutic properties. One area that has proved to be promising for such studies is the constituted by carbohydrates which participate in biological processes of paramount importance. On the one hand, carbohydrates have shown to be information carriers with similar, if not higher, importance than traditionally considered carriers as amino acids and nucleic acids. On the other hand, the knowledge on molecular recognition of sugars by lectins and other carbohydrate-binding proteins has been employed for the development of new biomedical strategies. Biophysical techniques such as X-Ray crystallography and NMR spectroscopy lead currently the investigation on this field. In this review, a description of traditional and novel NMR methodologies employed in the study of sugar-protein interactions is briefly presented in combination with a palette of NMR-based studies related to biologically and/or pharmaceutically relevant applications.

Simulation of carbohydrates, from molecular docking to dynamics in water

Modeling of carbohydrates is particularly challenging because of the variety of structures resulting for the high number of monosaccharides and possible linkages and also because of their intrinsic flexibility. The development of carbohydrate parameters for molecular modeling is still an active field. Nowadays, main carbohydrates force fields are GLYCAM06, CHARMM36, and GROMOS 45A4. GLYCAM06 includes the largest choice of compounds and is compatible with the AMBER force fields and associated. Furthermore, AMBER includes tools for the implementation of new parameters. When looking at protein-carbohydrate interaction, the choice of the starting structure is of importance. Such complex can be sometimes obtained from the Protein Data Bank-although the stereochemistry of sugars may require some corrections. When no experimental data is available, molecular docking simulation is generally used to the obtain protein-carbohydrate complex coordinates. As molecular docking parameters are not specifically dedicated to carbohydrates, inaccuracies should be expected, especially for the docking of polysaccharides. This issue can be addressed at least partially by combining molecular docking with molecular dynamics simulation in water.

Solvent structure improves docking prediction in lectin-carbohydrate complexes

Recognition and complex formation between proteins and carbohydrates is a key issue in many important biological processes. Determination of the three-dimensional structure of such complexes is thus most relevant, but particularly challenging because of their usually low binding affinity. In silico docking methods have a long-standing tradition in predicting protein-ligand complexes, and allow a potentially fast exploration of a number of possible protein-carbohydrate complex structures. However, determining which of these predicted complexes represents the correct structure is not always straightforward. In this work, we present a modification of the scoring function provided by AutoDock4, a widely used docking software, on the basis of analysis of the solvent structure adjacent to the protein surface, as derived from molecular dynamics simulations, that allows the definition and characterization of regions with higher water occupancy than the bulk solvent, called water sites. They mimic the interaction held between the carbohydrate -OH groups and the protein. We used this information for an improved docking method in relation to its capacity to correctly predict the protein-carbohydrate complexes for a number of tested proteins, whose ligands range in size from mono- to tetrasaccharide. Our results show that the presented method significantly improves the docking predictions. The resulting solvent-structure-biased docking protocol, therefore, appears as a powerful tool for the design and optimization of development of glycomimetic drugs, while providing new insights into protein-carbohydrate interactions. Moreover, the achieved improvement also underscores the relevance of the solvent structure to the protein carbohydrate recognition process.

Carbohydrate force fields

Carbohydrates present a special set of challenges to the generation of force fields. First, the tertiary structures of monosaccharides are complex merely by virtue of their exceptionally high number of chiral centers. In addition, their electronic characteristics lead to molecular geometries and electrostatic landscapes that can be challenging to predict and model. The monosaccharide units can also interconnect in many ways, resulting in a large number of possible oligosaccharides and polysaccharides, both linear and branched. These larger structures contain a number of rotatable bonds, meaning they potentially sample an enormous conformational space. This article briefly reviews the history of carbohydrate force fields, examining and comparing their challenges, forms, philosophies, and development strategies. Then it presents a survey of recent uses of these force fields, noting trends, strengths, deficiencies, and possible directions for future expansion.

The use of glycoinformatics in glycochemistry

Glycoinformatics is a small but growing branch of bioinformatics and chemoinformatics. Various resources are now available that can be of use to glycobiologists, but also to chemists who work on the synthesis or analysis of carbohydrates. This article gives an overview of existing glyco-specific databases and tools, with a focus on their application to glycochemistry: Databases can provide information on candidate glycan structures for synthesis, or on glyco-enzymes that can be used to synthesize carbohydrates. Statistical analyses of glycan databases help to plan glycan synthesis experiments. 3D-Structural data of protein-carbohydrate complexes are used in targeted drug design, and tools to support glycan structure analysis aid with quality control. Specific problems of glycoinformatics compared to bioinformatics for genomics or proteomics, especially concerning integration and long-term maintenance of the existing glycan databases, are also discussed.