Deciphering the glycan preference of bacterial lectins by glycan array and molecular docking with validation by microcalorimetry and crystallography

Comparative Bindings of Glycosaminoglycans with CXCL8 Monomer and Dimer: Insights from Conformational Dynamics and Kinetics of Hydrogen Bonds

GAGs bind to both the monomeric and dimeric forms of CXCL8, helping to form a concentration gradient of the chemokine that facilitates the recruitment of neutrophils to an injury site and supports other biological functions. In this study, atomistic molecular dynamics simulations were conducted to investigate the binding behavior of two hexameric GAGs sulfated at two different positions, chondroitin sulfate (CS) and heparan sulfate (HS), with the monomer (SIL8) and dimer (DIL8) forms of the CXCL8 protein. The results support that the conformational diversity of CS and HS appeared to be more when binding with monomer SIL8 than dimer DIL8. CS gained more configurational entropy from glycosidic linkage flexibility when bound to SIL8 than DIL8, with a higher energy barrier, whereas HS exhibited a lower energy barrier for configurational entropy when bound to SIL8 and DIL8. The monomer SIL8 exhibited more favorable and preferential binding with GAGs compared to DIL8. Formation of hydrogen bonds with the basic amino acids of SIL8 and GAG was more rigid and required higher activation energy to break than the other identified hydrogen bondings. Water molecules involved in hydrogen bonding with GAGs, excluding those with basic amino acids of DIL8, showed longer lifetimes and slower relaxation compared to SIL8. This suggests that water-mediated interactions also favor binding of DIL8 with GAGs. Despite having more basic amino acids, DIL8 did not display stronger binding than SIL8, indicating the significant role of basic residues in stabilizing the GAG-protein interactions in the monomers. The reason could be that the greater number of nonbasic amino acids in dimeric CXCL8 stabilizes the complex by forming water-mediated hydrogen bonds, reducing the conformational preferences for binding with GAGs. In contrast, the monomeric form of CXCL8 exhibits a higher conformational preference for protein-GAG interactions.

Conformational preferences of heparan sulfate to recognize the CXCL8 dimer in aqueous medium: degree of sulfation and hydrogen bonds

The sulfation pattern and epimerization of the long-chain sulfated polysaccharide heparan sulfate (HS) cause structural diversity and regulate various physiological and pathological processes when binding with proteins. In this work, we performed a series of molecular dynamics simulations of three variants of the octadecasaccharide HS with varying sulfation positions in aqueous medium in their free forms and in the presence of the chemokine CXCL8 dimer. The free energy of binding depicts the sulfation at the 6-O position of GlcNAc (HS6S), and both 3-O and 6-O positions of GlcNAc (HS3S6S) of HS variants are more likely to bind with the CXCL8 dimer than the triply sulfated HS2S3S6S, which is sulfated at the 2-O position of GlcUA additionally along with 3-O and 6-O positions of GlcNAc. Binding between HS and CXCL8 was driven by electrostatic and van der Waals interactions predominantly regardless of the sulfation pattern; however, unfavorable entropic contribution suppressed the interaction between HS and CXCL8. The contribution of different amino acid residues to the binding energetics suggested that basic amino acids line up the binding site of CXCL8. This study further acknowledges the role of interfacial water that is structured and bound with HS through hydrogen bonds, exhibiting differential hydrogen bond relaxation dynamics compared to when the HS molecules are free. Moreover, this study identifies that with the increase in sulfation, the HS-water hydrogen bond relaxation occurs faster with the complexation, while the reverse trend is followed in their free forms. Significant structural adaptation of the different sulfated HS molecules, as verified from the free energy landscapes generated from various reaction coordinates, root-mean-square-deviations, end-to-end distances, including ring pucker angles, dihedral flexibility, and the high conformational entropy cost arising from the glycosidic bonds, suggests that the different sulfated variants of HS undergo significant structural transformation to bind with CXCL8. The presence of a CXCL8 dimer imposes the bound forms of HS to adopt non-linear structures with skew-boat conformations. The atomistic details of the study would help in understanding the selectivity and conformational diversity, as well as the role of solvents in the recognition of CXCL8 by different sulfated variants of HS molecules.

Glycan microarrays from construction to applications

Through their specific interactions with proteins, cellular glycans play key roles in a wide range of physiological and pathological processes. One of the main goals of research in the areas of glycobiology and glycomedicine is to understand glycan-protein interactions at the molecular level. Over the past two decades, glycan microarrays have become powerful tools for the rapid evaluation of interactions between glycans and proteins. In this review, we briefly describe methods used for the preparation of glycan probes and the construction of glycan microarrays. Next, we highlight applications of glycan microarrays to rapid profiling of glycan-binding patterns of plant, animal and pathogenic lectins, as well as other proteins. Finally, we discuss other important uses of glycan microarrays, including the rapid analysis of substrate specificities of carbohydrate-active enzymes, the quantitative determination of glycan-protein interactions, discovering high-affinity or selective ligands for lectins, and identifying functional glycans within cells. We anticipate that this review will encourage researchers to employ glycan microarrays in diverse glycan-related studies.

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.

Advances in understanding N-glycosylation structure, function, and regulation in health and disease

N-linked glycosylation is a post-translational modification crucial for membrane protein folding, stability and other cellular functions. Alteration of membrane protein N-glycans is implicated in wide range of pathological conditions including cancer metastasis, chronic inflammatory diseases, and viral pathogenesis. Even though the roles of N-glycans have been studied extensively, our knowledge of their mechanisms remains unclear due to the lack of detailed structural analysis of the N-glycome. Mapping the N-glycome landscape will open new avenues to explore disease mechanisms and identify novel therapeutic targets. This review discusses the diverse structure of N-linked glycans, the function and regulation of N-glycosylation in health and disease, and ends with a focus on recent approaches to target N-glycans in rheumatoid arthritis and cancer metastasis.

Investigating Conformational Dynamics of Lewis Y Oligosaccharides and Elucidating Blood Group Dependency of Cholera Using Molecular Dynamics

Cholera is caused by Vibrio cholerae and is an example of a blood-group-dependent disease. Recent studies suggest that the receptor-binding B subunit of the cholera toxin (CT) binds histo-blood group antigens at a secondary binding site. Herein, we studied the conformational dynamics of Lewis Y (Le^Y) oligosaccharides, H-tetrasaccharides and A-pentasaccharides, in aqueous solution by conducting accelerated molecular dynamics (aMD) simulations. The flexible nature of both oligosaccharides was displayed in aMD simulations. Furthermore, aMD simulations revealed that for both oligosaccharides in the free form, ⁴C₁ and ¹C₄ puckers were sampled for all but GalNAc monosaccharides, while either the ⁴C₁ (GlcNAc, Gal, GalNAc) or ¹C₄ (Fuc2, Fuc3) pucker was sampled in the CT-bound forms. In aMD, the complete transition from the ⁴C₁ to ¹C₄ pucker was sampled for GlcNAc and Gal in both oligosaccharides. Further, we have observed a transition from the open to closed conformer in the case of A-pentasaccharide, while H-tetrasaccharide remains in the open conformation throughout the simulation. Both oligosaccharides adopted an open conformation in the CT binding site. Moreover, we have investigated the molecular basis of recognition of Le^Y oligosaccharides by the B subunit of the cholera toxin of classical and El Tor biotypes using the molecular mechanics generalized Born surface area (MM/GBSA) scheme. The O blood group determinant, H-tetrasaccharide, exhibits a stronger affinity to both biotypes compared to the A blood group determinant, A-pentasaccharide, which agrees with the experimental data. The difference in binding free energy between O and A blood group determinants mainly arises due to the increased entropic cost and desolvation energy in the case of A-pentasaccharide compared to that of H-tetrasaccharide. Our study also reveals that the terminal Fuc3 contributes most to the binding free energy compared to other carbohydrate residues as it forms multiple hydrogen bonds with CT. Overall, our study might help in designing glycomimetic drugs targeting the cholera toxin.

Ring Puckering Landscapes of Glycosaminoglycan-Related Monosaccharides from Molecular Dynamics Simulations

The conformational flexibility of the glycosaminoglycans (GAGs) is known to be key in their binding and biological function, for example in regulating coagulation and cell growth. In this work, we employ enhanced sampling molecular dynamics simulations to probe the ring conformations of GAG-related monosaccharides, including a range of acetylated and sulfated GAG residues. We first perform unbiased MD simulations of glucose anomers and the epimers glucuronate and iduronate. These calculations indicate that in some cases, an excess of 15 μs is required for adequate sampling of ring pucker due to the high energy barriers between states. However, by applying our recently developed msesMD simulation method (multidimensional swarm-enhanced sampling molecular dynamics), we were able to quantitatively and rapidly reproduce these ring pucker landscapes. From msesMD simulations, the puckering free energy profiles were then compared for 15 further monosaccharides related to GAGs; this includes to our knowledge the first simulation study of sulfation effects on β-GalNAc ring puckering. For the force field employed, we find that in general the calculated pucker free energy profiles for sulfated sugars were similar to the corresponding unsulfated profiles. This accords with recent experimental studies suggesting that variation in ring pucker of sulfated GAG residues is primarily dictated by interactions with surrounding residues rather than by intrinsic conformational preference. As an exception to this, however, we predict that 4-O-sulfation of β-GalNAc leads to reduced ring rigidity, with a significant lowering in energy of the ¹C₄ ring conformation; this observation may have implications for understanding the structural basis of the biological function of β-GalNAc-containing glycosaminoglycans such as dermatan sulfate.

Next-Generation Glycan Microarray Enabled by DNA-Coded Glycan Library and Next-Generation Sequencing Technology

Interactions of glycans with proteins, cells, and microorganisms play important roles in cell-cell adhesion and host-pathogen interaction. Glycan microarray technology, in which multiple glycan structures are immobilized on a single glass slide and interrogated with glycan-binding proteins (GBPs), has become an indispensable tool in the study of protein-glycan interactions. Despite its great success, the current format of the glycan microarray requires expensive, specialized instrumentation and labor-intensive assay and image processing procedures, which limit automation and possibilities for high-throughput analyses. Furthermore, the current microarray is not suitable for assaying interaction with intact cells due to their large size compared to the two-dimensional microarray surface. To address these limitations, we developed the next-generation glycan microarray (NGGM) based on artificial DNA coding of glycan structures. In this novel approach, a glycan library is presented as a mixture of glycans and glycoconjugates, each of which is coded with a unique oligonucleotide sequence (code). The glycan mixture is interrogated by GBPs followed by the separation of unbound coded glycans. The DNA sequences that identify individual bound glycans are quantitatively sequenced (decoded) by powerful next-generation sequencing (NGS) technology, and copied numbers of the DNA codes represent relative binding specificities of corresponding glycan structures to GBPs. We demonstrate that NGGM generates glycan-GBP binding data that are consistent with that generated in a slide-based glycan microarray. More importantly, the solution phase binding assay is directly applicable to identifying glycan binding to intact cells, which is often challenging using glass slide-based glycan microarrays.

Screening of a Library of Oligosaccharides Targeting Lectin LecB of Pseudomonas Aeruginosa and Synthesis of High Affinity Oligoglycoclusters

The Gram negative bacterium Pseudomonas aeruginosa (PA) is an opportunistic bacterium that causes severe and chronic infection of immune-depressed patients. It has the ability to form a biofilm that gives a selective advantage to the bacteria with respect to antibiotherapy and host defenses. Herein, we have focused on the tetrameric soluble lectin which is involved in bacterium adherence to host cells, biofilm formation, and cytotoxicity. It binds to l-fucose, d-mannose and glycan exposing terminal fucose or mannose. Using a competitive assay on microarray, 156 oligosaccharides and polysaccharides issued from fermentation or from the biomass were screened toward their affinity to LecB. Next, the five best ligands (Lewis^a, Lewis^b, Lewis^x, siayl-Lewis^x and 3-fucosyllactose) were derivatized with a propargyl aglycon allowing the synthesis of 25 trivalent, 25 tetravalent and 5 monovalent constructions thanks to copper catalyzed azide alkyne cycloaddition. The 55 clusters were immobilized by DNA Directed immobilization leading to the fabrication of a glycocluster microarray. Their binding to LecB was studied. Multivalency improved the binding to LecB. The binding structure relationship of the clusters is mainly influenced by the carbohydrate residues. Molecular simulations indicated that the simultaneous contact of both binding sites of monomer A and D seems to be energetically possible.

Assessing the Performance of MM/PBSA, MM/GBSA, and QM-MM/GBSA Approaches on Protein/Carbohydrate Complexes: Effect of Implicit Solvent Models, QM Methods, and Entropic Contributions

Rapid and accurate binding affinity prediction of protein-carbohydrate complexes is a major challenge in glycomimetics design. Among the existing computational techniques, end-point methods have received considerable interest because of their low computational cost. However, significant obstacles remain when such methods are applied to protein-glycan complexes. This article reports the performance of end-point free-energy calculation methods: molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA), MM/generalized Born surface area (MM/GBSA), and quantum mechanics-MM/GBSA (QM-MM/GBSA) on monosaccharides bound to RSL lectin from Ralstonia solanacearum. A careful investigation of the molecular dynamics simulation length, van der Waals radii sets, GB models, QM Hamiltonians, and entropic compensation has been made, and the results are compared with the experimental binding free energies from isothermal titration calorimetry/surface plasmon resonance measurements. The binding free energies using implicit solvent methods are found to be sensitive to the simulation length, radii set, GB model, and QM Hamiltonian. A simulation length of 10 ns using the radii set mbondi provides the best agreement with the experimental values ( r² = 0.96) by MM/PBSA. The GB^HCT model is in accord with the experimental values in MM/GBSA ( r² = 0.91) or in combination with parameterized model number 6 (PM6) ( r² = 0.98) in QM-MM/GBSA. Out of 12 QM Hamiltonians tested, PM6, density functional theory-based tight binding (DFTB), and their variants proved to be more efficient than other semiempirical methods. These methods perform equally well in predicting both absolute and relative binding free energies.

Influence of Trp flipping on carbohydrate binding in lectins. An example on Aleuria aurantia lectin AAL

Protein-carbohydrate interactions are very often mediated by the stacking CH-π interactions involving the side chains of aromatic amino acids such as tryptophan (Trp), tyrosine (Tyr) or phenylalanine (Phe). Especially suitable for stacking is the Trp residue. Analysis of the PDB database shows Trp stacking for 265 carbohydrate or carbohydrate like ligands in 5 208 Trp containing motives. An appropriate model system to study such an interaction is the AAL lectin family where the stacking interactions play a crucial role and are thought to be a driving force for carbohydrate binding. In this study we present data showing a novel finding in the stacking interaction of the AAL Trp side chain with the carbohydrate. High resolution X-ray structure of the AAL lectin from Aleuria aurantia with α-methyl-l-fucoside ligand shows two possible Trp side chain conformations with the same occupation in electron density. The in silico data shows that the conformation of the Trp side chain does not influence the interaction energy despite the fact that each conformation creates interactions with different carbohydrate CH groups. Moreover, the PDB data search shows that the conformations are almost equally distributed across all Trp-carbohydrate complexes, which would suggest no substantial preference for one conformation over another.

Molecular Simulations of Carbohydrates with a Fucose-Binding Burkholderia ambifaria Lectin Suggest Modulation by Surface Residues Outside the Fucose-Binding Pocket

Burkholderia ambifaria is an opportunistic respiratory pathogen belonging to the Burkholderia cepacia complex, a collection of species responsible for the rapidly fatal cepacia syndrome in cystic fibrosis patients. A fucose-binding lectin identified in the B. ambifaria genome, BambL, is able to adhere to lung tissue, and may play a role in respiratory infection. X-ray crystallography has revealed the bound complex structures for four fucosylated human blood group epitopes (blood group B, H type 1, H type 2, and Le^x determinants). The present study employed computational approaches, including docking and molecular dynamics (MD), to extend the structural analysis of BambL-oligosaccharide complexes to include four additional blood group saccharides (A, Le^a, Le^b, and Le^y) and a library of blood-group-related carbohydrates. Carbohydrate recognition is dominated by interactions with fucose via a hydrogen-bonding network involving Arg15, Glu26, Ala38, and Trp79 and a stacking interaction with Trp74. Additional hydrogen bonds to non-fucose residues are formed with Asp30, Tyr35, Thr36, and Trp74. BambL recognition is dominated by interactions with fucose, but also features interactions with other parts of the ligands that may modulate specificity or affinity. The detailed computational characterization of the BambL carbohydrate-binding site provides guidelines for the future design of lectin inhibitors.

Histo-blood group antigens as mediators of infections

The critical first step of a microbial infection is usually the attachment of pathogens to host cell glycans. Targets on host tissues are in particular the histo-blood group antigens (HBGAs), which are present in rich diversity in the mucus layer and on the underlying mucosa. Recent structural and functional studies have revealed significant new insight into the molecular mechanisms, explaining why individuals with certain blood groups are at increased risk of some infections. The most prominent example of blood-group-associated diseases is cholera, caused by infection with Vibrio cholerae. Many other microbial pathogens, for example Pseudomonas aeruginosa infecting the airways, and enterotoxigenic Escherichia coli (ETEC) causing traveler's diarrhea, also bind to histo-blood group antigens, but show a less clear correlation with blood group phenotype. Yet other pathogens, for example norovirus and Helicobacter pylori, recognize HBGAs differently depending on the strain. In all cases, milk oligosaccharides can aid the hosts' defenses, acting as natural receptor decoys, and anti-infectious therapy can be designed along similar strategies. In this review, we focus on important infections of humans, but the molecular mechanisms are of general relevance to a broad range of microbial infections of humans and animals.

The lectin-like protein 1 in Lactobacillus rhamnosus GR-1 mediates tissue-specific adherence to vaginal epithelium and inhibits urogenital pathogens

The probiotic Lactobacillus rhamnosus GR-1 has been documented to survive implantation onto the vaginal epithelium and interfere with urogenital pathogens. However, the molecular mechanisms involved are largely unknown. Here, we report for the first time the construction of dedicated knock-out mutants in L. rhamnosus GR-1 to enable the study of gene functions. In a search for genes responsible for the adherence capacity of L. rhamnosus GR-1, a genomic region encoding a protein with homology to lectin-like proteins was identified. Phenotypic analyses of the knock-out mutant of L. rhamnosus GR-1 revealed a two-fold decreased adhesion to the vaginal and ectocervical epithelial cell lines compared to wild-type. In contrast, the adhesion to gastro-intestinal epithelial (Caco2) and endocervical cell lines (Hela and End1/E6E7) was not drastically affected by the mutation, suggesting that the LGR-1_Llp1 lectins mediates tissue tropism. The purified LGR-1_Llp1 protein also inhibited biofilm formation and adhesion of uropathogenic Escherichia coli. For the first time, an important role for a novel lectin-like protein in the adhesion capacity and host cell-specific interaction of a vaginal probiotic Lactobacillus strain has been discovered, with an additional role in pathogen inhibition.

High mannose-specific lectin Msl mediates key interactions of the vaginal Lactobacillus plantarum isolate CMPG5300

To characterize the interaction potential of the human vaginal isolate Lactobacillus plantarum CMPG5300, its genome was mined for genes encoding lectin-like proteins. cmpg5300.05_29 was identified as the gene encoding a putative mannose-binding lectin. Phenotypic analysis of a gene knock-out mutant of cmpg5300.05_29 showed that expression of this gene is important for auto-aggregation, adhesion to the vaginal epithelial cells, biofilm formation and binding to mannosylated glycans. Purification of the predicted lectin domain of Cmpg5300.05_29 and characterization of its sugar binding capacity confirmed the specificity of the lectin for high- mannose glycans. Therefore, we renamed Cmpg5300.05_29 as a mannose-specific lectin (Msl). The purified lectin domain of Msl could efficiently bind to HIV-1 glycoprotein gp120 and Candida albicans, and showed an inhibitory activity against biofilm formation of uropathogenic Escherichia coli, Staphylococcus aureus and Salmonella Typhimurium. Thus, using a combination of molecular lectin characterization and functional assays, we could show that lectin-sugar interactions play a key role in host and pathogen interactions of a prototype isolate of the vaginal Lactobacillus microbiota.

Lectin-Like Molecules of Lactobacillus rhamnosus GG Inhibit Pathogenic Escherichia coli and Salmonella Biofilm Formation

Objectives

Increased antibiotic resistance has catalyzed the research on new antibacterial molecules and alternative strategies, such as the application of beneficial bacteria. Since lectin molecules have unique sugar-recognizing capacities, and pathogens are often decorated with sugars that affect their survival and infectivity, we explored whether lectins from the probiotic strain Lactobacillus rhamnosus GG have antipathogenic properties.

Methods

The genome sequence of L. rhamnosus GG was screened for the presence of lectin-like proteins. Two genes, LGG_RS02780 and LGG_RS02750, encoding for polypeptides with an N-terminal conserved L-type lectin domain were detected and designated Llp1 (lectin-like protein 1) and Llp2. The capacity of Llp1 and Llp2 to inhibit biofilm formation of various pathogens was investigated. Sugar specificity was determined by Sepharose beads assays and glycan array screening.

Results

The isolated lectin domains of Llp1 and Llp2 possess pronounced inhibitory activity against biofilm formation by various pathogens, including clinical Salmonella species and uropathogenic E. coli, with Llp2 being more active than Llp1. In addition, sugar binding assays with Llp1 and Llp2 indicate specificity for complex glycans. Both proteins are also involved in the adhesion capacity of L. rhamnosus GG to gastrointestinal and vaginal epithelial cells.

Conclusions

Lectins isolated from or expressed by beneficial lactobacilli could be considered promising bio-active ingredients for improved prophylaxis of urogenital and gastrointestinal infections.

The Hidden Conformation of Lewis x, a Human Histo-Blood Group Antigen, Is a Determinant for Recognition by Pathogen Lectins

Histo-blood group epitopes are fucosylated branched oligosaccharides with well-defined conformations in solution that are recognized by receptors, such as lectins from pathogens. We report here the results of a series of experimental and computational endeavors revealing the unusual distortion of histo-blood group antigens by bacterial and fungal lectins. The Lewis x trisaccharide adopts a rigid closed conformation in solution, while crystallography and molecular dynamics reveal several higher energy open conformations when bound to the Ralstonia solanacearum lectin, which is in agreement with thermodynamic and kinetic measurements. Extensive molecular dynamics simulations confirm rare transient Le(x) openings in solution, frequently assisted by distortion of the central N-acetyl-glucosamine ring. Additional directed molecular dynamic trajectories revealed the role of a conserved tryptophan residue in guiding the fucose into the binding site. Our findings show that conformational adaptation of oligosaccharides is of paramount importance in cell recognition and should be considered when designing anti-infective glyco-compounds.

Combining 3D structure with glycan array data provides insight into the origin of glycan specificity

Defining how a glycan-binding protein (GBP) specifically selects its cognate glycan from among the ensemble of glycans within the cellular glycome is an area of intense study. Powerful insight into recognition mechanisms can be gained from 3D structures of GBPs complexed to glycans; however, such structures remain difficult to obtain experimentally. Here an automated 3D structure generation technique, called computational carbohydrate grafting, is combined with the wealth of specificity information available from glycan array screening. Integration of the array data with modeling and crystallography allows generation of putative co-complex structures that can be objectively assessed and iteratively altered until a high level of agreement with experiment is achieved. Given an accurate model of the co-complexes, grafting is also able to discern which binding determinants are active when multiple potential determinants are present within a glycan. In some cases, induced fit in the protein or glycan was necessary to explain the observed specificity, while in other examples a revised definition of the minimal binding determinants was required. When applied to a collection of 10 GBP-glycan complexes, for which crystallographic and array data have been reported, grafting provided a structural rationalization for the binding specificity of >90% of 1223 arrayed glycans. A webtool that enables researchers to perform computational carbohydrate grafting is available at www.glycam.org/gr (accessed 03 March 2016).

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.

Evaluation of Selected Classical Force Fields for Alchemical Binding Free Energy Calculations of Protein-Carbohydrate Complexes

Protein-carbohydrate recognition is crucial in many vital biological processes including host-pathogen recognition, cell-signaling, and catalysis. Accordingly, computational prediction of protein-carbohydrate binding free energies is of enormous interest for drug design. However, the accuracy of current force fields (FFs) for predicting binding free energies of protein-carbohydrate complexes is not well understood owing to technical challenges such as the highly polar nature of the complexes, anomerization, and conformational flexibility of carbohydrates. The present study evaluated the performance of alchemical predictions of binding free energies with the GAFF1.7/AM1-BCC and GLYCAM06j force fields for modeling protein-carbohydrate complexes. Mean unsigned errors of 1.1 ± 0.06 (GLYCAM06j) and 2.6 ± 0.08 (GAFF1.7/AM1-BCC) kcal·mol(-1) are achieved for a large data set of monosaccharide ligands for Ralstonia solanacearum lectin (RSL). The level of accuracy provided by GLYCAM06j is sufficient to discriminate potent, moderate, and weak binders, a goal that has been difficult to achieve through other scoring approaches. Accordingly, the protocols presented here could find useful applications in carbohydrate-based drug and vaccine developments.

Lectin engineering, a molecular evolutionary approach to expanding the lectin utilities

In the post genomic era, glycomics—the systematic study of all glycan structures of a given cell or organism—has emerged as an indispensable technology in various fields of biology and medicine. Lectins are regarded as “decipherers of glycans”, being useful reagents for their structural analysis, and have been widely used in glycomic studies. However, the inconsistent activity and availability associated with the plant-derived lectins that comprise most of the commercially available lectins, and the limit in the range of glycan structures covered, have necessitated the development of innovative tools via engineering of lectins on existing scaffolds. This review will summarize the current state of the art of lectin engineering and highlight recent technological advances in this field. The key issues associated with the strategy of lectin engineering including selection of template lectin, construction of a mutagenesis library, and high-throughput screening methods are discussed.

Improvements, trends, and new ideas in molecular docking: 2012-2013 in review

Molecular docking is a computational method for predicting the placement of ligands in the binding sites of their receptor(s). In this review, we discuss the methodological developments that occurred in the docking field in 2012 and 2013, with a particular focus on the more difficult aspects of this computational discipline. The main challenges and therefore focal points for developments in docking, covered in this review, are receptor flexibility, solvation, scoring, and virtual screening. We specifically deal with such aspects of molecular docking and its applications as selection criteria for constructing receptor ensembles, target dependence of scoring functions, integration of higher-level theory into scoring, implicit and explicit handling of solvation in the binding process, and comparison and evaluation of docking and scoring methods.

Recent advances in employing molecular modelling to determine the specificity of glycan-binding proteins

Impressive improvements in docking performance can be achieved by applying energy bonuses to poses in which glycan hydroxyl groups occupy positions otherwise preferred by bound waters. In addition, inclusion of glycosidic conformational energies allows unlikely glycan conformations to be appropriately penalized. A method for predicting the binding specificity of glycan-binding proteins has been developed, which is based on grafting glycan branches onto a minimal binding determinant in the binding site. Grafting can be used either to screen virtual libraries of glycans, such as the known glycome, or to identify docked poses of minimal binding determinants that are consistent with specificity data. The reviewed advances allow accurate modelling of carbohydrate-protein 3D co-complexes, but challenges remain in ranking the affinity of congeners.

The challenge and promise of glycomics

Glycomics is a broad and emerging scientific discipline focused on defining the structures and functional roles of glycans in biological systems. The staggering complexity of the glycome, minimally defined as the repertoire of glycans expressed in a cell or organism, has resulted in many challenges that must be overcome; these are being addressed by new advances in mass spectrometry as well as by the expansion of genetic and cell biology studies. Conversely, identifying the specific glycan recognition determinants of glycan-binding proteins by employing the new technology of glycan microarrays is providing insights into how glycans function in recognition and signaling within an organism and with microbes and pathogens. The promises of a more complete knowledge of glycomes are immense in that glycan modifications of intracellular and extracellular proteins have critical functions in almost all biological pathways.