1
|
Fontana C, de Meirelles JL, Verli H. Theoretical models of staurosporine and analogs uncover detailed structural information in biological solution. J Mol Graph Model 2024; 126:108653. [PMID: 37922640 DOI: 10.1016/j.jmgm.2023.108653] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2023] [Revised: 10/12/2023] [Accepted: 10/12/2023] [Indexed: 11/07/2023]
Abstract
Staurosporine and its analogs (STA-analogs) are indolocarbazoles (ICZs) compounds able to inhibit kinase proteins in a non-specific way, while present antimicrobial and cytostatic properties. The knowledge of molecular features associated to the complexation, including the ligand shape in solution and thermodynamics of complexation, is substantial to the development of new bioactive ICZs with improved therapeutic properties. In this context, the empirical approach of GROMOS force field is able to accurately reproduce condensed phase physicochemical properties of molecular systems after parameterization. Hence, through parameterization under GROMOS force field and molecular simulations, we assessed STA-analogs dynamics in aqueous solution, as well as its interaction with water to probe conformational and structural features involved in complexation to therapeutic targets. The coexistence of multiple conformers observed in simulations, and confirmed by metadynamics calculations, expanding the conformational space knowledge of these ligands with potential implications in understanding the ligand conformational selection during complexation. Also, changes in availability to H-bonding concerning the different substituents and water can reflect on effects at complexation free energy due to variation at the desolvation energetic costs. Based on these results, we expect the obtained structural data provide systemic framework for rational chemical modification of STA-analogs.
Collapse
Affiliation(s)
- Crisciele Fontana
- Universidade Federal do Rio Grande do Sul, Centro de Biotecnologia, Av. Bento Gonçalves, 9500 (Caixa Postal 15005), Porto Alegre, CEP 91501-970, RS, Brazil
| | - João Luiz de Meirelles
- Universidade Federal do Rio Grande do Sul, Centro de Biotecnologia, Av. Bento Gonçalves, 9500 (Caixa Postal 15005), Porto Alegre, CEP 91501-970, RS, Brazil
| | - Hugo Verli
- Universidade Federal do Rio Grande do Sul, Centro de Biotecnologia, Av. Bento Gonçalves, 9500 (Caixa Postal 15005), Porto Alegre, CEP 91501-970, RS, Brazil.
| |
Collapse
|
2
|
Multiscale Simulations Examining Glycan Shield Effects on Drug Binding to Influenza Neuraminidase. Biophys J 2020; 119:2275-2289. [PMID: 33130120 DOI: 10.1016/j.bpj.2020.10.024] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2020] [Revised: 10/08/2020] [Accepted: 10/21/2020] [Indexed: 12/18/2022] Open
Abstract
Influenza neuraminidase is an important drug target. Glycans are present on neuraminidase and are generally considered to inhibit antibody binding via their glycan shield. In this work, we studied the effect of glycans on the binding kinetics of antiviral drugs to the influenza neuraminidase. We created all-atom in silico systems of influenza neuraminidase with experimentally derived glycoprofiles consisting of four systems with different glycan conformations and one system without glycans. Using Brownian dynamics simulations, we observe a two- to eightfold decrease in the rate of ligand binding to the primary binding site of neuraminidase due to the presence of glycans. These glycans are capable of covering much of the surface area of neuraminidase, and the ligand binding inhibition is derived from glycans sterically occluding the primary binding site on a neighboring monomer. Our work also indicates that drugs preferentially bind to the primary binding site (i.e., the active site) over the secondary binding site, and we propose a binding mechanism illustrating this. These results help illuminate the complex interplay between glycans and ligand binding on the influenza membrane protein neuraminidase.
Collapse
|
3
|
Hamark C, Pendrill R, Landström J, Dotson Fagerström A, Sandgren M, Ståhlberg J, Widmalm G. Enantioselective Binding of Propranolol and Analogues Thereof to Cellobiohydrolase Cel7A. Chemistry 2018; 24:17975-17985. [PMID: 30255965 DOI: 10.1002/chem.201803104] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2018] [Indexed: 12/28/2022]
Abstract
At the catalytic site for the hydrolysis of cellulose the enzyme cellobiohydrolase Cel7A binds the enantiomers of the adrenergic beta-blocker propranolol with different selectivity. Methyl-to-hydroxymethyl group modifications of propranolol, which result in higher affinity and improved selectivity, were herein studied by 1 H,1 H and 1 H,13 C scalar spin-spin coupling constants as well as utilizing the nuclear Overhauser effect (NOE) in conjunction with molecular dynamics simulations of the ligands per se, which showed the presence of all-antiperiplanar conformations, except for the one containing a vicinal oxygen-oxygen arrangement governed by the gauche effect. For the ligand-protein complexes investigated by NMR spectroscopy using, inter alia, transferred NOESY and saturation-transfer difference (STD) NMR experiments the S-isomers were shown to bind with a higher affinity and a conformation similar to that preferred in solution, in contrast to the R-isomer. The fact that the S-form of the propranolol enantiomer is pre-arranged for binding to the protein is also observed for a crystal structure of dihydroxy-(S)-propranolol and Cel7A presented herein. Whereas the binding of propranolol is entropy driven, the complexation with the dihydroxy analogue is anticipated to be favored also by an enthalpic term, such as for its enantiomer, that is, dihydroxy-(R)-propranolol, because hydrogen-bond donation replaces the corresponding bonding from hydroxyl groups in glucosyl residues of the natural substrate. In addition to a favorable entropy component, albeit lesser in magnitude, this represents an effect of enthalpy-to-entropy compensation in ligand-protein interactions.
Collapse
Affiliation(s)
- Christoffer Hamark
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, 10691, Stockholm, Sweden
| | - Robert Pendrill
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, 10691, Stockholm, Sweden
| | - Jens Landström
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, 10691, Stockholm, Sweden
| | | | - Mats Sandgren
- Department of Molecular Sciences, Swedish University of Agricultural Sciences, P.O. Box 7015, 75007, Uppsala, Sweden
| | - Jerry Ståhlberg
- Department of Molecular Sciences, Swedish University of Agricultural Sciences, P.O. Box 7015, 75007, Uppsala, Sweden
| | - Göran Widmalm
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, 10691, Stockholm, Sweden
| |
Collapse
|
4
|
Abstract
Complex carbohydrates are ubiquitous in nature, and together with proteins and nucleic acids they comprise the building blocks of life. But unlike proteins and nucleic acids, carbohydrates form nonlinear polymers, and they are not characterized by robust secondary or tertiary structures but rather by distributions of well-defined conformational states. Their molecular flexibility means that oligosaccharides are often refractory to crystallization, and nuclear magnetic resonance (NMR) spectroscopy augmented by molecular dynamics (MD) simulation is the leading method for their characterization in solution. The biological importance of carbohydrate-protein interactions, in organismal development as well as in disease, places urgency on the creation of innovative experimental and theoretical methods that can predict the specificity of such interactions and quantify their strengths. Additionally, the emerging realization that protein glycosylation impacts protein function and immunogenicity places the ability to define the mechanisms by which glycosylation impacts these features at the forefront of carbohydrate modeling. This review will discuss the relevant theoretical approaches to studying the three-dimensional structures of this fascinating class of molecules and interactions, with reference to the relevant experimental data and techniques that are key for validation of the theoretical predictions.
Collapse
Affiliation(s)
- Robert J Woods
- Complex Carbohydrate Research Center and Department of Biochemistry and Molecular Biology , University of Georgia , 315 Riverbend Road , Athens , Georgia 30602 , United States
| |
Collapse
|
5
|
Hackbusch S, Watson A, Franz AH. Synthesis and conformational analysis of d-gluco-pyranosyl-(6,6′)-d-gluco-pyranuronate, a model compound for the inter-glycan 6,6′-ester linkage. Carbohydr Res 2018; 458-459:1-12. [DOI: 10.1016/j.carres.2018.01.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2017] [Revised: 01/07/2018] [Accepted: 01/16/2018] [Indexed: 11/26/2022]
|
6
|
Angles d’Ortoli T, Sjöberg NA, Vasiljeva P, Lindman J, Widmalm G, Bergenstråhle-Wohlert M, Wohlert J. Temperature Dependence of Hydroxymethyl Group Rotamer Populations in Cellooligomers. J Phys Chem B 2015; 119:9559-70. [DOI: 10.1021/acs.jpcb.5b02866] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Thibault Angles d’Ortoli
- Department
of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-106
91 Stockholm, Sweden
| | - Nils A. Sjöberg
- Wallenberg
Wood Science Center, and the Department of Fibre and Polymer Technology, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden
| | - Polina Vasiljeva
- Wallenberg
Wood Science Center, and the Department of Fibre and Polymer Technology, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden
| | - Jonas Lindman
- Wallenberg
Wood Science Center, and the Department of Fibre and Polymer Technology, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden
| | - Göran Widmalm
- Department
of Organic Chemistry, Arrhenius Laboratory, Stockholm University, SE-106
91 Stockholm, Sweden
| | - Malin Bergenstråhle-Wohlert
- Wallenberg
Wood Science Center, and the Department of Fibre and Polymer Technology, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden
| | - Jakob Wohlert
- Wallenberg
Wood Science Center, and the Department of Fibre and Polymer Technology, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden
| |
Collapse
|
7
|
Kang Y, Barbirz S, Lipowsky R, Santer M. Conformational Diversity of O-Antigen Polysaccharides of the Gram-Negative Bacterium Shigella flexneri Serotype Y. J Phys Chem B 2014; 118:2523-34. [DOI: 10.1021/jp4111713] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Yu Kang
- Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany
| | - Stefanie Barbirz
- Physikalische
Biochemie, Universität Potsdam, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany
| | - Reinhard Lipowsky
- Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany
| | - Mark Santer
- Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany
| |
Collapse
|
8
|
Pendrill R, Säwén E, Widmalm G. Conformation and dynamics at a flexible glycosidic linkage revealed by NMR spectroscopy and molecular dynamics simulations: analysis of β-L-Fucp-(1→6)-α-D-Glcp-OMe in water solution. J Phys Chem B 2013; 117:14709-22. [PMID: 24175957 DOI: 10.1021/jp409985h] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
The intrinsic flexibility of carbohydrates facilitates different 3D structures in response to altered environments. At glycosidic (1→6)-linkages, three torsion angles are variable, and herein the conformation and dynamics of β-L-Fucp-(1→6)-α-D-Glcp-OMe are investigated using a combination of NMR spectroscopy and molecular dynamics (MD) simulations. The disaccharide shows evidence of conformational averaging for the ψ and ω torsion angles, best explained by a four-state conformational distribution. Notably, there is a significant population of conformations having ψ = 85° (clinal) in addition to those having ψ = 180° (antiperiplanar). Moderate differences in (13)C R1 relaxation rates are found to be best explained by axially symmetric tumbling in combination with minor differences in librational motion for the two residues, whereas the isomerization motions are occurring too slowly to be contributing significantly to the observed relaxation rates. The MD simulation was found to give a reasonably good agreement with experiment, especially with respect to diffusive properties, among which the rotational anisotropy, D∥/D⊥, is found to be 2.35. The force field employed showed too narrow ω torsion angles in the gauche-trans and gauche-gauche states as well as overestimating the population of the gauche-trans conformer. This information can subsequently be used in directing parameter developments and emphasizes the need for refinement of force fields for (1→6)-linked carbohydrates.
Collapse
Affiliation(s)
- Robert Pendrill
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University , S-106 91 Stockholm, Sweden
| | | | | |
Collapse
|
9
|
Sattelle BM, Almond A. Shaping up for structural glycomics: a predictive protocol for oligosaccharide conformational analysis applied to N-linked glycans. Carbohydr Res 2013; 383:34-42. [PMID: 24252626 PMCID: PMC3909462 DOI: 10.1016/j.carres.2013.10.011] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2013] [Revised: 10/17/2013] [Accepted: 10/19/2013] [Indexed: 11/28/2022]
Abstract
Aqueous 10 μs simulations of N-linked mannosyl cores and sialyl Lewis (sLe) antennae are validated. Sequence dependent glycosidic linkage and pyranose ring μs motions are implicated in bioactivity. Stacked pyranoses in sLea and sLex are predicted to be atypically rigid on μs timescales. In a 25 μs simulation of sLex, all known conformers were sampled within the initial 10 μs of dynamics. Unbiased 10 μs simulations are proposed as a route to systematic and accurate glycomic 3D-analysis.
The human glycome comprises a vast untapped repository of 3D-structural information that holds the key to glycan recognition and a new era of rationally designed mimetic chemical probes, drugs, and biomaterials. Toward routine prediction of oligosaccharide conformational populations and exchange rates at thermodynamic equilibrium, we apply hardware-accelerated aqueous molecular dynamics to model μs motions in N-glycans that underpin inflammation and immunity. In 10 μs simulations, conformational equilibria of mannosyl cores, sialyl Lewis (sLe) antennae, and constituent sub-sequences agreed with prior refinements (X-ray and NMR). Glycosidic linkage and pyranose ring flexing were affected by branching, linkage position, and secondary structure, implicating sequence dependent motions in glycomic functional diversity. Linkage and ring conformational transitions that have eluded precise quantification by experiment and conventional (ns) simulations were predicted to occur on μs timescales. All rings populated non-chair shapes and the stacked galactose and fucose pyranoses of sLea and sLex were rigidified, suggesting an exploitable 3D-signature of cell adhesion protein binding. Analyses of sLex dynamics over 25 μs revealed that only 10 μs were sufficient to explore all aqueous conformers. This simulation protocol, which yields conformational ensembles that are independent of initial 3D-structure, is proposed as a route to understanding oligosaccharide recognition and structure–activity relationships, toward development of carbohydrate-based novel chemical entities.
Collapse
Affiliation(s)
- Benedict M Sattelle
- Faculty of Life Sciences, The University of Manchester, Manchester Institute of Biotechnology, 131 Princess Street, Manchester M1 7DN, UK
| | - Andrew Almond
- Faculty of Life Sciences, The University of Manchester, Manchester Institute of Biotechnology, 131 Princess Street, Manchester M1 7DN, UK.
| |
Collapse
|
10
|
Widmalm G. A perspective on the primary and three-dimensional structures of carbohydrates. Carbohydr Res 2013; 378:123-32. [PMID: 23522728 DOI: 10.1016/j.carres.2013.02.005] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2013] [Accepted: 02/13/2013] [Indexed: 10/27/2022]
Abstract
Carbohydrates, in more biologically oriented areas referred to as glycans, constitute one of the four groups of biomolecules. The glycans, often present as glycoproteins or glycolipids, form highly complex structures. In mammals ten monosaccharides are utilized in building glycoconjugates in the form of oligo- (up to about a dozen monomers) and polysaccharides. Subsequent modifications and additions create a large number of different compounds. In bacteria, more than a hundred monosaccharides have been reported to be constituents of lipopolysaccharides, capsular polysaccharides, and exopolysaccharides. Thus, the number of polysaccharide structures possible to create is huge. NMR spectroscopy plays an essential part in elucidating the primary structure, that is, monosaccharide identity and ring size, anomeric configuration, linkage position, and sequence, of the sugar residues. The structural studies may also employ computational approaches for NMR chemical shift predictions (CASPER program). Once the components and sequence of sugar residues have been unraveled, the three-dimensional arrangement of the sugar residues relative to each other (conformation), their flexibility (transitions between and populations of conformational states), together with the dynamics (timescales) should be addressed. To shed light on these aspects we have utilized a combination of experimental liquid state NMR techniques together with molecular dynamics simulations. For the latter a molecular mechanics force field such as our CHARMM-based PARM22/SU01 has been used. The experimental NMR parameters acquired are typically (1)H,(1)H cross-relaxation rates (related to NOEs), (3)JCH and (3)JCCtrans-glycosidic coupling constants and (1)H,(13)C- and (1)H,(1)H-residual dipolar couplings. At a glycosidic linkage two torsion angles ϕ and ψ are defined and for 6-substituted residues also the ω torsion angle is required. Major conformers can be identified for which highly populated states are present. Thus, in many cases a well-defined albeit not rigid structure can be identified. However, on longer timescales, oligosaccharides must be considered as highly flexible molecules since also anti-conformations have been shown to exist with H-C-O-C torsion angles of ∼180°, compared to syn-conformations in which the protons at the carbon atoms forming the glycosidic linkage are in close proximity. The accessible conformational space governs possible interactions with proteins and both minor changes and significant alterations occur for the oligosaccharides in these interaction processes. Transferred NOE NMR experiments give information on the conformation of the glycan ligand when bound to the proteins whereas saturation transfer difference NMR experiments report on the carbohydrate part in contact with the protein. It is anticipated that the subtle differences in conformational preferences for glycan structures facilitate a means to regulate biochemical processes in different environments. Further developments in the analysis of glycan structure and in particular its role in interactions with other molecules, will lead to clarifications of the importance of structure in biochemical regulation processes essential to health and disease.
Collapse
Affiliation(s)
- Göran Widmalm
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, S-106 91 Stockholm, Sweden.
| |
Collapse
|
11
|
Toukach FV, Ananikov VP. Recent advances in computational predictions of NMR parameters for the structure elucidation of carbohydrates: methods and limitations. Chem Soc Rev 2013; 42:8376-415. [DOI: 10.1039/c3cs60073d] [Citation(s) in RCA: 91] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
|
12
|
Wehle M, Vilotijevic I, Lipowsky R, Seeberger PH, Varon Silva D, Santer M. Mechanical Compressibility of the Glycosylphosphatidylinositol (GPI) Anchor Backbone Governed by Independent Glycosidic Linkages. J Am Chem Soc 2012; 134:18964-72. [DOI: 10.1021/ja302803r] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Marko Wehle
- Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany
| | - Ivan Vilotijevic
- Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany
| | - Reinhard Lipowsky
- Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany
| | - Peter H. Seeberger
- Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany
- Department of Chemistry and
Biochemistry, Free University of Berlin, Arnimallee 22, 14195 Berlin, Germany
| | | | - Mark Santer
- Max Planck Institute of Colloids and Interfaces, 14424 Potsdam, Germany
| |
Collapse
|
13
|
Bestaoui-Berrekhchi-Berrahma N, Derreumaux P, Sekkal-Rahal M, Springborg M, Sayede A, Yousfi N, Kadoun AED. Density functional conformational study of 2-O-sulfated 3,6 anhydro-α-D-galactose and of neo-κ- and ι-carrabiose molecules in gas phase and water. J Mol Model 2012; 19:893-904. [PMID: 23086461 DOI: 10.1007/s00894-012-1621-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2012] [Accepted: 10/01/2012] [Indexed: 10/27/2022]
Abstract
We examined the conformational preferences of the 2-O-sulfated-3,6-α-D-anhydrogalactose (compound I) and two 1,3 linked disaccharides constituting-κ or ι-carrageenans using density functional and ab initio methods in gas phase and aqueous solution. Systematic modifications of two torsion angles leading to 324 and 144 starting geometries for the compound I and each disaccharide were used to generate adiabatic maps using B3LYP/6-31G(d). The lower energy conformers were then fully optimized using B3LYP, B3PW91 and MP2 with several basis sets. Overall, we discuss the impact of full relaxation on the energy and structure of the dominant conformations, present the performance comparison with previous molecular mechanics calculations if available, and determine whether our results are impacted, when polarization and diffuse functions are added to the 6-31G(d) basis set, or when the MP2 level of theory is used.
Collapse
|
14
|
Taha HA, Richards MR, Lowary TL. Conformational Analysis of Furanoside-Containing Mono- and Oligosaccharides. Chem Rev 2012; 113:1851-76. [DOI: 10.1021/cr300249c] [Citation(s) in RCA: 99] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Hashem A. Taha
- Alberta Glycomics Centre and Department of Chemistry, Gunning−Lemieux Chemistry Centre, University of Alberta, Edmonton, AB, Canada T6G 2G2
| | - Michele R. Richards
- Alberta Glycomics Centre and Department of Chemistry, Gunning−Lemieux Chemistry Centre, University of Alberta, Edmonton, AB, Canada T6G 2G2
| | - Todd L. Lowary
- Alberta Glycomics Centre and Department of Chemistry, Gunning−Lemieux Chemistry Centre, University of Alberta, Edmonton, AB, Canada T6G 2G2
| |
Collapse
|
15
|
Re S, Nishima W, Miyashita N, Sugita Y. Conformational flexibility of N-glycans in solution studied by REMD simulations. Biophys Rev 2012; 4:179-187. [PMID: 28510079 DOI: 10.1007/s12551-012-0090-y] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2012] [Accepted: 06/21/2012] [Indexed: 01/09/2023] Open
Abstract
Protein-glycan recognition regulates a wide range of biological and pathogenic processes. Conformational diversity of glycans in solution is apparently incompatible with specific binding to their receptor proteins. One possibility is that among the different conformational states of a glycan, only one conformer is utilized for specific binding to a protein. However, the labile nature of glycans makes characterizing their conformational states a challenging issue. All-atom molecular dynamics (MD) simulations provide the atomic details of glycan structures in solution, but fairly extensive sampling is required for simulating the transitions between rotameric states. This difficulty limits application of conventional MD simulations to small fragments like di- and tri-saccharides. Replica-exchange molecular dynamics (REMD) simulation, with extensive sampling of structures in solution, provides a valuable way to identify a family of glycan conformers. This article reviews recent REMD simulations of glycans carried out by us or other research groups and provides new insights into the conformational equilibria of N-glycans and their alteration by chemical modification. We also emphasize the importance of statistical averaging over the multiple conformers of glycans for comparing simulation results with experimental observables. The results support the concept of "conformer selection" in protein-glycan recognition.
Collapse
Affiliation(s)
- Suyong Re
- RIKEN Advanced Science Institute, 2-1 Hirosawa, Wako, Saitama, 351-0198, Japan
| | - Wataru Nishima
- RIKEN Advanced Science Institute, 2-1 Hirosawa, Wako, Saitama, 351-0198, Japan
| | - Naoyuki Miyashita
- RIKEN Quantitative Biology Center, IMDA 6F, 1-6-5 Minatojimaminamimachi, Chuo-ku, Kobe, Hyogo, 650-0047, Japan
| | - Yuji Sugita
- RIKEN Advanced Science Institute, 2-1 Hirosawa, Wako, Saitama, 351-0198, Japan. .,RIKEN Quantitative Biology Center, IMDA 6F, 1-6-5 Minatojimaminamimachi, Chuo-ku, Kobe, Hyogo, 650-0047, Japan. .,RIKEN Advanced Institute for Computational Science, 7-1-26 Minatojimaminamimachi, Chuo-ku, Kobe, Hyogo, 650-0047, Japan.
| |
Collapse
|
16
|
Yamamoto S, Zhang Y, Yamaguchi T, Kameda T, Kato K. Lanthanide-assisted NMR evaluation of a dynamic ensemble of oligosaccharide conformations. Chem Commun (Camb) 2012; 48:4752-4. [DOI: 10.1039/c2cc30353a] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
17
|
Borgogna M, Bellich B, Cesàro A. Marine polysaccharides in microencapsulation and application to aquaculture: "from sea to sea". Mar Drugs 2011; 9:2572-2604. [PMID: 22363241 PMCID: PMC3280570 DOI: 10.3390/md9122572] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2011] [Revised: 11/18/2011] [Accepted: 11/22/2011] [Indexed: 12/15/2022] Open
Abstract
This review's main objective is to discuss some physico-chemical features of polysaccharides as intrinsic determinants for the supramolecular structures that can efficiently provide encapsulation of drugs and other biological entities. Thus, the general characteristics of some basic polysaccharides are outlined in terms of their conformational, dynamic and thermodynamic properties. The analysis of some polysaccharide gelling properties is also provided, including the peculiarity of the charged polysaccharides. Then, the way the basic physical chemistry of polymer self-assembly is made in practice through the laboratory methods is highlighted. A description of the several literature procedures used to influence molecular interactions into the macroscopic goal of the encapsulation is given with an attempt at classification. Finally, a practical case study of specific interest, the use of marine polysaccharide matrices for encapsulation of vaccines in aquaculture, is reported.
Collapse
Affiliation(s)
| | | | - Attilio Cesàro
- Department of Life Sciences, University of Trieste, Via L. Giorgieri, 1-I-34127 Trieste, Italy; (M.B.); (B.B.)
| |
Collapse
|
18
|
Theillet FX, Frank M, Vulliez-Le Normand B, Simenel C, Hoos S, Chaffotte A, Bélot F, Guerreiro C, Nato F, Phalipon A, Mulard LA, Delepierre M. Dynamic aspects of antibody:oligosaccharide complexes characterized by molecular dynamics simulations and saturation transfer difference nuclear magnetic resonance. Glycobiology 2011; 21:1570-9. [PMID: 21610193 DOI: 10.1093/glycob/cwr059] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Carbohydrates are likely to maintain significant conformational flexibility in antibody (Ab):carbohydrate complexes. As demonstrated herein for the protective monoclonal Ab (mAb) F22-4 recognizing the Shigella flexneri 2a O-antigen (O-Ag) and numerous synthetic oligosaccharide fragments thereof, the combination of molecular dynamics simulations and nuclear magnetic resonance saturation transfer difference experiments, supported by physicochemical analysis, allows us to determine the binding epitope and its various contributions to affinity without using any modified oligosaccharides. Moreover, the methods used provide insights into ligand flexibility in the complex, thus enabling a better understanding of the Ab affinities observed for a representative set of synthetic O-Ag fragments. Additionally, these complementary pieces of information give evidence to the ability of the studied mAb to recognize internal as well as terminal epitopes of its cognate polysaccharide antigen. Hence, we show that an appropriate combination of computational and experimental methods provides a basis to explore carbohydrate functional mimicry and receptor binding. The strategy may facilitate the design of either ligands or carbohydrate recognition domains, according to needed improvements of the natural carbohydrate:receptor properties.
Collapse
|
19
|
Roldós V, Cañada FJ, Jiménez-Barbero J. Carbohydrate-Protein Interactions: A 3D View by NMR. Chembiochem 2011; 12:990-1005. [DOI: 10.1002/cbic.201000705] [Citation(s) in RCA: 72] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2010] [Indexed: 12/29/2022]
|
20
|
Sattelle BM, Hansen SU, Gardiner J, Almond A. Free energy landscapes of iduronic acid and related monosaccharides. J Am Chem Soc 2010; 132:13132-4. [PMID: 20809637 DOI: 10.1021/ja1054143] [Citation(s) in RCA: 79] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The pyranose ring of L-iduronic acid (IdoA), a major constituent of the anticoagulant heparin, is an equilibrium of multiple ring puckers that have evaded quantification by experiment or computation. In order to resolve this enigma, we have calculated the free energy landscape of IdoA and two related monosaccharides from extensive microsecond simulations. After establishing that the simulated puckers had reached equilibrium, hypotheses were confirmed that (a) IdoA (1)C(4)- and (4)C(1)-chair conformations exchange on the microsecond time scale, (b) C5 epimerization leads to a (4)C(1)-chair, and (c) IdoA 2-O-sulfation (IdoA2S) stabilizes the (1)C(4) conformer. The IdoA and IdoA2S (1)C(4) conformers were isoenergetic and computed to be 0.9 and 2.6 kcal mol(-1) lower in free energy than their respective (4)C(1)-chair conformations. The simulations also predicted that the IdoA (2)S(O)-skew-boat was less populated than previously thought. Novel chemical synthesis and ultra-high-field NMR supported these observations, but slight discrepancies in observed and predicted NMR vicinal couplings implied that the simulation overestimated the population of the IdoA (4)C(1)-chair with respect to (1)C(4)-chair due to small force field inaccuracies that only manifest in long simulations. These free-energy calculations drive improvements in computational methods and provide a novel route to carbohydrate mimetic biomaterials and pharmaceuticals.
Collapse
Affiliation(s)
- Benedict M Sattelle
- Manchester Interdisciplinary Biocentre, University of Manchester, 131 Princess Street, Manchester M1 7DN, UK
| | | | | | | |
Collapse
|
21
|
Pol-Fachin L, Serrato RV, Verli H. Solution conformation and dynamics of exopolysaccharides from Burkholderia species. Carbohydr Res 2010; 345:1922-31. [DOI: 10.1016/j.carres.2010.06.010] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2010] [Revised: 06/09/2010] [Accepted: 06/21/2010] [Indexed: 11/17/2022]
|
22
|
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]
|