1
|
Dhurua S, Jana M. Sulfation Effects of Chondroitin Sulfate to Bind a Chemokine in Aqueous Medium: Conformational Heterogeneity and Dynamics from Molecular Simulation. J Chem Inf Model 2023; 63:5660-5675. [PMID: 37611186 DOI: 10.1021/acs.jcim.3c00668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/25/2023]
Abstract
The sulfation patterns and degree of sulfation of chondroitin sulfate (CS), an important class of glycosaminoglycans (GAG), and their interactions with chemokines are accountable for various diseases. To realize the underlying mechanism of such complex biological phenomena at a molecular level and their application in rational drug design, a study on conformations and dynamics of CSs is necessary. To explore this, in this study, we performed a series of atomistic molecular dynamics (MD) simulations with different sulfated variants of octadecasaccharide CS, like CS-C, CS-E, and CS-T, in their free forms and when bound to the protein chemokine CXCL8 dimer in an aqueous medium. The calculated binding free energy of CSs with the CXCL8 dimer is favorable, and the degree of sulfation favors the complexation process further with prominent hydrophobic and hydrogen-bonded interactions. We find that the recognition is associated with the configurational entropy loss of the CS molecules as calculated from the Gaussian mixture approach, which supports that the degree of sulfation regulates the process. Cluster analysis through the k-means algorithm and end-to-end distance measurement revealed that although the free CS molecules adopted linear conformations, the nonlinear conformations during binding with protein were noted. Adaptation of nonlinear forms in the bound forms is noteworthy for the less-sulfated CS-C and CS-E. Apart from favorable 4C1 conformations, the occasional appearance of skew-boat forms from the free-energy map of ring pucker for the GlcUA unit was observed, which remains unaffected by the sulfation. We find that during recognition, the average relaxation time of intra-CS and inter-CS-CXCL8 hydrogen bonds (HBs) is about a magnitude lesser than that of CS-water HBs, most prominent on the involvement of higher sulfated CS-T analogues. The translational motion of surrounded water molecules in CSs exhibited sublinear diffusion, and the degree of sublinearity increases around the heavily sulfated molecules due to the hindrance created by them as well as the presence of the chemokine and exhibited markedly slow heterogeneous diffusion.
Collapse
Affiliation(s)
- Shakuntala Dhurua
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, Rourkela 769008, India
| | - Madhurima Jana
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, Rourkela 769008, India
| |
Collapse
|
2
|
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
|
3
|
Roy R, Jonniya NA, Kar P. Effect of Sulfation on the Conformational Dynamics of Dermatan Sulfate Glycosaminoglycan: A Gaussian Accelerated Molecular Dynamics Study. J Phys Chem B 2022; 126:3852-3866. [PMID: 35594147 DOI: 10.1021/acs.jpcb.2c01807] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Glycosaminoglycans (GAGs) are anionic biopolymers present on cell surfaces as a part of proteoglycans. The biological activities of GAGs depend on the sulfation pattern. In our study, we have considered three octadecasaccharide dermatan sulfate (DS) chains with increasing order of sulfation (dp6s, dp7s, and dp12s) to illuminate the role of sulfation on the GAG units and its chain conformation through 10 μs-long Gaussian accelerated molecular dynamics simulations. DS is composed of repeating disaccharide units of iduronic acid (IdoA) and N-acetylgalactosamine (N-GalNAc). Here, N-GalNAc is linked to IdoA via β(1-4), while IdoA is linked to N-GalNAc through α(1-3). With the increase in sulfation, the DS structure becomes more rigid and linear, as is evident from the distribution of root-mean-square deviations (RMSDs) and end-to-end distances. The tetrasaccharide linker region of the main chain shows a rigid conformation in terms of the glycosidic linkage. We have observed that upon sulfation (i.e., dp12s), the ring flip between two chair forms vanished for IdoA. The dynamic cross-correlation analysis reveals that the anticorrelation motions in dp12s are reduced significantly compared to dp6s or dp7s. An increase in sulfation generates relatively more stable hydrogen-bond networks, including water bridging with the neighboring monosaccharides. Despite the favorable linear structures of the GAG chains, our study also predicts few significant bendings related to the different puckering states, which may play a notable role in the function of the DS. The relation between the global conformation with the micro-level parameters such as puckering and water-mediated hydrogen bonds shapes the overall conformational space of GAGs. Overall, atomistic details of the DS chain provided in this study will help understand their functional and mechanical roles, besides developing new biomaterials.
Collapse
Affiliation(s)
- Rajarshi Roy
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore 453552, Madhya Pradesh, India
| | - Nisha Amarnath Jonniya
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore 453552, Madhya Pradesh, India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore 453552, Madhya Pradesh, India
| |
Collapse
|
4
|
Alibay I, Bryce RA. Ring Puckering Landscapes of Glycosaminoglycan-Related Monosaccharides from Molecular Dynamics Simulations. J Chem Inf Model 2019; 59:4729-4741. [PMID: 31609614 DOI: 10.1021/acs.jcim.9b00529] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The conformational flexibility of the glycosaminoglycans (GAGs) is known to be key in their binding and biological function, for example in regulating coagulation and cell growth. In this work, we employ enhanced sampling molecular dynamics simulations to probe the ring conformations of GAG-related monosaccharides, including a range of acetylated and sulfated GAG residues. We first perform unbiased MD simulations of glucose anomers and the epimers glucuronate and iduronate. These calculations indicate that in some cases, an excess of 15 μs is required for adequate sampling of ring pucker due to the high energy barriers between states. However, by applying our recently developed msesMD simulation method (multidimensional swarm-enhanced sampling molecular dynamics), we were able to quantitatively and rapidly reproduce these ring pucker landscapes. From msesMD simulations, the puckering free energy profiles were then compared for 15 further monosaccharides related to GAGs; this includes to our knowledge the first simulation study of sulfation effects on β-GalNAc ring puckering. For the force field employed, we find that in general the calculated pucker free energy profiles for sulfated sugars were similar to the corresponding unsulfated profiles. This accords with recent experimental studies suggesting that variation in ring pucker of sulfated GAG residues is primarily dictated by interactions with surrounding residues rather than by intrinsic conformational preference. As an exception to this, however, we predict that 4-O-sulfation of β-GalNAc leads to reduced ring rigidity, with a significant lowering in energy of the 1C4 ring conformation; this observation may have implications for understanding the structural basis of the biological function of β-GalNAc-containing glycosaminoglycans such as dermatan sulfate.
Collapse
Affiliation(s)
- Irfan Alibay
- Division of Pharmacy and Optometry, School of Health Sciences , University of Manchester , Oxford Road , Manchester M13 9PT , U.K.,Structural Bioinformatics and Computational Biochemistry Unit, Department of Biochemistry , University of Oxford , South Parks Road , Oxford OX1 3QU , U.K
| | - Richard A Bryce
- Division of Pharmacy and Optometry, School of Health Sciences , University of Manchester , Oxford Road , Manchester M13 9PT , U.K
| |
Collapse
|
5
|
Hosseini-Hashemi Z, Mirzaei M, Jafari A, Hosseinpour P, Yousefi M, Frontera A, Lari Dashtbayaz M, Shamsipur M, Ardalani M. Effects of N-oxidation on the molecular and crystal structures and properties of isocinchomeronic acid, its metal complexes and their supramolecular architectures: experimental, CSD survey, solution and theoretical approaches. RSC Adv 2019; 9:25382-25404. [PMID: 35530069 PMCID: PMC9070054 DOI: 10.1039/c9ra05143k] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2019] [Accepted: 08/01/2019] [Indexed: 12/17/2022] Open
Abstract
A series of new coordination complexes bearing H2pydco, an O-donor ligand, were synthesized and compared to corresponding complexes with H2pydc. The molecular and supramolecular structures were investigated using crystallography, solution phase and high level DFT studies.
Collapse
Affiliation(s)
| | - Masoud Mirzaei
- Department of Chemistry
- Faculty of Science
- Ferdowsi University of Mashhad
- Mashhad
- Iran
| | - Ameneh Jafari
- Department of Chemistry
- Faculty of Science
- Ferdowsi University of Mashhad
- Mashhad
- Iran
| | - Peyman Hosseinpour
- Department of Chemistry
- Faculty of Science
- Ferdowsi University of Mashhad
- Mashhad
- Iran
| | - Mohammad Yousefi
- Department of Chemistry
- Yadegar-e-Imam Khomeini (RAH) Shahr-e-Rey Branch
- Islamic Azad University
- Tehran
- Iran
| | - Antonio Frontera
- Department of Chemistry
- Universitat de les Illes Balears
- 07122 Palma de Mallorca (Baleares)
- Spain
| | - Mahmoud Lari Dashtbayaz
- Department of Economics and Administrative Sciences
- Ferdowsi University of Mashhad
- Mashhad
- Iran
| | | | - Mehdi Ardalani
- Department of Analytical Chemistry
- Razi University
- Kermanshah
- Iran
| |
Collapse
|
6
|
Lee J, Pothula KR, Kleinekathöfer U, Im W. Simulation Study of Occk5 Functional Properties in Pseudomonas aeruginosa Outer Membranes. J Phys Chem B 2018; 122:8185-8192. [DOI: 10.1021/acs.jpcb.8b07109] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Joonseong Lee
- Departments of Biological Sciences and Bioengineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| | - Karunakar R. Pothula
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Wonpil Im
- Departments of Biological Sciences and Bioengineering, Lehigh University, Bethlehem, Pennsylvania 18015, United States
| |
Collapse
|
7
|
Azurmendi HF, Battistel MD, Zarb J, Lichaa F, Negrete Virgen A, Shiloach J, Freedberg DI. The β-reducing end in α(2-8)-polysialic acid constitutes a unique structural motif. Glycobiology 2017; 27:900-911. [PMID: 28369425 PMCID: PMC6283323 DOI: 10.1093/glycob/cwx025] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2017] [Revised: 03/10/2017] [Accepted: 03/14/2017] [Indexed: 12/13/2022] Open
Abstract
Over the years, structural characterizations of α(2-8)-polysialic acid (polySia) in solution have produced inconclusive results. Efforts for obtaining detailed information in this important antigen have focused primarily on the α-linked residues and not on the distinctive characteristics of the terminal ones. The thermodynamically preferred anomeric configuration for the reducing end of sialic acids is β, which has the [I]CO2- group equatorial and the OH ([I]OH2) axial, while for all other residues the CO2- group is axial. We show that this purportedly minor difference has distinct consequences for the structure of α(2-8)-polySia near the reducing end, as the β configuration places the [I]OH2 in a favorable position for the formation of a hydrogen bond with the carboxylate group of the following residue ([II]CO2-). Molecular dynamics (MD) simulations predicted the hydrogen bond, which we subsequently directly detected by NMR. The combination of MD and residual dipolar couplings shows that the net result for the structure of Sia2-βOH is a stable conformation with well-defined hydration and charge patterns, and consistent with experimental NOE-based hydroxyl and aliphatic inter-proton distances. Moreover, we provide evidence that this distinct conformation is preserved on Sia oligosaccharides, thus constituting a motif that determines the structure and dynamics of α(2-8)-polySia for at least the first two residues of the polymer. We suggest the hypothesis that this structural motif sheds light on a longtime puzzling observation for the requirement of 10 residues of α(2-8)-polySia in order to bind effectively to specific antibodies, about four units more than for analogous cases.
Collapse
Affiliation(s)
- Hugo F Azurmendi
- Laboratory of Bacterial Polysaccharides, Center for Biologics Evaluation and Research, Food and Drug Administration, 10903 New Hampshire Avenue, Silver Spring, MD 20903, USA
| | - Marcos D Battistel
- Laboratory of Bacterial Polysaccharides, Center for Biologics Evaluation and Research, Food and Drug Administration, 10903 New Hampshire Avenue, Silver Spring, MD 20903, USA
| | - Jasmin Zarb
- Laboratory of Bacterial Polysaccharides, Center for Biologics Evaluation and Research, Food and Drug Administration, 10903 New Hampshire Avenue, Silver Spring, MD 20903, USA
| | - Flora Lichaa
- Laboratory of Bacterial Polysaccharides, Center for Biologics Evaluation and Research, Food and Drug Administration, 10903 New Hampshire Avenue, Silver Spring, MD 20903, USA
| | - Alejandro Negrete Virgen
- Laboratory of Bacterial Polysaccharides, Center for Biologics Evaluation and Research, Food and Drug Administration, 10903 New Hampshire Avenue, Silver Spring, MD 20903, USA
| | - Joseph Shiloach
- Biotechnology Unit, MSC 5522, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Darón I Freedberg
- Laboratory of Bacterial Polysaccharides, Center for Biologics Evaluation and Research, Food and Drug Administration, 10903 New Hampshire Avenue, Silver Spring, MD 20903, USA
| |
Collapse
|
8
|
Lee J, Patel DS, Kucharska I, Tamm LK, Im W. Refinement of OprH-LPS Interactions by Molecular Simulations. Biophys J 2017; 112:346-355. [PMID: 28122220 DOI: 10.1016/j.bpj.2016.12.006] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2016] [Revised: 11/27/2016] [Accepted: 12/05/2016] [Indexed: 12/16/2022] Open
Abstract
The outer membrane (OM) of Gram-negative bacteria is composed of lipopolysaccharide (LPS) in the outer leaflet and phospholipids in the inner leaflet. The outer membrane protein H (OprH) of Pseudomonas aeruginosa provides an increased stability to the OMs by directly interacting with LPS. Here we report the influence of various P. aeruginosa and, for comparison, Escherichia coli LPS environments on the physical properties of the OMs and OprH using all-atom molecular dynamics simulations. The simulations reveal that although the P. aeruginosa OMs are thinner hydrophobic bilayers than the E. coli OMs, which is expected from the difference in the acyl chain length of their lipid A, this effect is almost imperceptible around OprH due to a dynamically adjusted hydrophobic match between OprH and the OM. The structure and dynamics of the extracellular loops of OprH show distinct behaviors in different LPS environments. Including the O-antigen greatly reduces the flexibility of the OprH loops and increases the interactions between these loops and LPS. Furthermore, our study shows that the interactions between OprH and LPS mainly depend on the secondary structure of OprH and the chemical structure of LPS, resulting in distinctive patterns in different LPS environments.
Collapse
Affiliation(s)
- Joonseong Lee
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania; Bioengineering Program, Lehigh University, Bethlehem, Pennsylvania
| | - Dhilon S Patel
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania; Bioengineering Program, Lehigh University, Bethlehem, Pennsylvania
| | - Iga Kucharska
- Center for Membrane and Cell Physiology, University of Virginia School of Medicine, Charlottesville, Virginia; Department of Molecular Physiology and Biological Physics, University of Virginia School of Medicine, Charlottesville, Virginia
| | - Lukas K Tamm
- Center for Membrane and Cell Physiology, University of Virginia School of Medicine, Charlottesville, Virginia; Department of Molecular Physiology and Biological Physics, University of Virginia School of Medicine, Charlottesville, Virginia
| | - Wonpil Im
- Department of Biological Sciences, Lehigh University, Bethlehem, Pennsylvania; Bioengineering Program, Lehigh University, Bethlehem, Pennsylvania.
| |
Collapse
|
9
|
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]
|
10
|
Wu EL, Fleming PJ, Yeom MS, Widmalm G, Klauda JB, Fleming KG, Im W. E. coli outer membrane and interactions with OmpLA. Biophys J 2014; 106:2493-502. [PMID: 24896129 PMCID: PMC4052237 DOI: 10.1016/j.bpj.2014.04.024] [Citation(s) in RCA: 103] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2014] [Revised: 04/17/2014] [Accepted: 04/21/2014] [Indexed: 10/25/2022] Open
Abstract
The outer membrane of Gram-negative bacteria is a unique asymmetric lipid bilayer composed of phospholipids (PLs) in the inner leaflet and lipopolysaccharides (LPSs) in the outer leaflet. Its function as a selective barrier is crucial for the survival of bacteria in many distinct environments, and it also renders Gram-negative bacteria more resistant to antibiotics than their Gram-positive counterparts. Here, we report the structural properties of a model of the Escherichia coli outer membrane and its interaction with outer membrane phospholipase A (OmpLA) utilizing molecular dynamics simulations. Our results reveal that given the lipid composition used here, the hydrophobic thickness of the outer membrane is ∼3 Å thinner than the corresponding PL bilayer, mainly because of the thinner LPS leaflet. Further thinning in the vicinity of OmpLA is observed due to hydrophobic matching. The particular shape of the OmpLA barrel induces various interactions between LPS and PL leaflets, resulting in asymmetric thinning around the protein. The interaction between OmpLA extracellular loops and LPS (headgroups and core oligosaccharides) stabilizes the loop conformation with reduced dynamics, which leads to secondary structure variation and loop displacement compared to that in a DLPC bilayer. In addition, we demonstrate that the LPS/PL ratios in asymmetric bilayers can be reliably estimated by the per-lipid surface area of each lipid type, and there is no statistical difference in the overall membrane structure for the outer membranes with one more or less LPS in the outer leaflet, although individual lipid properties vary slightly.
Collapse
Affiliation(s)
- Emilia L Wu
- Department of Molecular Biosciences and Center for Bioinformatics, The University of Kansas, Lawrence, Kansas
| | - Patrick J Fleming
- T. C. Jenkins Department of Biophysics, John Hopkins University, Baltimore, Maryland
| | - Min Sun Yeom
- Korean Institute of Science and Technology Information, Daejeon, Korea
| | - Göran Widmalm
- Department of Organic Chemistry and Stockholm Center for Biomembrane Research, Arrhenius Laboratory, Stockholm University, Stockholm, Sweden
| | - Jeffery B Klauda
- Department of Chemical and Biomolecular Engineering, The University of Maryland, College Park, Maryland
| | - Karen G Fleming
- T. C. Jenkins Department of Biophysics, John Hopkins University, Baltimore, Maryland.
| | - Wonpil Im
- Department of Molecular Biosciences and Center for Bioinformatics, The University of Kansas, Lawrence, Kansas.
| |
Collapse
|
11
|
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: 17] [Impact Index Per Article: 1.7] [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
|
12
|
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
|
13
|
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
|
14
|
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
|
15
|
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
|
16
|
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
|
17
|
Sattelle BM, Bose-Basu B, Tessier M, Woods RJ, Serianni AS, Almond A. Dependence of pyranose ring puckering on anomeric configuration: methyl idopyranosides. J Phys Chem B 2012; 116:6380-6. [PMID: 22577942 DOI: 10.1021/jp303183y] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In the aldohexopyranose idose, the unique presence of three axial ring hydroxyl groups causes considerable conformational flexibility, rendering it challenging to study experimentally and an excellent model for rationalizing the relationship between puckering and anomeric configuration. Puckering in methyl α- and β-L-idopyranosides was predicted from kinetically rigorous 10 μs simulations using GLYCAM11 and three explicit water models (TIP3P, TIP4P, and TIP4P-EW). In each case, computed pyranose ring three-bond (vicinal) (1)H-(1)H spin couplings ((3)J(H,H)) trended with NMR measurements. These values, calculated puckering exchange rates and free energies, were independent of the water model. The α- and β-anomers were (1)C(4) chairs for 85 and >99% of their respective trajectories and underwent (1)C(4)→(4)C(1) exchange at rates of 20 μs(-1) and 1 μs(-1). Computed α-anomer (1)C(4)↔(4)C(1) puckering rates depended on the exocyclic C6 substituent, comparing hydroxymethyl with carboxyl from previous work. The slower kinetics and restricted pseudorotational profile of the β-anomer were caused by water occupying a cavity bounded by the anomeric 1-O-methyl and the C6 hydroxymethyl groups. This finding rationalizes the different methyl α- and β-L-idopyranoside (3)J(H,H) values. Identifying a relationship between idopyranose anomeric configuration, microsecond puckering, and water structure facilitates engineering of biologically and commercially important derivatives and underpins deciphering presently elusive structure-function relationships in the glycome.
Collapse
Affiliation(s)
- Benedict M Sattelle
- Faculty of Life Sciences, Manchester Interdisciplinary Biocentre, The University of Manchester, 131 Princess Street, Manchester M1 7DN, UK
| | | | | | | | | | | |
Collapse
|