1
|
Dorst KM, Widmalm G. Conformational Preferences at the Glycosidic Linkage of Saccharides in Solution as Deduced from NMR Experiments and MD Simulations: Comparison to Crystal Structures. Chemistry 2024; 30:e202304047. [PMID: 38180821 DOI: 10.1002/chem.202304047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2023] [Revised: 12/30/2023] [Accepted: 01/05/2024] [Indexed: 01/07/2024]
Abstract
Glycans are central to information content and regulation in biological systems. These carbohydrate molecules are active either as oligo- or polysaccharides, often in the form of glycoconjugates. The monosaccharide entities are joined by glycosidic linkages and stereochemical arrangements are of utmost importance in determining conformation and flexibility of saccharides. The conformational preferences and population distributions at the glycosidic torsion angles φ and ψ have been investigated for O-methyl glycosides of three disaccharides where the substitution takes place at a secondary alcohol, viz., in α-l-Fucp-(1→3)-β-d-Glcp-OMe, α-l-Fucp-(1→3)-α-d-Galp-OMe and α-d-Glcp-(1→4)-α-d-Galp-OMe, corresponding to disaccharide structural elements present in bacterial polysaccharides. Stereochemical differences at or adjacent to the glycosidic linkage were explored by solution state NMR spectroscopy using one-dimensional 1 H,1 H-NOESY NMR experiments to obtain transglycosidic proton-proton distances and one- and two-dimensional heteronuclear NMR experiments to obtain 3 JCH transglycosidic coupling constants related to torsion angles φ and ψ. Computed effective proton-proton distances from molecular dynamics (MD) simulations showed excellent agreement to experimentally derived distances for the α-(1→3)-linked disaccharides and revealed that for the bimodal distribution at the ψ torsion angle for the α-(1→4)-linked disaccharide experiment and simulation were at variance with each other, calling for further force field developments. The MD simulations disclosed a highly intricate inter-residue hydrogen bonding pattern for the α-(1→4)-linked disaccharide, including a nonconventional hydrogen bond between H5' in the glucosyl residue and O3 in the galactosyl residue, supported by a large downfield 1 H NMR chemical shift displacement compared to α-d-Glcp-OMe. Comparison of population distributions of the glycosidic torsion angles φ and ψ in the disaccharide entities to those of corresponding crystal structures highlighted the potential importance of solvation on the preferred conformation.
Collapse
Affiliation(s)
- Kevin M Dorst
- 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
|
2
|
Deng C, Zhang T, Zhang X, Gu T, Xu L, Yu Z, Zheng M, Zhou Y. Multiscale structure and precipitation mechanism of debranched starch precipitated by different alcohols. Int J Biol Macromol 2023; 241:124562. [PMID: 37088190 DOI: 10.1016/j.ijbiomac.2023.124562] [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: 12/30/2022] [Revised: 02/22/2023] [Accepted: 04/18/2023] [Indexed: 04/25/2023]
Abstract
Alcohol solution is a cheap, simple, and effective precipitating solvent frequently used for separating debranched starch (DBS), yet little is known about the precipitation mechanism of DBS by different alcohols. This study precipitated DBS from pullulanase-hydrolyzed starch using ethanol, n-butanol, and isopentanol. The multiscale structures of DBS were characterized, including chain length, single/double helix, and crystalline. The chain conformation and precipitation mechanism of DBS in different alcohols was investigated using molecular dynamics (MD) simulation. DBS precipitated by n-butanol contained the largest proportion of short chain (DP6-24, 83.2 %), the highest V-type crystallinity (21.1 %), and the largest single-helix content (24.7 %). A single helix conformation of DBS chain was determined in alcohols, where alcohol molecules entered the helix cavity. Intra/inter-molecular hydrogen bonds stabilized the helix, with a large number of hydrogen bonds leading to strong molecular interaction and stable helical structure. The solvent accessible surface area of DBS chain decreased by 7.88-19.32 % in alcohols, and the radial distribution function revealed that the first solvent layer of DBS chain at 0.29 nm was closely related to hydrogen bonding. This study provides a basis for the choice of precipitation solvent for preparing DBS with different chain lengths and physicochemical properties.
Collapse
Affiliation(s)
- Changyue Deng
- Key Laboratory of Jianghuai Agricultural Product Fine Processing and Resource Utilization of Ministry of Agriculture and Rural Affairs, College of Tea & Food Science and Technology, Anhui Agricultural University, Hefei 230036, China; Food Processing Research Institute, Anhui Agricultural University, Hefei 230036, China
| | - Tiantian Zhang
- Key Laboratory of Jianghuai Agricultural Product Fine Processing and Resource Utilization of Ministry of Agriculture and Rural Affairs, College of Tea & Food Science and Technology, Anhui Agricultural University, Hefei 230036, China; Food Processing Research Institute, Anhui Agricultural University, Hefei 230036, China
| | - Xiumei Zhang
- Key Laboratory of Jianghuai Agricultural Product Fine Processing and Resource Utilization of Ministry of Agriculture and Rural Affairs, College of Tea & Food Science and Technology, Anhui Agricultural University, Hefei 230036, China; Food Processing Research Institute, Anhui Agricultural University, Hefei 230036, China
| | - Tingting Gu
- Key Laboratory of Jianghuai Agricultural Product Fine Processing and Resource Utilization of Ministry of Agriculture and Rural Affairs, College of Tea & Food Science and Technology, Anhui Agricultural University, Hefei 230036, China; Food Processing Research Institute, Anhui Agricultural University, Hefei 230036, China
| | - Li Xu
- Key Laboratory of Jianghuai Agricultural Product Fine Processing and Resource Utilization of Ministry of Agriculture and Rural Affairs, College of Tea & Food Science and Technology, Anhui Agricultural University, Hefei 230036, China; Food Processing Research Institute, Anhui Agricultural University, Hefei 230036, China
| | - Zhenyu Yu
- Key Laboratory of Jianghuai Agricultural Product Fine Processing and Resource Utilization of Ministry of Agriculture and Rural Affairs, College of Tea & Food Science and Technology, Anhui Agricultural University, Hefei 230036, China; Food Processing Research Institute, Anhui Agricultural University, Hefei 230036, China
| | - Mingming Zheng
- Food Processing Research Institute, Anhui Agricultural University, Hefei 230036, China; Key Laboratory of Oilseeds Processing, Ministry of Agriculture, Hubei Key Laboratory of Lipid Chemistry and Nutrition, Oil Crops and Lipids Process Technology National & Local Joint Engineering Laboratory, Oil Crops Research Institute, Chinese Academy of Agricultural Sciences, Wuhan 430062, China
| | - Yibin Zhou
- Key Laboratory of Jianghuai Agricultural Product Fine Processing and Resource Utilization of Ministry of Agriculture and Rural Affairs, College of Tea & Food Science and Technology, Anhui Agricultural University, Hefei 230036, China; Food Processing Research Institute, Anhui Agricultural University, Hefei 230036, China.
| |
Collapse
|
3
|
Cao C, Deng C, Hu J, Zhou Y. Formation and molecular dynamics simulation of inclusion complex of large-ring cyclodextrin and 4-terpineol. J Food Sci 2022; 87:4609-4621. [PMID: 36086891 DOI: 10.1111/1750-3841.16303] [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: 04/11/2022] [Revised: 07/22/2022] [Accepted: 08/02/2022] [Indexed: 12/01/2022]
Abstract
In the present study, the formation and structure of the inclusion compound of large-ring cyclodextrin and 4-terpineol were obtained through different experiments and molecular dynamics (MD) simulation. The analysis of FTIR, 1 H-NMR, and thermodynamic results confirmed the formation of clathrates. Analysis of molecular structure (root-mean-square deviation and radius of gyration), solubility, and interaction energy (Coul, H bond) based on MD simulations further clarified the nature of the clathrate and the conformational changes caused by guest molecules as well as inclusion complexes process trends. The inclusion complex reportedly has a new crystal structure with improved thermal stability. PRACTICAL APPLICATION: This is the first work to demonstrate the complex formation between 4-terpineol and large-ring cyclodextrin by molecular dynamics simulation. Molecular dynamics simulation confirmed the formation of inclusion complexes theoretically. Conformational changes of the molecules and the formation of complexes with improved thermal stability were observed. Complexing with large-ring cyclodextrin can be used as an effective means to encapsulate the aroma/flavor compounds.
Collapse
Affiliation(s)
- Chuan Cao
- The College of Environment and Chemical Engineering, Anhui Vocational and Technical College, Hefei, China.,Anhui Engineering Laboratory of Agricultural Products Processing, Anhui Agricultural University, Hefei, China
| | - Changyue Deng
- Anhui Engineering Laboratory of Agricultural Products Processing, Anhui Agricultural University, Hefei, China
| | - Jinwei Hu
- Anhui Engineering Laboratory of Agricultural Products Processing, Anhui Agricultural University, Hefei, China
| | - Yibin Zhou
- Anhui Engineering Laboratory of Agricultural Products Processing, Anhui Agricultural University, Hefei, China
| |
Collapse
|
4
|
Shao L, Huo Z, Lei N, Yang M, He Z, Zhang Y, Wei Q, Chen C, Xiao M, Wang F, Gu G, Cai F. Reinvestigation of N, N-Diacetylimido-Protected 2-Aminothioglycosides in O-Glycosylation: Intermolecular Hydrogen Bonds Contributing to 1,2-Orthoamide Formation. J Org Chem 2021; 86:13212-13230. [PMID: 34533021 DOI: 10.1021/acs.joc.1c01009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
N,N-Diacetylimido protection of 2-aminoglycosides is an elegant strategy but has had limited applications due to unexpected side reactions in glycosylation. We found that high acid concentrations could diminish the side reactions. We observed intermolecular hydrogen bonding among alcohols and acids could disrupt. Assuming that intermolecular hydrogen bonding accelerates the formation of 1,2-orthoamides and disrupting intermolecular hydrogen bonds could turn to the desired glycosylation, we successfully employed sulfenyl triflate pre-activation in the glycosylation of a broad scope of alcohol acceptors, as well as in a one-pot synthesis of a protected human milk oligosaccharide, lacto-N-neotetraose.
Collapse
Affiliation(s)
- Liming Shao
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| | - Zhenni Huo
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| | - Na Lei
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| | - Min Yang
- Center for Analysis and Characterization, School of Physical Science and Technology, ShanghaiTech University, 393 Huaxia Middle Rd., Shanghai 201210, China
| | - Zehuan He
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| | - Yongliang Zhang
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| | - Qinlong Wei
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| | - Changsheng Chen
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| | - Mei Xiao
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| | - Fei Wang
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| | - Guofeng Gu
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| | - Feng Cai
- National Glycoengineering Research Center and Shandong Provincial Key Laboratory of Carbohydrate Chemistry and Glycobiology, Shandong University, 72 Binhai Road, Qingdao 266237, China
| |
Collapse
|
5
|
Qiu Y, Smith DGA, Stern CD, Feng M, Jang H, Wang LP. Driving torsion scans with wavefront propagation. J Chem Phys 2020; 152:244116. [PMID: 32610969 DOI: 10.1063/5.0009232] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
The parameterization of torsional/dihedral angle potential energy terms is a crucial part of developing molecular mechanics force fields. Quantum mechanical (QM) methods are often used to provide samples of the potential energy surface (PES) for fitting the empirical parameters in these force field terms. To ensure that the sampled molecular configurations are thermodynamically feasible, constrained QM geometry optimizations are typically carried out, which relax the orthogonal degrees of freedom while fixing the target torsion angle(s) on a grid of values. However, the quality of results and computational cost are affected by various factors on a non-trivial PES, such as dependence on the chosen scan direction and the lack of efficient approaches to integrate results started from multiple initial guesses. In this paper, we propose a systematic and versatile workflow called TorsionDrive to generate energy-minimized structures on a grid of torsion constraints by means of a recursive wavefront propagation algorithm, which resolves the deficiencies of conventional scanning approaches and generates higher quality QM data for force field development. The capabilities of our method are presented for multi-dimensional scans and multiple initial guess structures, and an integration with the MolSSI QCArchive distributed computing ecosystem is described. The method is implemented in an open-source software package that is compatible with many QM software packages and energy minimization codes.
Collapse
Affiliation(s)
- Yudong Qiu
- Department of Chemistry, UC Davis, Davis, California 95616, USA
| | - Daniel G A Smith
- The Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - Chaya D Stern
- Computational and Systems Biology Program, Sloan-Kettering Institute, New York, New York 10065, USA
| | - Mudong Feng
- Department of Chemistry and Biochemistry, UC San Diego, La Jolla, California 92093, USA
| | - Hyesu Jang
- Department of Chemistry, UC Davis, Davis, California 95616, USA
| | - Lee-Ping Wang
- Department of Chemistry, UC Davis, Davis, California 95616, USA
| |
Collapse
|
6
|
Cheng L, Zhu X, Hamaker BR, Zhang H, Campanella OH. Complexation process of amylose under different concentrations of linoleic acid using molecular dynamics simulation. Carbohydr Polym 2019; 216:157-166. [DOI: 10.1016/j.carbpol.2019.04.013] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2019] [Revised: 03/14/2019] [Accepted: 04/02/2019] [Indexed: 02/09/2023]
|
7
|
NMR solution geometry of saccharides containing the 6-O-(α-D-glucopyranosyl)-α/β-D-glucopyranose (isomaltose) or 6-O-(α-D-galactopyranosyl)-α/β-D-glucopyranose (melibiose) core. Carbohydr Res 2019; 473:18-35. [DOI: 10.1016/j.carres.2018.12.012] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2018] [Revised: 12/19/2018] [Accepted: 12/19/2018] [Indexed: 12/26/2022]
|
8
|
Vickman AE, Pohl NLB. Probing deoxysugar conformational preference: A comprehensive computational study investigating the effects of deoxygenation. Carbohydr Res 2018; 475:17-26. [PMID: 30771703 DOI: 10.1016/j.carres.2018.12.003] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2018] [Revised: 11/26/2018] [Accepted: 12/06/2018] [Indexed: 10/27/2022]
Abstract
Deoxysugars are intrinsic components in a number of antibiotics, antimicrobials, and therapeutic agents that often dictate receptor binding, improve efficacy, and provide a diverse toolbox in modifying glycoconjugate function due to an extensive number of unique isomers and inherent conformational flexibility. Hence, this work provides a comprehensive examination of the conformational effects associated with deoxygenation of the pyranose ring. Both the location and degree of deoxygenation were evaluated by interrogating the energetic landscape for a number of mono- and dideoxyhexopyranose derivatives using DFT methods (M05-2X/cc-pVTZ(-f)). Both anomeric forms and in some cases, the alternate chair form, have been investigated in the gas phase. As was documented in a preceding study, variation of the C-6 oxidation state has been shown to affect the anomeric preference of select glucose stereoisomers. Similar results were also observed for several deoxysugar isomers in this work, wherein the alternate anomer was favored upon reduction to the 6-deoxyhexose derivative or oxidation to the hexonic acid. Additionally, comparison of relative Gibbs free energies revealed C-3 deoxygenation imparts greater instability compared to C-2 or C-4 deoxygenation, as indicated by an increase in free energy for 3-deoxysugars. A polarizable continuum solvation model was also applied to empirically validate theoretical results for several deoxysugars, wherein good agreement with both carbon (σ = 1.6 ppm) and proton (σ = 0.20 ppm) NMR shifts was observed for the majority of isomers. Solvated and gas phase anomeric ratios were also calculated and compared favorably to reported literature values, although some discrepancies are noted.
Collapse
Affiliation(s)
- Alison E Vickman
- Department of Chemistry, Indiana University, Bloomington, IN, 47405, USA
| | - Nicola L B Pohl
- Department of Chemistry, Indiana University, Bloomington, IN, 47405, USA.
| |
Collapse
|
9
|
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
|
10
|
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
|