Practical tools for molecular modeling of complex carbohydrates and their interactions with proteins

Computer modeling has become a valuable component of studies of carbohydrate three-dimensional structures and their relationship to function and properties. In this paper we examine the methods required for conformational modeling of carbohydrates, and we present a series of tools that have been developed to this end. These tools can be integrated into three-dimensional real-time molecular modeling software. A data base of pre-optimized carbohydrate fragments has been established to be used further in the construction of much more complex molecules. In addition we describe some possible uses of a data base of three dimensional structures of the disaccharide fragments present in the glycan moiety ofN-glycoprotein. A molecular mechanical force field appropriate for the conformational analysis of oligosaccharides has been derived by the addition of new parameters to the Tripos force field and is compatible with protein simulations. The new parametrization has been assessed in three stages of increasing complexity: computations of potential energy surfaces, conformational refinement of relevant oligosaccharides, modeling at the atomic level of a protein/carbohydrate complex.     

The molecular mechanism of sulfated carbohydrate recognition by the cysteine-rich domain of mannose receptor

The mannose receptor (MR) binds foreign and host ligands through interactions with their carbohydrates. Two portions of MR have distinct carbohydrate recognition properties. One is conferred by the amino-terminal cysteine-rich domain (Cys-MR), which plays a critical role in binding sulfated glycoproteins including pituitary hormones. The other is achieved by tandemly arranged C-type lectin domains that facilitate carbohydrate-dependent uptake of infectious microorganisms. This dual carbohydrate binding specificity enables MR to bind ligands by interacting with both sulfated and non-sulfated polysaccharide chains. We previously determined crystal structures of Cys-MR complexed with 4-SO(4)-N-acetylglucosamine and with an unidentified ligand resembling Hepes (N-[2-hydroxyethyl]piperazine-N'-[2-ethanesulfonic acid]). In continued efforts to elucidate the mechanism of sulfated carbohydrate recognition by Cys-MR, we characterized the binding affinities between Cys-MR and potential carbohydrate ligands using a fluorescence-based assay. We find that Cys-MR binds sulfated carbohydrates with relatively high affinities (K(D)=0.1 mM to 1.0 mM) compared to the affinities of other lectins. Cys-MR also binds Hepes with a K(D) value of 3.9 mM, consistent with the suggestion that the ligand in the original Cys-MR crystal structure is Hepes. We also determined crystal structures of Cys-MR complexed with 3-SO(4)-Lewis(x), 3-SO(4)-Lewis(a), and 6-SO(4)-N-acetylglucosamine at 1.9 A, 2.2 A, and 2.5 A resolution, respectively, and the 2.0 A structure of Cys-MR that had been treated to remove Hepes. The conformation of the Cys-MR binding site is virtually identical in all Cys-MR crystal structures, suggesting that Cys-MR does not undergo conformational changes upon ligand binding. The structures are used to rationalize the binding affinities derived from the biochemical studies and to elucidate the molecular mechanism of sulfated carbohydrate recognition by Cys-MR.     

DNA polymorphism: a comparison of force fields for nucleic acids

The improvements of the force fields and the more accurate treatment of long-range interactions are providing more reliable molecular dynamics simulations of nucleic acids. The abilities of certain nucleic acid force fields to represent the structural and conformational properties of nucleic acids in solution are compared. The force fields are AMBER 4.1, BMS, CHARMM22, and CHARMM27; the comparison of the latter two is the primary focus of this paper. The performance of each force field is evaluated first on its ability to reproduce the B-DNA decamer d(CGATTAATCG)(2) in solution with simulations in which the long-range electrostatics were treated by the particle mesh Ewald method; the crystal structure determined by Quintana et al. (1992) is used as the starting point for all simulations. A detailed analysis of the structural and solvation properties shows how well the different force fields can reproduce sequence-specific features. The results are compared with data from experimental and previous theoretical studies.     

Structure,dynamics, and interactions of jacalin. Insights from molecular dynamics simulations examined in conjunction with results of X‐ray studies

Molecular dynamics simulations have been carried out on all the jacalin–carbohydrate complexes of known structure, models of unliganded molecules derived from the complexes and also models of relevant complexes where X‐ray structures are not available. Results of the simulations and the available crystal structures involving jacalin permit delineation of the relatively rigid and flexible regions of the molecule and the dynamical variability of the hydrogen bonds involved in stabilizing the structure. Local flexibility appears to be related to solvent accessibility. Hydrogen bonds involving side chains and water bridges involving buried water molecules appear to be important in the stabilization of loop structures. The lectin–carbohydrate interactions observed in crystal structures, the average parameters pertaining to them derived from simulations, energetic contribution of the stacking residue estimated from quantum mechanical calculations, and the scatter of the locations of carbohydrate and carbohydrate‐binding residues are consistent with the known thermodynamic parameters of jacalin–carbohydrate interactions. The simulations, along with X‐ray results, provide a fuller picture of carbohydrate binding by jacalin than provided by crystallographic analysis alone. The simulations confirm that in the unliganded structures water molecules tend to occupy the positions occupied by carbohydrate oxygens in the lectin–carbohydrate complexes. Population distributions in simulations of the free lectin, the ligands, and the complexes indicate a combination of conformational selection and induced fit. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

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.     

Structural monitoring of oligosaccharides through 13C enrichment and NMR observation of acetyl groups

Structural characterization of biomolecules by NMR methods frequently requires the enrichment of magnetically active isotopes at particular molecular sites. Introduction is usually achieved biosynthetically through the use of bacterial cultures grown on isotopically enriched media, but for certain types of molecules-cell-surface carbohydrates of mammalian origin, for example-this is not practical. Here we explore a means of introducing 13C-enriched sites, postisolation in natural carbohydrate products, and illustrate an ability to acquire sufficient information to select appropriate conformational models from among energetically allowed sets. The application presented involves replacement of native N-acetyl groups with 13C-labeled acetyl groups in a simple disaccharide derivative, (GlcNAc)2-OBu, or O-butyl-chitobiose. The assignment of the two acetyl groups introduced is based on a novel combination of NMR and mass spectrometry data. Structural information is obtained from chemical shift anisotropy offsets of 13C carbonyl resonances and 13C-13C dipolar couplings between the labeled methyl and carbonyl carbons of the acetyl groups. Although the application is to a relatively simple system, it lays the groundwork for application to biologically important complex carbohydrate systems.     

Role of glycosylation on the ensemble of conformations in the MUC1 immunodominant epitope

MUC1 is a membrane glycoprotein, which in adenocarninomas is overexpressed and exhibits truncated O‐glycosylation. Overexpression and altered glycosylation make MUC1 into a candidate for immunotherapy. Monoclonal antibodies directed against MUC1 frequently bind an immunodominant epitope that contains a single site for O‐glycosylation. Glycosylation with tumor carbohydrate antigens such as the Tn‐antigen (GalNAc‐O‐Ser/Thr) results in antibodies binding with higher affinity. One proposed model to explain the enhanced affinity of antibodies for the glycosylated antigen is that the addition of a carbohydrate alters the conformational properties, favoring a binding‐competent state. The conformational effects associated with Tn glycosylation of the MUC1 antigen was investigated using solution‐state NMR and molecular dynamics. NMR experiments revealed distinct substructures of the glycosylated MUC1 peptides compared with the unglycosylated peptide. Molecular dynamics simulations of the MUC1 glycopeptide and peptide revealed distinguishing differences in their conformational preferences. Furthermore, the glycopeptide displayed a smaller conformational sampling compared with the peptide, suggesting that the glycopeptide sampled a narrower conformational space and is less dynamic. A comparison of the computed ensemble of conformations assuming random distribution, NMR models, and molecular dynamics simulations indicated that the MUC1 glycopeptide and aglycosylated peptide sampled structurally distinctly ensembles and that these ensembles were different from that of the random coil. Together, these data support the hypothesis that that conformational pre‐selection could be an essential feature of these peptides that dictates the binding affinities to MUC1 specific antibodies.     

Conformational study of cyclolinopeptide A. A distance geometry and molecular dynamics approach

The conformation of cyclolinopeptide A, c(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val), a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of distance geometry calculations and molecular dynamics simulations. The starting points for all the calculations were an X-ray structure and other structures obtained from distance geometry calculations based on NMR data. Restrained and unrestrained molecular dynamics simulations are reported in vacuo and in CCl4. Structural and dynamic properties are investigated and compared with those experimentally determined. The conformation obtained from the MD simulations which best reproduces the NMR parameters is at the same time one of the most stable ones and is also fairly similar to the crystal structure. An explanation for the occurrence of multiple conformations in solution at room temperature is given.     

Analysis of protein-carbohydrate interaction at the lower size limit of the protein part (15-mer peptide) by NMR spectroscopy,electrospray ionization mass spectrometry,and molecular modeling

Structural analysis of minimally sized lectins will offer insights into fundamentals of intermolecular recognition and potential for biomedical applications. We thus moved significantly beyond the natural limit of lectin size to determine the structure of synthetic mini-lectins in solution, their carbohydrate selectivity and the impact of ligand binding on their conformational behavior. Using three disaccharide (Thomsen-Friedenreich antigen; Gal beta 1,3GalNAc alpha 1,R)-binding pentadecapeptides without internal disulfide bridges as role models, we successfully tested a combined strategy with different techniques of NMR spectroscopy, electrospray ionization mass spectrometry, and molecular modeling. In solution, the peptides invariably displayed flexibility with rather limited restrictions, shown by NMR experiments including nearly complete resonance assignments and molecular dynamics simulations. The occurrence of aromatic/nonpolar amino acids in the sequence did not lead to formation of a hydrophobic core known from microbial chitinase modules. Selectivity of disaccharide binding was independently observed by mass spectrometry and NMR analysis. Specific ligand interaction yielded characteristic NMR signal alterations but failed to reduce conformational flexibility significantly. We have thereby proven effectiveness of our approach to analyze even low-affinity interactions (not restricted to carbohydrates as ligands). It will be useful to evaluate the impact of rational manipulation of lead peptide sequences.     

糖类结构的光谱分析的特点

光谱分析方法是糖生物学化学方面研究及开发的关键技术手段。由于糖类物质种类繁多、结构多样、构效关系极为复杂,如何合理高效使用波谱分析这一便捷方法,对于糖类物质的准确分析尤为重要,文章重点分析了现代波谱技术,如红外光谱、核磁共振、质谱等在糖生物学研究和化学化工领域的应用特点,特别指出糖分析中需要注意的一些特殊性,以促进谱学分析技术在糖生物医学、化学等领域的有效应用。     

The structure of DC-SIGNR with a portion of its repeat domain lends insights to modeling of the receptor tetramer

The dendritic cell-specific ICAM-3 non-integrin (DC-SIGN) and its close relative DC-SIGNR recognize various glycoproteins, both pathogenic and cellular, through the receptor lectin domain-mediated carbohydrate recognition. While the carbohydrate-recognition domains (CRD) exist as monomers and bind individual carbohydrates with low affinity and are permissive in nature, the full-length receptors form tetramers through their repeat domain and recognize specific ligands with high affinity. To understand the tetramer-based ligand binding avidity, we determined the crystal structure of DC-SIGNR with its last repeat region. Compared to the carbohydrate-bound CRD structure, the structure revealed conformational changes in the calcium and carbohydrate coordination loops of CRD, an additional disulfide bond between the N and the C termini of the CRD, and a helical conformation for the last repeat. On the basis of the current crystal structure and other published structures with sequence homology to the repeat domain, we generated a tetramer model for DC-SIGN/R using homology modeling and propose a ligand-recognition index to identify potential receptor ligands.     

Structure and behaviour of biomolecules from Raman optical activity

Raman optical activity, which can be measured as a small circularly polarized component in Raman-scattered light from chiral molecules, holds much promise for studying a large range of biomolecules in aqueous solution. Among other things, it provides information about motif and fold, as well as secondary structure, of proteins; the solution structure of carbohydrates; and the structure of the polypeptide and carbohydrate components of intact glycoproteins. In addition, new insights into the structural elements present in unfolded protein sequences, and the structure of the protein and nucleic acid components of intact viruses can be obtained. Ab initio quantum-chemical simulations of observed Raman optical activity spectra provide the complete three-dimensional structure of small biomolecules. Raman optical activity measurements are now routine thanks to the availability of a commercial instrument based on a novel design.     

Effect of molecular structure on thermodynamic properties of carbohydrates. A calorimetric study of aqueous di- and oligosaccharides at subzero temperatures

For aqueous solutions of di- and oligosaccharides thermodynamic properties have been investigated at subzero temperatures using differential scanning calorimetry. The amount of unfrozen water observed is found to increase linearly with the glass transition temperatures of anhydrous carbohydrates. Furthermore, the amount of unfrozen water shows a linear relationship with known solution properties of aqueous carbohydrates, such as partial molar compressibility and heat of solution. The different effectiveness among various di- and oligosaccharides to avoid ice formation is associated with the combination of constitutive monosaccharides and attendant molecular structure features including the position and type of the glycosidic linkage between the constituent units. More unfrozen water is induced in the presence of a carbohydrate having a poorer compatibility with the three-dimensional hydrogen-bond network of water. A series of these results obtained imply that there is a common key of carbohydrate stereochemistry governing several different thermodynamic amounts of a given system involving carbohydrates. In this context, a modified stereospecific-hydration model can be used to interpret the present results in terms of stereochemical effects of carbohydrates.     

Structural basis for ligand recognition in a mushroom lectin: solvent structure as specificity predictor

Lectins are able to recognize specific carbohydrate structures through their carbohydrate recognition domain (CRD). The lectin from the mushroom Agaricus bisporus (ABL) has the remarkable ability of selectively recognizing the TF-antigen, composed of Galβ1-3GalNAc, Ser/Thr linked to proteins, specifically exposed in neoplastic tissues. Strikingly, the recently solved crystal structure of tetrameric ABL in the presence of TF-antigen and other carbohydrates showed that each monomer has two CRDs, each being able to bind specifically to different monosaccharides that differ only in the configuration of a single hydroxyl, like N-acetyl-d-galactosamine (GalNAc) and N-acetyl-d-glucosamine (GlcNAc). Understanding how lectin CRDs bind and discriminate mono and/or (poly)-saccharides is an important issue in glycobiology, with potential impact in the design of better and selective lectin inhibitors with potential therapeutic properties. In this work, and based on the unusual monosaccharide epimeric specificity of the ABL CRDs, we have performed molecular dynamics simulations of the natural (crystallographic) and inverted (changing GalNAc for GlcNAc and vice-versa) ABL–monosaccharide complexes in order to understand the selective ligand recognition properties of each CRD. We also performed a detailed analysis of the CRD local solvent structure, using previously developed methodology, and related it with the recognition mechanism. Our results provide a detailed picture of each ABL CRD specificity, allowing a better understanding of the carbohydrate selective recognition process in this particular lectin.     

Effect of glycosylation on structure and dynamics of MHC class I glycoprotein: a molecular dynamics study

Complex carbohydrates linked to glycoproteins are recently being implicated to play a variety of biological roles. The lack of well-resolved crystallographic coordinates of the carbohydrates makes it difficult to assess the contributions of the glycan chain on protein structure and dynamics. We have modeled two different oligosaccharides NeuNAc2Gal3Man3GlcNAc5Fuc and Man3GlcNAc4 to generate two glycosylation variants of major histocompatibility complex (MHC) class I glycoprotein. Molecular dynamics simulations of the isolated fourteen- and seven-residue oligosaccharides have been done in vacuo and in solution. The dynamics of the two glycoforms of MHC class I protein have been simulated in solution in the free as well as in the peptide-bound form. Good agreement between the calculated solution conformations of the oligosaccharides in isolated and conjugated forms and the average conformations obtained from x-ray or NMR data was observed for most of the glycosidic linkages. These molecular dynamics simulations of the isolated glycan chains and the glycoconjugates reveal the details of the conformational flexibility of the glycan chains; they also provide atomic level details of protein-carbohydrate interactions and the effect of the ligand binding on the carbohydrate structure and dynamics. It was found that though there is some flexibility in some of the glycosidic linkages in the isolated oligosaccharides, in the protein-conjugated form the linkages adopt more restricted conformations. The glycan chains protrude out into the solvent and might hinder the lateral association of the proteins. The presence of the bulky glycan chains does not affect the average backbone fold of the protein but induces local changes in protein structure and dynamics. It has been noted that the extent of the changes depends upon the nature of the attached glycan chain. The glycan chains do not appear to influence the peptide binding property of the protein directly, but may stabilize the protein residues that are involved in ligand binding.     

Structure and internal dynamics of the bovine pancreatic trypsin inhibitor in aqueous solution from long-time molecular dynamics simulations

Structural and dynamic properties of bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution are investigated using two molecular dynamics (MD) simulations: one of 1.4 ns length and one of 0.8 ns length in which atom-atom distance bounds derived from NMR spectroscopy are included in the potential energy function to make the trajectory satisfy these experimental data more closely. The simulated properties of BPTI are compared with crystal and solution structures of BPTI, and found to be in agreement with the available experimental data. The best agreement with experiment was obtained when atom-atom distance restraints were applied in a time-averaged manner in the simulation. The polypeptide segments found to be most flexible in the MD simulations coincide closely with those showing differences between the crystal and solution structures of BPTI. © 1995 Wiley-Liss, Inc.     

Solution and Crystal Molecular Dynamics Simulation Study of m4-Cyanovirin-N Mutants Complexed with Di-Mannose

Cyanovirin-N (CVN) is a highly potent anti-HIV carbohydrate-binding agent that establishes its microbicide activity through interaction with mannose-rich glycoprotein gp120 on the virion surface. The m4-CVN and P51G-m4-CVN mutants represent simple models for studying the high-affinity binding site, B^M. A recently determined 1.35 Å high-resolution structure of P51G-m4-CVN provided details on the di-mannose binding mechanism, and suggested that the Arg-76 and Glu-41 residues are critical components of high mannose specificity and affinity. We performed molecular-dynamics simulations in solution and a crystal environment to study the role of Arg-76. Network analysis and clustering were used to characterize the dynamics of Arg-76. The results of our explicit solvent solution and crystal simulations showed a significant correlation with conformations of Arg-76 proposed from x-ray crystallographic studies. However, the crystal simulation showed that the crystal environment strongly biases conformational sampling of the Arg-76 residue. The solution simulations demonstrated no conformational preferences for Arg-76, which would support its critical role as the residue that locks the ligand in the bound state. Instead, a comparative analysis of trajectories from >50 ns of simulation for two mutants revealed the existence of a very stable eight-hydrogen-bond network between the di-mannose ligand and predominantly main-chain atoms. This network may play a key role in the specific recognition and strong binding of mannose oligomers in CVN and its homologs.     

A CHARMm Based Force Field for Carbohydrates Using the CHEAT Approach: Carbohydrate Hydroxyl Groups Represented by Extended Atoms

Abstract

Computational studies of carbohydrates that do not consider explicit solvent molecules suffer from the strong tendency of the carbohydrate pendant hydroxyl groups to form intramolecular hydrogen bonds that are unlikely to be present in protic media. In this paper a novel approach towards molecular modelling of carbohydrates is described. The average effect of intra- and intermolecular hydrogen bonding is introduced into the potential energy function by adding a new (extended) atom type representing a carbohydrate hydroxyl group to the CHARMm force field; we coin the name CHEAT (Carbohydrate Hydroxyls represented by Extended AToms) for the resulting force field. As a training set for the parametrisation of CHEAT we used ethylene glycol, 10 cyclohexanols, 5 inositols, and 12 glycopyranoses for which in total 64 conformational energy differences were estimated using a set of steric interaction energies between hydroxyl and/or methyl groups on six-membered ring compounds as derived by Angyal (Angew. Chem., 8, 172-182, (1969)). The root-mean-square deviation between the estimated energy differences and the corresponding values obtained by our CHEAT approach amounts to 0.37 kcal/mol (n = 64). The CHEAT approach, which is claimed to calculate aqueous state compatible energetical and conformational properties of carbohydrates, is computationally very efficient and facilitates the calculation of nanosecond range MD trajectories as well as systematic conformational searches of oligosaccharides.