1
|
Lutsyk V, Plazinski W. Exploring Ring Conformation in Uronate Monosaccharides: Insights from Ab Initio Calculations and Classical Molecular Dynamics Simulations. J Phys Chem B 2024; 128:472-491. [PMID: 38170925 DOI: 10.1021/acs.jpcb.3c06556] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
The study focuses on the conformational properties of biologically relevant monosaccharides belonging to the group of uronates: α-l-iduronate, O2-sulfated-α-l-iduronate, and O2-sulfated-α-l-guluronate, either unfunctionalized or O1-methylated. We applied the previously proposed two-step methodology, combining classical MD simulations and subsequent ab initio (QM) calculations, performed on a rationally subsampled set of molecular configurations. We found that, regardless of the number of molecular configurations considered, the level of theory, and the weighting scheme applied, none of the QM approaches is capable of predicting the correct conformational equilibrium of sulfated iduronates as long as the tight counterion binding is not considered. Multicenter, ring-shape-specific binding of either Na+ or Ca2+ ions drastically shifts the conformational equilibrium of the pyranose ring in sulfated iduronates toward 1C4 but does not significantly affect the conformation of non-sulfated compounds. A similar shift is observed upon the protonation of carboxyl groups in all iduronates. In addition, we report a set of average J-coupling constant values related to vicinal protons bound to the pyranose ring of iduronates and corresponding to each of the three main groups of ring conformers, i.e., 4C1, B/S (boat/skew boat), and 1C4. In combination with the conformational energies or with the experimental data, these values allowed the relative proportions of the ring conformers to be estimated and the Karplus-type equations linking the 3JHH-coupling constants to the torsion angles within the pyranose ring to be refined.
Collapse
Affiliation(s)
- Valery Lutsyk
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, 30-239 Krakow, Poland
| | - Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, 30-239 Krakow, Poland
- Department of Biopharmacy, Medical University of Lublin, Chodzki 4a, 20-093 Lublin, Poland
| |
Collapse
|
2
|
Dhurua S, Jana M. Understanding Conformational Properties and Role of Hydrogen Bonds in Glycosaminoglycans-Interleukin8 Complexes in Aqueous Medium by Molecular Dynamics Simulation. Chemphyschem 2023; 24:e202200440. [PMID: 36239153 DOI: 10.1002/cphc.202200440] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2022] [Revised: 10/13/2022] [Indexed: 11/11/2022]
Abstract
Atomistic molecular dynamics simulations were performed under ambient conditions to explore the conformational features and binding affinities of hexameric glycosaminoglycans (GAGs) with chemokine Interleukin8 (IL8) in an aqueous medium. We tried to understand the role of hydrogen bonds (HBs) involving conserved water in mediating the interactions. The Luzar-Chandler model was adopted to study the kinetics of HB breaking and formation concerning different water-mediated HBs. The conformational flexibilities of bound GAGs are due to the flexible glycosidic linkages than the occasional/rare ring pucker conformation. The free energy landscape constructed with ϕ, and ψ, depicted that different conformational minima associated with the glycosidic linkage flexibility of the GAGs in bound states are separated by energy barriers. The binding affinities of IL8 towards GAGs are favored through the electrostatic and non-polar solvation interactions. 4-different types of conserved water were explored in the solvent-mediated binding of GAGs with IL8. The average lifetime of the IL8-GAG direct HB pairs was ∼ten times less than the IL8-GAG-shared water HBs. This is due to the rapid establishment of HB breaking and reformation kinetics involving water of a shared layer. We find that despite the highly negatively charged surface of GAGs, the IL8 surface populated by non-cationic amino acids could serve as a promising binding site in addition to the cationic surface of the protein.
Collapse
Affiliation(s)
- Shakuntala Dhurua
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, 769008, Rourkela, India
| | - Madhurima Jana
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, 769008, Rourkela, India
| |
Collapse
|
3
|
Roy R, Poddar S, Kar P. Comparison of the conformational dynamics of an N-glycan in implicit and explicit solvents. Carbohydr Res 2022; 522:108700. [DOI: 10.1016/j.carres.2022.108700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2022] [Revised: 09/30/2022] [Accepted: 10/03/2022] [Indexed: 11/28/2022]
|
4
|
Kogut MM, Danielsson A, Ricard-Blum S, Samsonov SA. Impact of calcium ions on the structural and dynamic properties of heparin oligosaccharides by computational analysis. Comput Biol Chem 2022; 99:107727. [PMID: 35841830 DOI: 10.1016/j.compbiolchem.2022.107727] [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: 05/12/2022] [Revised: 06/24/2022] [Accepted: 07/07/2022] [Indexed: 11/13/2022]
Abstract
Heparin (HP) belongs to glycosaminoglycans (GAGs), anionic linear polysaccharides composed of repetitive disaccharide units. They are key players in many biological processes occurring in the extracellular matrix and at the cell surface. GAGs are challenging molecules for computational research due to their high chemical heterogeneity, flexibility, periodicity, pseudosymmetry, predominantly electrostatics-driven nature of interactions with their protein partners and potential multipose binding. The molecular mechanisms underlying GAG interactions mediated by divalent ions, which are important for GAG binding to several proteins, are not well understood. The goal of this study was to characterize the binding of Ca2+ to two HP oligosaccharides of different lengths (dp10 and dp18, dp: degree of polymerization) and their impact on HP conformational space and their dynamic behavior with the use of molecular dynamics (MD)-based approaches with two Ca2+ parameter sets. MD data suggested that the flexibility of the monosaccharides, the glycosidic linkages and ring puckering were not affected by the presence of Ca2+, in contrast to H-bond propensities and the calculated Rg for a fraction of the oligosaccharide populations in both dp10 and dp18. Moreover, the essential differences in the data obtained by using two Ca2+ parameter sets were reported.
Collapse
Affiliation(s)
- Małgorzata M Kogut
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Annemarie Danielsson
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, Gdańsk 80-308, Poland
| | - Sylvie Ricard-Blum
- Univ Lyon, University Claude Bernard Lyon 1, CNRS, Institute of Molecular and Supramolecular Chemistry and Biochemistry, UMR 5246, Villeurbanne CEDEX F-69622, France
| | - Sergey A Samsonov
- Faculty of Chemistry, University of Gdańsk, ul. Wita Stwosza 63, Gdańsk 80-308, Poland.
| |
Collapse
|
5
|
Homayonia S, Ling CC. Highly Efficient and Stereoselective Synthesis of 6,7‐Dideoxy‐β‐D‐ido‐octopyranuronates. European J Org Chem 2022. [DOI: 10.1002/ejoc.202200175] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Saba Homayonia
- University of Calgary Department of Chemistry 2500 University Drive NW T2N 1N4 Calgary CANADA
| | - Chang-Chun Ling
- University of Calgary Department of Chemistry 2500 University Drive NW T2N 1N4 Calgary CANADA
| |
Collapse
|
6
|
Janke JJ, Yu Y, Pomin VH, Zhao J, Wang C, Linhardt RJ, García AE. Characterization of Heparin's Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy. J Chem Theory Comput 2022; 18:1894-1904. [PMID: 35108013 PMCID: PMC9027489 DOI: 10.1021/acs.jctc.1c00760] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin's solution conformational dynamics and structure-function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule's underlying conformational ensemble, which then provides important information necessary for understanding structure-function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin's conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. The anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.
Collapse
Affiliation(s)
- J Joel Janke
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, New York 12180, United States
| | - Yanlei Yu
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, New York 12180, United States
| | - Vitor H Pomin
- Department of BioMolecular Sciences, University of Mississippi, Oxford, Mississippi 38677, United States
| | - Jing Zhao
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, New York 12180, United States
| | - Chunyu Wang
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, New York 12180, United States
| | - Robert J Linhardt
- Center for Biotechnology and Interdisciplinary Studies, Rensselaer Polytechnic Institute, Troy, New York 12180, United States
| | - Angel E García
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
7
|
Feng X, Li F, Ding M, Zhang R, Shi T, Lu Y, Jiang W. Molecular dynamic simulation: Study on the recognition mechanism of linear β-(1 → 3)-D-glucan by Dectin-1. Carbohydr Polym 2022; 286:119276. [DOI: 10.1016/j.carbpol.2022.119276] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2021] [Revised: 02/13/2022] [Accepted: 02/18/2022] [Indexed: 12/26/2022]
|
8
|
Guvench O, Martin D, Greene M. Pyranose Ring Puckering Thermodynamics for Glycan Monosaccharides Associated with Vertebrate Proteins. Int J Mol Sci 2021; 23:473. [PMID: 35008898 PMCID: PMC8745717 DOI: 10.3390/ijms23010473] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2021] [Revised: 12/28/2021] [Accepted: 12/28/2021] [Indexed: 12/22/2022] Open
Abstract
The conformational properties of carbohydrates can contribute to protein structure directly through covalent conjugation in the cases of glycoproteins and proteoglycans and indirectly in the case of transmembrane proteins embedded in glycolipid-containing bilayers. However, there continue to be significant challenges associated with experimental structural biology of such carbohydrate-containing systems. All-atom explicit-solvent molecular dynamics simulations provide a direct atomic resolution view of biomolecular dynamics and thermodynamics, but the accuracy of the results depends on the quality of the force field parametrization used in the simulations. A key determinant of the conformational properties of carbohydrates is ring puckering. Here, we applied extended system adaptive biasing force (eABF) all-atom explicit-solvent molecular dynamics simulations to characterize the ring puckering thermodynamics of the ten common pyranose monosaccharides found in vertebrate biology (as represented by the CHARMM carbohydrate force field). The results, along with those for idose, demonstrate that the CHARMM force field reliably models ring puckering across this diverse set of molecules, including accurately capturing the subtle balance between 4C1 and 1C4 chair conformations in the cases of iduronate and of idose. This suggests the broad applicability of the force field for accurate modeling of carbohydrate-containing vertebrate biomolecules such as glycoproteins, proteoglycans, and glycolipids.
Collapse
Affiliation(s)
- Olgun Guvench
- Department of Pharmaceutical Sciences and Administration, School of Pharmacy, University of New England, 716 Stevens Avenue, Portland, ME 04103, USA; (D.M.); (M.G.)
- Graduate School of Biomedical Science and Engineering, University of Maine, 5775 Stodder Hall, Orono, ME 04469, USA
| | - Devon Martin
- Department of Pharmaceutical Sciences and Administration, School of Pharmacy, University of New England, 716 Stevens Avenue, Portland, ME 04103, USA; (D.M.); (M.G.)
- Graduate School of Biomedical Science and Engineering, University of Maine, 5775 Stodder Hall, Orono, ME 04469, USA
| | - Megan Greene
- Department of Pharmaceutical Sciences and Administration, School of Pharmacy, University of New England, 716 Stevens Avenue, Portland, ME 04103, USA; (D.M.); (M.G.)
| |
Collapse
|
9
|
Feng X, Li F, Ding M, Zhang R, Shi T, Jiang W. Molecular dynamic simulation: Structural insights of multi-stranded curdlan in aqueous solution. Carbohydr Polym 2021; 261:117844. [PMID: 33766340 DOI: 10.1016/j.carbpol.2021.117844] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2020] [Revised: 02/01/2021] [Accepted: 02/18/2021] [Indexed: 12/28/2022]
Abstract
In this work, by using molecular dynamic simulation we provide microscale structure information which helps to reveal the molecular mechanisms concerning the multi-chain conformational behavior of short curdlan. Through simulations starting with different conformations of curldan dodecasaccharides, it is found that the right-handed triple helix is thermodynamically the most stable conformation in aqueous solutions, which is well maintained and stabilized by an inter-strand hydrogen bonding network of the C2 hydroxyls. Unlike any predicted forms, the inter-strand hydrogen bonds exhibit a left-handed double helix pattern with preferred global orientations. Temperature REMD results suggest that the formation of triple helix is temperature sensitive, but the already formed triple helix is not. Investigation of curdlan with numbers of repeating units from 3 to 12 captures a critical value of 6, which in a way elucidates the relationship between the formation of triple helix and the chain length.
Collapse
Affiliation(s)
- Xuan Feng
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China; School of Applied Chemistry and Engineering, University of Science and Technology of China, Hefei 230026, PR China
| | - Fan Li
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China
| | - Mingming Ding
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China
| | - Ran Zhang
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China.
| | - Tongfei Shi
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China; School of Applied Chemistry and Engineering, University of Science and Technology of China, Hefei 230026, PR China.
| | - Wei Jiang
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China; School of Applied Chemistry and Engineering, University of Science and Technology of China, Hefei 230026, PR China
| |
Collapse
|
10
|
Feng X, Li F, Ding M, Zhang R, Shi T. Molecular dynamic simulation: Conformational properties of single-stranded curdlan in aqueous solution. Carbohydr Polym 2020; 250:116906. [PMID: 33049882 DOI: 10.1016/j.carbpol.2020.116906] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2020] [Revised: 08/04/2020] [Accepted: 08/04/2020] [Indexed: 01/08/2023]
Abstract
Recently, molecular dynamic simulation technique has been proved to be a powerful tool providing structural insights for better understanding the functionality of carbohydrates. Here, by using molecular dynamic simulation method we investigated the detailed conformational properties of the single-stranded curdlan with 12 glucose units. The results showed that the right-handed 6/1 helix structure was thermodynamically the most stable conformation in solution. The formation of the helix conformation was determined by many factors such as the glycosidic linkage, explicit water solvation and hydrogen bonds. When temperature was increased, the representative helix conformation was found becoming unstable giving rise to metastable conformations because when water mobility was accelerated with temperature, the hydrogen bonding strength between the curdlan chain and water went down, breaking the continuity of the hydrogen bonding network of water and hydroxyl groups. When the number of repeating glucose units varied from 6 to 24, the major helix conformation remained, but the conformational properties of longer chains were more apparently aff ;ected by chain flexibility.
Collapse
Affiliation(s)
- Xuan Feng
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China; School of Applied Chemistry and Engineering, University of Science and Technology of China, Hefei 230026, PR China
| | - Fan Li
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China.
| | - Mingming Ding
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China
| | - Ran Zhang
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China.
| | - Tongfei Shi
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, PR China; School of Applied Chemistry and Engineering, University of Science and Technology of China, Hefei 230026, PR China
| |
Collapse
|
11
|
The conformation of the idopyranose ring revisited: How subtle O-substituent induced changes can be deduced from vicinal 1H-NMR coupling constants. Carbohydr Res 2020; 496:108052. [PMID: 32738719 DOI: 10.1016/j.carres.2020.108052] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2020] [Revised: 05/17/2020] [Accepted: 05/27/2020] [Indexed: 01/18/2023]
Abstract
The idopyranose ring plays a pivotal role in the conformational, dynamical, and intermolecular binding aspects of glycosaminoglycans like heparin and dermatan sulfate and it was early on assigned a role in the Sugar Code governing biological recognition processes. There is consensus that next to the two canonical 1C4 and 4C1 chair conformations, the conformational space accessible to the idopyranose ring entails a 2SO skew-boat conformation, but the equilibrium between these three ring puckers has evaded satisfactory quantification. In this study a meta-analysis of X-ray solid-state data and vicinal NMR coupling constants is presented, based on the Truncated Fourier Puckering (TFP) formalism and the generalized Karplus (CAGPLUS) equation. This approach yields a model-free, granular and consistent reckoning of 159 idopyranose solution puckering equilibria studied by NMR and allows us to reproduce the involved 636 NMR vicinal couplings with an overall residual RMS(Jobs-Jcalc) of 0.184 Hz. Our analyses show that for all ring systems examined, the idopyranosyl chair conformations take up the same ring pucker irrespective of the ring substituent pattern or a vast variety in experimental conditions. Instead, it is the (skew-)boat conformation that adapts to the substitution pattern of the idopyranose ring or a specific sulfation pattern of neighboring saccharides. All idopyranose rings are involved in conformational equilibria that subsume the aforementioned conformers which turn out to differ only a few kJ/mole in conformational energy. Thus, the plasticity and flexibility of idopyranose remains intact under practically all circumstances and, as the glycosidic linkages in heparin are considered to be relatively stiff, the iduronic moiety functions as the linchpin of heparin flexibility thereby being rather a "space(r)" than a "letter" in the alleged Sugar Code alphabet.
Collapse
|
12
|
Pandey P, Aytenfisu AH, MacKerell AD, Mallajosyula SS. Drude Polarizable Force Field Parametrization of Carboxylate and N-Acetyl Amine Carbohydrate Derivatives. J Chem Theory Comput 2019; 15:4982-5000. [PMID: 31411469 PMCID: PMC6852669 DOI: 10.1021/acs.jctc.9b00327] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
In this work, we report the development of Drude polarizable force field parameters for the carboxylate and N-acetyl amine derivatives, extending the functionality of the existing Drude polarizable carbohydrate force field. The force field parameters have been developed in a hierarchical manner, reproducing the quantum mechanical gas-phase properties of small model compounds representing the key functional group in the carbohydrate derivatives, including optimization of the electrostatic and bonded parameters. The optimized parameters were then used to generate the models for carboxylate and N-acetyl amine carbohydrate derivatives. The transferred parameters were further tested and optimized to reproduce crystal geometries and J-coupling data from nuclear magnetic resonance experiments. The parameter development resulted in the incorporation of d-glucuronate, l-iduronate, N-acetyl-d-glucosamine (GlcNAc), and N-acetyl-d-galactosamine (GalNAc) sugars into the Drude polarizable force field. The parameters developed in this study were then applied to study the conformational properties of glycosaminoglycan polymer hyaluronan, composed of d-glucuronate and N-acetyl-d-glucosamine, in aqueous solution. Upon comparing the results from the additive and polarizable simulations, it was found that the inclusion of polarization improved the description of the electrostatic interactions observed in hyaluronan, resulting in enhanced conformational flexibility. The developed Drude polarizable force field parameters in conjunction with the remainder of the Drude polarizable force field parameters can be used for future studies involving carbohydrates and their conjugates in complex, heterogeneous systems.
Collapse
Affiliation(s)
| | - Asaminew H Aytenfisu
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , 20 Penn Street , Baltimore , Maryland 21201 , United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , 20 Penn Street , Baltimore , Maryland 21201 , United States
| | | |
Collapse
|
13
|
Nagarajan B, Sankaranarayanan NV, Desai UR. Perspective on computational simulations of glycosaminoglycans. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2019; 9:e1388. [PMID: 31080520 PMCID: PMC6504973 DOI: 10.1002/wcms.1388] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/03/2018] [Accepted: 07/07/2018] [Indexed: 01/06/2023]
Abstract
Glycosaminoglycans (GAGs) represent a formidable frontier for chemists, biochemists, biologists, medicinal chemists and drug delivery specialists because of massive structural complexity. GAGs are arguably the most complex, natural linear biopolymers with theoretical diversity orders of magnitude higher than proteins and nucleic acids. Yet, this diversity remains generally untapped. Computational approaches offer major routes to understand GAG structure and dynamics so as to enable novel applications of these biopolymers. In fact, computational algorithms, softwares, online tools and techniques have reached a level of sophistication that help understand atomistic details of conformational variation and protein recognition of individual GAG sequences. This review describes current approaches and challenges in computational study of GAGs. It presents a history of major findings since the earliest mention of GAGs (the 1960s), the development of parameters and force fields specific for GAGs, and the application of these tools in understanding GAG structure-function relationship. This review also presents a section on how to perform simulation of GAGs, which is directed toward researchers interested in entering this promising field with potential to impact therapy.
Collapse
Affiliation(s)
- Balaji Nagarajan
- Institute for Structural Biology, Drug Discovery and Development, Virginia Commonwealth University, Richmond,
VA 23298, USA
- Department of Medicinal Chemistry, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Nehru Viji Sankaranarayanan
- Institute for Structural Biology, Drug Discovery and Development, Virginia Commonwealth University, Richmond,
VA 23298, USA
- Department of Medicinal Chemistry, Virginia Commonwealth University, Richmond, VA 23298, USA
| | - Umesh R. Desai
- Institute for Structural Biology, Drug Discovery and Development, Virginia Commonwealth University, Richmond,
VA 23298, USA
- Department of Medicinal Chemistry, Virginia Commonwealth University, Richmond, VA 23298, USA
| |
Collapse
|
14
|
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
|
15
|
Panczyk K, Gaweda K, Drach M, Plazinski W. Extension of the GROMOS 56a6 CARBO/CARBO_R Force Field for Charged, Protonated, and Esterified Uronates. J Phys Chem B 2018; 122:3696-3710. [PMID: 29558620 DOI: 10.1021/acs.jpcb.7b11548] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
An extension of the GROMOS 56a6CARBO/CARBO_R force field for hexopyranose-based carbohydrates is presented. The additional parameters describe the conformational properties of uronate residues. The three distinct chemical states of the carboxyl group are considered: deprotonated (negatively charged), protonated (neutral), and esterified (neutral). The parametrization procedure was based on quantum-chemical calculations, and the resulting parameters were tested in the context of (i) flexibility of the pyranose rings under different pH conditions, (ii) conformation of the glycosidic linkage of the (1 → 4)-type for uronates with different chemical states of carboxyl moieties, (iii) conformation of the exocyclic (i.e., carboxylate and lactol) moieties, and (iv) structure of the Ca2+-linked chain-chain complexes of uronates. The presently proposed parameters in combination with the 56a6CARBO/CARBO_R set can be used to describe the naturally occurring polyuronates, composed either of homogeneous (e.g., glucuronans) or heterogeneous (e.g., pectins, alginates) segments. The results of simulations relying on the new set of parameters indicate that the conformation of glycosidic linkage is nearly unaffected by the chemical state of the carboxyl group, in contrary to the ring conformational equilibria. The calculations for the poly(α-d-galacturonate)-Ca2+ and poly(α-l-guluronate)-Ca2+ complexes show that both parallel and anitiparallel arrangements of uronate chains are possible but differ in several structural aspects.
Collapse
Affiliation(s)
- Karina Panczyk
- Jerzy Haber Institute of Catalysis and Surface Chemistry , Polish Academy of Sciences , ul. Niezapominajek 8 , 30-239 Cracow , Poland
| | - Karolina Gaweda
- Jerzy Haber Institute of Catalysis and Surface Chemistry , Polish Academy of Sciences , ul. Niezapominajek 8 , 30-239 Cracow , Poland
| | - Mateusz Drach
- Department of Theoretical Chemistry, Faculty of Chemistry , M. Curie-Sklodowska University , pl. M. Curie-Sklodowskiej 3 , 20-031 Lublin , Poland
| | - Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry , Polish Academy of Sciences , ul. Niezapominajek 8 , 30-239 Cracow , Poland
| |
Collapse
|
16
|
Almond A. Multiscale modeling of glycosaminoglycan structure and dynamics: current methods and challenges. Curr Opin Struct Biol 2017; 50:58-64. [PMID: 29253714 DOI: 10.1016/j.sbi.2017.11.008] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2017] [Revised: 10/24/2017] [Accepted: 11/26/2017] [Indexed: 01/24/2023]
Abstract
Glycosaminoglycans are long unbranched and complex polysaccharides that are an essential component of mammalian extracellular matrices. Characterization of their molecular structure, dynamics and interactions are essential to understand important biological phenomena in health and disease, and will lead to novel therapeutics and medical devices. However, this has proven to be a challenge experimentally and theoretical techniques are needed to develop new hypotheses, and interpret experiments. This review aims to examine the current theoretical (rather than experimental) methods used by researchers to investigate glycosaminoglycan structure, dynamics and interactions, from the monosaccharide to the macromolecular scale. It will consider techniques such as quantum mechanics, molecular mechanics, molecular dynamics, coarse graining and docking.
Collapse
Affiliation(s)
- Andrew Almond
- School of Chemistry, The University of Manchester, Manchester Institute of Biotechnology, 131 Princess Street, Manchester M1 7DN, UK.
| |
Collapse
|
17
|
Komarova BS, Gerbst AG, Finogenova AM, Dmitrenok AS, Tsvetkov YE, Nifantiev NE. 1,3-syn-Diaxial Repulsion of Typical Protecting Groups Used in Carbohydrate Chemistry in 3-O-Substituted Derivatives of Isopropyl d-Idopyranosides. J Org Chem 2017; 82:8897-8908. [DOI: 10.1021/acs.joc.7b01167] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Bozhena S. Komarova
- Laboratory of Glycoconjugate
Chemistry, N.D. Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, Leninsky Prospect 47, 119991 Moscow, Russia
| | - Alexey G. Gerbst
- Laboratory of Glycoconjugate
Chemistry, N.D. Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, Leninsky Prospect 47, 119991 Moscow, Russia
| | - Anastasiia M. Finogenova
- Laboratory of Glycoconjugate
Chemistry, N.D. Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, Leninsky Prospect 47, 119991 Moscow, Russia
| | - Andrey S. Dmitrenok
- Laboratory of Glycoconjugate
Chemistry, N.D. Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, Leninsky Prospect 47, 119991 Moscow, Russia
| | - Yury E. Tsvetkov
- Laboratory of Glycoconjugate
Chemistry, N.D. Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, Leninsky Prospect 47, 119991 Moscow, Russia
| | - Nikolay E. Nifantiev
- Laboratory of Glycoconjugate
Chemistry, N.D. Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, Leninsky Prospect 47, 119991 Moscow, Russia
| |
Collapse
|
18
|
Matzner D, Schüller A, Seitz T, Wittmann V, Mayer G. Fluoro-Carba-Sugars are Glycomimetic Activators of the glmS Ribozyme. Chemistry 2017; 23:12604-12612. [PMID: 28661578 DOI: 10.1002/chem.201702371] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2017] [Indexed: 11/05/2022]
Abstract
The glmS ribozyme is a bacterial gene-regulating riboswitch that controls cell wall synthesis, depending on glucosamine-6-phosphate as a cofactor. Due to the presence of this ribozyme in several human pathogen bacteria (e.g., MRSA, VRSA), the glmS ribozyme represents an attractive target for the development of artificial cofactors. The substitution of the ring oxygen in carbohydrates by functionalized methylene groups leads to a new generation of glycomimetics that exploits distinct interaction possibilities with their target structure in biological systems. Herein, we describe the synthesis of mono-fluoro-modified carba variants of α-d-glucosamine and β-l-idosamine. (5aR)-Fluoro-carba-α-d-glucosamine-6-phosphate is a synthetic mimic of the natural ligand of the glmS ribozyme and is capable of effectively addressing its unique self-cleavage mechanism. However, in contrast to what was expected, the activity is significantly decreased compared to its non-fluorinated analog. By combining self-cleavage assays with the Bacillus subtilis and Staphylococcus aureus glmS ribozyme and molecular docking studies, we provide a structure-activity relationship for fluorinated carba-sugars.
Collapse
Affiliation(s)
- Daniel Matzner
- Life & Medical Sciences Institute, University of Bonn, Gerhard-Domagk-Str. 1, 53121, Bonn, Germany
| | - Anna Schüller
- Life & Medical Sciences Institute, University of Bonn, Gerhard-Domagk-Str. 1, 53121, Bonn, Germany
| | - Torben Seitz
- Department of Chemistry, University of Konstanz, Universitätsstrasse 10, 78457, Konstanz, Germany.,Current address: Cilag AG, Schaffhausen, Switzerland
| | - Valentin Wittmann
- Department of Chemistry, University of Konstanz, Universitätsstrasse 10, 78457, Konstanz, Germany
| | - Günter Mayer
- Life & Medical Sciences Institute, University of Bonn, Gerhard-Domagk-Str. 1, 53121, Bonn, Germany
| |
Collapse
|
19
|
Bose-Basu B, Zhang W, Kennedy JLW, Hadad MJ, Carmichael I, Serianni AS. 13C-Labeled Idohexopyranosyl Rings: Effects of Methyl Glycosidation and C6 Oxidation on Ring Conformational Equilibria. J Org Chem 2017; 82:1356-1370. [DOI: 10.1021/acs.joc.6b02399] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Bidisha Bose-Basu
- Department
of Chemistry and Physics, Fayetteville State University, Fayetteville, North Carolina 28301, United States
| | | | | | | | | | | |
Collapse
|
20
|
Gaweda K, Plazinski W. Pyranose ring conformations in mono- and oligosaccharides: a combined MD and DFT approach. Phys Chem Chem Phys 2017; 19:20760-20772. [DOI: 10.1039/c7cp02920a] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
A two-step computational protocol is proposed to efficiently study the conformational properties of hexopyranoses with a special emphasis on their ring-inversion-properties. By applying it, the errors resulting from overestimating the contribution of the hydrogen bond-rich, low-energy structures that are not abundant in aqueous solutions are avoided.
Collapse
Affiliation(s)
- Karolina Gaweda
- Jerzy Haber Institute of Catalysis and Surface Chemistry
- Polish Academy of Sciences
- 30-239 Cracow
- Poland
| | - Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry
- Polish Academy of Sciences
- 30-239 Cracow
- Poland
| |
Collapse
|
21
|
Singh A, Tessier MB, Pederson K, Wang X, Venot AP, Boons GJ, Prestegard JH, Woods RJ. Extension and validation of the GLYCAM force field parameters for modeling glycosaminoglycans. CAN J CHEM 2016; 94:927-935. [PMID: 28603292 PMCID: PMC5464424 DOI: 10.1139/cjc-2015-0606] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Glycosaminoglycans (GAGs) are an important class of carbohydrates that serve critical roles in blood clotting, tissue repair, cell migration and adhesion, and lubrication. The variable sulfation pattern and iduronate ring conformations in GAGs influence their polymeric structure and nature of interaction. This study characterizes several heparin-like GAG disaccharides and tetrasaccharides using NMR and molecular dynamics simulations to assist in the development of parameters for GAGs within the GLYCAM06 force field. The force field additions include parameters and charges for a transferable sulfate group for O- and N-sulfation, neutral (COOH) forms of iduronic and glucuronic acid, and Δ4,5-unsaturated uronate (ΔUA) residues. ΔUA residues frequently arise from the enzymatic digestion of heparin and heparin sulfate. Simulations of disaccharides containing ΔUA reveal that the presence of sulfation on this residue alters the relative populations of 1H2 and 2H1 ring conformations. Simulations of heparin tetrasaccharides containing N-sulfation in place of N-acetylation on glucosamine residues influence the ring conformations of adjacent iduronate residues.
Collapse
Affiliation(s)
- Arunima Singh
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, GA 30602, USA
| | - Matthew B Tessier
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, GA 30602, USA
| | - Kari Pederson
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, GA 30602, USA
| | - Xiaocong Wang
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, GA 30602, USA
| | - Andre P Venot
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, GA 30602, USA
| | - Geert-Jan Boons
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, GA 30602, USA
| | - James H Prestegard
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, GA 30602, USA
| | - Robert J Woods
- Complex Carbohydrate Research Center, University of Georgia, 315 Riverbend Road, Athens, GA 30602, USA
| |
Collapse
|
22
|
Hsieh PH, Thieker DF, Guerrini M, Woods RJ, Liu J. Uncovering the Relationship between Sulphation Patterns and Conformation of Iduronic Acid in Heparan Sulphate. Sci Rep 2016; 6:29602. [PMID: 27412370 PMCID: PMC4944151 DOI: 10.1038/srep29602] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2016] [Accepted: 06/21/2016] [Indexed: 12/19/2022] Open
Abstract
The L-iduronic acid (IdoA) residue is a critically important structural component in heparan sulphate polysaccharide for the biological functions. The pyranose ring of IdoA is present in (1)C4-chair, (2)SO-skew boat, and less frequently, in (4)C1-chair conformations. Here, we analyzed the conformation of IdoA residue in eight hexasaccharides by NMR. The data demonstrate a correlation between the conformation of IdoA and sulphations in the surrounding saccharide residues. For the 2-O-sulpho IdoA residue, a high degree of sulphation on neighboring residues drives ring dynamics towards the (2)SO-skew boat conformer. In contrast, the nonsulphated IdoA residue is pushed towards the (1)C4-chair conformer when the neighboring residues are highly sulphated. Our data suggest that the conformation of IdoA is regulated by the sulphation pattern of nearby saccharides that is genetically controlled by the heparan sulphate biosynthetic pathway.
Collapse
Affiliation(s)
- Po-Hung Hsieh
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| | - David F. Thieker
- Complex Carbohydrate Research Center, University of Georgia, Athens, Georgia 30602, USA
| | - Marco Guerrini
- Istituto di Ricerche Chimiche e Biochimiche ‘G. Ronzoni’, via G. Colombo 81, 20133 Milan, Italy
| | - Robert J. Woods
- Complex Carbohydrate Research Center, University of Georgia, Athens, Georgia 30602, USA
| | - Jian Liu
- Division of Chemical Biology and Medicinal Chemistry, Eshelman School of Pharmacy, University of North Carolina, Chapel Hill, North Carolina 27599, USA
| |
Collapse
|
23
|
Feng T, Zhu X, Campanella O. Molecular modeling tools to characterize the structure and complexation behavior of carbohydrates. Curr Opin Food Sci 2016. [DOI: 10.1016/j.cofs.2016.08.003] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
24
|
Ring inversion properties of 1→2, 1→3 and 1→6-linked hexopyranoses and their correlation with the conformation of glycosidic linkages. Carbohydr Res 2016; 423:43-8. [DOI: 10.1016/j.carres.2016.01.007] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2015] [Revised: 01/18/2016] [Accepted: 01/19/2016] [Indexed: 10/22/2022]
|
25
|
Plazinski W, Lonardi A, Hünenberger PH. Revision of the GROMOS 56A6(CARBO) force field: Improving the description of ring-conformational equilibria in hexopyranose-based carbohydrates chains. J Comput Chem 2015; 37:354-65. [PMID: 26525424 DOI: 10.1002/jcc.24229] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2015] [Revised: 10/03/2015] [Accepted: 10/06/2015] [Indexed: 01/19/2023]
Abstract
This article describes a revised version 56A6(CARBO_R) of the GROMOS 56A6(CARBO) force field for hexopyranose-based carbohydrates. The simulated properties of unfunctionalized hexopyranoses are unaltered with respect to 56A6CARBO . In the context of both O1 -alkylated hexopyranoses and oligosaccharides, the revision stabilizes the regular (4) C1 chair for α-anomers, with the opposite effect for β-anomers. As a result, spurious ring inversions observed in α(1→4)-linked chains when using the original 56A6(CARBO) force field are alleviated. The (4) C1 chair is now the most stable conformation for all d-hexopyranose residues, irrespective of the linkage type and anomery, and of the position of the residue along the chain. The methylation of a d-hexopyranose leads to a systematic shift in the ring-inversion free energy ((4) C1 to (1) C4 ) by 7-8 kJ mol(-1), positive for the α-anomers and negative for the β-anomers, which is qualitatively compatible with the expected enhancement of the anomeric effect upon methylation at O1. The ring-inversion free energies for residues within chains are typically smaller in magnitude compared to those of the monomers, and correlate rather poorly with the latter. This suggests that the crowding of ring substituents upon chain formation alters the ring flexibility in a nonsystematic fashion. In general, the description of carbohydrate chains afforded by 56A6(CARBO_R) suggests a significant extent of ring flexibility, i.e., small but often non-negligible equilibrium populations of inverted chairs, and challenges the "textbook" picture of conformationally locked carbohydrate rings.
Collapse
Affiliation(s)
- Wojciech Plazinski
- Laboratory of Physical Chemistry, ETH Hönggerberg, HCI, Zürich, CH-8093, Switzerland.,J. Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Ul. Niezapominajek 8, Cracow, 30-239, Poland
| | - Alice Lonardi
- Laboratory of Physical Chemistry, ETH Hönggerberg, HCI, Zürich, CH-8093, Switzerland
| | | |
Collapse
|
26
|
Unione L, Xu B, Díaz D, Martín-Santamaría S, Poveda A, Sardinha J, Rauter AP, Blériot Y, Zhang Y, Cañada FJ, Sollogoub M, Jiménez-Barbero J. Conformational Plasticity in Glycomimetics: Fluorocarbamethyl-L-idopyranosides Mimic the Intrinsic Dynamic Behaviour of Natural Idose Rings. Chemistry 2015; 21:10513-21. [PMID: 26096911 DOI: 10.1002/chem.201501249] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2015] [Revised: 05/08/2015] [Indexed: 11/06/2022]
Abstract
Sugar function, structure and dynamics are intricately correlated. Ring flexibility is intrinsically related to biological activity; actually plasticity in L-iduronic rings modulates their interactions with biological receptors. However, the access to the experimental values of the energy barriers and free-energy difference for conformer interconversion in water solution has been elusive. Here, a new generation of fluorine-containing glycomimetics is presented. We have applied a combination of organic synthesis, NMR spectroscopy and computational methods to investigate the conformational behaviour of idose- and glucose-like rings. We have used low-temperature NMR spectroscopic experiments to slow down the conformational exchange of the idose-like rings. Under these conditions, the exchange rate becomes slow in the (19) F NMR spectroscopic chemical shift timescale and allows shedding light on the thermodynamic and kinetic features of the equilibrium. Despite the minimal structural differences between these compounds, a remarkable difference in their dynamic behaviour indeed occurs. The importance of introducing fluorine atoms in these sugars mimics is also highlighted. Only the use of (19) F NMR spectroscopic experiments has permitted the unveiling of key features of the conformational equilibrium that would have otherwise remained unobserved.
Collapse
Affiliation(s)
- Luca Unione
- Chemical and Physical Biology Department, Centro de Investigaciones Biológicas, CIB-CSIC, Ramiro de Maeztu 9, 28040 Madrid (Spain).,Infectious Disease Programme, CIC bioGUNE, 48160 Derio, Bizkaia (Spain)
| | - Bixue Xu
- Sorbonne Universités, UPMC Univ Paris 06, Institut Universitaire de France, UMR CNRS 8232, IPCM, 4, place Jussieu, 75005 Paris (France).,Present address: The Key Laboratory of Chemistry for Natural Products of Guizhou Province and Chinese Academy of Sciences, 202 Shachong South Road, Guiyang, 550002 (P. R. China)
| | - Dolores Díaz
- Chemical and Physical Biology Department, Centro de Investigaciones Biológicas, CIB-CSIC, Ramiro de Maeztu 9, 28040 Madrid (Spain)
| | - Sonsoles Martín-Santamaría
- Chemical and Physical Biology Department, Centro de Investigaciones Biológicas, CIB-CSIC, Ramiro de Maeztu 9, 28040 Madrid (Spain)
| | - Ana Poveda
- Infectious Disease Programme, CIC bioGUNE, 48160 Derio, Bizkaia (Spain)
| | - João Sardinha
- Sorbonne Universités, UPMC Univ Paris 06, Institut Universitaire de France, UMR CNRS 8232, IPCM, 4, place Jussieu, 75005 Paris (France).,Centro de Química e Bioquímica, Faculdade de Ciências da Universidade de Lisboa, Edifício C8, 1749-016 Lisboa (Portugal)
| | - Amelia Pilar Rauter
- Centro de Química e Bioquímica, Faculdade de Ciências da Universidade de Lisboa, Edifício C8, 1749-016 Lisboa (Portugal)
| | - Yves Blériot
- Université de Poitiers, UMR CNRS 7285, IC2MP, Equipe Synthèse organique, Groupe Glycochimie, 4, avenue Michel Brunet, 86022 Poitiers Cedex (France)
| | - Yongmin Zhang
- Sorbonne Universités, UPMC Univ Paris 06, Institut Universitaire de France, UMR CNRS 8232, IPCM, 4, place Jussieu, 75005 Paris (France)
| | - F Javier Cañada
- Chemical and Physical Biology Department, Centro de Investigaciones Biológicas, CIB-CSIC, Ramiro de Maeztu 9, 28040 Madrid (Spain)
| | - Matthieu Sollogoub
- Sorbonne Universités, UPMC Univ Paris 06, Institut Universitaire de France, UMR CNRS 8232, IPCM, 4, place Jussieu, 75005 Paris (France).
| | - Jesus Jiménez-Barbero
- Chemical and Physical Biology Department, Centro de Investigaciones Biológicas, CIB-CSIC, Ramiro de Maeztu 9, 28040 Madrid (Spain). .,Infectious Disease Programme, CIC bioGUNE, 48160 Derio, Bizkaia (Spain). .,IKERBASQUE, Basque Foundation for Science, 48011 Bilbao (Spain).
| |
Collapse
|
27
|
Force fields and scoring functions for carbohydrate simulation. Carbohydr Res 2015; 401:73-81. [DOI: 10.1016/j.carres.2014.10.028] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2014] [Revised: 10/28/2014] [Accepted: 10/30/2014] [Indexed: 12/31/2022]
|
28
|
Makeneni S, Foley BL, Woods RJ. BFMP: a method for discretizing and visualizing pyranose conformations. J Chem Inf Model 2014; 54:2744-50. [PMID: 25289680 PMCID: PMC4210171 DOI: 10.1021/ci500325b] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report a new classification method for pyranose ring conformations called Best-fit, Four-Membered Plane (BFMP), which describes pyranose ring conformations based on reference planes defined by four atoms. The method is able to characterize all asymmetrical and symmetrical shapes of a pyran ring, is readily automated, easy to interpret, and maps trivially to IUPAC definitions. It also provides a qualitative measurement of the distortion of the ring. Example applications include the analysis of data from crystal structures and molecular dynamics simulations.
Collapse
Affiliation(s)
- Spandana Makeneni
- Complex Carbohydrate Research Center, University of Georgia , 315 Riverbend Road, Athens, Georgia 30602, United States
| | | | | |
Collapse
|
29
|
Sattelle BM, Almond A. Microsecond kinetics in model single- and double-stranded amylose polymers. Phys Chem Chem Phys 2014; 16:8119-26. [PMID: 24652085 PMCID: PMC4006424 DOI: 10.1039/c4cp00570h] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2014] [Accepted: 03/14/2014] [Indexed: 02/04/2023]
Abstract
Amylose, a component of starch with increasing biotechnological significance, is a linear glucose polysaccharide that self-organizes into single- and double-helical assemblies. Starch granule packing, gelation and inclusion-complex formation result from finely balanced macromolecular kinetics that have eluded precise experimental quantification. Here, graphics processing unit (GPU) accelerated multi-microsecond aqueous simulations are employed to explore conformational kinetics in model single- and double-stranded amylose. The all-atom dynamics concur with prior X-ray and NMR data while surprising and previously overlooked microsecond helix-coil, glycosidic linkage and pyranose ring exchange are hypothesized. In a dodecasaccharide, single-helical collapse was correlated with linkages and rings transitioning from their expected syn and (4)C1 chair conformers. The associated microsecond exchange rates were dependent on proximity to the termini and chain length (comparing hexa- and trisaccharides), while kinetic features of dodecasaccharide linkage and ring flexing are proposed to be a good model for polymers. Similar length double-helices were stable on microsecond timescales but the parallel configuration was sturdier than the antiparallel equivalent. In both, tertiary organization restricted local chain dynamics, implying that simulations of single amylose strands cannot be extrapolated to dimers. Unbiased multi-microsecond simulations of amylose are proposed as a valuable route to probing macromolecular kinetics in water, assessing the impact of chemical modifications on helical stability and accelerating the development of new biotechnologies.
Collapse
Affiliation(s)
- Benedict M. Sattelle
- Faculty of Life Sciences , Manchester Institute of Biotechnology , The University of Manchester , 131 Princess Street , Manchester , M1 7DN , UK . ; Tel: +44 (0)161 306 4199
| | - Andrew Almond
- Faculty of Life Sciences , Manchester Institute of Biotechnology , The University of Manchester , 131 Princess Street , Manchester , M1 7DN , UK . ; Tel: +44 (0)161 306 4199
| |
Collapse
|
30
|
Mayes HB, Broadbelt LJ, Beckham GT. How Sugars Pucker: Electronic Structure Calculations Map the Kinetic Landscape of Five Biologically Paramount Monosaccharides and Their Implications for Enzymatic Catalysis. J Am Chem Soc 2014; 136:1008-22. [DOI: 10.1021/ja410264d] [Citation(s) in RCA: 117] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Heather B. Mayes
- Department
of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, United States
- National
Bioenergy Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Linda J. Broadbelt
- Department
of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, United States
| | - Gregg T. Beckham
- National
Bioenergy Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| |
Collapse
|
31
|
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
|
32
|
Sattelle BM, Shakeri J, Almond A. Does Microsecond Sugar Ring Flexing Encode 3D-Shape and Bioactivity in the Heparanome? Biomacromolecules 2013; 14:1149-59. [DOI: 10.1021/bm400067g] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Benedict M. Sattelle
- Faculty of Life Sciences, University of Manchester, Manchester
Institute of Biotechnology,
131 Princess Street, Manchester, M1 7DN, United Kingdom
| | - Javad Shakeri
- Faculty of Life Sciences, University of Manchester, Manchester
Institute of Biotechnology,
131 Princess Street, Manchester, M1 7DN, United Kingdom
| | - Andrew Almond
- Faculty of Life Sciences, University of Manchester, Manchester
Institute of Biotechnology,
131 Princess Street, Manchester, M1 7DN, United Kingdom
| |
Collapse
|
33
|
Oborský P, Tvaroška I, Králová B, Spiwok V. Toward an accurate conformational modeling of iduronic acid. J Phys Chem B 2013; 117:1003-9. [PMID: 23286518 DOI: 10.1021/jp3100552] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Iduronic acid (IdoA), unlike most other monosaccharides, can adopt different ring conformations, depending on the context of the molecular structure. Accurate modeling of this building block is essential for understanding the role of glycosaminoglycans and other glycoconjugates. Here, we use metadynamics to predict equilibria of (1)C(4), (4)C(1) and (2)S(O) conformations of α-L-IdoA-OMe and α-L-IdoA2S-OMe. Different schemes of scaling of atoms separated by three bonds (1-4 interaction) were tested. It was found that scaling (reduction) of 1-4 electrostatic interactions significantly changes conformational preferences toward the (4)C(1) conformation. More interestingly, scaling of 1-4 van der Waals interaction favors skew-boat conformations. This shows that a minor modification of noncovalent 1-4 interactions parameters can provide a good agreement between populations of conformers of iduronic acid in water from simulations and experiments.
Collapse
Affiliation(s)
- Pavel Oborský
- Department of Biochemistry and Microbiology, Institute of Chemical Technology, Prague, Technická 3, Prague 6 166 28, Czech Republic
| | | | | | | |
Collapse
|
34
|
Kirschner KN, Lins RD, Maass A, Soares TA. A Glycam-Based Force Field for Simulations of Lipopolysaccharide Membranes: Parametrization and Validation. J Chem Theory Comput 2012; 8:4719-31. [DOI: 10.1021/ct300534j] [Citation(s) in RCA: 83] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Karl N. Kirschner
- Fraunhofer-Institute for Algorithms
and Scientific Computing (SCAI), Department of Bioinformatics, Schloss
Birlinghoven, 53754 Sankt Augustin, Germany
| | - Roberto D. Lins
- Department of Fundamental Chemistry,
Federal University of Pernambuco, Cidade Universitária, Recife,
PE 50740-560 Brazil
| | - Astrid Maass
- Fraunhofer-Institute for Algorithms
and Scientific Computing (SCAI), Department of Simulation Engineering,
Schloss Birlinghoven, 53754 Sankt Augustin, Germany
| | - Thereza A. Soares
- Department of Fundamental Chemistry,
Federal University of Pernambuco, Cidade Universitária, Recife,
PE 50740-560 Brazil
| |
Collapse
|