Aromatic-carbohydrate interactions: an NMR and computational study of model systems

Three-dimensional potential energy surface of selected carbohydrates' CH/π dispersion interactions calculated by high-level quantum mechanical methods

In this study we present the first systematic computational three-dimensional scan of carbohydrate hydrophobic patches for the ability to interact through CH/π dispersion interactions. The carbohydrates β-d-glucopyranose, β-d-mannopyranose and α-l-fucopyranose were studied in a complex with a benzene molecule, which served as a model of the CH/π interaction in carbohydrate/protein complexes. The 3D relaxed scans were performed at the SCC-DFTB-D level with 3 757 grid points for both carbohydrate hydrophobic sides. The interaction energy of all grid points was recalculated at the DFT-D BP/def2-TZVPP level. The results obtained clearly show highly delimited and separated areas around each CH group, with an interaction energy up to -5.40 kcal mol(-1) . The results also show that with increasing H⋅⋅⋅π distance these delimited areas merge and form one larger region, which covers all hydrogen atoms on that specific carbohydrate side. Simultaneously, the interaction becomes weaker with an energy of -2.5 kcal mol(-1) . All local energy minima were optimized at the DFT-D BP/def2-TZVPP level and the interaction energies of these complexes were refined by use of the high-level ab initio computation at the CCSD(T)/CBS level. Results obtained from the optimization suggest that the CH group hydrogen atoms are not equivalent and the interaction energy at the CCSD(T)/CBS level range from -3.54 to -5.40 kcal mol(-1) . These results also reveal that the optimal H⋅⋅⋅π distance for the CH/π dispersion interaction is approximately (2.310±0.030) Å, and the angle defined as carbon-hydrogen-benzene geometrical centre is (180±30)°. These results reveal that whereas the dispersion interactions with the lowest interaction energies are quite strictly located in space, the slightly higher interaction energy regions adopt a much larger space.

Quantification of binding affinities of essential sugars with a tryptophan analogue and the ubiquitous role of C-H···π interactions

The role of noncovalent interactions in carbohydrate recognition by aromatic amino acids has long been reported. To develop a molecular understanding of noncovalent interactions in the recognition process, we have examined a series of binary complexes between 3-methylindole (3-MeIn) and sugars. In particular, the geometries and binding affinities of 3-MeIn with α/β-D-glucose, β-D-galactose, α-D-mannose and α/β-L-fucose are obtained using the MP2(full)/6-31G(d,p) and the M06/TZV2D//MP2/6-31G(d,p) level of theories. The conventional hydrogen bonding such as N-H···O and C-H···O as well as nonconventional O-H···π and C-H···π type of interactions is, in general, identified as responsible for the moderately strong interaction energies. Large variations in the position-orientations of 3-MeIn with respect to saccharide are noticed, within the same sugar family, as well as across different sugar series. Furthermore, complexes with large differences in their geometries are recognized as capable of exhibiting very similar interaction energies, underscoring the significance of exhaustive conformation sampling, as carried out in the present study. These observations are readily attributed to the differences in the efficiency of the type of interactions enlisted above. The highest and lowest interaction energies, upon inclusion of 50% BSSE correction, are found to be -16.02 and -6.22 kcal mol(-1), respectively, for α-D-glucose (1a) and α-L-fucose (5j). While more number of prominent conventional hydrogen bonding contacts remains as a characteristic feature of the strongly bound complexes, the lower end of the interaction energy spectrum is dominated by multiple C-H···π interactions. The complexes exhibiting as many as four C-H···π contacts are identified in the case of α/β-D-glucose, β-D-galactose, and α/β-L-fucose with an interaction energy hovering around -8 kcal mol(-1). The presence of effective C-H···π interactions is found to be dependent on the saccharide configuration as well as the area of the apolar patch constituted by the C-H groups. The study offers a comprehensive set of binary complexes, across different saccharides, which serves as an illustration of the significance and ubiquitous nature of C-H···π interactions in carbohydrate binding in saccharide-protein complexes.

Free energy surfaces for the interaction of D-glucose with planar aromatic groups in aqueous solution

Multidimensional potentials of mean force for the interactions in aqueous solution of both anomers of D-glucopyranose with two planar aromatic molecules, indole and para-methyl-phenol, have been calculated using molecular dynamics simulations with umbrella sampling and were subsequently used to estimate binding free energies. Indole and para-methyl-phenol serve as models for the side chains of the amino acids tryptophan and tyrosine, respectively. In all cases, a weak affinity between the glucose molecules and the flat aromatic surfaces was found. The global minimum for these interactions was found to be for the case when the pseudoplanar face of β-D-glucopyranose is stacked against the planar surfaces of the aromatic residues. The calculated binding free energies are in good agreement with both experiment and previous simulations. The multidimensional free energy maps suggest a mechanism that could lend kinetic stability to the complexes formed by sugars bound to sugar-binding proteins.

Structural motif screening reveals a novel, conserved carbohydrate-binding surface in the pathogenesis-related protein PR-5d

BACKGROUND

Aromatic amino acids play a critical role in protein-glycan interactions. Clusters of surface aromatic residues and their features may therefore be useful in distinguishing glycan-binding sites as well as predicting novel glycan-binding proteins. In this work, a structural bioinformatics approach was used to screen the Protein Data Bank (PDB) for coplanar aromatic motifs similar to those found in known glycan-binding proteins.

RESULTS

The proteins identified in the screen were significantly associated with carbohydrate-related functions according to gene ontology (GO) enrichment analysis, and predicted motifs were found frequently within novel folds and glycan-binding sites not included in the training set. In addition to numerous binding sites predicted in structural genomics proteins of unknown function, one novel prediction was a surface motif (W34/W36/W192) in the tobacco pathogenesis-related protein, PR-5d. Phylogenetic analysis revealed that the surface motif is exclusive to a subfamily of PR-5 proteins from the Solanaceae family of plants, and is absent completely in more distant homologs. To confirm PR-5d's insoluble-polysaccharide binding activity, a cellulose-pulldown assay of tobacco proteins was performed and PR-5d was identified in the cellulose-binding fraction by mass spectrometry.

CONCLUSIONS

Based on the combined results, we propose that the putative binding site in PR-5d may be an evolutionary adaptation of Solanaceae plants including potato, tomato, and tobacco, towards defense against cellulose-containing pathogens such as species of the deadly oomycete genus, Phytophthora. More generally, the results demonstrate that coplanar aromatic clusters on protein surfaces are a structural signature of glycan-binding proteins, and can be used to computationally predict novel glycan-binding proteins from 3 D structure.

Bioinformatics and molecular modeling in glycobiology

The field of glycobiology is concerned with the study of the structure, properties, and biological functions of the family of biomolecules called carbohydrates. Bioinformatics for glycobiology is a particularly challenging field, because carbohydrates exhibit a high structural diversity and their chains are often branched. Significant improvements in experimental analytical methods over recent years have led to a tremendous increase in the amount of carbohydrate structure data generated. Consequently, the availability of databases and tools to store, retrieve and analyze these data in an efficient way is of fundamental importance to progress in glycobiology. In this review, the various graphical representations and sequence formats of carbohydrates are introduced, and an overview of newly developed databases, the latest developments in sequence alignment and data mining, and tools to support experimental glycan analysis are presented. Finally, the field of structural glycoinformatics and molecular modeling of carbohydrates, glycoproteins, and protein-carbohydrate interaction are reviewed.

Tryptophan as a molecular shovel in the glycosyl transfer activity of Trypanosoma cruzi trans-sialidase

Molecular dynamics investigations into active site plasticity of Trypanosoma cruzi trans-sialidase, a protein implicated in Chagas disease, suggest that movement of the Trp(312) loop plays an important role in the enzyme's sialic acid transfer mechanism. The observed Trp(312) flexibility equates to a molecular shovel action, which leads to the expulsion of the donor aglycone leaving group from the catalytic site. These computational simulations provide detailed structural insights into sialyl transfer by the trans-sialidase and may aid the design of inhibitors effective against this neglected tropical disease.

Enthalpic nature of the CH/pi interaction involved in the recognition of carbohydrates by aromatic compounds, confirmed by a novel interplay of NMR, calorimetry, and theoretical calculations

Specific interactions between molecules, including those produced by a given solute, and the surrounding solvent are essential to drive molecular recognition processes. A simple molecule such as benzene is capable of recognizing and differentiating among very similar entities, such as methyl 2,3,4,6-tetra-O-methyl-alpha-D-galactopyranoside (alpha-Me(5)Gal), methyl 2,3,4,6-tetra-O-methyl-beta-D-galactopyranoside (beta-Me(5)Gal), 1,2,3,4,6-penta-O-acetyl-beta-D-galactopyranose (beta-Ac(5)Gal), and methyl 2,3,4,6-tetra-O-methyl-alpha-D-mannopyranoside (alpha-Me(5)Man). In order to determine if these complexes are formed, the interaction energy between benzene and the different carbohydrates was determined, using Calvet microcalorimetry, as the enthalpy of solvation. These enthalpy values were -89.0 +/- 2.0, -88.7 +/- 5.5, -132.5 +/- 6.2, and -78.8 +/- 3.9 kJ mol(-1) for the four complexes, respectively. Characterization of the different complexes was completed by establishing the molecular region where the interaction takes place using NMR. It was determined that beta-Me(5)Gal is stabilized by the CH/pi interaction produced by the nonpolar region of the carbohydrate on the alpha face. In contrast, alpha-Me(5)Man is not specifically solvated by benzene and does not present any stacking interaction. Although alpha-Me(5)Gal has a geometry similar to that of its epimer, the obtained NMR data seem to indicate that the axial methoxy group at the anomeric position increases the distance of the benzene molecules from the pyranose ring. Substitution of the methoxy groups by acetate moieties, as in beta-Ac(5)Gal, precludes the approach of benzene to produce the CH/pi interaction. In fact, the elevated stabilization energy of beta-Ac(5)Gal is probably due to the interaction between benzene and the methyl groups of the acetyls. Therefore, methoxy and acetyl substituents have different effects on the protons of the pyranose ring.

Progress in biomimetic carbohydrate recognition

The importance of carbohydrate recognition in biology, and the unusual challenges involved, have lead to great interest in mimicking saccharide-binding proteins such as lectins. In this review, we discuss the design of artificial carbohydrate receptors, focusing on those which work under natural (i.e. aqueous) conditions. The problem is intrinsically difficult because of the similarity between substrate (carbohydrate) and solvent (water) and, accordingly, progress has been slow. However, recent developments suggest that solutions can be found. In particular, the "temple" family of carbohydrate receptors show good affinities and excellent selectivities for certain all-equatorial substrates. One example is selective for O-linked beta-N-acetylglucosamine (GlcNAc, as in the O-GlcNAc protein modification), while another is specific for beta-cellobiosyl and closely related disaccharides. Both show roughly millimolar affinities, matching the strength of some lectin-carbohydrate interactions.

Carbohydrate-pi interactions: what are they worth?

Protein-carbohydrate interactions play an important role in many biologically important processes. The recognition is mediated by a number of noncovalent interactions, including an interaction between the alpha-face of the carbohydrate and the aromatic side chain of the protein. To elucidate this interaction, it has been studied in the context of a beta-hairpin in aqueous solution, in which the interaction can be investigated in the absence of other cooperative noncovalent interactions. In this beta-hairpin system, both the aromatic side chain and the carbohydrate were varied in an effort to gain greater insight into the driving force and magnitude of the carbohydrate-pi interaction. The magnitude of the interaction was found to vary from -0.5 to -0.8 kcal/mol, depending on the nature of the aromatic ring and the carbohydrate. Replacement of the aromatic ring with an aliphatic group resulted in a decrease in interaction energy to -0.1 kcal/mol, providing evidence for the contribution of CH-pi interactions to the driving force. These findings demonstrate the significance of carbohydrate-pi interactions within biological systems and also their utility as a molecular recognition element in designed systems.