GLYCAM06: a generalizable biomolecular force field. Carbohydrates

Computational glycoscience: characterizing the spatial and temporal properties of glycans and glycan-protein complexes

Modern computational methods offer the tools to provide insight into the structural and dynamic properties of carbohydrate-protein complexes, beyond that provided by experimental structural biology. Dynamic properties such as the fluctuation of inter-molecular hydrogen bonds, the residency times of bound water molecules, side chain motions and ligand flexibility may be readily determined computationally. When taken with respect to the unliganded states, these calculations can also provide insight into the entropic and enthalpic changes in free energy associated with glycan binding. In addition, virtual ligand screening may be employed to predict the three dimensional (3D) structures of carbohydrate-protein complexes, given 3D structures for the components. In principle, the 3D structure of the protein may itself be derived by modeling, leading to the exciting--albeit high risk--realm of virtual structure prediction. This latter approach is appealing, given the difficulties associated with generating experimental 3D structures for some classes of glycan binding proteins; however, it is also the least robust. An unexpected outcome of the development of algorithms for modeling carbohydrate-protein interactions has been the discovery of errors in reported experimental 3D structures and a heightened awareness of the need for carbohydrate-specific computational tools for assisting in the refinement and curation of carbohydrate-containing crystal structures. Here we present a summary of the basic strategies associated with employing classical force field based modeling approaches to problems in glycoscience, with a focus on identifying typical pitfalls and limitations. This is not an exhaustive review of the current literature, but hopefully will provide a guide for the glycoscientist interested in modeling carbohydrates and carbohydrate-protein complexes, as well as the computational chemist contemplating such tasks.

Plant O-hydroxyproline arabinogalactans are composed of repeating trigalactosyl subunits with short bifurcated side chains

Classical arabinogalactan proteins partially defined by type II O-Hyp-linked arabinogalactans (Hyp-AGs) are structural components of the plant extracellular matrix. Recently we described the structure of a small Hyp-AG putatively based on repetitive trigalactosyl subunits and suggested that AGs are less complex and varied than generally supposed. Here we describe three additional AGs with similar subunits. The Hyp-AGs were isolated from two different arabinogalactan protein fusion glycoproteins expressed in tobacco cells; that is, a 22-residue Hyp-AG and a 20-residue Hyp-AG, both isolated from interferon alpha2b-(Ser-Hyp)(20), and a 14-residue Hyp-AG isolated from (Ala-Hyp)(51)-green fluorescent protein. We used NMR spectroscopy to establish the molecular structure of these Hyp-AGs, which share common features: (i) a galactan main chain composed of two 1-->3 beta-linked trigalactosyl blocks linked by a beta-1-->6 bond; (ii) bifurcated side chains with Ara, Rha, GlcUA, and a Gal 6-linked to Gal-1 and Gal-2 of the main-chain trigalactosyl repeats; (iii) a common side chain structure composed of up to six residues, the largest consisting of an alpha-L-Araf-(1-->5)-alpha-L-Araf-(1-->3)-alpha-L-Araf-(1-->3- unit and an alpha-L-Rhap-(1-->4)-beta-D-GlcUAp-(1-->6)-unit, both linked to Gal. The conformational ensemble obtained by using nuclear Overhauser effect data in structure calculations revealed a galactan main chain with a reverse turn involving the beta-1-->6 link between the trigalactosyl blocks, yielding a moderately compact structure stabilized by H-bonds.

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.

The R.E.D. tools: advances in RESP and ESP charge derivation and force field library building

Deriving atomic charges and building a force field library for a new molecule are key steps when developing a force field required for conducting structural and energy-based analysis using molecular mechanics. Derivation of popular RESP charges for a set of residues is a complex and error prone procedure because it depends on numerous input parameters. To overcome these problems, the R.E.D. Tools (RESP and ESP charge Derive, ) have been developed to perform charge derivation in an automatic and straightforward way. The R.E.D. program handles chemical elements up to bromine in the periodic table. It interfaces different quantum mechanical programs employed for geometry optimization and computing molecular electrostatic potential(s), and performs charge fitting using the RESP program. By defining tight optimization criteria and by controlling the molecular orientation of each optimized geometry, charge values are reproduced at any computer platform with an accuracy of 0.0001 e. The charges can be fitted using multiple conformations, making them suitable for molecular dynamics simulations. R.E.D. allows also for defining charge constraints during multiple molecule charge fitting, which are used to derive charges for molecular fragments. Finally, R.E.D. incorporates charges into a force field library, readily usable in molecular dynamics computer packages. For complex cases, such as a set of homologous molecules belonging to a common family, an entire force field topology database is generated. Currently, the atomic charges and force field libraries have been developed for more than fifty model systems and stored in the RESP ESP charge DDataBase. Selected results related to non-polarizable charge models are presented and discussed.

Mimicking chitin: chemical synthesis, conformational analysis, and molecular recognition of the beta(1-->3) N-acetylchitopentaose analogue

Mimicking Nature by using synthetic molecules that resemble natural products may open avenues to key knowledge that is difficult to access by using substances from natural sources. In this context, a novel N-acetylchitooligosaccharide analogue, beta-1,3-N-acetamido-gluco-pentasaccharide, has been designed and synthesized by using aminoglucose as the starting material. A phthalic group has been employed as the protecting group of the amine moiety, whereas a thioalkyl was used as the leaving group on the reducing end. The conformational properties of this new molecule have been explored and compared to those of the its chito analogue, with the beta-1,3 linkages, by a combined NMR spectroscopic/molecular modeling approach. Furthermore, the study of its molecular recognition properties towards two proteins, a lectin (wheat germ agglutinin) and one enzyme (a chitinase) have also been performed by using NMR spectroscopy and docking protocols. There are subtle differences in the conformational behavior of the mimetic versus the natural chitooligosaccharide, whereas this mimetic is still recognized by these two proteins and can act as a moderate inhibitor of chitin hydrolysis.

Role of water molecules in structure and energetics of Pseudomonas aeruginosa lectin I interacting with disaccharides

Calcium-dependent lectin I from Pseudomonas aeruginosa (PA-IL) binds specifically to oligosaccharides presenting an alpha-galactose residue at their nonreducing end, such as the disaccharides alphaGal1-2betaGalOMe, alphaGal1-3betaGalOMe, and alphaGal1-4betaGalOMe. This provides a unique model for studying the effect of the glycosidic linkage of the ligands on structure and thermodynamics of the complexes by means of experimental and theoretical tools. The structural features of PA-IL in complex with the three disaccharides were established by docking and molecular dynamics simulations and compared with those observed in available crystal structures, including PA-IL.alphaGal1-2betaGalOMe complex, which was solved at 2.4 A resolution and reported herein. The role of a structural bridge water molecule in the binding site of PA-IL was also elucidated through molecular dynamics simulations and free energy calculations. This water molecule establishes three very stable hydrogen bonds with O6 of nonreducing galactose, oxygen from Pro-51 main chain, and nitrogen from Gln-53 main chain of the lectin binding site. Binding free energies for PA-IL in complex with the three disaccharides were investigated, and the results were compared with the experimental data determined by titration microcalorimetry. When the bridge water molecule was included in the free energy calculations, the simulations predicted the correct binding affinity trends with the 1-2-linked disaccharide presenting three times stronger affinity ligand than the other two. These results highlight the role of the water molecule in the binding site of PA-IL and indicate that it should be taken into account when designing glycoderivatives active against P. aeruginosa adhesion.

Population distribution of flexible molecules from maximum entropy analysis using different priors as background information: application to the Φ, Ψ-conformational space of the α-(1-->2)-linked mannose disaccharide present in N- and O-linked glycoproteins

The conformational space available to the flexible molecule α-D-Manp-(1-->2)-α-D-Manp-OMe, a model for the α-(1-->2)-linked mannose disaccharide in N- or O-linked glycoproteins, is determined using experimental data and molecular simulation combined with a maximum entropy approach that leads to a converged population distribution utilizing different input information. A database survey of the Protein Data Bank where structures having the constituent disaccharide were retrieved resulted in an ensemble with >200 structures. Subsequent filtering removed erroneous structures and gave the database (DB) ensemble having three classes of mannose-containing compounds, viz., N- and O-linked structures, and ligands to proteins. A molecular dynamics (MD) simulation of the disaccharide revealed a two-state equilibrium with a major and a minor conformational state, i.e., the MD ensemble. These two different conformation ensembles of the disaccharide were compared to measured experimental spectroscopic data for the molecule in water solution. However, neither of the two populations were compatible with experimental data from optical rotation, NMR (1)H,(1)H cross-relaxation rates as well as homo- and heteronuclear (3)J couplings. The conformational distributions were subsequently used as background information to generate priors that were used in a maximum entropy analysis. The resulting posteriors, i.e., the population distributions after the application of the maximum entropy analysis, still showed notable deviations that were not anticipated based on the prior information. Therefore, reparameterization of homo- and heteronuclear Karplus relationships for the glycosidic torsion angles Φ and Ψ were carried out in which the importance of electronegative substituents on the coupling pathway was deemed essential resulting in four derived equations, two (3)J(COCC) and two (3)J(COCH) being different for the Φ and Ψ torsions, respectively. These Karplus relationships are denoted JCX/SU09. Reapplication of the maximum entropy analysis gave excellent agreement between the MD- and DB-posteriors. The information entropies show that the current reparametrization of the Karplus relationships constitutes a significant improvement. The Φ(H) torsion angle of the disaccharide is governed by the exo-anomeric effect and for the dominating conformation Φ(H) = -40 degrees and Ψ(H) = 33 degrees. The minor conformational state has a negative Ψ(H) torsion angle; the relative populations of the major and the minor states are approximately 3 : 1. It is anticipated that application of the methodology will be useful to flexible molecules ranging from small organic molecules to large biomolecules.

Molecular simulations of carbohydrates and protein-carbohydrate interactions: motivation, issues and prospects

The characterization of the 3D structure of oligosaccharides, their conjugates and analogs is particularly challenging for traditional experimental methods. Molecular simulation methods provide a basis for interpreting sparse experimental data and for independently predicting conformational and dynamic properties of glycans. Here, we summarize and analyze the issues associated with modeling carbohydrates, with a detailed discussion of four of the most recently developed carbohydrate force fields, reviewed in terms of applicability to natural glycans, carbohydrate-protein complexes and the emerging area of glycomimetic drugs. In addition, we discuss prospectives and new applications of carbohydrate modeling in drug discovery.

Common sialylated glycan in Actinobacillus suis

A sialylated oligosaccharide was identified in four representative strains of the Gram-negative swine pathogen, Actinobacillus suis. As characterized, the glycan consists of a free oligosaccharide with a N-acetyl-lactosamine-like backbone decorated with sialic acid, phosphoethanolamine (PEA) and O-acetyl units: 9-O-Ac-Neu5Ac-(2-->6)-beta-d-Galp-(1-->4)-beta-d-6-O-Ac-GlcpNAc-(1-->3)-[PEA-->6]-beta-d-Galp-(1-->3)-beta-d-GlcpNAc-(1-->2)-[9-O-Ac-Neu5Ac-(2-->6)]-beta-d-Galp-(1-->4)-beta-d-6-O-Ac-GlcpNAc-(1-->3)-[PEA-->6]-beta-d-Galp-(1-->3)-d-GlcpNAc. The ubiquitous expression of this sialylated glycan suggests that this carbohydrate may play an important role in the survival of A. suis in the host.

Computational studies on the interaction of ABO-active saccharides with the norovirus VA387 capsid protein can explain experimental binding data

Norovirus strains are known to cause recurring epidemics of winter vomiting disease. The crystal structure of the capsid protein of VA387, a representative of the clinically important GII.4 genocluster, was recently solved in complex with histo-blood group A- and B-trisaccharides. However, the VA387 strain is known to bind also to other natural carbohydrates for which detailed structural information of the complexes is not available. In this study we have computationally explored the fit of the VA387 with a set of naturally occurring carbohydrate ligands containing a terminal alpha1,2-linked fucose. MD simulations both with explicit and implicit solvent models indicate that type 1 and 3 extensions of the ABO-determinant including ALe(b) and BLe(b) pentasaccharides can be well accommodated in the site. Scoring with Glide XP indicates that the downstream extensions of the ABO-determinants give an increase in binding strength, although the alpha1,2-linked fucose is the single strongest interacting residue. An error was discovered in the geometry of the GalNAc-Gal moiety of the published crystal structure of the A-trisaccharide/VA387 complex. The present modeling of the complexes with histo-blood group A-active structures shows some contacts which provide insight into mutational data, explaining the involvement of I389 and Q331. Our results can be applicable in structure-based design of adhesion inhibitors of noroviruses.

Presentation of membrane-anchored glycosphingolipids determined from molecular dynamics simulations and NMR paramagnetic relaxation rate enhancement

Challenges for structural characterization of membrane-bound glycosphingolipids include their high internal dynamic motions and their physical proximity to membrane surfaces. Here we demonstrate that NMR paramagnetic relaxation enhancement can be used, alongside independent molecular dynamics simulations and an outer-sphere relaxation model, to quantitatively characterize the presentation (insertion depth and orientation relative to a membrane surface) of ganglioside GM1 in biologically relevant membrane environments. Longitudinal and transverse paramagnetic relaxation enhancement effects were measured for GM1, anchored to phospholipid bicelles, using both water-soluble and membrane-anchored paramagnetic probes, respectively. A method was developed to rapidly calculate paramagnetic relaxation enhancement effects from thousands of structures taken from a simulation of GM1 in a phospholipid bilayer. The combined computational and experimental approach yielded experimentally verified atomic-resolution 3D models of a highly plastic membrane-bound biomolecule.

Structure and function of glycosylated tandem repeats from Candida albicans Als adhesins

Tandem repeat (TR) regions are common in yeast adhesins, but their structures are unknown, and their activities are poorly understood. TR regions in Candida albicans Als proteins are conserved glycosylated 36-residue sequences with cell-cell aggregation activity (J. M. Rauceo, R. De Armond, H. Otoo, P. C. Kahn, S. A. Klotz, N. K. Gaur, and P. N. Lipke, Eukaryot. Cell 5:1664-1673, 2006). Ab initio modeling with either Rosetta or LINUS generated consistent structures of three-stranded antiparallel beta-sheet domains, whereas randomly shuffled sequences with the same composition generated various structures with consistently higher energies. O- and N-glycosylation patterns showed that each TR domain had exposed hydrophobic surfaces surrounded by glycosylation sites. These structures are consistent with domain dimensions and stability measurements by atomic force microscopy (D. Alsteen, V. Dupres, S. A. Klotz, N. K. Gaur, P. N. Lipke, and Y. F. Dufrene, ACS Nano 3:1677-1682, 2009) and with circular dichroism determination of secondary structure and thermal stability. Functional assays showed that the hydrophobic surfaces of TR domains supported binding to polystyrene surfaces and other TR domains, leading to nonsaturable homophilic binding. The domain structures are like "classic" subunit interaction surfaces and can explain previously observed patterns of promiscuous interactions between TR domains in any Als proteins or between TR domains and surfaces of other proteins. Together, the modeling techniques and the supporting data lead to an approach that relates structure and function in many kinds of repeat domains in fungal adhesins.

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 A 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.

NMR and MD investigations of human galectin-1/oligosaccharide complexes

The specific recognition of carbohydrates by lectins plays a major role in many cellular processes. Galectin-1 belongs to a family of 15 structurally related beta-galactoside binding proteins that are able to control a variety of cellular events, including cell cycle regulation, adhesion, proliferation, and apoptosis. The three-dimensional structure of galectin-1 has been solved by x-ray crystallography in the free form and in complex with various carbohydrate ligands. In this work, we used a combination of two-dimensional NMR titration experiments and molecular-dynamics simulations with explicit solvent to study the mode of interaction between human galectin-1 and five galactose-containing ligands. Isothermal titration calorimetry measurements were performed to determine their affinities for galectin-1. The contribution of the different hexopyranose units in the protein-carbohydrate interaction was given particular consideration. Although the galactose moiety of each oligosaccharide is necessary for binding, it is not sufficient by itself. The nature of both the reducing sugar in the disaccharide and the interglycosidic linkage play essential roles in the binding to human galectin-1.

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.

Barley xyloglucan xyloglucosyl transferases bind xyloglucan-derived oligosaccharides in their acceptor-binding regions in multiple conformational states

Three barley xyloglucan endotransglycosylases (HvXETs), known as xyloglucan xyloglucosyl transferases (EC 2.4.1.207), were subjected to kinetic and computational docking studies. The k(cat) x K(m)(-1) values with the reduced [3H]-labelled XXXG, XXLG/XLXG and XLLG acceptor substrates were 0.02 x 10(-2), 0.1 x 10(-2) and 3.2 x 10(-2) s(-1) microM(-1), while the K(m) constants were 10.6, 8.6 and 5.3 mM, obtained for HvXET3, HvXET4 and HvXET6, respectively. Docking of XLLG in acceptor-binding regions revealed that at least two conformational states were likely to participate in all isoforms. The assessments of kinetic and computational data indicated that the disposition of aromatic residues at the entrance to the active sites and the flexibility of proximal COOH-terminal loops could orient acceptors more or less favourably during binding, thus leading to tighter or weaker K(m) constants. The data suggested that binding of acceptors in HvXETs is guided by contributions from the conserved residues in the active sites and by the of neighbouring loops.

NMR characterization of the interactions between lyso-GM1 aqueous micelles and amyloid beta

Gangliosides are targets for a variety of pathologically relevant proteins, including amyloid beta (Abeta), an important component implicated in Alzheimer's disease (AD). To provide a structural basis for this pathogenic interaction associated with AD, we conducted NMR analyses of the Abeta interactions with gangliosides using lyso-GM1 micelles as a model system. Our NMR data revealed that the sugar-lipid interface is primarily perturbed upon binding of Abeta to the micelles, underscoring the importance of the inner part of the ganglioside cluster for accommodating Abeta in comparison with the outer carbohydrate branches that provide microbial toxin- and virus-binding sites.

Context-specific target definition in influenza a virus hemagglutinin-glycan receptor interactions

Protein-glycan interactions are important regulators of a variety of biological processes, ranging from immune recognition to anticoagulation. An important area of active research is directed toward understanding the role of host cell surface glycans as recognition sites for pathogen protein receptors. Recognition of cell surface glycans is a widely employed strategy for a variety of pathogens, including bacteria, parasites, and viruses. We present here a representative example of such an interaction: the binding of influenza A hemagglutinin (HA) to specific sialylated glycans on the cell surface of human upper airway epithelial cells, which initiates the infection cycle. We detail a generalizable strategy to understand the nature of protein-glycan interactions both structurally and biochemically, using HA as a model system. This strategy combines a top-down approach using available structural information to define important contacts between glycans and HA, with a bottom-up approach using data-mining and informatics approaches to identify the common motifs that distinguish glycan binders from nonbinders. By probing protein-glycan interactions simultaneously through top-down and bottom-up approaches, we can scientifically validate a series of observations. This in turn provides additional confidence and surmounts known challenges in the study of protein-glycan interactions, such as accounting for multivalency, and thus truly defines concepts such as specificity, affinity, and avidity. With the advent of new technologies for glycomics-including glycan arrays, data-mining solutions, and robust algorithms to model protein-glycan interactions-we anticipate that such combination approaches will become tractable for a wide variety of protein-glycan interactions.

Nuclear magnetic resonance structural characterization of substrates bound to the alpha-2,6-sialyltransferase, ST6Gal-I

The alpha-2,6-sialyltransferase (ST6Gal-I) is a key enzyme that regulates the distribution of sialic acid-containing molecules on mammalian cell surfaces. However, the fact that its native form is membrane-bound and glycosylated has made structural characterization by X-ray crystallography of this eukaryotic protein difficult. Its large size ( approximately 40 kDa for just the catalytic domain) also poses a challenge for complete structure determination by nuclear magnetic resonance (NMR). However, even without complete structure determination, there are NMR strategies that can return targeted information about select regions of the protein, including information about the active site as seen from the perspective of its bound ligands. Here, in a continuation of a previous study, a spin-labeled mimic of a glycan acceptor ligand is used to identify additional amino acids located in the protein active site. In addition, the spin-labeled donor is used to characterize the relative placement of the two bound ligands. The ligand conformation and protein-ligand contact surfaces are studied by transferred nuclear Overhauser effects (trNOEs) and saturation transfer difference (STD) experiments. The data afforded by the methods mentioned above lead to a geometric model of the bound substrates that in many ways carries an imprint of the ST6Gal-I binding site.

Mechanism of glycan receptor recognition and specificity switch for avian, swine, and human adapted influenza virus hemagglutinins: a molecular dynamics perspective

Hemagglutinins (HA's) from duck, swine, and human influenza viruses have previously been shown to prefer avian and human glycan receptor analogues with distinct topological profiles, pentasaccharides LSTa (alpha-2,3 linkage) and LSTc (alpha-2,6 linkage), in comparative molecular dynamics studies. On the basis of detailed analyses of the dynamic motions of the receptor binding domains (RBDs) and interaction energy profiles with individual glycan residues, we have identified approximately 30 residue positions in the RBD that present distinct profiles with the receptor analogues. Glycan binding constrained the conformational space sampling by the HA. Electrostatic steering appeared to play a key role in glycan binding specificity. The complex dynamic behaviors of the major SSE and trimeric interfaces with or without bound glycans suggested that networks of interactions might account for species specificity in these low affinity and high avidity (multivalent) interactions between different HA and glycans. Contact frequency, energetic decomposition, and H-bond analyses revealed species-specific differences in HA-glycan interaction profiles, not readily discernible from crystal structures alone. Interaction energy profiles indicated that mutation events at the set of residues such as 145, 156, 158, and 222 would favor human or avian receptor analogues, often through interactions with distal asialo-residues. These results correlate well with existing experimental evidence, and suggest new opportunities for simulation-based vaccine and drug development.

A 3D-structural model of unsulfated chondroitin from high-field NMR: 4-sulfation has little effect on backbone conformation

The glycosaminoglycan chondroitin sulfate is essential in human health and disease but exactly how sulfation dictates its 3D-structure at the atomic level is unclear. To address this, we have purified homogenous oligosaccharides of unsulfated chondroitin (with and without (15)N-enrichment) and analysed them by high-field NMR to make a comparison published chondroitin sulfate and hyaluronan 3D-structures. The result is the first full assignment of the tetrasaccharide and an experimental 3D-model of the hexasaccharide (PDB code 2KQO). In common with hyaluronan, we confirm that the amide proton is not involved in strong, persistent inter-residue hydrogen bonds. However, in contrast to hyaluronan, a hydrogen bond is not inferred between the hexosamine OH-4 and the glucuronic acid O5 atoms across the beta(1-->3) glycosidic linkage. The unsulfated chondroitin bond geometry differs slightly from hyaluronan by rotation about the beta(1-->3) psi dihedral (as previously predicted by simulation), while the beta(1-->4) linkage is unaffected. Furthermore, comparison shows that this glycosidic linkage geometry is similar in chondroitin-4-sulfate. We therefore hypothesise that both hexosamine OH-4 and OH-6 atoms are solvent exposed in chondroitin, explaining why it is amenable to sulfation and hyaluronan is not, and also that 4-sulfation has little effect on backbone conformation. Our conclusions exemplify the value of the 3D-model presented here and progress our understanding of glycosaminoglycan molecular properties.

Conformational free energy surface of alpha-N-acetylneuraminic acid: an interplay between hydrogen bonding and solvation

The conformational free energy surface of alpha-N-acetylneuraminic acid (Neu5Ac, sialic acid) in the space of ring-puckering coordinates was calculated using the metadynamics method. Free energy surfaces in vacuum and with an explicit solvent were calculated in GLYCAM 06 force field. In vacuum three structures are almost equivalently populated, namely, the (2)C(5) chair and the B(3,6)/(2)S(6) and (O)S(3) boat/skew-boat conformations. The B(3,6)/(2)S(6) structure is stabilized by an ionic hydrogen bond between the amide N-H bond and the carboxylic group. However, this structure is unfavorable in a water environment in which the experimentally observed (2)C(5) chair conformation is predicted to be more stable than the other structures. These results indicate that environment significantly influences conformation of Neu5Ac and that Neu5Ac-processing enzymes might modify a conformation of their substrates solely by a changing polarity of the environment. The structure of Neu5Ac bound in influenza neuraminidase ((4)S(2)/B(2,5)) belongs to conformations preferred in a water environment. The free energy penalty of this conformational change was calculated (relative to (2)C(5)) as 10.2 +/- 2.0 and 17.3 +/- 2.0 kJ/mol for (4,O)B/(O)S(3) and (4)S(2), respectively. This result indicates that mimicking of the enzyme-bound conformation is likely to be a viable strategy for the design of neuraminidase inhibitors.

Conformational entropy changes upon lactose binding to the carbohydrate recognition domain of galectin-3

The conformational entropy of proteins can make significant contributions to the free energy of ligand binding. NMR spin relaxation enables site-specific investigation of conformational entropy, via order parameters that parameterize local reorientational fluctuations of rank-2 tensors. Here we have probed the conformational entropy of lactose binding to the carbohydrate recognition domain of galectin-3 (Gal3), a protein that plays an important role in cell growth, cell differentiation, cell cycle regulation, and apoptosis, making it a potential target for therapeutic intervention in inflammation and cancer. We used (15)N spin relaxation experiments and molecular dynamics simulations to monitor the backbone amides and secondary amines of the tryptophan and arginine side chains in the ligand-free and lactose-bound states of Gal3. Overall, we observe good agreement between the experimental and computed order parameters of the ligand-free and lactose-bound states. Thus, the (15)N spin relaxation data indicate that the molecular dynamics simulations provide reliable information on the conformational entropy of the binding process. The molecular dynamics simulations reveal a correlation between the simulated order parameters and residue-specific backbone entropy, re-emphasizing that order parameters provide useful estimates of local conformational entropy. The present results show that the protein backbone exhibits an increase in conformational entropy upon binding lactose, without any accompanying structural changes.

Comparison of different force fields for the study of disaccharides

Eighteen empirical force fields and the semi-empirical quantum method PM3CARB-1 were compared for studying beta-cellobiose, alpha-maltose, and alpha-galabiose [alpha-D-Galp-(1-->4)-alpha-D-Galp]. For each disaccharide, the energies of 54 conformers with differing hydroxymethyl, hydroxyl, and glycosidic linkage orientations were minimized by the different methods, some at two dielectric constants. By comparing these results and the available crystal structure data and/or higher level density functional theory results, it was concluded that the newer parameterizations for force fields (GROMOS, GLYCAM06, OPLS-2005 and CSFF) give results that are reasonably similar to each other, whereas the older parameterizations for Amber, CHARMM or OPLS were more divergent. However, MM3, an older force field, gave energy and geometry values comparable to those of the newer parameterizations, but with less sensitivity to dielectric constant values. These systems worked better than MM2 variants, which were still acceptable. PM3CARB-1 also gave adequate results in terms of linkage and exocyclic torsion angles. GROMOS, GLYCAM06, and MM3 appear to be the best choices, closely followed by MM4, CSFF, and OPLS-2005. With GLYCAM06 and to a lesser extent, CSFF, and OPLS-2005, a number of the conformers that were stable with MM3 changed to other forms.

CHARMM Additive All-Atom Force Field for Glycosidic Linkages between Hexopyranoses

We present an extension of the CHARMM hexopyranose monosaccharide additive all-atom force field to enable modeling of glycosidic-linked hexopyranose polysaccharides. The new force field parameters encompass 1→1, 1→2, 1→3, 1→4, and 1→6 hexopyranose glycosidic linkages, as well as O-methylation at the C(1) anomeric carbon, and are developed to be consistent with the CHARMM all-atom biomolecular force fields for proteins, nucleic acids, and lipids. The parameters are developed in a hierarchical fashion using model compounds containing the key atoms in the full carbohydrates, in particular O-methyl-tetrahydropyran and glycosidic-linked dimers consisting of two molecules of tetrahyropyran or one of tetrahydropyran and one of cyclohexane. Target data for parameter optimization include full two-dimensional energy surfaces defined by the Φ/Ψ glycosidic dihedral angles in the disaccharide analogs as determined by quantum mechanical MP2/cc-pVTZ single point energies on MP2/6-31G(d) optimized structures (MP2/cc-pVTZ//MP2/6-31G(d)). In order to achieve balanced, transferable dihedral parameters for the Φ/Ψ glycosidic dihedral angles, surfaces for all possible chiralities at the ring carbon atoms involved in the glycosidic linkages are considered, resulting in over 5000 MP2/cc-pVTZ//MP2/6-31G(d) conformational energies. Also included as target data are vibrational frequencies, pair interaction energies and distances with water molecules, and intramolecular geometries including distortion of the glycosidic valence angle as a function of the glycosidic dihedral angles. The model-compound optimized force field parameters are validated on full disaccharides through comparison of molecular dynamics results to available experimental data. Good agreement is achieved with experiment for a variety of properties including crystal cell parameters and intramolecular geometries, aqueous densities, and aqueous NMR coupling constants associated with the glycosidic linkage. The newly-developed parameters allow for the modeling of linear, branched, and cyclic hexopyranose glycosides both alone and in heterogenous systems including proteins, nucleic acids and/or lipids when combined with existing CHARMM biomolecular force fields.

7th HUPO World Congress: the human disease glycomics/proteomics initiative (HGPI) session 17 August 2008, Amsterdam, The Netherlands

The Human Disease Glycomics/Proteomics Initiative (HGPI) Session was held on August 17, 2008, at the HUPO World Congress in Amsterdam. Reports were made on the progress of the first and second analytical pilot studies to profile N- and O-linked glycan structures of standard glycoproteins utilizing laboratories from around the world. In addition, recent advances in glycan structural analyses were presented, including the use of O-linked glycan libraries of standards, use of negative mode nano-LC-MS for O-linked glycan analysis, and identification of aberrant O-glycosylation of IgA1 as a cause of IgA nephropathy. A report was made of a newly discovered lectin, malectin, which appears to function in the folding/quality control of glycoproteins with N-linked glycans and may regulate several human disorders whose etiology involves protein quality control in the ER. The major glycan ligand for malectin was identified using a novel printed glycan microarray. Advances in the analysis of the genes that are associated with glycan expression and recognition--the glycotranscriptome--were described, as well as technologies to determine the relative quantitation of N- and O-linked glycans from as few as 2 x 10(6) cells. These technologies are being applied to identify potential biomarkers of stem and cancer cells.

Gangliosides in cell recognition and membrane protein regulation

Gangliosides, sialic acid-bearing glycosphingolipids, are expressed on all vertebrate cells, and are the major glycans on nerve cells. They are anchored to the plasma membrane through their ceramide lipids with their varied glycans extending into the extracellular space. Through sugar-specific interactions with glycan-binding proteins on apposing cells, gangliosides function as receptors in cell-cell recognition, regulating natural killer cell cytotoxicity via Siglec-7, myelin-axon interactions via Siglec-4 (myelin-associated glycoprotein), and inflammation via E-selectin. Gangliosides also interact laterally in their own membranes, regulating the responsiveness of signaling proteins including the insulin, epidermal growth factor, and vascular endothelial growth factor receptors. In these ways, gangliosides act as regulatory elements in the immune system, in the nervous system, in metabolic regulation, and in cancer progression.

Metadynamics modelling of the solvent effect on primary hydroxyl rotamer equilibria in hexopyranosides

Accurate modelling of rotamer equilibria for the primary hydroxyl groups of monosaccharides continues to be a great challenge of computational glycochemistry. The metadynamics technique was applied to study the conformational free energy surfaces of methyl alpha-D-glucopyranoside and methyl alpha-D-galactopyranoside, employing the GLYCAM06 force field. For both molecules, seven to eight conformational free-energy minima, differing in the omega (O-5-C-5-C-6-O-6) and chi (C-3-C-4-O-4-HO-4) dihedral angles, were identified in vacuum or in a water environment. The calculated rotamer equilibrium of the primary hydroxyl group is significantly different in vacuum than in water. The major effect of a water environment is the destabilisation of a hydrogen bond between O-4-HO-4 and O-6-HO-6 groups. It was possible to calculate the free-energy differences of individual rotamers with an accuracy of better than 2 kJ/mol. The calculated gg, gt and tg rotamer populations in water are in close agreement with experimental measurements, and therefore support the theoretical background of metadynamics.

Force-field parameters of the Psi and Phi around glycosidic bonds to oxygen and sulfur atoms

The Psi and Phi torsion angles around glycosidic bonds in a glycoside chain are the most important determinants of the conformation of a glycoside chain. We determined force-field parameters for Psi and Phi torsion angles around a glycosidic bond bridged by a sulfur atom, as well as a bond bridged by an oxygen atom as a preparation for the next study, i.e., molecular dynamics free energy calculations for protein-sugar and protein-inhibitor complexes. First, we extracted the Psi or Phi torsion energy component from a quantum mechanics (QM) total energy by subtracting all the molecular mechanics (MM) force-field components except for the Psi or Phi torsion angle. The Psi and Phi energy components extracted (hereafter called "the remaining energy components") were calculated for simple sugar models and plotted as functions of the Psi and Phi angles. The remaining energy component curves of Psi and Phi were well represented by the torsion force-field functions consisting of four and three cosine functions, respectively. To confirm the reliability of the force-field parameters and to confirm its compatibility with other force-fields, we calculated adiabatic potential curves as functions of Psi and Phi for the model glycosides by adopting the Psi and Phi force-field parameters obtained and by energetically optimizing other degrees of freedom. The MM potential energy curves obtained for Psi and Phi well represented the QM adiabatic curves and also these curves' differences with regard to the glycosidic oxygen and sulfur atoms. Our Psi and Phi force-fields of glycosidic oxygen gave MM potential energy curves that more closely represented the respective QM curves than did those of the recently developed GLYCAM force-field.

CHARMM Additive All-Atom Force Field for Acyclic Polyalcohols, Acyclic Carbohydrates and Inositol

Parametrization of the additive all-atom CHARMM force field for acyclic polyalcohols, acyclic carbohydrates and inositol is conducted. Initial parameters were transferred from the alkanes and hexopyranose carbohydrates, with subsequent development and optimization of parameters unique to the molecules considered in this study. Using the model compounds acetone and acetaldehyde, nonbonded parameters for carbonyls were optimized targeting quantum mechanical interaction data for solute-water pairs and pure solvent thermodynamic data. Bond and angle parameters were adjusted by comparing optimized geometries to small molecule crystal survey data and by performing vibrational analyses on acetone, acetaldehyde and glycerol. C-C-C-C, C-C-C-O, C-C-OH and O-C-C-O torsional parameters for polyol chains were fit to quantum mechanical dihedral potential energy scans comprising over 1500 RIMP2/cc-pVTZ//MP2/6-31G(d) conformations using an automated Monte Carlo simulated annealing procedure. Comparison of computed condensed-phase data, including crystal lattice parameters and densities, NMR proton-proton couplings, densities and diffusion coefficients of aqueous solutions, to experimental data validated the optimized parameters. Parameter development for these compounds proved particularly challenging because of the flexibility of the acyclic sugars and polyalcohols as well as the intramolecular hydrogen bonding between vicinal hydroxyls for all of the compounds. The newly optimized additive CHARMM force field parameters are anticipated to be of utility for atomic level of detail simulations of acyclic polyalcohols, acyclic carbohydrates and inositol in solution.

Computational approaches for the design of peptides with anti-breast cancer properties

Background: Breast cancer is the most common cancer among women. Tamoxifen is the preferred drug for estrogen receptor-positive breast cancer treatment, yet many of these cancers are intrinsically resistant to tamoxifen or acquire resistance during treatment. Therefore, scientists are searching for breast cancer drugs that have different molecular targets. Methodology: Recently, a computational approach was used to successfully design peptides that are new lead compounds against breast cancer. We used replica exchange molecular dynamics to predict the structure and dynamics of active peptides, leading to the discovery of smaller bioactive peptides. Conclusions: These analogs inhibit estrogen-dependent cell growth in a mouse uterine growth assay, a test showing reliable correlation with human breast cancer inhibition. We outline the computational methods that were tried and used along with the experimental information that led to the successful completion of this research.

The conformational properties of methyl alpha-(2,8)-di/trisialosides and their N-acyl analogues: implications for anti-Neisseria meningitidis B vaccine design

The conformational properties of di- and trisaccharide fragments of the polysialic acid O-antigen capsular polysaccharide (CPS) of Neisseria meningitidis B (NmB) have been investigated by a combination of solution phase NMR spectroscopy and explicit-solvent molecular dynamics (MD) simulations. Simulations employing 100 ns of conventional MD, as well as 160 ns of replica exchange MD (REMD), with the GLYCAM06 force field were shown to be in agreement with experimental NMR scalar J-coupling and NOE values. The presence of conformational families has been determined by monitoring interglycosidic torsion angles, by comparing structural superimpositions, as well as via a Bayesian statistical analysis of the torsional data. Attempts to augment the immunogenicity of NmB CPS often involve chemical modifications of the N-acetyl moiety. Here the effects of these chemical group modifications on the conformational properties of the trisialoside have been probed via REMD simulations of the N-glycolyl, N-propionyl, N-propyl and N-butanoyl analogues. Although there were conformational families unique to each non-native analogue, the chemical modifications resulted in largely equivalent overall conformational phase-spaces compared to the native trisialoside. On the basis of the conformational distributions, these shared conformational properties suggest that a recurrent global conformational epitope may be present in both the native and chemically modified CPS fragments. Explanations are therefore provided for monoclonal antibody cross-reactivity, in terms of recognition of a shared global CPS conformation, as well as for lack of cross-reactivity, in terms of fine structural differences associated with the N-acyl groups, which may be dominant in highly matured antibody responses.