NMR, molecular modeling, and crystallographic studies of lentil lectin-sucrose interaction

Molecular Level Sucrose Quantification: A Critical Review

Sucrose is a primary metabolite in plants, a source of energy, a source of carbon atoms for growth and development, and a regulator of biochemical processes. Most of the traditional analytical chemistry methods for sucrose quantification in plants require sample treatment (with consequent tissue destruction) and complex facilities, that do not allow real-time sucrose quantification at ultra-low concentrations (nM to pM range) under in vivo conditions, limiting our understanding of sucrose roles in plant physiology across different plant tissues and cellular compartments. Some of the above-mentioned problems may be circumvented with the use of bio-compatible ligands for molecular recognition of sucrose. Nevertheless, problems such as the signal-noise ratio, stability, and selectivity are some of the main challenges limiting the use of molecular recognition methods for the in vivo quantification of sucrose. In this review, we provide a critical analysis of the existing analytical chemistry tools, biosensors, and synthetic ligands, for sucrose quantification and discuss the most promising paths to improve upon its limits of detection. Our goal is to highlight the criteria design need for real-time, in vivo, highly sensitive and selective sucrose sensing capabilities to enable further our understanding of living organisms, the development of new plant breeding strategies for increased crop productivity and sustainability, and ultimately to contribute to the overarching need for food security.

Machine learning recognition of protein secondary structures based on two-dimensional spectroscopic descriptors

Protein secondary structure discrimination is crucial for understanding their biological function. It is not generally possible to invert spectroscopic data to yield the structure. We present a machine learning protocol which uses two-dimensional UV (2DUV) spectra as pattern recognition descriptors, aiming at automated protein secondary structure determination from spectroscopic features. Accurate secondary structure recognition is obtained for homologous (97%) and nonhomologous (91%) protein segments, randomly selected from simulated model datasets. The advantage of 2DUV descriptors over one-dimensional linear absorption and circular dichroism spectra lies in the cross-peak information that reflects interactions between local regions of the protein. Thanks to their ultrafast (∼200 fs) nature, 2DUV measurements can be used in the future to probe conformational variations in the course of protein dynamics.

A review on the antimicrobial properties of lectins

Lectins are biologically versatile biomolecules with remarkable antimicrobial effects, notably against bacteria, fungi and protozoa, in addition to modulating host immunity. For this, the lectins bind to carbohydrates on the surface of the pathogen, which can cause damage to the cell wall and prevent the attachment of microorganisms to host cells. Thus, this study intends to review the biological activities of lectins, with an emphasis on antimicrobial activity. Lectins of plant stood out for its antimicrobial effects, demonstrating that they act against a variety of strains, where in vitro were able to inhibit their development and affect their morphology. In vivo, they modulated host immunity, signaling and activating defense cells. Some of these lectins were capable to modulate the action of antibiotics, indicating their potential to minimize the antibiotic resistance. The results suggest that lectins have antimicrobial activity with potential to be used in drug development.

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.

A review of Vicieae lectins studies: End of the book or a story in the writing?

Vicieae tribe, Leguminosae family (Fabaceae), has been extensively studied. In particular, the study of lectins. The purification, physicochemical and structural characterizations of the various purified lectins and the analysis of their relevant biological activities are ongoing. In this review, several works already published about Vicieae lectins are addressed. Initially, we presented the purification protocols and the physicochemical aspects, such as specificity for carbohydrates, optimal activity in the face of variations in temperature and pH, as well metals-dependence. Following, structural characterization studies are highlighted and, finally, various biological activities already reported are summarized. Studies on lectins in almost all genera (Lathyrus, Lens, Pisum and Vicia) are considered, with the exception of Vavilovia which studies of lectins have not yet been reported. Like other leguminous lectins, Vicieae lectins present heterogeneous profiles of agglutination profiles for erythrocytes and other cells of the immune system, and glycoproteins. Most Vicieae lectins consist of two subunits, α and β, products of a single precursor protein derived from a single gene. The differences between the isoforms result from varying degrees of proteolytic processing. Along with the identification of these molecules and their characteristics, biological activities become very relevant and robust for both basic and applied research.

Are Dietary Lectins Relevant Allergens in Plant Food Allergy?

Lectins or carbohydrate-binding proteins are widely distributed in seeds and vegetative parts of edible plant species. A few lectins from different fruits and vegetables have been identified as potential food allergens, including wheat agglutinin, hevein (Hev b 6.02) from the rubber tree and chitinases containing a hevein domain from different fruits and vegetables. However, other well-known lectins from legumes have been demonstrated to behave as potential food allergens taking into account their ability to specifically bind IgE from allergic patients, trigger the degranulation of sensitized basophils, and to elicit interleukin secretion in sensitized people. These allergens include members from the different families of higher plant lectins, including legume lectins, type II ribosome-inactivating proteins (RIP-II), wheat germ agglutinin (WGA), jacalin-related lectins, GNA (Galanthus nivalis agglutinin)-like lectins, and Nictaba-related lectins. Most of these potentially active lectin allergens belong to the group of seed storage proteins (legume lectins), pathogenesis-related protein family PR-3 comprising hevein and class I, II, IV, V, VI, and VII chitinases containing a hevein domain, and type II ribosome-inactivating proteins containing a ricin B-chain domain (RIP-II). In the present review, we present an exhaustive survey of both the structural organization and structural features responsible for the allergenic potency of lectins, with special reference to lectins from dietary plant species/tissues consumed in Western countries.

The role of hydration effects in 5-fluorouridine binding to SOD1: insight from a new 3D-RISM-KH based protocol for including structural water in docking simulations

Misfolded Cu/Zn superoxide dismutase enzyme (SOD1) shows prion-like propagation in neuronal cells leading to neurotoxic aggregates that are implicated in amyotrophic lateral sclerosis (ALS). Tryptophan-32 (W32) in SOD1 is part of a potential site for templated conversion of wild type SOD1. This W32 binding site is located on a convex, solvent exposed surface of the SOD1 suggesting that hydration effects can play an important role in ligand recognition and binding. A recent X-ray crystal structure has revealed that 5-Fluorouridine (5-FUrd) binds at the W32 binding site and can act as a pharmacophore scaffold for the development of anti-ALS drugs. In this study, a new protocol is developed to account for structural (non-displaceable) water molecules in docking simulations and successfully applied to predict the correct docked conformation binding modes of 5-FUrd at the W32 binding site. The docked configuration is within 0.58 Å (RMSD) of the observed configuration. The docking protocol involved calculating a hydration structure around SOD1 using molecular theory of solvation (3D-RISM-KH, 3D-Reference Interaction Site Model-Kovalenko-Hirata) whereby, non-displaceable water molecules are identified for docking simulations. This protocol was also used to analyze the hydrated structure of the W32 binding site and to explain the role of solvation in ligand recognition and binding to SOD1. Structural water molecules mediate hydrogen bonds between 5-FUrd and the receptor, and create an environment favoring optimal placement of 5-FUrd in the W32 binding site.

Functional similarities between the dictyostelium protein AprA and the human protein dipeptidyl-peptidase IV

Autocrine proliferation repressor protein A (AprA) is a protein secreted by Dictyostelium discoideum cells. Although there is very little sequence similarity between AprA and any human protein, AprA has a predicted structural similarity to the human protein dipeptidyl peptidase IV (DPPIV). AprA is a chemorepellent for Dictyostelium cells, and DPPIV is a chemorepellent for neutrophils. This led us to investigate if AprA and DPPIV have additional functional similarities. We find that like AprA, DPPIV is a chemorepellent for, and inhibits the proliferation of, D. discoideum cells, and that AprA binds some DPPIV binding partners such as fibronectin. Conversely, rAprA has DPPIV-like protease activity. These results indicate a functional similarity between two eukaryotic chemorepellent proteins with very little sequence similarity, and emphasize the usefulness of using a predicted protein structure to search a protein structure database, in addition to searching for proteins with similar sequences.

Coronavirus receptor switch explained from the stereochemistry of protein-carbohydrate interactions and a single mutation

Hemagglutinin-esterases (HEs) are bimodular envelope proteins of orthomyxoviruses, toroviruses, and coronaviruses with a carbohydrate-binding "lectin" domain appended to a receptor-destroying sialate-O-acetylesterase ("esterase"). In concert, these domains facilitate dynamic virion attachment to cell-surface sialoglycans. Most HEs (type I) target 9-O-acetylated sialic acids (9-O-Ac-Sias), but one group of coronaviruses switched to using 4-O-Ac-Sias instead (type II). This specificity shift required quasisynchronous adaptations in the Sia-binding sites of both lectin and esterase domains. Previously, a partially disordered crystal structure of a type II HE revealed how the shift in lectin ligand specificity was achieved. How the switch in esterase substrate specificity was realized remained unresolved, however. Here, we present a complete structure of a type II HE with a receptor analog in the catalytic site and identify the mutations underlying the 9-O- to 4-O-Ac-Sia substrate switch. We show that (i) common principles pertaining to the stereochemistry of protein-carbohydrate interactions were at the core of the transition in lectin ligand and esterase substrate specificity; (ii) in consequence, the switch in O-Ac-Sia specificity could be readily accomplished via convergent intramolecular coevolution with only modest architectural changes in lectin and esterase domains; and (iii) a single, inconspicuous Ala-to-Ser substitution in the catalytic site was key to the emergence of the type II HEs. Our findings provide fundamental insights into how proteins "see" sugars and how this affects protein and virus evolution.

Introducing DInaMo: A Package for Calculating Protein Circular Dichroism Using Classical Electromagnetic Theory

The dipole interaction model is a classical electromagnetic theory for calculating circular dichroism (CD) resulting from the π-π* transitions of amides. The theoretical model, pioneered by J. Applequist, is assembled into a package, DInaMo, written in Fortran allowing for treatment of proteins. DInaMo reads Protein Data Bank formatted files of structures generated by molecular mechanics or reconstructed secondary structures. Crystal structures cannot be used directly with DInaMo; they either need to be rebuilt with idealized bond angles and lengths, or they need to be energy minimized to adjust bond lengths and bond angles because it is common for crystal structure geometries to have slightly short bond lengths, and DInaMo is sensitive to this. DInaMo reduces all the amide chromophores to points with anisotropic polarizability and all nonchromophoric aliphatic atoms including hydrogens to points with isotropic polarizability; all other atoms are ignored. By determining the interactions among the chromophoric and nonchromophoric parts of the molecule using empirically derived polarizabilities, the rotational and dipole strengths are determined leading to the calculation of CD. Furthermore, ignoring hydrogens bound to methyl groups is initially explored and proves to be a good approximation. Theoretical calculations on 24 proteins agree with experiment showing bands with similar morphology and maxima.

Comprehensive list of lectins: origins, natures, and carbohydrate specificities

More than 100 years have passed since the first lectin ricin was discovered. Since then, a wide variety of lectins (lect means "select" in Latin) have been isolated from plants, animals, fungi, bacteria, as well as viruses, and their structures and properties have been characterized. At present, as many as 48 protein scaffolds have been identified as functional lectins from the viewpoint of three-dimensional structures as described in this chapter. In this chapter, representative 53 lectins are selected, and their major properties that include hemagglutinating activity, mitogen activity, blood group specificity, molecular weight, metal requirement, and sugar specificities are summarized as a comprehensive table. The list will provide a practically useful, comprehensive list for not only experienced lectin users but also many other non-expert researchers, who are not familiar to lectins and, therefore, have no access to advanced lectin biotechnologies described in other chapters.

Involvement of water in carbohydrate-protein binding: concanavalin A revisited

Ordered water molecules bound to protein surfaces, or in protein-ligand interfaces, are frequently observed by crystallography. The investigation of the impact of such conserved water molecules on protein stability and ligand affinity requires detailed structural, dynamic, and thermodynamic analyses. Several crystal structures of the legume lectin concanavalin A (Con A) bound to closely related carbohydrate ligands show the presence of a conserved water molecule that mediates ligand binding. Experimental thermodynamic and theoretical studies have examined the role of this conserved water in the complexation of Con A with a synthetic analog of the natural trisaccharide, in which a hydroxyethyl side chain replaces the hydroxyl group at the C-2 position in the central mannosyl residue. Molecular modeling earlier indicated (Clarke, C.; Woods, R. J.; Glushka, J.; Cooper, A.; Nutley, M. A.; Boons, G.-J. J. Am. Chem. Soc. 2001, 123, 12238-12247) that the hydroxyl group in this synthetic side chain could occupy a position equivalent to that of the conserved water, and thus might displace it. An interpretation of the experimental thermodynamic data, which was consistent with the displacement of the conserved water, was also presented. The current work reports the crystal structure of Con A with this synthetic ligand and shows that even though the position and interactions of the conserved water are distorted, this key water is not displaced by the hydroxyethyl moiety. This new structural data provides a firm basis for molecular dynamics simulations and thermodynamic integration calculations whose results indicate that differences in van der Waals contacts (insertion energy), rather than electrostatic interactions (charging energy) are fundamentally responsible for the lower affinity of the synthetic ligand. When combined with the new crystallographic data, this study provides a straightforward interpretation for the lower affinity of the synthetic analog; specifically, that it arises primarily from weaker interactions with the protein via the positionally perturbed conserved water. This interpretation is fully consistent with the experimental observations that the free energy of binding is enthalpy driven, that there is both less enthalpic gain and less entropic penalty for binding the synthetic ligand, relative to the natural trisaccharide, and that the entropic component does not arise from releasing an ordered water molecule from the protein surface to the bulk solvent.

Specificity analysis of lectins and antibodies using remodeled glycoproteins

Due to their ability to bind specifically to certain carbohydrate sequences, lectins are a frequently used tool in cytology, histology, and glycan analysis but also offer new options for drug targeting and drug delivery systems. For these and other potential applications, it is necessary to be certain as to the carbohydrate structures interacting with the lectin. Therefore, we used glycoproteins remodeled with glycosyltransferases and glycosidases for testing specificities of lectins from Aleuria aurantia (AAL), Erythrina cristagalli (ECL), Griffonia simplicifolia (GSL I-B(4)), Helix pomatia agglutinin (HPA), Lens culinaris (LCA), Lotus tetragonolobus (LTA), peanut (Arachis hypogaeae) (PNA), Ricinus communis (RCA I), Sambucus nigra (SNA), Vicia villosa (VVA), and wheat germ (Triticum vulgaris) (WGA) as well as reactivities of anti-carbohydrate antibodies (anti-bee venom, anti-horseradish peroxidase [anti-HRP], and anti-Lewis(x)). After enzymatic remodeling, the resulting neoglycoforms display defined carbohydrate sequences and can be used, when spotted on nitrocellulose or in enzyme-linked lectinosorbent assays, to identify the sugar moieties bound by the lectins. Transferrin with its two biantennary complex N-glycans was used as scaffold for gaining diverse N-glycosidic structures, whereas fetuin was modified using glycosidases to test the specificities of lectins toward both N- and O-glycans. In addition, alpha(1)-acid glycoprotein and Schistosoma mansoni egg extract were chosen as controls for lectin interactions with fucosylated glycans (Lewis(x) and core alpha1,3-fucose). Our data complement and expand the existing knowledge about the binding specificity of a range of commercially available lectins.

Structural glycobiology: a game of snakes and ladders

Oligo- and polysaccharides are infamous for being extremely flexible molecules, populating a series of well-defined rotational isomeric states under physiological conditions. Characterization of this heterogeneous conformational ensemble has been a major obstacle impeding high-resolution structure determination of carbohydrates and acting as a bottleneck in the effort to understand the relationship between the carbohydrate structure and function. This challenge has compelled the field to develop and apply theoretical and experimental methods that can explore conformational ensembles by both capturing and deconvoluting the structural and dynamic properties of carbohydrates. This review focuses on computational approaches that have been successfully used in combination with experiment to detail the three-dimensional structure of carbohydrates in a solution and in a complex with proteins. In addition, emerging experimental techniques for three-dimensional structural characterization of carbohydrate-protein complexes and future challenges in the field of structural glycobiology are discussed. The review is divided into five sections: (1) The complexity and plasticity of carbohydrates, (2) Predicting carbohydrate-protein interactions, (3) Calculating relative and absolute binding free energies for carbohydrate-protein complexes, (4) Emerging and evolving techniques for experimental characterization of carbohydrate-protein structures, and (5) Current challenges in structural glycoscience.

Ligand specificity of CS-35, a monoclonal antibody that recognizes mycobacterial lipoarabinomannan: a model system for oligofuranoside-protein recognition

The CS-35 antibody is widely used in the characterization of glycans containing D-arabinofuranose residues, in particular polysaccharides present in the mycobacterial cell wall. A detailed understanding of the combining site of this antibody and the measurement of its binding to different ligands is of interest as this knowledge will have implications in the characterization of arabinofuranose-containing glycoconjugates that are increasingly recognized as important biological molecules. Of even greater significance is that an in-depth study of this carbohydrate-protein interaction will provide insights into the mechanisms by which oligosaccharides containing furanose rings are bound by proteins, an area that has, to date, received little attention. This system has been refractory to X-ray crystallography, and thus we report here a study of the interaction of CS-35 with its ligands using a combination of chemical synthesis, mass spectrometry, titration microcalorimetry, and NMR spectroscopy. Through these investigations we have established that the binding pocket recognizes, as a minimum epitope, a linear tetrasaccharide motif and that the residues at the reducing and non-reducing end of the oligosaccharide are essential for tight binding. The residue at the non-reducing end appears to be bound in an aliphatic pocket, whereas the rest of the tetrasaccharide interacts more strongly with aromatic amino acids.

Amphiphilic N-Glycosyl-thiocarbamoyl Cyclodextrins: Synthesis, Self-Assembly, and Fluorimetry of Recognition by Lens culinaris Lectin

Amphiphilic beta-cyclodextrins have been synthesized bearing hexylthio, dodecylthio, and hexadecylthio chains at the 6-positions and glycosylthiocarbamoyl-oligo(ethylene glycol) units at the 2-positions. The glycosyl residues (alpha-D-mannosyl and beta-L-fucosyl) are intended for cell-targeting. Self-assembly of these new amphiphilic glycosylated cyclodextrins in water to form vesicles was investigated by dynamic light scattering and transmission electron microscopy. Selective binding of the hexylthio assemblies to a protein receptor (Lens culinaris lectin) was confirmed by fluorescence spectroscopy.

The utility of residual dipolar couplings in detecting motion in carbohydrates: application to sucrose

The solution structure and dynamics of sucrose are examined using a combination of NMR residual dipolar coupling and molecular mechanics force fields. It is found that the alignment tensors of the individual rings are different, and that fitting 35 measured residual dipolar couplings to structures with specific phi, psi values indicates the presence of three major conformations: phi, psi=(120 degrees ,270 degrees), (45 degrees, 300 degrees) and (90 degrees ,180 degrees). Furthermore, fitting two structures simultaneously to the 35 residual dipolar couplings results in a substantial improvement in the fits. The existence of multiple conformations having similar stabilities is a strong indication of motion, due to the interconversion among these states. Results from four molecular mechanics force fields are in general agreement with the experimental results. However, there are major disagreements between force fields. Because fits of residual dipolar couplings to structures are dependent on the force field used to calculate the structures, multiple force fields were used to interpret NMR data. It is demonstrated that the pucker of the fructofuranosyl ring affects the calculated potential energy surface, and the fit to the residual dipolar couplings data. Previously published 13C nuclear relaxation results suggesting that sucrose is rigid are not inconsistent with the present results when motional timescales are considered.

Studies of Disaccharide Solvation—Molecular Dynamics versus HPLC Retention

Molecular dynamics simulations have been used to assess the conformational behaviour of seven disaccharides in aqueous solution. Solvation decreased the overall conformational fluctuations of the sugars, compared to in vacuo simulations using a high dielectric constant. The most significant finding was a linear correlation between the experimental chromatographic retention parameter K´ and a molecular modelling parameter based on the next-nearest oxygen–oxygen distances in the disaccharides. The results support previous proposals for a stereospecific hydration model for carbohydrates and demonstrate the utility of a combined experimental/molecular modelling approach to its study.

Crystal structure of Pterocarpus angolensis lectin in complex with glucose, sucrose, and turanose

The crystal structure of the Man/Glc-specific seed lectin from Pterocarpus angolensis was determined in complex with methyl-alpha-d-glucose, sucrose, and turanose. The carbohydrate binding site contains a classic Man/Glc type specificity loop. Its metal binding loop on the other hand is of the long type, different from what is observed in other Man/Glc-specific legume lectins. Glucose binding in the primary binding site is reminiscent of the glucose complexes of concanavalin A and lentil lectin. Sucrose is found to be bound in a conformation similar as seen in the binding site of lentil lectin. A direct hydrogen bond between Ser-137(OG) to Fru(O2) in Pterocarpus angolensis lectin replaces a water-mediated interaction in the equivalent complex of lentil lectin. In the turanose complex, the binding site of the first molecule in the asymmetric unit contains the alphaGlc1-3betaFruf form of furanose while the second molecule contains the alphaGlc1-3betaFrup form in its binding site.

A novel modeling protocol for protein receptors guided by bound-ligand conformation

A novel protocol for protein homology modeling is described, in which uncertainty in protein structure is resolved by applying the criterion that the protein must bind to a complementary ligand. A diverse library of protein models is created and then screened by docking with ligands of known conformation. The more accurate protein models form higher-quality docked complexes, and the quality of the fit is used to select the best models. The effectiveness of this technique with both natural and unnatural ligands is demonstrated by modeling the Fv fragment of an antibody and comparing the results to known crystal structures.

Conformation of glycomimetics in the free and protein-bound state: structural and binding features of the C-glycosyl analogue of the core trisaccharide alpha-D-Man-(1 --> 3)-[alpha-D-Man-(1 --> 6)]-D-Man

The conformational properties of the C-glycosyl analogue of the core trisaccharide alpha-D-Man-(1 --> 3)-[alpha-D-Man-(1 --> 6)]-D-Man in solution have been carefully analyzed by a combination of NMR spectroscopy and time-averaged restrained molecular dynamics. It has been found that both the alpha-1,3- and the alpha-1,6-glycosidic linkages show a major conformational averaging. Unusual Phi ca. 60 degrees orientations for both Phi torsion angles are found. Moreover, a major conformational distinction between the natural compound and the glycomimetic affects to the behavior of the omega(16) torsion angle around the alpha-1 --> 6-linkage. Despite this increased flexibility, the C-glycosyl analogue is recognized by three mannose binding lectins, as shown by NMR (line broadening, TR-NOE, and STD) and surface plasmon resonance (SPR) methods. Moreover, a process of conformational selection takes place, so that these lectins probably bind the glycomimetic similarly to the way they recognize the natural analogue. Depending upon the architecture and extension of the binding site of the lectin, loss or gain of binding affinity with respect to the natural analogue is found.

Crystal structure of the carbohydrate recognition domain of p58/ERGIC-53, a protein involved in glycoprotein export from the endoplasmic reticulum

p58/ERGIC-53 is an animal calcium-dependent lectin that cycles between the endoplasmic reticulum (ER) and the Golgi complex and appears to act as a cargo receptor for a subset of soluble glycoproteins exported from the ER. We have determined the crystal structure of the carbohydrate recognition domain (CRD) of p58, the rat homologue of human ERGIC-53, to 1.46 A resolution. The fold and ligand binding site are most similar to those of leguminous lectins. The structure also resembles that of the CRD of the ER folding chaperone calnexin and the neurexins, a family of non-lectin proteins expressed on neurons. The CRD comprises one concave and one convex beta-sheet packed into a beta-sandwich. The ligand binding site resides in a negatively charged cleft formed by conserved residues. A large surface patch of conserved residues with a putative role in protein-protein interactions and oligomerization lies on the opposite side of the ligand binding site. Together with previous functional data, the structure defines a new and expanding class of calcium-dependent animal lectins and provides a starting point for the understanding of glycoprotein sorting between the ER and the Golgi.

The diluted aqueous solvation of carbohydrates as inferred from molecular dynamics simulations and NMR spectroscopy

The purpose of this paper is to review our understanding of the dilute hydration (aqueous solvation) behaviour of disaccharide compounds. To this end we discuss and scrutinize the results that have been obtained for the three model disaccharides: maltose, sucrose and trehalose from experimental NMR studies and from theoretical molecular dynamics studies in explicit aqueous solutions. The focus is on the description of molecular hydration features that will influence macroscopic entities such as diffusion and relaxation: residence times of hydration waters, hydration numbers and hydration densities. The principles of molecular dynamics simulation are briefly outlined while a detailed presentation is given of the key features that characterise hydration: the solvation of the glycosidic linkage, the radial hydration of the solute, the water density anisotropy around the solute, the residential behaviour of water molecules in the periphery of the solute, and rotational and translational diffusion coefficients. With respect to the use of NMR in characterising the structure and dynamics of the hydration, the hydrodynamic theory of rotational and translational diffusion of biomolecules as well as the use of pulse field gradient spin echo experiments are briefly presented. The NMR-defined rotational diffusion coefficients (D(r)) and the experimentally determined translational diffusion (D(t)) coefficients are reported for 4% (w/w) solutions of sucrose, trehalose and maltose. These results are compared with theoretical data obtained from molecular dynamics simulations of sucrose, trehalose and maltose under identical conditions (concentration, temperature, etc.). With our present level of knowledge we can propose that although carbohydrates share a number of hydration characteristics, evidence is accumulating in support of the notion that it is not the amount or overall hydration but rather the detailed individual carbohydrate-water interaction that is likely to determine carbohydrate structure and functionality.

Stereochemical metrics of lectin-carbohydrate interactions: comparison with protein-protein interfaces

A global census of stereochemical metrics including interface size, hydropathy, amino acid propensities, packing and hydrogen bonding was carried out on 32 x-ray-elucidated structures of lectin-carbohydrate complexes covering eight different lectin families. It is shown that the interactions at primary binding subsites are more efficient than at other subsites. Another salient behavior found for primary subsites was a marked negative correlation between the interface size and the polar surface content. It is noteworthy that this demographic rule is delineated by lectins with unrelated phylogenetic origin, indicating that independent interface architectures have evolved through common optimization paths. The structural properties of lectin-carbohydrate interfaces were compared with those characterizing a set of 32 protein homodimers. Overall, the analysis shows that the stereochemical bases of lectin-carbohydrate and protein-protein interfaces differ drastically from each other. In comparison with protein-protein complexes, lectin-carbohydrate interfaces have superior packing efficiency, better hydrogen bonding stereochemistry, and higher interaction cooperativity. A similar conclusion holds in the comparison with protein-protein heterocomplexes. We propose that the energetic consequence of this better interaction geometry is a larger decrease in free energy per unit of area buried, feature that enables lectins and carbohydrates to form stable complexes with relatively small interface areas. These observations lend support to the emerging notion that systems differing from each other in their stereochemical metrics may rely on different energetic bases.

Structural basis of carbohydrate recognition by lectin II from Ulex europaeus, a protein with a promiscuous carbohydrate-binding site

Protein-carbohydrate interactions are the language of choice for inter- cellular communication. The legume lectins form a large family of homologous proteins that exhibit a wide variety of carbohydrate specificities. The legume lectin family is therefore highly suitable as a model system to study the structural principles of protein-carbohydrate recognition. Until now, structural data are only available for two specificity families: Man/Glc and Gal/GalNAc. No structural data are available for any of the fucose or chitobiose specific lectins. The crystal structure of Ulex europaeus (UEA-II) is the first of a legume lectin belonging to the chitobiose specificity group. The complexes with N-acetylglucosamine, galactose and fucosylgalactose show a promiscuous primary binding site capable of accommodating both N-acetylglucos amine or galactose in the primary binding site. The hydrogen bonding network in these complexes can be considered suboptimal, in agreement with the low affinities of these sugars. In the complexes with chitobiose, lactose and fucosyllactose this suboptimal hydrogen bonding network is compensated by extensive hydrophobic interactions in a Glc/GlcNAc binding subsite. UEA-II thus forms the first example of a legume lectin with a promiscuous binding site and illustrates the importance of hydrophobic interactions in protein-carbohydrate complexes. Together with other known legume lectin crystal structures, it shows how different specificities can be grafted upon a conserved structural framework.

Structure/thermodynamics relationships of lectin-saccharide complexes: the Erythrina corallodendron case

Molecular dynamics (MD) simulations of Erythrina corallodendron lectin binding to a monosaccharide, alpha-galactose, and a disaccharide, N-acetyl lactosamine, have been performed in order to investigate the relationship between structure and thermodynamics. A simulated annealing protocol has been used to generate ensembles of structures for the two complexes, from which both qualitative and quantitative information on binding dynamics have been extracted. The ensembled averaged lectin-saccharide interaction enthalpy is equivalent for both sugars, whereas the calculation based on the X-ray structures does show a difference. Within large statistical errors, the calculated 'binding enthalpy' is also the same for the two systems. These errors arise largely from terms involving solvent and are a typical limitation of current MD simulations. Significant qualitative differences in binding between the two complexes are, however, observed over the ensembles. These could be important for unraveling the structure/thermodynamic relationship. Stated simply, there are a greater number of binding options available to the disaccharide compared to the monosaccharide. The implications of alternative binding states on thermodynamic parameters and the 'breaking of enthalpy-entropy compensation' are discussed. The role of solvent in lectin-saccharide complex formation is suggested to be significant.

An unusual carbohydrate binding site revealed by the structures of two Maackia amurensis lectins complexed with sialic acid-containing oligosaccharides

Seeds from the legume tree Maackia amurensis contain two lectins that can agglutinate different blood cell types. Their specificity toward sialylated oligosaccharides is unique among legume lectins; the leukoagglutinin preferentially binds to sialyllactosamine (alphaNeuAc(2-3)betaGal(1-4)betaGlcNAc), whereas the hemagglutinin displays higher affinity for a disialylated tetrasaccharide (alphaNeuAc(2-3)betaGal(1-3)[alphaNeuAc(2-6)]alphaG alNAc). The three-dimensional structure of the complex between M. amurensis leukoagglutinin and sialyllactose has been determined at 2.75-A resolution using x-ray crystallography. The carbohydrate binding site consists of a deep cleft that accommodates the three carbohydrate residues of the sialyllactose. The central galactose sits in the primary binding site in an orientation that has not been observed previously in other legume lectins. The carboxyl group of sialic acid establishes a salt bridge with a lysine side chain. The glucose residue is very efficiently docked between two tyrosine aromatic rings. The complex between M. amurensis hemagglutinin and a disialylated tetrasaccharide could be modeled from the leukoagglutinin/sialyllactose crystal structure. The substitution of one tyrosine by an alanine residue is responsible for the difference in fine specificity between the two isolectins. Comparison with other legume lectins indicates that oligosaccharide specificity within this family is achieved by the recycling of structural loops in different combinations.

The mean hydration of carbohydrates as studied by normalized two-dimensional radial pair distributions

The hydration of carbohydrates plays a key role in many biological processes. Molecular dynamics simulations provide an effective tool for investigating the hydration of complex solutes such as carbohydrates. In this article we devise an algorithm for the calculation of two-dimensional radial pair distributions describing the probability of finding a water molecule in a site defined by two reference atoms. The normalized 2D radial pair distribution is proposed as an effective tool for investigating and comparing localized or ordered water sites around flexible molecules such as carbohydrates when analyzing molecular dynamics simulations and the utility of 2D radial pair distributions is demonstrated using sucrose as an example. In this relatively simple structure, 2D radial pair distributions were able to characterize and quantify the importance of two unique interresidue hydration sites in which a water molecule is forming a bridge between the glycopyranosyl and fructofuranosyl residues. The approach is proposed to be a valuable tool for comparing and understanding the hydration of flexible biomolecules such as carbohydrates.

The structure and function of a foot-and-mouth disease virus-oligosaccharide receptor complex

Heparan sulfate has an important role in cell entry by foot-and-mouth disease virus (FMDV). We find that subtype O1 FMDV binds this glycosaminoglycan with a high affinity by immobilizing a specific highly abundant motif of sulfated sugars. The binding site is a shallow depression on the virion surface, located at the junction of the three major capsid proteins, VP1, VP2 and VP3. Two pre-formed sulfate-binding sites control receptor specificity. Residue 56 of VP3, an arginine in this virus, is critical to this recognition, forming a key component of both sites. This residue is a histidine in field isolates of the virus, switching to an arginine in adaptation to tissue culture, forming the high affinity heparan sulfate-binding site. We postulate that this site is a conserved feature of FMDVs, such that in the infected animal there is a biological advantage to low affinity, or more selective, interactions with glycosaminoglycan receptors.

Three dimensional structure of the soybean agglutinin Gal/GalNAc complexes by homology modeling

Complexes of soybean agglutinin (SBA) with galactose (Gal) and N-acetyl galactosamine (GalNAc) have been modeled based on its homology to erythrina corallodendron (EcorL) lectin. The three dimensional structure of SBA-Gal modeled with homology techniques agrees well with SBA-(beta-LacNAc)2Gal-R complex determined by X-ray crystallographic techniques at the beta-sheet regions and the regions where Ca2+ and Mn2+ ions bind. However, significant deviations have been observed between the modeled and the X-ray structures, particularly at the loop regions where the polypeptide chain could not be unequivocally traced in the X-ray structure. The hydrogen bonding scheme, predicted from the homology model, shows that the invariant residues i.e. Asp, Gly, Asn, and aromatic residues (Phe) found in all other legume lectins, bind Gal, slightly in a different way than reported in X-ray structure of SBA-pentasaccharide complex. The higher binding affinity of GalNAc over Gal to SBA is due to additional hydrophobic interactions with Tyr107 rather than a hydrogen bond between N-acetamide group of the sugar and the side chain of Asp88 as suggested from X-ray crystal structure studies. Our modeling also suggest that the variation in the length of the loop D observed among galactose binding legume lectins may not have any effect on the binding of sugar at the monosaccharide specific site of the lectins. Soybean agglutinin (SBA) is a member of the leguminous family of lectins. They generally possess a single carbohydrate binding site, besides the tightly bound Ca2+ and Mn2+ ions which are required for their carbohydrate binding activity. They possess a high degree of sequence homology and about 50% of the amino acid residues are invariant. Some of these invariant amino acid residues are involved in the binding of sugar moieties and in metal ion coordination. X-ray crystallographic studies showed that their three-dimensional structures are very similar, though they differ in their carbohydrate binding specificity (1-6). Three of the invariant residues Asp, Gly, and Asn, besides an aromatic residue (Phe or Tyr), are involved in carbohydrate binding. Independent of their sugar specificity, these four residues in legume lectins provide the basic frame for the sugar to bind.

Legume lectin structure

The legume lectins are a large family of homologous carbohydrate binding proteins that are found mainly in the seeds of most legume plants. Despite their strong similarity on the level of their amino acid sequences and tertiary structures, their carbohydrate specificities and quaternary structures vary widely. In this review we will focus on the structural features of legume lectins and their complexes with carbohydrates. These will be discussed in the light of recent mutagenesis results when appropriate. Monosaccharide specificity seems to be achieved by the use of a conserved core of residues that hydrogen bond to the sugar, and a variable loop that determines the exact shape of the monosaccharide binding site. The higher affinity for particular oligosaccharides and monosaccharides containing a hydrophobic aglycon results mainly from a few distinct subsites next to the monosaccharide binding site. These subsites consist of a small number of variable residues and are found in both the mannose and galactose specificity groups. The quaternary structures of these proteins form the basis of a higher level of specificity, where the spacing between individual epitopes of multivalent carbohydrates becomes important. This results in homogeneous cross-linked lattices even in mixed precipitation systems, and is of relevance for their effects on the biological activities of cells such as mitogenic responses. Quaternary structure is also thought to play an important role in the high affinity interaction between some legume lectins and adenine and a series of adenine-derived plant hormones. The molecular basis of the variation in quaternary structure in this group of proteins is poorly understood.

Computational carbohydrate chemistry: what theoretical methods can tell us

Computational methods have had a long history of application to carbohydrate systems and their development in this regard is discussed. The conformational analysis of carbohydrates differs in several ways from that of other biomolecules. Many glycans appear to exhibit numerous conformations coexisting in solution at room temperature and a conformational analysis of a carbohydrate must address both spatial and temporal properties. When solution nuclear magnetic resonance data are used for comparison, the simulation must give rise to ensemble-averaged properties. In contrast, when comparing to experimental data obtained from crystal structures a simulation of a crystal lattice, rather than of an isolated molecule, is appropriate. Molecular dynamics simulations are well suited for such condensed phase modeling. Interactions between carbohydrates and other biological macromolecules are also amenable to computational approaches. Having obtained a three-dimensional structure of the receptor protein, it is possible to model with accuracy the conformation of the carbohydrate in the complex. An example of the application of free energy perturbation simulations to the prediction of carbohydrate-protein binding energies is presented.

The how and why of protein-carbohydrate interaction: a primer to the theoretical concept and a guide to application in drug design

The common principles of molecular recognition with cooperative or bidentate hydrogen bonds, dispersion forces and hydrophobic packing govern the specificity of protein-carbohydrate interaction. Enthalpy/entropy-compensation is also valid, maintaining KD-values in the range of 30 mM to 200 nM. The individual contributions of the enthalpic and entropic factors which originate from the receptor, the ligand and/or the solvent to the overall free energy change can at least be estimated by a combination of computer-assisted molecular modeling, NMR spectroscopy of the reactants before and after complex formation and thermodynamic measurements. The delineation of adaptable parameters such as ligand or receptor side chain flexibility points to a route to practicable guidelines for a rational design of mimetics in glycosciences.

Structure of the sucrose-specific porin ScrY from Salmonella typhimurium and its complex with sucrose

The X-ray structure of a sucrose-specific porin (ScrY) from Salmonella typhimurium has been determined by multiple isomorphous replacement at 2.4 A resolution both in its uncomplexed form and with bound sucrose. ScrY is a noncrystallographic trimer of identical subunits, each with 413 structurally well-defined amino acids. A monomer is built up of 18 anti-parallel beta-strands surrounding a hydrophilic pore, with a topology closely similar to that of maltoporin. Two non-overlapping sucrose-binding sites were identified in difference Fourier maps. The higher permeability for sucrose of ScrY as compared to maltoporin is mainly accounted for by differences in their pore-lining residues.

Internal motions and hydration of sucrose in a diluted water solution

This paper represents a synthesis of our most recent work on the hydration, internal and overall dynamics of sucrose in a diluted water solution. The studies were carried out as a total ensemble of 1.2 nanosecond condensed phase molecular dynamics trajectories. In this study the focus is on a 500 ps trajectory starting with the solute in the crystalline conformation. The presence of water was found to significantly alter the accessible conformational space of the solute. All potential intra-molecular hydrogen bonds were found to be exchanged to surrounding water molecules and the simulations suggest that the sucrose conformation is stabilized by the dynamic presence of two interring bridging water molecules: O-2g...Ow...O-3f and O-2g...Ow...O-lf. The overall shift in conformation of the solute induced by the presence of water was found to improve the theoretical models of experimental traits. It is demonstrated that the hydration structure and the internal and overall motions of sucrose compare extremely well with NMR data such as glycosidic heteronuclear coupling constants and the molecular tumbling time, with X-ray data of two partially hydrated sucrose structures in a protein complex and with translational diffusion coefficients and hydration numbers established from experimental studies.

Conformational analysis of biantennary glycans and molecular modeling of their complexes with lentil lectin

Some mannose-binding legume lectins show higher affinity for fucosylated glycans than for glycans without fucose. These lectins possess a secondary binding site. Owing to the possibility of additional fucose binding, oligosaccharides adopt different conformations depending on whether they contain fucose or not. To study these conformational differences, complexes of fucosylated and unfucosylated glycans with Lens culinaris lectin have been modeled. Starting points were X-ray structures of lentil lectin and complexes of the homologous Lathyrus ochrus lectin. The SYBYL molecular modeling package with the TRIPOS force field was used. Two different models were built, displaying in both a network of hydrogen bonds between the saccharide and the binding site. Furthermore, to compare the free and bound ligand, conformational analysis in the free state has been performed. A complete analysis of all possible disaccharide fragments has been performed using the MM3 force field. A CICADA analysis employing the same force field was carried out to study the complete oligosaccharide. Low-energy conformers found by CICADA were clustered in conformational families and analyzed in terms of flexibility and rotational barriers. All values of glycosidic torsion angles are in the range as calculated by MM3 for the disaccharides.

Applications of nuclear magnetic resonance spectroscopy and molecular modeling to the study of protein-carbohydrate interactions

This work provides an overview of the applications of NMR to the study of protein-carbohydrate interactions. The use of TR-NOE experiments in this context is given. In particular, the study of Ricin/lactose and Hevein/chitobiose complexes is detailed.

Knowledge-based modeling of a legume lectin and docking of the carbohydrate ligand: the Ulex europaeus lectin I and its interaction with fucose

Ulex europaeus isolectin I is specific for fucose-containing oligosaccharide such as H type 2 trisaccharide alpha-L-Fuc (1-->2) beta-D-Gal (1-->4) beta-D-GlcNAc. Several legume lectins have been crystallized and modeled, but no structural data are available concerning such fucose-binding lectin. The three-dimensional structure of Ulex europaeus isolectin I has been constructed using seven legume lectins for which high-resolution crystal structures were available. Some conserved water molecules, as well as the structural cations, were taken into account for building the model. In the predicted binding site, the most probable locations of the secondary hydroxyl groups were determined using the GRID method. Several possible orientations could be determined for a fucose residue. All of the four possible conformations compatible with energy calculations display several hydrogen bonds with Asp-87 and Ser-132 and a stacking interaction with Tyr-220 and Phe-136. In two orientations, the O-3 and O-4 hydroxyl groups of fucose are the most buried ones, whereas two other, the O-2 and O-3 hydroxyl groups are at the bottom of the site. Possible docking modes are also studied by analysis of the hydrophobic and hydrophilic surfaces for both the ligand and the protein. The SCORE method allows for a quantitative evaluation of the complementarity of these surfaces, on the basis of molecular lipophilicity calculations. The predictions presented here are compared with known biochemical data.

A structure of the complex between concanavalin A and methyl-3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside reveals two binding modes

The structure of concanavalin A in complex with the trimannoside methyl-3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside has been determined in a novel space group. In three of the four subunits of the concanavalin A tetramer, the interactions between the protein and the bound saccharide are essentially identical to those reported previously by other authors (Naismith, J. H., and Field, R. A. (1996) J. Biol. Chem. 271, 972-976). In the fourth subunit, however, the alpha1-->3 linkage has a different conformation, resulting in a different part of the alpha1-->3-linked mannose interacting with essentially the same surface of the protein. Furthermore, significant differences are observed in the quaternary associations of the subunits compared with the saccharide-free structures and other carbohydrate complexes, suggesting that the concanavalin A tetramer is a rather flexible entity.