Combined NMR and molecular modeling study of an iduronic acid-containing trisaccharide related to antithrombotic heparin fragments

An iduronic acid-containing trisaccharide, methyl-O-(4-O-methyl-2,3,6-tri-O-sulfo-alpha-D-glucopyranosyl-(1-->4)-O- (2-O-sulfo-alpha-L-idopyranosyluronic acid)-(1-->4)-O-2,6-di-O-sulfo-alpha-D-glucopyranoside, related to antithrombotic heparin fragments has been subjected to a combined NMR and molecular modeling investigation. The conformational behavior of the two constituting disaccharide segments was investigated using a systematic grid search approach with the MM3 force field along with the proper parameters for the sulfate ester group. The exploration of the potential energy surfaces of the trisaccharide was performed through the use of the CICADA methods interfaced with the MM3 force field. In all cases, the 2-O-sulfo-alpha-L-iduronate moiety was given the three favored ring conformations (1)C4, (4)C1, and (2)S0. Conformations were clustered into families, four of which are likely to exhibit significant occupancy in solution. The different low-energy conformational families display different orientations at the glycosidic linkages and/or different ring shapes for the iduronate ring. The (2)S0 conformation is the major one for the 2-O-sulfo-alpha-L-iduronate but is still in equilibrium with the (1)C4 ring shape. The occurrence of such a conformational equilibrium in solution was probed via high-resolution NMR spectroscopy through measurements of coupling constants and NOE build-up. These results are in keeping with the observation that 2-O-sulfated pentasaccharides display a similar affinity for antithrombin III as their 2-N-sulfated counterparts.

Transferred nuclear Overhauser enhancement (NOE) and rotating-frame NOE experiments reflect the size of the bound segment of the Forssman pentasaccharide in the binding site of Dolichos biflorus lectin

A complex between the Forssman pentasaccharide alpha-D-GalNAc-(1-->3)-beta-D-GalNAc-(1-->3)-alpha-D-Gal-(1-->4)-beta-D- Gal-(1-->4)-D-Glc and the seed lectin from Dolichos biflorus was studied using transfer-NOESY and transfer rotating frame NOE spectroscopy (ROESY) experiments. The evolution of transferred NOEs and ROEs as a function of the pentasaccharide/lectin ratio was different for the non-reducing disaccharide moiety alpha-D-GalNAc-(1-->3)-beta-D-GalNac compared to the rest of the molecule, which reflects distinct relaxation properties and effects of exchange broadening of the corresponding ligand resonances. Significantly, several intermolecular transferred NOEs were observed between protons of the nonreducing disaccharide moiety alpha-D-GalNAc-(1-->3)-beta-D-GalNAc and aliphatic as well as aromatic amino acid side chain protons in the binding pocket of the lectin. It is concluded that the non-reducing disaccharide fragment is buried in the lectin-binding pocket, whereas the reducing trisaccharide portion alpha-D-Gal-(1-->4)-beta-D-Gal-(1-->4)-D-Glc has no immediate contacts with the protein. The experimental transfer NOE data were qualitatively compared to theoretical proton-proton distances from a model that was based on a previous homology modeling study of a complex between the disaccharide fragment alpha-D-GalNAc-(1-->3)-beta-D-GalNAc and D. biflorus lectin. It was found that all intermolecular transferred NOEs matched short interatomic distances between ligand protons and aliphatic or aromatic amino acid side chain protons predicted by the theoretical model.

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.

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.

How do antibodies and lectins recognize histo-blood group antigens? A 3D-QSAR study by comparative molecular field analysis (CoMFA)

The cross-reaction patterns of nine antibodies and three lectins against 12 H type 2 related oligosaccharides have been analysed by means of 3D-QSAR study. Three-dimensional descriptors of the molecular properties have been used in comparative molecular field analysis (CoMFA). Three different alignments were considered for the oligosaccharides. One, based on the superimposition of the oligosaccharide core, could be correlated to most of the antibody activities. A second alignment, based on a superimposition of the fucose residue, had to be taken into account for explaining the binding properties of Ulex europaeus isolectin I. Analysis of the QSAR data gives indications on the carbohydrate epitopes essential for antibody recognition and yields some insights about the nature of the molecular recognition. This study complements previous biochemical estimates of the H type 2 related oligosaccharide binding areas (Mollicone, R.; Cailleau, A.; Imberty, A.; Gane, P.; Pérez, S.; Oriol, R. Glycoconj. J. 1996, 13, 263-271).

Solution conformations of pectin polysaccharides: determination of chain characteristics by small angle neutron scattering, viscometry, and molecular modeling

The solution behavior of pectin polysaccharides has been investigated by small angle neutron scattering (SANS), viscosimetric, and molecular modeling studies. The samples used in the experimental study were obtained from apple and citrus and had degrees of methylation ranging from 28 to 73%, with a rhamnose content lying between 0.6 and 2.2%. Persistence lengths, derived from intrinsic viscosity measurements, ranged from 59 to 126 A, whereas those derived by SANS were between 45 and 75 A. These values correspond to 10-17 monomer units. The modeling simulations were performed for both homogalacturonan itself and homogalacturonan carrying various degrees of rhamnose inserts (rhamnogalacturonan). This required the evaluation of the accessible conformational space for the eight disaccharides that represent the constituent repeating segments of the homogalacturonan and rhamnogalacturonan polysaccharides. For each dimer, complete conformational analysis was accomplished using the flexible residue method of the MM3 molecular mechanics procedure and the results used to access the configurational statistics of representative pectic polysaccharide chains. For homogalacturonan, an extended chain conformation having a persistence length of 135 A (corresponding to 30 monomers) was predicted. The inclusion of varying amounts of rhamnose units (5-25%) in the model in strict alternating sequence with galacturonate residues (equivalent to the rhamnogalacturonan "hairy region" chains) only slightly reduced the calculated persistence length. The extended overall chain conformation remained relatively unchanged as a consequence of the self-cancellation of the kinking effects of successive paired rhamnose units.

Conformational analysis of blood group A trisaccharide in solution and in the binding site of Dolichos biflorus lectin using transient and transferred nuclear Overhauser enhancement (NOE) and rotating-frame NOE experiments

The present study is concerned with the elucidation of the conformation of the blood group A trisaccharide (alpha-D-GalNAc(1-3)[alpha-L-Fuc(1-->2)] beta-D-Gal-O-R) in the combining site of Dolichos biflorus seed lectin by use of 400-MHz and 600-MHz NMR spectroscopy. D. biflorus lectin displays a unique specificity for GalNAc residues. It occurs in solution as a tetrameric assembly having a molecular mass of 110 kDa, with two carbohydrate-binding sites per molecule. First, NOE build-up curves were obtained for the free blood group A trisaccharide from one-dimensional transient NOE experiments. Simulated NOE build-up curves were constructed from an ensemble of low-energy conformers derived from previous investigations. The comparison of theoretical and experimental data indicates that an equilibrium between two families of low-energy conformers most likely reflects the solution behavior of the trisaccharide in solution. Two-dimensional transferred NOE and rotating-frame enhancements (ROE) were subsequently measured for the trisaccharide complexed with the D. biflorus seed lectin. In addition to the NOEs observed for the free trisaccharide, the transferred NOESY spectrum showed several new NOEs that were identified as spin diffusion using a rotating-frame NOESY (ROESY) experiment. Experimental interglycosidic transferred nuclear Overhauser effect (TRNOE) build-up curves were compared to theoretical curves calculated for both low-energy conformers located in the D. biflorus lectin-binding site. Calculations of theoretical TRNOE were performed using a combination of the full relaxation matrix and the protein-ligand exchange matrix. Comparison between experimental and simulated TRNOE volumes leads to the conclusion that one conformation of blood group A trisaccharide is selected upon binding by D. biflorus lectin.

Predicting helical structures of the exopolysaccharide produced by Lactobacillus sake 0- 1

The viscous exopolysaccharide (EPS) produced by Lactobacillus sake 0- 1 is a high molecular mass polymer (Mm 6 x 10(6) Da) consisting of pentasaccharide repeating units with a composition of D-glucose, L-rhamnose, and sn-glycerol 3-phosphate in molar ratios of 3:2:1. One of the rhamnose residues in the repeating unit is partially 2-O-acetylated. The O-deacetylated, deglycerophosphorylated EPS has been investigated by molecular mechanics calculations. A complete conformational analysis of each of the constituent disaccharide fragments has been performed using the flexible residue approach with the MM3(92) force field. Furthermore, using the same force field, CICADA analyses were accomplished on hexa- and octasaccharide substructure of the polysaccharide. Based on these analyses, insight was obtained into nine conformational minima for the polysaccharide. The low energy conformations found by CICADA were extrapolated to regular polysaccharide structures using a polysaccharide builder program. The generated helices exhibit either 2-fold or 3- or 4-fold right-handed chiralities, and in each case the helices are highly extended.

Crystal and molecular structure of a histo-blood group antigen involved in cell adhesion: the Lewis x trisaccharide

This work describes the first crystal structure ever reported of a histo-blood group carbohydrate antigen: Le(x). This study provides a detailed description of the conformation of two crystallographic independent molecules in a highly hydrated environment along with their hydrogen bonding properties and packing features. Some interactions observed between adjacent trisaccharides can provide the basis for involvement of Le(x)-Le(x) interactions in cell-cell adhesion.

Recognition of the blood group H type 2 trisaccharide epitope by 28 monoclonal antibodies and three lectins

The patterns of cross-reaction of 30 monoclonal antibodies and three lectins were determined by ELISA with 21 ABH, Ii or Lewis related synthetic oligosaccharides coupled to bovine serum albumin. At least seven main groups of cross-reactive patterns were identified among the antibodies, plus several isolated antibodies which had intermediate patterns between two of the main antibody groups. The three lectins had different cross-reaction patterns, Galactia tenuiflora was different from all the antibodies, Ulex europaeus lectin 1 and Lotus tetragonolobus were similar, but not identical to groups III and V of antibodies respectively. The anti-H antibodies cross-reacting with A type 2 gave similar agglutination scores with all the normal ABO erythrocytes, while the anti-H antibodies not cross-reacting with A type 2 reacted with different scores: O > A2 > A2B > B > A1 > A1B > O(h), suggesting that these antibodies react better with the free H epitopes and do not recognize the H in A or B epitopes. Based on the ELISA and agglutination results and the lowest energy conformations of each oligosaccharide obtained by computer modelling, the most probable oligosaccharide surface areas recognized by each antibody main group are illustrated.

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

The conformational features of sucrose in the combining site of lentil lectin have been characterized through elucidation of a crystalline complex at 1.9-A resolution, transferred nuclear Overhauser effect experiments performed at 600 Mhz, and molecular modeling. In the crystal, the lentil lectin dimer binds one sucrose molecule per monomer. The locations of 229 water molecules have been identified. NMR experiments have provided 11 transferred NOEs. In parallel, the docking study and conformational analysis of sucrose in the combining site of lentil lectin indicate that three different conformations can be accommodated. Of these, the orientation with lowest energy is identical with the one observed in the crystalline complex and provides good agreement with the observed transferred NOEs. These structural investigations indicate that the bound sucrose has a unique conformation for the glycosidic linkage, close to the one observed in crystalline sucrose, whereas the fructofuranose ring remains relatively flexible and does not exhibit any strong interaction with the protein. Major differences in the hydrogen bonding network of sucrose are found. None of the two inter-residue hydrogen bonds in crystalline sucrose are conserved in the complex with the lectin. Instead, a water molecule bridges hydroxyl groups O2-g and O3-f of sucrose.

Molecular modelling of the interaction between the catalytic site of pig pancreatic alpha-amylase and amylose fragments

A stereo chemical refinement of the crystalline complex between porcine pancreatic alpha-amylase and a pseudopentasaccharide from the amylostatin family has been performed through molecular mechanics calculations, using a set of parameters appropriate for protein and protein-carbohydrate interactions. The refinement provided a starting point for docking a maltopentaose moiety within the catalytic site, in the absence of water. A thorough exploration of the different orientations and conformations of maltopentaose established the sense of binding of the amylosic substrate in the amylase cleft. After optimising the geometry of the binding site, the conformations adopted by the four contiguous linkages could be rationalised by considering the environment, either hydrophobic or hydrophilic, of the different glucose moieties. Seemingly, details of the non-bonded interactions (hydrogen bonds, van der Waals and stacking interactions) that underlie this molecular recognition have been established. In particular, it was confirmed that the three acidic amino acids of the catalytic site (Asp197, Asp300 and Glu233) are close to their glucosidic target, and that there is no steric reason to propose an alteration of the 4C1 conformation of the glucose residue prior to hydrolysis. However, in the absence of water molecules, it is difficult to elucidate the details of the catalysis. Additional macroscopic information has been gained, such as the impossibility to fit a double-helical arrangement of amylose chains in the amylasic cleft. This explains why some native starches containing such motifs resist amylolytic enzymes. Tentative models involving longer amylosic chains have been elaborated, which extend our knowledge of the interaction and orientation of starch fragments in the vicinity of the hydrolytic sites.

Stereochemistry of the N-glycosylation sites in glycoproteins

The stereochemical features displayed by the N-glycosidic linkage in crystalline N-linked glycoproteins are analyzed. From the statistical analysis of 44 different glycosylation sites belonging to 26 glycoproteins of the Brookhaven Protein Data Bank, a mean standard geometry for the GlcNAc moiety, along with a rationalization of its conformational behavior, can be proposed. As for the glycopeptide linkage, the distribution of observed conformations has been analyzed on the basis of molecular mechanics calculations. The rotamer distribution of the Asn side chains conforms to that observed on non-glycosylated structures, and it agrees with the pattern of flexible conformations gathered from NMR measurements. In characterizing the protein-glycan interactions, some hydrogen bonds occur. Stacking between the amphiphilic moiety of the glycan and some surrounding aromatic, or at least hydrophobic, amino acid residues is also found. When looking at the secondary structure of the glycosylated peptide, only 25% of the glycosylation sites correspond to situations where Asn is located at the top of a beta-turn. Other types of secondary structure exist which fulfill the spatial requirement of having the glycan exposed at the surface of the protein. These data can be compared with the most recent studies on the peptide conformation which would be required for glycosylation.

Computer simulation of histo-blood group oligosaccharides: energy maps of all constituting disaccharides and potential energy surfaces of 14 ABH and Lewis carbohydrate antigens

The three-dimensional structures of fourteen histo-blood groups carbohydrate antigens have been established through a combination of molecular mechanics and conformational searching methods. The conformational space available for each disaccharide, constituents of these determinants, has been throroughly characterized. The results have been organized in a data bank fashion. Larger relatives, i.e. 14 tri- and tetrasaccharides of histo-blood group antigens, have been modelled using a different method for exploring the complex potential energy surface. This approach is aimed at establishing all the possible families of conformations, along with the conformational pathways. Different conformational behaviours are exhibited by these oligosaccharides. Some of them, i.e. Le(x) and Le(y) tri and tetrasaccharides, are very rigid; 99% of their populations belong to the same conformational family. Others, like H type 1, H type 2 or H type 6 oligosaccharides, are essentially rigid, but a secondary conformational family, corresponding to 3-4% of the total population, can arise. Finally, the H types 3 and 4 trisaccharides, and the A type 1 and A type 2 tetrasaccharides are predicted to behave rather flexibly. The information gathered in the present investigation has been used to analyse the body of experimental evidence, either physical or biological, available for this series of carbohydrate antigens. Of special interest are the several different alignments that can be proposed for these molecules. They yield a realistic definition of the three-dimensional features of the epitopes thereby providing essential information about how carbohydrate antigens are recognized by proteins.

The monosaccharide binding site of lentil lectin: an X-ray and molecular modelling study

The X-ray crystal structure of lentil lectin in complex with alpha-D-glucopyranose has been determined by molecular replacement and refined to an R-value of 0.20 at 3.0 A resolution. The glucose interacts with the protein in a manner similar to that found in the mannose complexes of concanavalin A, pea lectin and isolectin I from Lathyrus ochrus. The complex is stabilized by a network of hydrogen bonds involving the carbohydrate oxygens O6, O4, O3 and O5. In addition, the alpha-D-glucopyranose residue makes van der Waals contacts with the protein, involving the phenyl ring of Phe123 beta. The overall structure of lentil lectin, at this resolution, does not differ significantly from the highly refined structures of the uncomplexed lectin. Molecular docking studies were performed with mannose and its 2-O and 3-O-m-nitro-benzyl derivatives to explain their high affinity binding. The interactions of the modelled mannose with lentil lectin agree well with those observed experimentally for the protein-carbohydrate complex. The highly flexible Me-2-O-(m-nitro-benzyl)-alpha-D-mannopyranoside and Me-3-O-(m-nitro-benzyl)-alpha-D-mannopyranoside become conformationally restricted upon binding to lentil lectin. For best orientations of the two substrates in the combining site, the loss of entropy is accompanied by the formation of a strong hydrogen bond between the nitro group and one amino acid, Gly97 beta and Asn125 beta, respectively, along with the establishment of van der Waals interactions between the benzyl group and the aromatic amino acids Tyr100 beta and Trp128 beta.

Molecular modelling of the Dolichos biflorus seed lectin and its specific interactions with carbohydrates: alpha-D-N-acetyl-galactosamine, Forssman disaccharide and blood group A trisaccharide

The three-dimensional structure of Dolichos biflorus seed lectin has been constructed using five legume lectins for which high resolution crystal structures were available. The validity of the resulting model has been thoroughly investigated. Final structure optimization was conducted for the lectin complexed with alpha GalNAc, providing thereby the first three-dimensional structure of lectin/GalNAc complex. The role of the N-acetyl group was clearly evidenced by the occurrence of a strong hydrogen bond between the protein and the carbonyl oxygen of the carbohydrate and by hydrophobic interaction between the methyl group and aromatic amino acids. Since the lectin specificity is maximum for the Forssman disaccharide alpha GalNAc(1-3) beta GalNAc-O-Me and the blood group A trisaccharide alpha GalNAc(1-3)[alpha Fuc(1-2)] beta Gal-O-Me, the complexes with these oligosaccharides have been also modelled.

Molecular modelling of protein-carbohydrate interactions. Understanding the specificities of two legume lectins towards oligosaccharides

By means of a series of new molecular modelling tools, the conformational behaviour of mannose-containing di- and trisaccharides bound to either concanavalin A or Lathyrus ochrus isolectin I (LOLI) has been assessed. Tools for estimating and analysing either the 'rigid' or the 'relaxed' potential energy surfaces, representing the conformational space available for carbohydrates once interacting with lectins, are reported for the first time. Restrictions of conformational space are predicted to occur with different magnitudes, depending on the nature of the glycosidic linkages, as well as the size of the carbohydrates. Results from these molecular modelling studies are compared to existing structural data. Not only could the observed conformations and orientations of carbohydrates in crystalline lectin-oligosaccharides complexes be reproduced, but several other likely situations were also predicted to occur. Entropy calculations have been performed for comparison with experimental thermodynamics data. The results of the stimulation can also help giving an explanation of some observed affinity constants at the molecular level.