1551
|
Woods RJ, Tessier MB. Computational glycoscience: characterizing the spatial and temporal properties of glycans and glycan-protein complexes. Curr Opin Struct Biol 2010; 20:575-83. [PMID: 20708922 DOI: 10.1016/j.sbi.2010.07.005] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2010] [Revised: 07/08/2010] [Accepted: 07/19/2010] [Indexed: 01/09/2023]
Abstract
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.
Collapse
Affiliation(s)
- Robert J Woods
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, GA 30602, USA.
| | | |
Collapse
|
1552
|
Tan L, Varnai P, Lamport DTA, Yuan C, Xu J, Qiu F, Kieliszewski MJ. Plant O-hydroxyproline arabinogalactans are composed of repeating trigalactosyl subunits with short bifurcated side chains. J Biol Chem 2010; 285:24575-83. [PMID: 20489210 PMCID: PMC2915693 DOI: 10.1074/jbc.m109.100149] [Citation(s) in RCA: 83] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2010] [Revised: 04/14/2010] [Indexed: 11/06/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Li Tan
- From the Department of Chemistry and Biochemistry, Biochemistry Research Facility, Ohio University, Athens, Ohio 45701
| | - Peter Varnai
- the Department of Chemistry and Biochemistry and
| | - Derek T. A. Lamport
- the School of Life Sciences, University of Sussex, Falmer, Brighton BN1 9QG, United Kingdom
| | - Chunhua Yuan
- the Campus Chemical Instrument Center, Ohio State University, Columbus, Ohio 43210, and
| | - Jianfeng Xu
- From the Department of Chemistry and Biochemistry, Biochemistry Research Facility, Ohio University, Athens, Ohio 45701
| | - Feng Qiu
- the Department of Chemistry, Temple University, Philadelphia, Pennsylvaina 19122
| | - Marcia J. Kieliszewski
- From the Department of Chemistry and Biochemistry, Biochemistry Research Facility, Ohio University, Athens, Ohio 45701
| |
Collapse
|
1553
|
Frank M, Schloissnig S. Bioinformatics and molecular modeling in glycobiology. Cell Mol Life Sci 2010; 67:2749-72. [PMID: 20364395 PMCID: PMC2912727 DOI: 10.1007/s00018-010-0352-4] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2009] [Revised: 03/08/2010] [Accepted: 03/11/2010] [Indexed: 12/11/2022]
Abstract
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.
Collapse
Affiliation(s)
- Martin Frank
- Molecular Structure Analysis Core Facility-W160, Deutsches Krebsforschungszentrum (German Cancer Research Centre), 69120 Heidelberg, Germany.
| | | |
Collapse
|
1554
|
Dupradeau FY, Pigache A, Zaffran T, Savineau C, Lelong R, Grivel N, Lelong D, Rosanski W, Cieplak P. The R.E.D. tools: advances in RESP and ESP charge derivation and force field library building. Phys Chem Chem Phys 2010; 12:7821-39. [PMID: 20574571 PMCID: PMC2918240 DOI: 10.1039/c0cp00111b] [Citation(s) in RCA: 677] [Impact Index Per Article: 48.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- François-Yves Dupradeau
- CNRS UMR 6219 & Faculté de Pharmacie, Université de Picardie - Jules Verne, 1, rue des Louvels, F-80037 Amiens cedex 1, France
| | - Adrien Pigache
- CNRS UMR 6219 & Faculté de Pharmacie, Université de Picardie - Jules Verne, 1, rue des Louvels, F-80037 Amiens cedex 1, France
| | - Thomas Zaffran
- CNRS UMR 6219 & Faculté de Pharmacie, Université de Picardie - Jules Verne, 1, rue des Louvels, F-80037 Amiens cedex 1, France
| | - Corentin Savineau
- CNRS UMR 6219 & Faculté de Pharmacie, Université de Picardie - Jules Verne, 1, rue des Louvels, F-80037 Amiens cedex 1, France
| | - Rodolphe Lelong
- CNRS UMR 6219 & Faculté de Pharmacie, Université de Picardie - Jules Verne, 1, rue des Louvels, F-80037 Amiens cedex 1, France
| | - Nicolas Grivel
- CNRS UMR 6219 & Faculté de Pharmacie, Université de Picardie - Jules Verne, 1, rue des Louvels, F-80037 Amiens cedex 1, France
| | - Dimitri Lelong
- CNRS UMR 6219 & Faculté de Pharmacie, Université de Picardie - Jules Verne, 1, rue des Louvels, F-80037 Amiens cedex 1, France
| | - Wilfried Rosanski
- CNRS UMR 6219 & Faculté de Pharmacie, Université de Picardie - Jules Verne, 1, rue des Louvels, F-80037 Amiens cedex 1, France
| | - Piotr Cieplak
- Sanford|Burnham Institute for Medical Research, La Jolla, California, 92037, USA
| |
Collapse
|
1555
|
Hernández-Gay JJ, Ardá A, Eller S, Mezzato S, Leeflang BR, Unverzagt C, Cañada FJ, Jiménez-Barbero J. Insights into the Dynamics and Molecular Recognition Features of Glycopeptides by Protein Receptors: The 3D Solution Structure of Hevein Bound to the Trisaccharide Core of N-Glycoproteins. Chemistry 2010; 16:10715-26. [DOI: 10.1002/chem.201000939] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
1556
|
Morando M, Yao Y, Martín-Santamaría S, Zhu Z, Xu T, Cañada FJ, Zhang Y, Jiménez-Barbero J. Mimicking chitin: chemical synthesis, conformational analysis, and molecular recognition of the beta(1-->3) N-acetylchitopentaose analogue. Chemistry 2010; 16:4239-49. [PMID: 20229531 DOI: 10.1002/chem.200902860] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
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.
Collapse
Affiliation(s)
- Maria Morando
- Chemical and Physical Biology, Centro de Investigacione, Biológicas, CSIC, Ramiro de Maeztu 9, 28040 Madrid, Spain
| | | | | | | | | | | | | | | |
Collapse
|
1557
|
Nurisso A, Blanchard B, Audfray A, Rydner L, Oscarson S, Varrot A, Imberty A. Role of water molecules in structure and energetics of Pseudomonas aeruginosa lectin I interacting with disaccharides. J Biol Chem 2010; 285:20316-27. [PMID: 20410292 PMCID: PMC2888444 DOI: 10.1074/jbc.m110.108340] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2010] [Revised: 04/11/2010] [Indexed: 11/06/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Alessandra Nurisso
- From the Centre de Rechèrche sur les Macromolécules Végétales-CNRS (affiliated with Université Joseph Fourier and Institut de Chimie Moléculaire de Grenoble), BP 53, 38041 Grenoble Cedex 9, France and
| | - Bertrand Blanchard
- From the Centre de Rechèrche sur les Macromolécules Végétales-CNRS (affiliated with Université Joseph Fourier and Institut de Chimie Moléculaire de Grenoble), BP 53, 38041 Grenoble Cedex 9, France and
| | - Aymeric Audfray
- From the Centre de Rechèrche sur les Macromolécules Végétales-CNRS (affiliated with Université Joseph Fourier and Institut de Chimie Moléculaire de Grenoble), BP 53, 38041 Grenoble Cedex 9, France and
| | - Lina Rydner
- the Centre for Synthesis and Chemical Biology, School of Chemistry and Chemical Biology, University College Dublin, Belfield, Dublin 4, Ireland
| | - Stefan Oscarson
- the Centre for Synthesis and Chemical Biology, School of Chemistry and Chemical Biology, University College Dublin, Belfield, Dublin 4, Ireland
| | - Annabelle Varrot
- From the Centre de Rechèrche sur les Macromolécules Végétales-CNRS (affiliated with Université Joseph Fourier and Institut de Chimie Moléculaire de Grenoble), BP 53, 38041 Grenoble Cedex 9, France and
| | - Anne Imberty
- From the Centre de Rechèrche sur les Macromolécules Végétales-CNRS (affiliated with Université Joseph Fourier and Institut de Chimie Moléculaire de Grenoble), BP 53, 38041 Grenoble Cedex 9, France and
| |
Collapse
|
1558
|
Säwén E, Massad T, Landersjö C, Damberg P, Widmalm G. 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. Org Biomol Chem 2010; 8:3684-95. [PMID: 20574564 DOI: 10.1039/c003958f] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- Elin Säwén
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, S-106 91, Stockholm, Sweden
| | | | | | | | | |
Collapse
|
1559
|
Less is more when simulating unsulfated glycosaminoglycan 3D-structure: Comparison of GLYCAM06/TIP3P, PM3-CARB1/TIP3P, and SCC-DFTB-D/TIP3P predictions with experiment. J Comput Chem 2010; 31:2932-47. [DOI: 10.1002/jcc.21589] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
1560
|
Fadda E, Woods RJ. Molecular simulations of carbohydrates and protein-carbohydrate interactions: motivation, issues and prospects. Drug Discov Today 2010; 15:596-609. [PMID: 20594934 DOI: 10.1016/j.drudis.2010.06.001] [Citation(s) in RCA: 148] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2009] [Revised: 04/20/2010] [Accepted: 06/01/2010] [Indexed: 11/16/2022]
Abstract
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.
Collapse
Affiliation(s)
- Elisa Fadda
- School of Chemistry, National University of Ireland, Galway, Ireland
| | | |
Collapse
|
1561
|
Petersen L, Ardèvol A, Rovira C, Reilly PJ. Molecular Mechanism of the Glycosylation Step Catalyzed by Golgi α-Mannosidase II: A QM/MM Metadynamics Investigation. J Am Chem Soc 2010; 132:8291-300. [DOI: 10.1021/ja909249u] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Luis Petersen
- Department of Chemical and Biological Engineering, Iowa State University, Ames, Iowa 50011, Computer Simulation and Modeling Laboratory (CoSMoLab), Parc Científic de Barcelona, 08028 Barcelona, Spain, Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain, and Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Spain
| | - Albert Ardèvol
- Department of Chemical and Biological Engineering, Iowa State University, Ames, Iowa 50011, Computer Simulation and Modeling Laboratory (CoSMoLab), Parc Científic de Barcelona, 08028 Barcelona, Spain, Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain, and Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Spain
| | - Carme Rovira
- Department of Chemical and Biological Engineering, Iowa State University, Ames, Iowa 50011, Computer Simulation and Modeling Laboratory (CoSMoLab), Parc Científic de Barcelona, 08028 Barcelona, Spain, Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain, and Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Spain
| | - Peter J. Reilly
- Department of Chemical and Biological Engineering, Iowa State University, Ames, Iowa 50011, Computer Simulation and Modeling Laboratory (CoSMoLab), Parc Científic de Barcelona, 08028 Barcelona, Spain, Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain, and Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Spain
| |
Collapse
|
1562
|
Deutschmann R, Boncheff AG, Daraban L, MacInnes JI, Monteiro MA. Common sialylated glycan in Actinobacillus suis. Glycobiology 2010; 20:1227-32. [PMID: 20501522 DOI: 10.1093/glycob/cwq079] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
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.
Collapse
|
1563
|
Koppisetty CAK, Nasir W, Strino F, Rydell GE, Larson G, Nyholm PG. Computational studies on the interaction of ABO-active saccharides with the norovirus VA387 capsid protein can explain experimental binding data. J Comput Aided Mol Des 2010; 24:423-31. [PMID: 20407802 DOI: 10.1007/s10822-010-9353-5] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2009] [Accepted: 03/29/2010] [Indexed: 01/08/2023]
Abstract
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.
Collapse
|
1564
|
Demarco ML, Woods RJ, Prestegard JH, Tian F. Presentation of membrane-anchored glycosphingolipids determined from molecular dynamics simulations and NMR paramagnetic relaxation rate enhancement. J Am Chem Soc 2010; 132:1334-8. [PMID: 20058858 DOI: 10.1021/ja907518x] [Citation(s) in RCA: 53] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Mari L Demarco
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, Georgia 30602, USA
| | | | | | | |
Collapse
|
1565
|
Babin V, Sagui C. Conformational free energies of methyl-α-L-iduronic and methyl-β-D-glucuronic acids in water. J Chem Phys 2010; 132:104108. [DOI: 10.1063/1.3355621] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
1566
|
Frank AT, Ramsook CB, Otoo HN, Tan C, Soybelman G, Rauceo JM, Gaur NK, Klotz SA, Lipke PN. Structure and function of glycosylated tandem repeats from Candida albicans Als adhesins. EUKARYOTIC CELL 2010; 9:405-14. [PMID: 19820118 PMCID: PMC2837987 DOI: 10.1128/ec.00235-09] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/16/2009] [Accepted: 10/02/2009] [Indexed: 01/09/2023]
Abstract
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.
Collapse
Affiliation(s)
- Aaron T. Frank
- Department of Biophysics, University of Michigan, Ann Arbor, Michigan 48109
| | - Caleen B. Ramsook
- Department of Biology, Brooklyn College of CUNY, Brooklyn, New York 11210
| | - Henry N. Otoo
- Department of Biology, Brooklyn College of CUNY, Brooklyn, New York 11210
| | - Cho Tan
- Department of Biology, Brooklyn College of CUNY, Brooklyn, New York 11210
| | - Gregory Soybelman
- Department of Biology, Brooklyn College of CUNY, Brooklyn, New York 11210
| | - Jason M. Rauceo
- Department of Biology, John Jay College of Criminal Justice of CUNY, New York, New York 10019; and
| | - Nand K. Gaur
- University of Arizona and Southern Arizona VA Health Care System, Tucson, Arizona 85724
| | - Stephen A. Klotz
- University of Arizona and Southern Arizona VA Health Care System, Tucson, Arizona 85724
| | - Peter N. Lipke
- Department of Biology, Brooklyn College of CUNY, Brooklyn, New York 11210
| |
Collapse
|
1567
|
Kozlovskaya L, Osolodkin D, Shevtsova A, Romanova L, Rogova Y, Dzhivanian T, Lyapustin V, Pivanova G, Gmyl A, Palyulin V, Karganova G. GAG-binding variants of tick-borne encephalitis virus. Virology 2010; 398:262-72. [DOI: 10.1016/j.virol.2009.12.012] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2009] [Revised: 11/25/2009] [Accepted: 12/10/2009] [Indexed: 11/30/2022]
|
1568
|
Gandhi NS, Mancera RL. Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study. Carbohydr Res 2010; 345:689-95. [DOI: 10.1016/j.carres.2009.12.020] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2009] [Revised: 11/30/2009] [Accepted: 12/20/2009] [Indexed: 10/20/2022]
|
1569
|
Vorontsov II, Miyashita O. Solution and crystal molecular dynamics simulation study of m4-cyanovirin-N mutants complexed with di-mannose. Biophys J 2010; 97:2532-40. [PMID: 19883596 DOI: 10.1016/j.bpj.2009.08.011] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2009] [Revised: 08/03/2009] [Accepted: 08/04/2009] [Indexed: 11/19/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Ivan I Vorontsov
- Department of Chemistry and Biochemistry, University of Arizona, Tucson, Arizona, USA
| | | |
Collapse
|
1570
|
Meynier C, Feracci M, Espeli M, Chaspoul F, Gallice P, Schiff C, Guerlesquin F, Roche P. NMR and MD investigations of human galectin-1/oligosaccharide complexes. Biophys J 2010; 97:3168-77. [PMID: 20006954 DOI: 10.1016/j.bpj.2009.09.026] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2008] [Revised: 08/31/2009] [Accepted: 09/10/2009] [Indexed: 12/18/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Christophe Meynier
- Unité Interactions et Modulateurs de Réponses, Institut de Microbiologie de la Méditerrannée, Centre National de la Recherche Scientifique, Marseille, France
| | | | | | | | | | | | | | | |
Collapse
|
1571
|
Sharma A, Sekar K, Vijayan M. Structure, dynamics, and interactions of jacalin. Insights from molecular dynamics simulations examined in conjunction with results of X-ray studies. Proteins 2010; 77:760-77. [PMID: 19544573 DOI: 10.1002/prot.22486] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
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.
Collapse
Affiliation(s)
- Alok Sharma
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India 560 012
| | | | | |
Collapse
|
1572
|
Modelling of β-d-glucopyranose ring distortion in different force fields: a metadynamics study. Carbohydr Res 2010; 345:530-7. [DOI: 10.1016/j.carres.2009.12.011] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2009] [Revised: 12/02/2009] [Accepted: 12/11/2009] [Indexed: 12/16/2022]
|
1573
|
Vaaje-Kolstad G, Farkas V, Fincher GB, Hrmova M. Barley xyloglucan xyloglucosyl transferases bind xyloglucan-derived oligosaccharides in their acceptor-binding regions in multiple conformational states. Arch Biochem Biophys 2010; 496:61-8. [PMID: 20117073 DOI: 10.1016/j.abb.2010.01.011] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2009] [Revised: 01/21/2010] [Accepted: 01/24/2010] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Gustav Vaaje-Kolstad
- Department of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, 1432 As, Norway
| | | | | | | |
Collapse
|
1574
|
Yagi-Utsumi M, Kameda T, Yamaguchi Y, Kato K. NMR characterization of the interactions between lyso-GM1 aqueous micelles and amyloid beta. FEBS Lett 2010; 584:831-6. [PMID: 20074569 DOI: 10.1016/j.febslet.2010.01.005] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2009] [Revised: 12/27/2009] [Accepted: 01/05/2010] [Indexed: 11/27/2022]
Abstract
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.
Collapse
Affiliation(s)
- Maho Yagi-Utsumi
- Graduate School of Pharmaceutical Sciences, Nagoya City University, Mizuho-ku, Nagoya, Japan
| | | | | | | |
Collapse
|
1575
|
Raju RK, Hillier IH, Burton NA, Vincent MA, Doudou S, Bryce RA. The effects of perfluorination on carbohydrate–π interactions: computational studies of the interaction of benzene and hexafluorobenzene with fucose and cyclodextrin. Phys Chem Chem Phys 2010; 12:7959-67. [DOI: 10.1039/c002058c] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
1576
|
Kövér KE, Szilágyi L, Batta G, Uhrín D, Jiménez-Barbero J. Biomolecular Recognition by Oligosaccharides and Glycopeptides: The NMR Point of View. COMPREHENSIVE NATURAL PRODUCTS II 2010:197-246. [DOI: 10.1016/b978-008045382-8.00193-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/27/2023]
|
1577
|
Chlubnová I, Filipp D, Spiwok V, Dvořáková H, Daniellou R, Nugier-Chauvin C, Králová B, Ferrières V. Enzymatic synthesis of oligo-d-galactofuranosides and l-arabinofuranosides: from molecular dynamics to immunological assays. Org Biomol Chem 2010; 8:2092-102. [DOI: 10.1039/b926988f] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
1578
|
Shriver Z, Raman R, Viswanathan K, Sasisekharan R. Context-specific target definition in influenza a virus hemagglutinin-glycan receptor interactions. ACTA ACUST UNITED AC 2009; 16:803-14. [PMID: 19716471 DOI: 10.1016/j.chembiol.2009.08.002] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2009] [Revised: 07/28/2009] [Accepted: 08/03/2009] [Indexed: 12/20/2022]
Abstract
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.
Collapse
Affiliation(s)
- Zachary Shriver
- Koch Institute for Integrative Cancer Research, Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, 02139, USA
| | | | | | | |
Collapse
|
1579
|
Liu S, Meng L, Moremen KW, Prestegard JH. Nuclear magnetic resonance structural characterization of substrates bound to the alpha-2,6-sialyltransferase, ST6Gal-I. Biochemistry 2009; 48:11211-9. [PMID: 19845399 DOI: 10.1021/bi9015154] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Shan Liu
- Complex Carbohydrate Research Center, University of Georgia, Athens, Georgia 30602, USA
| | | | | | | |
Collapse
|
1580
|
Newhouse EI, Xu D, Markwick PRL, Amaro RE, Pao HC, Wu KJ, Alam M, McCammon JA, Li WW. Mechanism of glycan receptor recognition and specificity switch for avian, swine, and human adapted influenza virus hemagglutinins: a molecular dynamics perspective. J Am Chem Soc 2009; 131:17430-42. [PMID: 19891427 PMCID: PMC2782351 DOI: 10.1021/ja904052q] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2009] [Indexed: 11/28/2022]
Abstract
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.
Collapse
|
1581
|
Sattelle BM, Shakeri J, Roberts IS, Almond A. A 3D-structural model of unsulfated chondroitin from high-field NMR: 4-sulfation has little effect on backbone conformation. Carbohydr Res 2009; 345:291-302. [PMID: 20022001 PMCID: PMC3098369 DOI: 10.1016/j.carres.2009.11.013] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2009] [Revised: 11/06/2009] [Accepted: 11/10/2009] [Indexed: 11/26/2022]
Abstract
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.
Collapse
Affiliation(s)
- Benedict M Sattelle
- Manchester Interdisciplinary Biocentre, 131 Princess Street, Manchester, M1 7DN, UK
| | | | | | | |
Collapse
|
1582
|
Conformational Analysis of Trimannoside and Bisected Trimannoside Using Aqueous Molecular Dynamics Simulations. B KOREAN CHEM SOC 2009. [DOI: 10.5012/bkcs.2009.30.11.2723] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
1583
|
López CA, Rzepiela AJ, de Vries AH, Dijkhuizen L, Hünenberger PH, Marrink SJ. Martini Coarse-Grained Force Field: Extension to Carbohydrates. J Chem Theory Comput 2009; 5:3195-210. [DOI: 10.1021/ct900313w] [Citation(s) in RCA: 307] [Impact Index Per Article: 20.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Cesar A. López
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), Centre for Carbohydrate Bioprocessing, University of Groningen, Kerklaan 30, 9751 NN Haren, The Netherlands, and Laboratorium für Physikalische Chemie, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Andrzej J. Rzepiela
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), Centre for Carbohydrate Bioprocessing, University of Groningen, Kerklaan 30, 9751 NN Haren, The Netherlands, and Laboratorium für Physikalische Chemie, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Alex H. de Vries
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), Centre for Carbohydrate Bioprocessing, University of Groningen, Kerklaan 30, 9751 NN Haren, The Netherlands, and Laboratorium für Physikalische Chemie, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Lubbert Dijkhuizen
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), Centre for Carbohydrate Bioprocessing, University of Groningen, Kerklaan 30, 9751 NN Haren, The Netherlands, and Laboratorium für Physikalische Chemie, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Philippe H. Hünenberger
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), Centre for Carbohydrate Bioprocessing, University of Groningen, Kerklaan 30, 9751 NN Haren, The Netherlands, and Laboratorium für Physikalische Chemie, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Siewert J. Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute & Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands, Groningen Biomolecular Sciences and Biotechnology Institute (GBB), Centre for Carbohydrate Bioprocessing, University of Groningen, Kerklaan 30, 9751 NN Haren, The Netherlands, and Laboratorium für Physikalische Chemie, ETH Zürich, CH-8093 Zürich, Switzerland
| |
Collapse
|
1584
|
Xia J, Margulis CJ. Computational Study of the Conformational Structures of Saccharides in Solution Based on J Couplings and the “Fast Sugar Structure Prediction Software”. Biomacromolecules 2009; 10:3081-8. [DOI: 10.1021/bm900756q] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Junchao Xia
- Department of Chemistry, University of Iowa, Iowa City, Iowa 52242
| | | |
Collapse
|
1585
|
Spiwok V, Tvaroska I. Conformational free energy surface of alpha-N-acetylneuraminic acid: an interplay between hydrogen bonding and solvation. J Phys Chem B 2009; 113:9589-94. [PMID: 19374421 DOI: 10.1021/jp8113495] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Vojtech Spiwok
- Department of Structure and Function of Saccharides, Centre for Glycomics, Slovak Academy of Sciences, 84538 Bratislava, Slovak Republic.
| | | |
Collapse
|
1586
|
Diehl C, Genheden S, Modig K, Ryde U, Akke M. Conformational entropy changes upon lactose binding to the carbohydrate recognition domain of galectin-3. JOURNAL OF BIOMOLECULAR NMR 2009; 45:157-169. [PMID: 19641853 DOI: 10.1007/s10858-009-9356-5] [Citation(s) in RCA: 62] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/21/2009] [Accepted: 06/29/2009] [Indexed: 05/28/2023]
Abstract
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.
Collapse
Affiliation(s)
- Carl Diehl
- Center for Molecular Protein Science, Biophysical Chemistry, Lund University, Lund, Sweden
| | | | | | | | | |
Collapse
|
1587
|
Stortz CA, Johnson GP, French AD, Csonka GI. Comparison of different force fields for the study of disaccharides. Carbohydr Res 2009; 344:2217-28. [PMID: 19758584 DOI: 10.1016/j.carres.2009.08.019] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2009] [Revised: 08/13/2009] [Accepted: 08/18/2009] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Carlos A Stortz
- Departamento de Química Orgánica-CIHIDECAR, FCEyN-Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina.
| | | | | | | |
Collapse
|
1588
|
Guvench O, Hatcher ER, Venable RM, Pastor RW, MacKerell AD. CHARMM Additive All-Atom Force Field for Glycosidic Linkages between Hexopyranoses. J Chem Theory Comput 2009; 5:2353-2370. [PMID: 20161005 PMCID: PMC2757763 DOI: 10.1021/ct900242e] [Citation(s) in RCA: 489] [Impact Index Per Article: 32.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Olgun Guvench
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20 Penn Street HSF II, Baltimore MD 21201
| | - Elizabeth R. Hatcher
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20 Penn Street HSF II, Baltimore MD 21201
| | - Richard M. Venable
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892
| | - Richard W. Pastor
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20892
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, University of Maryland School of Pharmacy, 20 Penn Street HSF II, Baltimore MD 21201
| |
Collapse
|
1589
|
Pierce JM, Taniguchi N. 7th HUPO World Congress: the human disease glycomics/proteomics initiative (HGPI) session 17 August 2008, Amsterdam, The Netherlands. Proteomics 2009; 9:1738-41. [PMID: 19294704 DOI: 10.1002/pmic.200900007] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- J Michael Pierce
- Department of Biochemistry and Molecular Biology, Complex Carbohydrate Research Center, UGA Cancer Center, University of Georgia, Athens, GA, USA
| | | |
Collapse
|
1590
|
Lopez PHH, Schnaar RL. Gangliosides in cell recognition and membrane protein regulation. Curr Opin Struct Biol 2009; 19:549-57. [PMID: 19608407 DOI: 10.1016/j.sbi.2009.06.001] [Citation(s) in RCA: 218] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2009] [Accepted: 06/10/2009] [Indexed: 12/16/2022]
Abstract
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.
Collapse
Affiliation(s)
- Pablo H H Lopez
- Department of Pharmacology and Molecular Sciences, The Johns Hopkins School of Medicine, 725 N. Wolfe Street, Baltimore, MD 21205, USA
| | | |
Collapse
|
1591
|
Spiwok V, Tvaroska I. Metadynamics modelling of the solvent effect on primary hydroxyl rotamer equilibria in hexopyranosides. Carbohydr Res 2009; 344:1575-81. [PMID: 19539272 DOI: 10.1016/j.carres.2009.05.019] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2009] [Revised: 04/29/2009] [Accepted: 05/13/2009] [Indexed: 11/16/2022]
Abstract
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.
Collapse
Affiliation(s)
- Vojtech Spiwok
- Department of Biochemistry and Microbiology, Institute of Chemical Technology in Prague, Technická 3, 166 28 Prague 6, Czech Republic.
| | | |
Collapse
|
1592
|
Salisburg AM, Deline AL, Lexa KW, Shields GC, Kirschner KN. Ramachandran-type plots for glycosidic linkages: Examples from molecular dynamic simulations using the Glycam06 force field. J Comput Chem 2009; 30:910-21. [DOI: 10.1002/jcc.21099] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
1593
|
Petersen L, Ardèvol A, Rovira C, Reilly PJ. Mechanism of Cellulose Hydrolysis by Inverting GH8 Endoglucanases: A QM/MM Metadynamics Study. J Phys Chem B 2009; 113:7331-9. [DOI: 10.1021/jp811470d] [Citation(s) in RCA: 85] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Luis Petersen
- Department of Chemical and Biological Engineering, Iowa State University, Ames, Iowa 50011, Computer Simulation and Modeling Laboratory (CoSMoLab), Parc Científic de Barcelona, 08028 Barcelona, Spain, Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain, and Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Spain
| | - Albert Ardèvol
- Department of Chemical and Biological Engineering, Iowa State University, Ames, Iowa 50011, Computer Simulation and Modeling Laboratory (CoSMoLab), Parc Científic de Barcelona, 08028 Barcelona, Spain, Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain, and Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Spain
| | - Carme Rovira
- Department of Chemical and Biological Engineering, Iowa State University, Ames, Iowa 50011, Computer Simulation and Modeling Laboratory (CoSMoLab), Parc Científic de Barcelona, 08028 Barcelona, Spain, Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain, and Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Spain
| | - Peter J. Reilly
- Department of Chemical and Biological Engineering, Iowa State University, Ames, Iowa 50011, Computer Simulation and Modeling Laboratory (CoSMoLab), Parc Científic de Barcelona, 08028 Barcelona, Spain, Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08028 Barcelona, Spain, and Institució Catalana de Recerca i Estudis Avançats (ICREA), 08010 Barcelona, Spain
| |
Collapse
|
1594
|
Saito M, Okazaki I. Force-field parameters of the Psi and Phi around glycosidic bonds to oxygen and sulfur atoms. J Comput Chem 2009; 30:2656-65. [PMID: 19396813 DOI: 10.1002/jcc.21245] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
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.
Collapse
Affiliation(s)
- Minoru Saito
- Graduate School of Science and Technology, Hirosaki University, 3 Bunkyo-cho, Hirosaki, Aomori 036-8561, Japan.
| | | |
Collapse
|
1595
|
Hatcher E, Guvench O, MacKerell AD. CHARMM Additive All-Atom Force Field for Acyclic Polyalcohols, Acyclic Carbohydrates and Inositol. J Chem Theory Comput 2009; 5:1315-1327. [PMID: 20160980 PMCID: PMC2760998 DOI: 10.1021/ct9000608] [Citation(s) in RCA: 126] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Elizabeth Hatcher
- Department of Pharmaceutical Sciences, 20 Penn Street HSF II, University of Maryland, Baltimore, Maryland 21201
| | - Olgun Guvench
- Department of Pharmaceutical Sciences, 20 Penn Street HSF II, University of Maryland, Baltimore, Maryland 21201
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, 20 Penn Street HSF II, University of Maryland, Baltimore, Maryland 21201
| |
Collapse
|
1596
|
Computational approaches for the design of peptides with anti-breast cancer properties. Future Med Chem 2009; 1:201-12. [DOI: 10.4155/fmc.09.13] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
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.
Collapse
|
1597
|
Yongye AB, Gonzalez-Outeiriño J, Glushka J, Schultheis V, Woods RJ. The conformational properties of methyl alpha-(2,8)-di/trisialosides and their N-acyl analogues: implications for anti-Neisseria meningitidis B vaccine design. Biochemistry 2009; 47:12493-514. [PMID: 18954144 DOI: 10.1021/bi800431c] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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.
Collapse
Affiliation(s)
- Austin B Yongye
- Complex Carbohydrate Research Center, 315 Riverbend Road, Athens, Georgia 30602, USA
| | | | | | | | | |
Collapse
|
1598
|
Csonka GI, French AD, Johnson GP, Stortz CA. Evaluation of Density Functionals and Basis Sets for Carbohydrates. J Chem Theory Comput 2009; 5:679-92. [DOI: 10.1021/ct8004479] [Citation(s) in RCA: 165] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Gábor I. Csonka
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology, Szent Gellért tér 4, Budapest, H-1521 Hungary, Southern Regional Research Center, U.S. Department of Agriculture, 1100 Robert E. Lee Boulevard, New Orleans, Louisiana 70124, and Departamento de Química Orgánica-CIHIDECAR, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Alfred D. French
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology, Szent Gellért tér 4, Budapest, H-1521 Hungary, Southern Regional Research Center, U.S. Department of Agriculture, 1100 Robert E. Lee Boulevard, New Orleans, Louisiana 70124, and Departamento de Química Orgánica-CIHIDECAR, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Glenn P. Johnson
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology, Szent Gellért tér 4, Budapest, H-1521 Hungary, Southern Regional Research Center, U.S. Department of Agriculture, 1100 Robert E. Lee Boulevard, New Orleans, Louisiana 70124, and Departamento de Química Orgánica-CIHIDECAR, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Carlos A. Stortz
- Department of Inorganic and Analytical Chemistry, Budapest University of Technology, Szent Gellért tér 4, Budapest, H-1521 Hungary, Southern Regional Research Center, U.S. Department of Agriculture, 1100 Robert E. Lee Boulevard, New Orleans, Louisiana 70124, and Departamento de Química Orgánica-CIHIDECAR, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| |
Collapse
|
1599
|
A short guide for molecular dynamics simulations of RNA systems. Methods 2009; 47:187-97. [DOI: 10.1016/j.ymeth.2008.09.020] [Citation(s) in RCA: 64] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2008] [Accepted: 09/19/2008] [Indexed: 01/13/2023] Open
|
1600
|
Pol-Fachin L, Fernandes CL, Verli H. GROMOS96 43a1 performance on the characterization of glycoprotein conformational ensembles through molecular dynamics simulations. Carbohydr Res 2009; 344:491-500. [DOI: 10.1016/j.carres.2008.12.025] [Citation(s) in RCA: 60] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2008] [Revised: 12/03/2008] [Accepted: 12/26/2008] [Indexed: 10/21/2022]
|