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
|
Plazinski W, Angles d'Ortoli T, Widmalm G. Conformational flexibility of the disaccharide β-L-Fuc p-(1→4)-α-D-Glc p-OMe as deduced from NMR spectroscopy experiments and computer simulations. Org Biomol Chem 2023; 21:6979-6994. [PMID: 37584331 DOI: 10.1039/d3ob01153d] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/17/2023]
Abstract
Carbohydrates in biological systems are referred to as glycans and modification of their structures is a hallmark indicator of disease. Analysis of the three-dimensional structure forms the basis for further insight into how they function and comparison of crystal structure with solution-state conformation(s) is particularly relevant, which has been performed for the disaccharide β-L-Fucp-(1→4)-α-D-Glcp-OMe. In water solution the conformational space at the glycosidic linkage between the two sugar residues is identified from molecular dynamics (MD) simulations as having a low-energy exo-syn conformation, deviating somewhat from the solid-state conformation, and two anti-conformational states, i.e., anti-ϕ and anti-ψ, indicating flexibility at the glycosidic linkage. NMR data were obtained from 1D 1H,1H-NOESY and STEP-NOESY experiments, measurement of transglycosidic 3JCH coupling constants and NMR spin-simulation. The free energy profile of the ω torsion angle computed from MD simulation was in excellent agreement with the rotamer distribution from NMR experiment being for gt:gg:tg 38 : 53 : 9, respectively, with a proposed inter-residue O5'⋯HO6 hydrogen bond being predominant in the gg rotamer. Quantum mechanics methodology was used to calculate transglycosidic NMR 3JCH coupling constants, averaged over a conformational ensemble of structures representing various rotamers of exocyclic groups, in good to excellent agreement with Karplus-type relationships previously developed. Furthermore, 1H and 13C NMR chemical shifts were calculated using the same methodology and were found to be in excellent agreement with experimental data.
Collapse
Affiliation(s)
- Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, 30-239 Krakow, Poland
- Department of Biopharmacy, Medical University of Lublin, 20-093 Lublin, Poland
| | - Thibault Angles d'Ortoli
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, S-106 91 Stockholm, Sweden.
| | - Göran Widmalm
- Department of Organic Chemistry, Arrhenius Laboratory, Stockholm University, S-106 91 Stockholm, Sweden.
| |
Collapse
|
3
|
Tsou PK, Huynh HT, Phan HT, Kuo JL. A self-adapting first-principles exploration on the dissociation mechanism in sodiated aldohexose pyranoses assisted with neural network potentials. Phys Chem Chem Phys 2023; 25:3332-3342. [PMID: 36633012 DOI: 10.1039/d2cp04421h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Understanding the mechanism of collision-induced dissociation (CID) in mono-saccharides with density functional theory (DFT) is challenging because of many possible reaction paths that originate from their high structural diversity. To search for the transition state (TS) from the huge number of conformers, we propose a three-step search scheme with the assistance of neural network potential (NNP). The search starts from a cross-checking of sugars, to a global search of all possible channels, and in the end, an exhaustive exploration around the low-lying channels. The cross-checking step quickly adapts the NNP from the studied molecules to the target ones. The other two steps utilize the adapted NNP to find the available pathways via random sampling of the structures. The study of the CID reactions in all eight types of aldohexose pyranoses was applied using the search scheme. The DFT calculations on AH-0 (Glc, Gal, and Man) in the previous study were utilized to construct an NNP and provide the TS structure database for searching AH-1 (All, Alt, Gul, Ido, and Tal). In total, we identified around 5200 TSs in AH-0 and AH-1, and the final NNP covers an energy range of more than 500 kJ mol-1 with a mean absolute error of energy less than 4 kJ mol-1. The search scheme is useful not only for saccharides but also for highly flexible bio-molecules.
Collapse
Affiliation(s)
- Pei-Kang Tsou
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, 10617, Taiwan.
| | - Hai Thi Huynh
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, 10617, Taiwan. .,Molecular Science and Technology Program, Taiwan International Graduate Program, Academia Sinica, Taipei, 11529, Taiwan.,Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan
| | - Huu Trong Phan
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, 10617, Taiwan. .,Molecular Science and Technology Program, Taiwan International Graduate Program, Academia Sinica, Taipei, 11529, Taiwan.,Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei, 10617, Taiwan. .,Molecular Science and Technology Program, Taiwan International Graduate Program, Academia Sinica, Taipei, 11529, Taiwan.,Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan.,International Graduate Program of Molecular Science and Technology (NTU-MST), National Taiwan University, Taipei 10617, Taiwan
| |
Collapse
|
4
|
J N C, Mallajosyula SS. Impact of Polarization on the Ring Puckering Dynamics of Hexose Monosaccharides. J Chem Inf Model 2023; 63:208-223. [PMID: 36475659 DOI: 10.1021/acs.jcim.2c01286] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Analysis of crystal structures of hexose monosaccharides α-d-mannose (α-MAN), β-d-mannose (β-MAN), α-d-glucose (α-GLC), β-d-glucose (β-GLC), α-d-galactose (α-GAL), β-d-galactose (β-GAL), α-d-altrose (α-ALT), β-d-altrose (β-ALT), α-d-idose (α-IDO), and β-d-idose (β-IDO) reveals that the monosaccharide ring adopts multiple ring conformations. These ring conformations can be broadly classified as chair, half-chair, envelope, boat, and skew-boat conformations. The ability of the monosaccharide ring to adopt multiple conformations has been closely tied with their bioactivity. However, it has been difficult to capture the dynamic information of these conformations from experimental studies. Even from simulations, capturing these different conformations is challenging because of the energy barriers involved in the transitions between the stable 4C1 and 1C4 chair forms. In this study, we analyze the influence of the polarizable force field on the ring dynamics of five major types of unsubstituted aldohexoses─glucose, mannose, galactose, altrose, and idose─and their anomers. We simulate microsecond trajectories to capture the influence of the CHARMM36 additive and polarizable carbohydrate force fields on the ring dynamics. The microsecond trajectories allow us to comment on the issues associated with equilibrium molecular dynamics simulations. Further, we use the extended system adaptive biasing force (eABF) method to compare the conformational sampling efficiencies of the additive and polarizable force fields. Our studies reveal that inclusion of polarization enhances the sampling of ring conformations and lowers the energy barriers between the 4C1 and 1C4 conformations. Overall, the CHARMM36 additive force field is observed to be rigid and favor the 4C1 conformations. Although the inclusion of polarizability results in enhancing ring flexibility, we observe sampling that does not agree with experimental results, warranting a revision of the polarizable Drude parameters.
Collapse
Affiliation(s)
- Chythra J N
- Discipline of Chemistry, Indian Institute of Technology Gandhinagar, Gandhinagar, Gujarat382355, India
| | - Sairam S Mallajosyula
- Discipline of Chemistry, Indian Institute of Technology Gandhinagar, Gandhinagar, Gujarat382355, India
| |
Collapse
|
5
|
Jaeschke SO, Lindhorst TK, Auer A. Between Two Chairs: Combination of Theory and Experiment for the Determination of the Conformational Dynamics of Xylosides. Chemistry 2022; 28:e202201544. [PMID: 35754398 PMCID: PMC9804333 DOI: 10.1002/chem.202201544] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2022] [Indexed: 01/05/2023]
Abstract
The conformational properties of monosaccharides constitute fundamental features of oligosaccharides. While the energy landscape of monosaccharides can be altered by a specific biochemical environment or by chemical modifications, the analysis of resulting dynamic conformational equilibria is not feasible by experimental means alone. In this work, a series of β-d-xylopyranosides is used to outline how a combination of experimental NMR parameters and computed molecular properties can be used to determine conformers and quantify the composition of conformational equilibria. We demonstrate that identifying the most stable conformers using energy calculations is challenging and computing of NMR shieldings is typically not sensitive enough. On the other hand, computed spin-spin coupling constants for the xyloside ring can be used to unambiguously assign experimental NMR data of dynamic conformational equilibria and quantify the ratio of different conformers in the mixture. As a proof of principle, this procedure allowed to analyze a hitherto unknown dynamic equilibrium of a diamino-xyloside as a precursor of a molecular switch.
Collapse
Affiliation(s)
- Sven Ole Jaeschke
- Otto Diels Institute of Organic ChemistryChristiana Albertina University of KielOtto-Hahn-Platz 3–424118KielGermany
| | - Thisbe K. Lindhorst
- Otto Diels Institute of Organic ChemistryChristiana Albertina University of KielOtto-Hahn-Platz 3–424118KielGermany
| | - Alexander Auer
- Max-Planck-Institut für KohlenforschungKaiser-Wilhelm-Platz 145470Mülheim an der RuhrGermany
| |
Collapse
|
6
|
Computational quantitation of the aldehyde forms of aldohexoses and disaccharides composed of d-glucose: Predictions of their reactivities in the Maillard reaction. COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113605] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
7
|
Lutsyk V, Plazinski W. Conformational Properties of Glycosaminoglycan Disaccharides: A Molecular Dynamics Study. J Phys Chem B 2021; 125:10900-10916. [PMID: 34550710 DOI: 10.1021/acs.jpcb.1c04860] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
The structure and conformation of glycosaminoglycans (GAGs) are of central importance to understand the mechanisms behind their functions in biological systems. Due to the inherent chemical and structural heterogeneity of GAGs, focusing on longer, naturally existing GAG chains hinders drawing conclusions on the influence of the chemical functionalization on the basic conformational degree of freedom, that is, the dynamic shape of glycosidic linkage present in the particular disaccharide repeating unit. In the present study, we have considered the complete set of 106 GAG-related disaccharides, being potential building blocks for longer GAG chains (including hyaluronan, chondroitin, keratan, dermatan, and heparan). Both the unfunctionalized units and all possible combinations of either partially or fully sulfated derivatives contribute to this number. The unbiased and enhanced sampling molecular dynamics simulations provide a link to understand the influence of sulfation on the conformational properties of GAG glycosidic linkages. Residue-residue hydrogen bonding is not significant for either the glycosidic linkage conformation or its flexibility. It was found that in the majority of cases, the dominating conformation of the linkage is weakly affected by sulfation and the main role is played by the steric and stereoelectronic effects. However, there exist numerous cases where sulfation increases the contribution of alternative conformations to a nonnegligible extent and, in some rare cases (restricted to disaccharides building heparan), leads to the reorientation of the glycosidic linkage. The identified sulfation sites, being the most important in this context, are C6 and C3 at the GlcNAc residue. Finally, the full set of free energy maps relying on the glycosidic dihedral angle values for diverse GAG disaccharides are provided; they may be used for further studies, focused on longer GAG chains.
Collapse
Affiliation(s)
- Valery Lutsyk
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, Krakow 30-239, Poland
| | - Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, Krakow 30-239, Poland
| |
Collapse
|
8
|
Plazinski W, Roslund MU, Säwén E, Engström O, Tähtinen P, Widmalm G. Tautomers of N-acetyl-d-allosamine: an NMR and computational chemistry study. Org Biomol Chem 2021; 19:7190-7201. [PMID: 34382051 DOI: 10.1039/d1ob01139a] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
d-Allosamine is a rare sugar in Nature but its pyranoid form has been found α-linked in the core region of the lipopolysaccharide from the Gram-negative bacterium Porphyromonas gingivalis and in the chitanase inhibitor allosamidin, then β-linked and N-acetylated. In water solution the monosaccharide N-acetyl-d-allosamine (d-AllNAc) shows a significant presence of four tautomers arising from pyranoid and furanoid ring forms and anomeric configurations. The furanoid ring forms both showed 3JH1,H2≈ 4.85 Hz and to differentiate the anomeric configurations a series of chemical shift anisotropy/dipole-dipole cross-correlated relaxation NMR experiments was performed in which the α-anomeric form showed notable different relaxation rates for its components of the H1 doublet, thereby making it possible to elucidate the anomeric configuration of each of the furanoses. The conformational preferences of the different forms of d-AllNAc were investigated by 3JHH, 2JCH and 3JCH coupling constants from NMR experiments, molecular dynamics simulations and density functional theory calculations. The pyranose form resides in the 4C1 conformation and the furanose ring form has the majority of its conformers located on the South-East region of the pseudorotation wheel, with a small population in the Northern hemisphere. The tautomeric equilibrium was quite sensitive to changes in temperature, where the β-anomer of the pyranoid ring form decreased upon a temperature increase while the other forms increased.
Collapse
Affiliation(s)
- Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, 30-239 Krakow, Poland
| | | | | | | | | | | |
Collapse
|
9
|
Nester K, Plazinski W. Deciphering the conformational preferences of furanosides. A molecular dynamics study. J Biomol Struct Dyn 2020; 38:3359-3370. [DOI: 10.1080/07391102.2019.1656670] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Karina Nester
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Krakow, Poland
| | - Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Krakow, Poland
| |
Collapse
|
10
|
Farkas V, Nagy A, Menyhárd DK, Perczel A. Assignment of Vibrational Circular Dichroism Cross-Referenced Electronic Circular Dichroism Spectra of Flexible Foldamer Building Blocks: Towards Assigning Pure Chiroptical Properties of Foldamers. Chemistry 2019; 25:14890-14900. [PMID: 31464009 PMCID: PMC6899845 DOI: 10.1002/chem.201903023] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2019] [Revised: 08/15/2019] [Indexed: 01/09/2023]
Abstract
Assignment of the most established electronic circular dichroism (ECD) spectra of polypeptides and foldamers is either “evidence based” or relies on the 3D structures of longer oligomers of limited internal dynamics, which are derived from NMR spectroscopy (or X‐ray) data. Critics warn that the use of NMR spectroscopy and ECD side by side has severe limitations for flexible molecules because explicit knowledge of conformational ensembles is a challenge. Herein, an old–new method of comparing ab initio computed and measured vibrational circular dichroism (VCD) data is presented to validate both the structures (conf(i)) and their relative weights (c(i)) that make up the conformational ensemble. Based on the array of {conf(i), c(i)}, the pure ECD spectra, g(i)conf(i), can be ab initio calculated. The reconstructed spectrum Σc(i)g(i)conf(i) can thus help to assign any experimental ECD counterparts. Herein, such a protocol is successfully applied to flexible foldamer building blocks of sugar β‐amino acid diamides. The epimeric pair of the model system was selected because these molecules were conformationally tunable by simple chemical modification, and thus, the robustness of the current approach could be probed. The initial hydrogen bond (NH⋅⋅⋅O) eliminated by N‐methylation reorients the amide plain, which influences the chiroptical properties of the foldamer building block; this structural change is successfully monitored by changes to the VCD and ECD transitions, which are now assigned to pure conformers. The current method seems to be general and effective without requiring extensive CPU and spectroscopic resources.
Collapse
Affiliation(s)
- Viktor Farkas
- MTA-ELTE Protein Modelling Research Group, Institute of Chemistry, Eötvös Loránd University, Pázmány P. stny. 1/A, Budapest, 1117, Hungary
| | - Adrienn Nagy
- Laboratory of Structural Chemistry and Biology, Institute of Chemistry, Eötvös Loránd University, Pázmány P. stny. 1/A, Budapest, 1117, Hungary
| | - Dóra K Menyhárd
- MTA-ELTE Protein Modelling Research Group, Institute of Chemistry, Eötvös Loránd University, Pázmány P. stny. 1/A, Budapest, 1117, Hungary
| | - András Perczel
- MTA-ELTE Protein Modelling Research Group, Institute of Chemistry, Eötvös Loránd University, Pázmány P. stny. 1/A, Budapest, 1117, Hungary.,Laboratory of Structural Chemistry and Biology, Institute of Chemistry, Eötvös Loránd University, Pázmány P. stny. 1/A, Budapest, 1117, Hungary
| |
Collapse
|
11
|
Tautomeric and epimeric equilibria of aldo- and ketohexoses studied by the MD simulations and QM calculations. Carbohydr Res 2019; 474:8-15. [DOI: 10.1016/j.carres.2019.01.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2018] [Revised: 01/04/2019] [Accepted: 01/09/2019] [Indexed: 11/17/2022]
|
12
|
Nester K, Gaweda K, Plazinski W. A GROMOS Force Field for Furanose-Based Carbohydrates. J Chem Theory Comput 2019; 15:1168-1186. [DOI: 10.1021/acs.jctc.8b00838] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Karina Nester
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek Str., 8, 30-239 Cracow, Poland
| | - Karolina Gaweda
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek Str., 8, 30-239 Cracow, Poland
| | - Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek Str., 8, 30-239 Cracow, Poland
| |
Collapse
|
13
|
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
|
14
|
Gebhardt J, Kleist C, Jakobtorweihen S, Hansen N. Validation and Comparison of Force Fields for Native Cyclodextrins in Aqueous Solution. J Phys Chem B 2018; 122:1608-1626. [PMID: 29287148 DOI: 10.1021/acs.jpcb.7b11808] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Molecular dynamics simulations of native α-, β-, and γ-cyclodextrin in aqueous solution have been conducted with the goal to investigate the performance of the CHARMM36 force field, the AMBER-compatible q4md-CD force field, and five variants of the GROMOS force field. The properties analyzed are structural parameters derived from X-ray diffraction and NMR experiments as well as hydrogen bonds and hydration patterns, including hydration free enthalpies. Recent revisions of the torsional-angle parameters for carbohydrate systems within the GROMOS family of force fields lead to a significant improvement of the agreement between simulated and experimental NMR data. Therefore, we recommend using the variant 53A6GLYC instead of 53A6 and 56A6CARBO_R or 2016H66 instead of 56A6CARBO to simulate cyclodextrins in solution. The CHARMM36 and q4md-CD force fields show a similar performance as the three recommended GROMOS parameter sets. A significant difference is the more flexible nature of the cyclodextrins modeled with the CHARMM36 and q4md-CD force fields compared to the three recommended GROMOS parameter sets.
Collapse
Affiliation(s)
- Julia Gebhardt
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart , D-70569 Stuttgart, Germany
| | - Catharina Kleist
- Institute of Thermal Separation Processes, Hamburg University of Technology , D-21073 Hamburg, Germany
| | - Sven Jakobtorweihen
- Institute of Thermal Separation Processes, Hamburg University of Technology , D-21073 Hamburg, Germany
| | - Niels Hansen
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart , D-70569 Stuttgart, Germany
| |
Collapse
|
15
|
Panczyk K, Plazinski W. Pyranose ring puckering in aldopentoses, ketohexoses and deoxyaldohexoses. A molecular dynamics study. Carbohydr Res 2017; 455:62-70. [PMID: 29175656 DOI: 10.1016/j.carres.2017.11.011] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2017] [Revised: 11/14/2017] [Accepted: 11/17/2017] [Indexed: 01/03/2023]
Abstract
Conformation of monosaccharides, including the ring shape, has for years been the subject of intensive research. Although d-aldohexopyranoses are the most extensively studied pyranoses, there also exist other groups of saccharides that contain analogous chemical system of the six-membered ring. Here we describe in details the results of the molecular dynamics-based conformational analysis concerning a series of pyranoses, namely: d-aldopentoses, d-ketohexoses as well as deoxy- (d-quinovose, l-fucose, l-rhamnose) and dideoxy- (abequose, paratose, tyvelose, digitoxose) derivatives of aldohexoses. By using the carbohydrate-dedicated GROMOS 56a6CARBO force field, we determined the conformational properties of both the lactol and hydroxymethyl groups as well as the anomeric populations for all considered compounds. The orientation of the lactol group follows the trend expected on the basis of the exo-anomeric effect for all compounds whereas the conformation of the hydroxymethyl group in d-ketohexoses is represented by the two gauche (with respect to the ring oxygen atom) rotamers. The special emphasis is put on the ring-inversion properties studied in the context of both the full chair-chair inversion and the chair-boat/skew-boat rearrangement. The calculated ring-distortion energies, compared with those obtained for regular d-aldohexopyranoses allowed for estimating the influence of particular substituents on the ring flexibility. Overall, such influence is correlated with the dimension of the substituent and its orientation but is limited to the case of the chair-chair inversion whereas the chair-to-boat/skew-boat rearrangement exhibits roughly the same properties for all pyranoses. For all d-aldopyranoses the α anomers exhibit lower ring-inversion free energies in comparison to the β anomers whereas this trend is inverted in the case of d-ketohexopyranoses.
Collapse
Affiliation(s)
- Karina Panczyk
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, ul. Niezapominajek 8, 30-239, Cracow, Poland
| | - Wojciech Plazinski
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, ul. Niezapominajek 8, 30-239, Cracow, Poland.
| |
Collapse
|