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

NMR analysis of the binding mode of two fungal endo-β-1,4-mannanases from GH5 and GH26 families

A different mode of action among two endo-β-1,4 mannanases from Podospora anserina has been revealed through an accurate NMR analysis.

Carbohydrate-protein interactions: molecular modeling insights

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

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

BACKGROUND

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

SCOPE OF REVIEW

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

MAJOR CONCLUSIONS

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

GENERAL SIGNIFICANCE

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

Reduced hematopoietic toxicity, enhanced cellular uptake and altered pharmacokinetics of azidothymidine loaded galactosylated liposomes

In order to target liposomes to the lectin receptors present on macrophages, galactosylated liposomes were prepared and characterized in vitro. O-palmitoylgalactose (OPG) for liposomal coating was synthesized by esterification of galactose with palmitoyl chloride. The galactose binding Ricinus communis lectin was employed as a model system for the determination of in vitro ligand binding capacity. Cellular drug uptake studies were performed using alveolar macrophages. Hematological changes, bone marrow toxicity, plasma and tissue distribution study of free, uncoated plain liposomal and galactosylated liposomal encapsulated azidothymidine (AZT) were determined following a bolus intravenous injection in Sprague-Dawley rats. Lectin (R. communis) carbohydrate interaction has been utilized for the effective delivery of AZT entrapped in galactosylated vesicles. Aggregation of galactosylated liposomes increased as lectin concentration was increased from 5 to 30 microg/ml. Cellular uptake of galactosylated liposomal formulation was maximum. No hematological toxicity was observed even after 10 days in case of galactosylated vesicle entrapped AZT. This formulation maintained a significant level of AZT in tissues rich in galactose specific receptors and had a prolonged residence in the body resulting in enhanced half-life of AZT. Conclusively, galactosylated liposomes are the potential candidate for targeted drug delivery and are anticipated to be promising in the treatment of AIDS6.

Identification of common structural features of binding sites in galactose-specific proteins

Galactose-binding proteins characterize an important subgroup of sugar-binding proteins that are involved in a variety of biological processes. Structural studies have shown that the Gal-specific proteins encompass a diverse range of primary and tertiary structures. The binding sites for galactose also seem to vary in different protein-galactose complexes. No common binding site features that are shared by the Gal-specific proteins to achieve ligand specificity are so far known. With the assumption that common recognition principles will exist for common substrate recognition, the present study was undertaken to identify and characterize any unique galactose-binding site signature by analyzing the three-dimensional (3D) structures of 18 protein-galactose complexes. These proteins belong to 7 nonhomologous families; thus, there is no sequence or structural similarity across the families. Within each family, the binding site residues and their relative distances were well conserved, but there were no similarities across families. A novel, yet simple, approach was adopted to characterize the binding site residues by representing their relative spatial dispositions in polar coordinates. A combination of the deduced geometrical features with the structural characteristics, such as solvent accessibility and secondary structure type, furnished a potential galactose-binding site signature. The signature was evaluated by incorporation into the program COTRAN to search for potential galactose-binding sites in proteins that share the same fold as the known galactose-binding proteins. COTRAN is able to detect galactose-binding sites with a very high specificity and sensitivity. The deduced galactose-binding site signature is strongly validated and can be used to search for galactose-binding sites in proteins. PROSITE-type signature sequences have also been inferred for galectin and C-type animal lectin-like fold families of Gal-binding proteins.

The conformation of the T-antigen disaccharide bound to Maclura pomifera agglutinin in aqueous solution

The complex of Maclura pomifera agglutinin with the T-antigen disaccharide (beta-d-Gal-(1-->3)-alpha-d-GalNAc-(1-->O)-Me) was investigated by NMR spectroscopy in aqueous solution. Intramolecular transferred nuclear Overhauser enhancement (NOE) effects between the monosaccharide moieties were used to derive the ligand conformation in the lectin-bound state. Ligand protons in contact with the protein were identified by saturation transfer difference experiments and intermolecular transferred NOE effects. It is demonstrated that structural differences exist for the ligand-lectin complex in aqueous solution as compared with the previously published crystal structure (Lee, X., Thompson, A., Zhiming, Z., Ton-that, H., Biesterfeldt, J., Ogata, C., Xu, L., Johnston, R. A. Z. , and Young, N. M. (1998) J. Biol. Chem. 273, 6312-6318). In order to accommodate the O-methyl group of the disaccharide, the amino acid side chain of Tyr-122 has to rotate from its position in the crystal. The NMR data are in accord with two conformational families at the beta-(1-->3)glycosidic linkage in the solution complex with interglycosidic angles phi/psi = 45/-65 degrees and -65/-18 degrees. These differ from the bound conformation of the ligand in the crystal (phi/psi = 39/-8 degrees ) and are not highly populated by the ligand in the free state. The reason for the structural differences at the beta-(1-->3)glycosidic linkage are hydrogen bonds that stabilize the relative orientation of the monosaccharide units in the crystal. Our results demonstrate that the crystallization of a protein-carbohydrate complex can interfere with the delicate process of carbohydrate recognition in solution.

Complexes of glucoamylase with maltoside heteroanalogues: bound ligand conformations by use of transferred NOE experiments and molecular modeling

Transferred nuclear Overhauser effect (trNOE) experiments have been performed to investigate the conformations of the competitive inhibitors, methyl 5'-thio-4-N-alpha-maltoside 3a and methyl 5'-thio-4-S-alpha-maltoside 4 when bound to the catalytic subunit of the enzyme glucoamylase. These NMR data suggest that, although each of the free ligands populates two conformational families, both heteroanalogues are bound by the enzyme in conformations in the area of the global energy minimum. These conformations have been used as initial points for docking into the active site of the enzyme taken from a X-ray crystal structure of the related glucoamylase-D-gluco-dihydroacarbose 2 complex. Minimization of the resulting complexes has yielded structures for the bound complexes. Corroboration of the structures is provided by fast T(1)(rho)-relaxation effects for certain ligand protons as a result of close contacts with protons in the enzyme active site. The results auger well for the combined use of transferred NOE spectroscopy and molecular modeling based on X-ray crystal structures of complexes of suitable congeners for the rapid analysis of ligand-receptor interactions.

Structure and conformation of complex carbohydrates of glycoproteins, glycolipids, and bacterial polysaccharides

For nuclear magnetic resonance determinations of the conformation of oligosaccharides in solution, simple molecular mechanics calculations and nuclear Overhauser enhancement measurements are adequate for small oligosaccharides that adopt single, relatively rigid conformations. Polysaccharides and larger or more flexible oligosaccharides generally require additional types of data, such as scalar and dipolar coupling constants, which are most conveniently measured in 13C-enriched samples. Nuclear magnetic resonance relaxation data provide information on the dynamics of oligosaccharides, which involves several different types of internal motion. Oligosaccharides complexed with lectins and antibodies have been successfully studied both by X-ray crystallography and by nuclear magnetic resonance spectroscopy. The complexes have been shown to be stabilized by a combination of polar hydrogen bonding interactions and van der Waals attractions. Although theoretical calculations of the conformation and stability of free oligosaccharides and of complexes with proteins can be carried out by molecular mechanics methods, the role of solvent water for these highly polar molecules continues to present computational problems.

Structure of heparin-derived tetrasaccharide complexed to the plasma protein antithrombin derived from NOEs, J-couplings and chemical shifts

A complex of the synthetic tetrasaccharide AGA*IM [GlcN, 6-SO3-alpha(1-4)-GlcA-beta(1-4)-GlcN,3, 6-SO3-alpha(1-4)-IdoA-alphaOMe] and the plasma protein antithrombin has been studied by NMR spectroscopy. 1H and 13C chemical shifts, three-bond proton-proton (3JH-H) and one-bond proton-carbon coupling constants (1JC-H) as well as transferred NOEs and rotating frame Overhauser effects (ROEs) were monitored as a function of the protein : ligand molar ratio and temperature. Considerable changes were observed at both 20 : 1 and 10 : 1 ratios (AGA*IM : antithrombin) in 1H as well as 13C chemical shifts. The largest changes in 1H chemical shifts, and the linewidths, were found for proton resonances (A1, A2, A6, A6', A1*, A2*, A3*, A4*) in GlcN, 6-SO3 and GlcN,3,6-SO3 units, indicating that both glucosamine residues are strongly involved in the binding process. The changes in the linewidths in the IdoA residue were considerably smaller than those in other residues, suggesting that the IdoA unit experienced different internal dynamics during the binding process. This observation was supported by measurements of 3JH-H and 1JC-H. The magnitude of the three-bond proton-proton couplings (3JH1-H2 = 2.51 Hz and 3JH4-H5 = 2.23 Hz) indicate that in the free state an equilibrium exists between 1C4 and 2S0 conformers in the ratio of approximately 75 : 25. The chair form appears the more favourable in the presence of antithrombin, as inferred from the magnitude of the coupling constants. In addition, two-dimensional NOESY and ROESY experiments in the free ligand, as well as transferred NOESY and ROESY spectra of the complex, were measured and interpreted using full relaxation and conformational exchange matrix analysis. The theoretical NOEs were computed using the geometry of the tetrasaccharide found in a Monte Carlo conformational search, and the three-dimensional structures of AGA*IM in both free and bound forms were derived. All monitored NMR variables, 1H and 13C chemical shifts, 1JC-H couplings and transferred NOEs, indicated that the changes in conformation at the glycosidic linkage GlcN, 6-SO3-alpha(1-4)-GlcA were induced by the presence of antithrombin and suggested that the receptor selected a conformer different from that in the free state. Such changes are compatible with the two-step model [Desai, U.R., Petitou, M., Bjork, I. & Olson, S. (1998) J. Biol. Chem. 273, 7478-7487] for the interaction of heparin-derived oligosaccharides with antithrombin, but with a minor extension: in the first step a low-affinity recognition complex between ligand and receptor is formed, accompanied by a conformational change in the tetrasaccharide, possibly creating a complementary three-dimensional structure to fit the protein-binding site. During the second step, as observed in a structurally similar pentasaccharide [Skinner, R., Abrahams, J.-P., Whisstock, J.C., Lesk, A.M., Carrell, R.W. & Wardell, M.R. (1997) J. Mol. Biol. 266, 601-609; Jin, L., Abrahams, J.-P., Skinner, R., Petitou, M., Pike, R. N. & Carrell, R.W. (1997) Proc. Natl Acad. Sci. USA 94, 14683-14688], conformational changes in the binding site of the protein result in a latent conformation.