1
|
Sandilya AA, Priya MH. The counteracting influence of 2-hydroxypropyl substitution and the presence of a guest molecule on the shape and size of the β-cyclodextrin cavity. Phys Chem Chem Phys 2024; 26:11531-11544. [PMID: 38323885 DOI: 10.1039/d3cp05354g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2024]
Abstract
The aqueous solubility of β-cyclodextrin (β-CD), a cyclic carbohydrate comprising seven α-D-glucose molecules, is enhanced by 2-hydroxypropyl (2-HP) substitution of the hydroxyl groups at the CD rims. Our thorough analysis of the structural and solvation properties with different degrees of 2-hydroxypropyl substitution on β-CD using molecular dynamics simulations reveals that the solubility is enhanced at the cost of the structural distortion of the CD cyclic structure. Substitution at the secondary rim predominantly enhances the favourable interactions between CD and water by decreasing CD-CD hydrogen bonding and promoting CD-water hydrogen bonding. However, the effect of substitution at the primary rim on the CD-water interactions is minimal; the hydrogen bonds between water and the primary hydroxyl group in native CD merely get replaced by those between water and 2-HP, since the substitution makes the primary hydroxyl oxygen (O6 atom) inaccessible to water. In contrast, substitution at the primary rim maintains the structural integrity of CD, while substitution at the secondary rim results in structural distortion due to the disruption of the intramolecular hydrogen bond belt, even leading to cavity closure. Certain strategic substitutions of the primary hydroxyl groups can help in the reduction of structural distortion, depending upon the degree of substitution at the secondary hydroxyl rim. A detailed inspection of the simulation trajectory revealed that the tilting of glucose units with the primary hydroxyl oxygen (O6) pointing inward is the primary driver for cavity closure. Even though the dynamics of glucose tilting can influence the kinetics of host-guest complex formation, once the guest is well incorporated into the cavity, glucose tilting is inhibited and the cavity opens up as in native β-CD.
Collapse
Affiliation(s)
- Avilasha A Sandilya
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai-600036, India.
| | - M Hamsa Priya
- Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai-600036, India.
| |
Collapse
|
2
|
Lebedinskiy K, Barvík I, Tošner Z, Císařová I, Jindřich J, Hrdina R. Spatial arrangements of cyclodextrin host-guest complexes in solution studied by 13C NMR and molecular modelling. Beilstein J Org Chem 2024; 20:331-335. [PMID: 38410781 PMCID: PMC10896224 DOI: 10.3762/bjoc.20.33] [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: 12/20/2023] [Accepted: 02/09/2024] [Indexed: 02/28/2024] Open
Abstract
13C NMR spectroscopic analyses of Cs symmetric guest molecules in the cyclodextrin host cavity, combined with molecular modelling and solid-state X-ray analysis, provides a detailed description of the spatial arrangement of cyclodextrin host-guest complexes in solution. The chiral cavity of the cyclodextrin molecule creates an anisotropic environment for the guest molecule resulting in a splitting of its prochiral carbon signals in 13C NMR spectra. This signal split can be correlated to the distance of the guest atoms from the wall of the host cavity and to the spatial separation of binding sites preferred by pairs of prochiral carbon atoms. These measurements complement traditional solid-state analyses, which rely on the crystallization of host-guest complexes and their crystallographic analysis.
Collapse
Affiliation(s)
- Konstantin Lebedinskiy
- Department of Organic Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 43 Praha, Czech Republic
| | - Ivan Barvík
- Institute of Physics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 2026/5, 121 16 Praha, Czech Republic
| | - Zdeněk Tošner
- Department of Organic Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 43 Praha, Czech Republic
| | - Ivana Císařová
- Department of Inorganic Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 43 Praha, Czech Republic
| | - Jindřich Jindřich
- Department of Organic Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 43 Praha, Czech Republic
| | - Radim Hrdina
- Department of Organic Chemistry, Faculty of Science, Charles University, Hlavova 8, 128 43 Praha, Czech Republic
| |
Collapse
|
3
|
Kou X, Su D, Pan F, Xu X, Meng Q, Ke Q. Molecular dynamics simulation techniques and their application to aroma compounds/cyclodextrin inclusion complexes: A review. Carbohydr Polym 2024; 324:121524. [PMID: 37985058 DOI: 10.1016/j.carbpol.2023.121524] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2023] [Revised: 10/16/2023] [Accepted: 10/22/2023] [Indexed: 11/22/2023]
Abstract
Homeostatic technologies play a crucial role in maintaining the quality and extending the service life of aroma compounds (ACs). Commercial cyclodextrins (CDs) are commonly used to form inclusion complexes (ICs) with ACs to enhance their solubility, stability, and morphology. The selection of suitable CDs and ACs is of paramount importance in this process. Molecular dynamics (MD) simulations provide an in-depth understanding of the interactions between ACs and CDs, aiding researchers in optimising the properties and effects of ICs. This review offers a systematic discussion of the application of MD simulations in ACs/CDs ICs, covering the establishment of the simulation process, parameter selection, model evaluation, and various application cases, along with their advantages and disadvantages. Additionally, this review summarises the major achievements and challenges of this method while identifying areas that require further exploration. These findings may contribute to a comprehensive understanding of the formation and stabilization mechanisms of ACs/CDs ICs and offer guidance for the selection and computational characterisation of CDs in the AC steady state.
Collapse
Affiliation(s)
- Xingran Kou
- Collaborative Innovation Center of Fragrance Flavour and Cosmetics, School of Perfume and Aroma Technology (Shanghai Research Institute of Fragrance & Flavour Industry), Shanghai Institute of Technology, Shanghai, China; Key Laboratory of Textile Science & Technology, Ministry of Education, College of Textiles, Donghua University, Shanghai, China
| | - Dongdong Su
- Collaborative Innovation Center of Fragrance Flavour and Cosmetics, School of Perfume and Aroma Technology (Shanghai Research Institute of Fragrance & Flavour Industry), Shanghai Institute of Technology, Shanghai, China
| | - Fei Pan
- State Key Laboratory of Resource Insects, Institute of Apicultural Research, Chinese Academy of Agricultural Sciences, Beijing 100093, China
| | - Xiwei Xu
- Collaborative Innovation Center of Fragrance Flavour and Cosmetics, School of Perfume and Aroma Technology (Shanghai Research Institute of Fragrance & Flavour Industry), Shanghai Institute of Technology, Shanghai, China
| | - Qingran Meng
- Collaborative Innovation Center of Fragrance Flavour and Cosmetics, School of Perfume and Aroma Technology (Shanghai Research Institute of Fragrance & Flavour Industry), Shanghai Institute of Technology, Shanghai, China.
| | - Qinfei Ke
- Collaborative Innovation Center of Fragrance Flavour and Cosmetics, School of Perfume and Aroma Technology (Shanghai Research Institute of Fragrance & Flavour Industry), Shanghai Institute of Technology, Shanghai, China; Key Laboratory of Textile Science & Technology, Ministry of Education, College of Textiles, Donghua University, Shanghai, China.
| |
Collapse
|
4
|
Alvira E. Theoretical study of structures and charge distributions of 2-, 3- and 6-hydroxypropyl-β-cyclodextrin with different degrees of substitution. J Biomol Struct Dyn 2023; 41:10582-10590. [PMID: 36524475 DOI: 10.1080/07391102.2022.2155700] [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: 07/13/2022] [Accepted: 11/29/2022] [Indexed: 12/23/2022]
Abstract
Hydroxypropyl-β-cyclodextrin (HPβCD) is a derivative of cyclodextrin extensively used in the pharmaceutical industry, since it improves the solubility of drugs, and widens their oral bioavailability and safety profile. Theoretical studies about HPβCD configurations are important so as to simulate by molecular mechanics and dynamics, the complex formation and enantiomeric separations of these cyclodextrins with other molecules. Twelve structures and charge distributions of 2-, 3- and 6-hydroxypropyl-β-cyclodextrin (2-, 3-, 6-HPβCD) with different degrees of substitution were obtained using ab initio methods. The atomic positions of glucose units, dihedral angles of hydroxypropyl groups (HPs), radius of gyration, and H-bonds formed were analysed in the structures. The 3-HPβCD configurations showed the greatest variations in their atomic positions, their HPs groups being slightly rotated towards the interior of the rim of cyclodextrin and hence narrowing its opening. The structures of 2-HPβCD are the least influenced by the degree of substitution and include fewer H-bonds. Different values for the partial atomic charges for each glucose unit in a molecule were obtained, as well as for all HPβCD configurations. This result, consequence of the ab initio methods used, is a new contribution that can be important to simulate processes in which small energy differences decisively influence the results. The electric potential due to the charge distribution of HPβCDs was calculated inside and outside the cavity. It was found that the changes at each position were mainly due to atomic configurations, since the differences in partial atomic charges are one or two orders of magnitude smaller than in atomic positions.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Elena Alvira
- Department of Physics, University of La Laguna, La Laguna, Spain
| |
Collapse
|
5
|
Wang X, Huai Z, Sun Z. Host Dynamics under General-Purpose Force Fields. Molecules 2023; 28:5940. [PMID: 37630194 PMCID: PMC10458655 DOI: 10.3390/molecules28165940] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2023] [Revised: 08/02/2023] [Accepted: 08/07/2023] [Indexed: 08/27/2023] Open
Abstract
Macrocyclic hosts as prototypical receptors to gaseous and drug-like guests are crucial components in pharmaceutical research. The external guests are often coordinated at the center of these macromolecular containers. The formation of host-guest coordination is accompanied by the broken of host-water and host-ion interactions and sometimes also involves some conformational rearrangements of the host. A balanced description of various components of interacting terms is indispensable. However, up to now, the modeling community still lacks a general yet detailed understanding of commonly employed general-purpose force fields and the host dynamics produced by these popular selections. To fill this critical gap, in this paper, we profile the energetics and dynamics of four types of popular macrocycles, including cucurbiturils, pillararenes, cyclodextrins, and octa acids. The presented investigations of force field definitions, refitting, and evaluations are unprecedently detailed. Based on the valuable observations and insightful explanations, we finally summarize some general guidelines on force field parametrization and selection in host-guest modeling.
Collapse
Affiliation(s)
- Xiaohui Wang
- Beijing Leto Laboratories Co., Ltd., Beijing 100083, China
- College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Zhe Huai
- XtalPi—AI Research Center, 7F, Tower A, Dongsheng Building, No. 8, Zhongguancun East Road, Beijing 100083, China
| | - Zhaoxi Sun
- College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
6
|
Triolo A, Lo Celso F, Fourmentin S, Russina O. Liquid Structure Scenario of the Archetypal Supramolecular Deep Eutectic Solvent: Heptakis(2,6-di- O-methyl)-β-cyclodextrin/levulinic Acid. ACS SUSTAINABLE CHEMISTRY & ENGINEERING 2023; 11:9103-9110. [PMID: 37351462 PMCID: PMC10283020 DOI: 10.1021/acssuschemeng.3c01858] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/29/2023] [Revised: 05/24/2023] [Indexed: 06/24/2023]
Abstract
The concept of supramolecular solvents has been recently introduced, and the extended liquid-state window accessible for mixtures of functionalized cyclodextrins (CDs) with hydrogen bond (HB) donor species, e.g., levulinic acid, led to the debut of supramolecular deep eutectic solvents (SUPRA-DES). These solvents retain CD's inclusion ability and complement it with enhanced solvation effectiveness due to an extended HB network. However, so far, these promising features were not rationalized in terms of a microscopic description, thus hindering a more complete capitalization. This is the first joint experimental and computational study on the archetypal SUPRA-DES: heptakis(2,6-di-O-methyl)-β-CD/levulinic acid (1:27). We used X-ray scattering to probe CD's aggregation level and molecular dynamics simulation to determine the nature of interactions between SUPRA-DES components. We discover that CDs are homogeneously distributed in bulk and that HB interactions, together with the electrostatic ones, play a major role in determining mutual interaction between components. However, dispersive forces act in synergy with HB to accomplish a fundamental task in hindering hydrophobic interactions between neighbor CDs and maintaining the system homogeneity. The mechanism of mutual solvation of CD and levulinic acid is fully described, providing fundamental indications on how to extend the spectrum of SUPRA-DES combinations. Overall, this study provides the key to interpreting structural organization and solvation tunability in SUPRA-DES to extend the range of sustainable applications for these new, unique solvents.
Collapse
Affiliation(s)
- Alessandro Triolo
- Laboratorio
Liquidi Ionici, Istituto Struttura della
Materia-Consiglio Nazionale delle Ricerche (ISM-CNR), Rome 00133, Italy
| | - Fabrizio Lo Celso
- Laboratorio
Liquidi Ionici, Istituto Struttura della
Materia-Consiglio Nazionale delle Ricerche (ISM-CNR), Rome 00133, Italy
- Department
of Physics and Chemistry, Università
di Palermo, Palermo 90133, Italy
| | - Sophie Fourmentin
- Unité
de Chimie Environnementale et Interactions sur le Vivant (UCEIV, UR
4492), Université du Littoral Côte
d’Opale (ULCO), 59140 Dunkerque, France
| | - Olga Russina
- Laboratorio
Liquidi Ionici, Istituto Struttura della
Materia-Consiglio Nazionale delle Ricerche (ISM-CNR), Rome 00133, Italy
- Department
of Chemistry, Sapienza University of Rome, Rome 00185, Italy
| |
Collapse
|
7
|
Ganjali Koli M, Eshaghi Malekshah R, Hajiabadi H. Insights from molecular dynamics and DFT calculations into the interaction of 1,4-benzodiazepines with 2-hydroxypropyl-βCD in a theoretical study. Sci Rep 2023; 13:9866. [PMID: 37332009 DOI: 10.1038/s41598-023-36385-w] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Accepted: 06/02/2023] [Indexed: 06/20/2023] Open
Abstract
This study delves into the interaction between benzodiazepine (BZD) drugs and 2-hydroxypropyl-β-cyclodextrin (2HPβCD), a cyclodextrin (CD) known to improve drug delivery and enhance therapeutic outcomes. We find that the 2HPβCD's atoms become more rigid in the presence of chlordiazepoxide (CDP), clonazepam (CLZ), and diazepam (DZM), whereas they become more flexible in the presence of nordazepam (NDM) and nitrazepam (NZP). We also investigated the structure of 2HPβCD and found that loading these drugs increases both the area and volume of the 2HPβCD cavity, making it more suitable for drug delivery. Moreover, this research found that all drugs exhibited negative values for the binding free energy, indicating thermodynamic favorability and improved solubility. The binding free energy order of the BZDs was consistent in both molecular dynamics and Monte Carlo methods, with CDP and DZM having the highest affinity for binding. We also analyzed the contribution of different interaction energies in binding between the carrier and the drugs and found that Van der Waals energy is the primary component. Our results indicate that the number of hydrogen bonds between 2HPβCD/water slightly decreases in the presence of BZDs, but the hydrogen bond's quality remains constant.
Collapse
Affiliation(s)
- Mokhtar Ganjali Koli
- InSilicoSci Computational Research Centre, Nikopardazesh Ltd., Karaj, Iran
- Department of Chemistry, University of Kurdistan, Sanandaj, Iran
| | | | - Hossein Hajiabadi
- InSilicoSci Computational Research Centre, Nikopardazesh Ltd., Karaj, Iran
| |
Collapse
|
8
|
Agles AA, Bourg IC. Structure-Thermodynamic Relationship of a Polysaccharide Gel (Alginate) as a Function of Water Content and Counterion Type (Na vs Ca). J Phys Chem B 2023; 127:1828-1841. [PMID: 36791328 PMCID: PMC10159261 DOI: 10.1021/acs.jpcb.2c07129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
Biofilms are the predominant mode of microbial life on Earth, and so a deep understanding of microbial communities─and their impacts on environmental processes─requires a firm understanding of biofilm properties. Because of the importance of biofilms to their microbial inhabitants, microbes have evolved different ways of engineering and reconfiguring the matrix of extracellular polymeric substances (EPS) that constitute the main non-living component of biofilms. This ability makes it difficult to distinguish between the biotic and abiotic origins of biofilm properties. An important route toward establishing this distinction has been the study of simplified models of the EPS matrix. This study builds on such efforts by using atomistic simulations to predict the nanoscale (≤10 nm scale) structure of a model EPS matrix and the sensitivity of this structure to interpolymer interactions and water content. To accomplish this, we use replica exchange molecular dynamics (REMD) simulations to generate all-atom configurations of ten 3.4 kDa alginate polymers at a range of water contents and Ca-Na ratios. Simulated systems are solvated with explicitly modeled water molecules, which allows us to capture the discrete structure of the hydrating water and to examine the thermodynamic stability of water in the gels as they are progressively dehydrated. Our primary findings are that (i) the structure of the hydrogels is highly sensitive to the identity of the charge-compensating cations, (ii) the thermodynamics of water within the gels (specific enthalpy and free energy) are, surprisingly, only weakly sensitive to cation identity, and (iii) predictions of the differential enthalpy and free energy of hydration include a short-ranged enthalpic term that promotes hydration and a longer-ranged (presumably entropic) term that promotes dehydration, where short and long ranges refer to distances shorter or longer than ∼0.6 nm between alginate strands.
Collapse
Affiliation(s)
- Avery A Agles
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Ian C Bourg
- Department of Civil and Environmental Engineering and High Meadows Environmental Institute, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
9
|
Boroushaki T, Dekamin MG. Interactions between β-cyclodextrin as a carrier for anti-cancer drug delivery: a molecular dynamics simulation study. J Biomol Struct Dyn 2023; 41:11551-11563. [PMID: 36597916 DOI: 10.1080/07391102.2022.2164058] [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: 09/28/2022] [Accepted: 12/24/2022] [Indexed: 01/05/2023]
Abstract
A series of molecular dynamics simulations were performed on 5-fluorouracil (5-Fu), Alendronate (Ald), and Temozolomide (TMZ) anticancer drugs in the presence and absence of β-cyclodextrin (βCD) as a carrier. Thermodynamic investigations showed that the van der Waals interaction energy was dominant in loading all drugs inside the βCD cavity. The sum of the interaction energies illustrated that the highest affinity was related to Ald (-136.5 kJ/mol), which in turn was due to the presence of bulky and charged atoms of phosphorus and oxygen, although TMZ (-115.92 kJ/mol) showed a very high affinity as well. At the same time, the hydrogen bond analysis also represented that Ald had the most hydrogen bond (1.97) with the highest half-life (3.13 ps) with βCD. Investigation of the root mean fluctuation (RMSF) indicated that all the drugs had a relatively rigid structure and maintain this rigidity during loading in the βCD cavity, and in the meantime, Ald was slightly more flexible than 5-Fu and TMZ. The area of the primary hydroxyl rim decreased in all drug-containing systems, which in turn was caused by the attractive interaction of drugs with oxygens in the primary hydroxyl rim. Especially for those drugs that were able to penetrate to the end of the primary hydroxyl rim of the βCD, that means TMZ and 5-Fu. Meanwhile, due to the lack of Ald penetration to the end of the primary hydroxyl rim, the area change in the Ald-containing system was less than in the two others.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Tahereh Boroushaki
- Pharmaceutical and Heterocyclic Compounds Research Laboratory, Department of Chemistry, Iran University of Science and Technology, Tehran, Iran
| | - Mohammad G Dekamin
- Pharmaceutical and Heterocyclic Compounds Research Laboratory, Department of Chemistry, Iran University of Science and Technology, Tehran, Iran
| |
Collapse
|
10
|
Liu X, Zheng L, Cong Y, Gong Z, Yin Z, Zhang JZH, Liu Z, Sun Z. Comprehensive evaluation of end-point free energy techniques in carboxylated-pillar[6]arene host-guest binding: II. regression and dielectric constant. J Comput Aided Mol Des 2022; 36:879-894. [PMID: 36394776 DOI: 10.1007/s10822-022-00487-w] [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: 08/26/2022] [Accepted: 10/29/2022] [Indexed: 11/18/2022]
Abstract
End-point free energy calculations as a powerful tool have been widely applied in protein-ligand and protein-protein interactions. It is often recognized that these end-point techniques serve as an option of intermediate accuracy and computational cost compared with more rigorous statistical mechanic models (e.g., alchemical transformation) and coarser molecular docking. However, it is observed that this intermediate level of accuracy does not hold in relatively simple and prototypical host-guest systems. Specifically, in our previous work investigating a set of carboxylated-pillar[6]arene host-guest complexes, end-point methods provide free energy estimates deviating significantly from the experimental reference, and the rank of binding affinities is also incorrectly computed. These observations suggest the unsuitability and inapplicability of standard end-point free energy techniques in host-guest systems, and alteration and development are required to make them practically usable. In this work, we consider two ways to improve the performance of end-point techniques. The first one is the PBSA_E regression that varies the weights of different free energy terms in the end-point calculation procedure, while the second one is considering the interior dielectric constant as an additional variable in the end-point equation. By detailed investigation of the calculation procedure and the simulation outcome, we prove that these two treatments (i.e., regression and dielectric constant) are manipulating the end-point equation in a somehow similar way, i.e., weakening the electrostatic contribution and strengthening the non-polar terms, although there are still many detailed differences between these two methods. With the trained end-point scheme, the RMSE of the computed affinities is improved from the standard ~ 12 kcal/mol to ~ 2.4 kcal/mol, which is comparable to another altered end-point method (ELIE) trained with system-specific data. By tuning PBSA_E weighting factors with the host-specific data, it is possible to further decrease the prediction error to ~ 2.1 kcal/mol. These observations along with the extremely efficient optimized-structure computation procedure suggest the regression (i.e., PBSA_E as well as its GBSA_E extension) as a practically applicable solution that brings end-point methods back into the library of usable tools for host-guest binding. However, the dielectric-constant-variable scheme cannot effectively minimize the experiment-calculation discrepancy for absolute binding affinities, but is able to improve the calculation of affinity ranks. This phenomenon is somehow different from the protein-ligand case and suggests the difference between host-guest and biomacromolecular (protein-ligand and protein-protein) systems. Therefore, the spectrum of tools usable for protein-ligand complexes could be unsuitable for host-guest binding, and numerical validations are necessary to screen out really workable solutions in these 'prototypical' situations.
Collapse
Affiliation(s)
- Xiao Liu
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China.
| | - Lei Zheng
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| | - Yalong Cong
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Zhihao Gong
- School of Micro-Nano Electronics, Zhejiang University, Hangzhou, 310027, China.,Hangzhou Global Scientific and Technological Innovation Center, Zhejiang University, Hangzhou, 310027, China
| | - Zhixiang Yin
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China
| | - John Z H Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China. .,School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China. .,Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong, China. .,Department of Chemistry, New York University, NY, NY, 10003, USA.
| | - Zhirong Liu
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China
| | - Zhaoxi Sun
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China.
| |
Collapse
|
11
|
Paulino M, Pérez-Juste I, Cid MM, Da Silva JP, Pereira MMA, Basílio N. 2-Hydroxychalcone-β-Cyclodextrin Conjugate with pH-Modulated Photoresponsive Binding Properties. J Org Chem 2022; 87:14422-14432. [PMID: 36242558 PMCID: PMC9776619 DOI: 10.1021/acs.joc.2c01875] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Stimuli-responsive supramolecular receptors are important building blocks for the construction of self-assembled functional materials. We report the design and synthesis of a pH- and light-responsive 2-hydroxychalcone-β-cyclodextrin conjugate (1-Ct) and its characterization by spectroscopic and computational methods. 1-Ct follows the typical reaction network of trans-chalcone-flavylium photoswitches. Upon light irradiation, 1-Ct can be photochemically converted into the cis-chalcone/hemiketal forms (1-Cc/1-B) under neutral pH conditions or to the flavylium cation (1-AH+) at acidic pH values. This stimuli-responsive β-cyclodextrin host, 1-Ct, was found to form stronger intramolecular self-inclusion complexes (Kintra = 14) than 1-AH+ (Kintra = 3) and weaker than 1-Cc/1-B (overall Kintra = 179), allowing control over their stability and binding properties by combinations of pH and light stimuli.
Collapse
Affiliation(s)
- Micael Paulino
- Laboratório
Associado para a Química Verde (LAQV), Rede de Química
e Tecnologia (REQUIMTE), Departamento de Química, Faculdade
de Ciências e Tecnologia, Universidade
NOVA de Lisboa, 2829-516Caparica, Portugal
| | - Ignacio Pérez-Juste
- Facultade
de Química, Edificio de Ciencias Experimentais, Universidade de Vigo, Campus Lagoas-Marcosende, 36310Vigo, Spain
| | - María Magdalena Cid
- Facultade
de Química, Edificio de Ciencias Experimentais, Universidade de Vigo, Campus Lagoas-Marcosende, 36310Vigo, Spain,
| | - José P. Da Silva
- Centre
of Marine Sciences (CCMAR/CIMAR LA), University
of Algarve, Campus de
Gambelas, 8005-139Faro, Portugal
| | - M. Manuela A. Pereira
- Laboratório
Associado para a Química Verde (LAQV), Rede de Química
e Tecnologia (REQUIMTE), Departamento de Química, Faculdade
de Ciências e Tecnologia, Universidade
NOVA de Lisboa, 2829-516Caparica, Portugal,
| | - Nuno Basílio
- Laboratório
Associado para a Química Verde (LAQV), Rede de Química
e Tecnologia (REQUIMTE), Departamento de Química, Faculdade
de Ciências e Tecnologia, Universidade
NOVA de Lisboa, 2829-516Caparica, Portugal,
| |
Collapse
|
12
|
Triolo A, Celso FL, Perez J, Russina O. Solubility and solvation features of native cyclodextrins in 1-ethyl-3-methylimidazolium acetate. Carbohydr Polym 2022; 291:119622. [DOI: 10.1016/j.carbpol.2022.119622] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Revised: 04/29/2022] [Accepted: 05/11/2022] [Indexed: 11/02/2022]
|
13
|
Sun Z, Zheng L, Wang K, Huai Z, Liu Z. Primary vs secondary: Directionalized guest coordination in β-cyclodextrin derivatives. Carbohydr Polym 2022; 297:120050. [DOI: 10.1016/j.carbpol.2022.120050] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 08/17/2022] [Accepted: 08/25/2022] [Indexed: 02/01/2023]
|
14
|
Abstract
The intrinsic nature of macrocyclic molecules to preferentially absorb a specific solute has been opening up supramolecular chemistry. Nevertheless, the determinant factor with molecular perspectives in promoting host-guest complexations remains inconclusive, due to the lack of rigorous thermodynamic examination on the guest solubility inside the host. Here, we quantify the solute-solvent energetic and entropic contributions between the end states and on the docking route during inclusion of noble gases in cucurbit[5]uril, cucurbit[6]uril, and α-cyclodextrin, using molecular dynamics simulations in combination with the potential distribution theorem. Results show that in all of the pairs examined both the solute-solvent energy and entropy favor the inclusion, while the former is rather dominant. The frequency of interior drying, which pertains to the entropic contribution, differs between the hosts and is controlled by the existence of lid water at portal and the flexibility of host framework. Moreover, the hosts exhibit various types of absorption manners, involving non-, single-, and double-free-energy barriers.
Collapse
Affiliation(s)
- Yifeng Yao
- Department of Chemistry, Zhejiang University, Hangzhou, 310028, P. R. China
| | - Xuan Zhang
- Department of Chemistry, Zhejiang University, Hangzhou, 310028, P. R. China
| | - Kenji Mochizuki
- Department of Chemistry, Zhejiang University, Hangzhou, 310028, P. R. China
| |
Collapse
|
15
|
Kraus H, Hansen N. An atomistic view on the uptake of aromatic compounds by cyclodextrin immobilized on mesoporous silica. ADSORPTION 2022. [DOI: 10.1007/s10450-022-00356-w] [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
AbstractThe effect of immobilized $$\upbeta$$
β
-cyclodextrin (bCD) molecules inside a mesoporous silica support on the uptake of benzene and p-nitrophenol from aqueous solution was investigated using all-atom molecular dynamics (MD) simulations. The calculated adsorption isotherms are discussed with respect to the free energies of binding for a 1:1 complex of bCD and the aromatic guest molecule. The adsorption capacity of the bCD-containing material significantly exceeds the amount corresponding to a 1:1 binding scenario, in agreement with experimental observations. Beside the formation of 1:2 and, to a lesser extent, 1:3 host:guest complexes, also host–host interactions on the surface as well as more unspecific host–guest interactions govern the adsorption process. The demonstrated feasibility of classical all-atom MD simulations to calculate liquid phase adsorption isotherms paves the way to a molecular interpretation of experimental data that are too complex to be described by empirical models.
Collapse
|
16
|
Lazar RD, Akher FB, Ravenscroft N, Kuttel MM. Carbohydrate Force Fields: The Role of Small Partial Atomic Charges in Preventing Conformational Collapse. J Chem Theory Comput 2022; 18:1156-1172. [PMID: 35015958 DOI: 10.1021/acs.jctc.1c00534] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Although the quality of current additive all-atom force fields for carbohydrates has been demonstrated in many applications, occasional significant differences reported for the hydrodynamic behavior of specific polysaccharides modeled with different force fields is a cause for concern. In particular, irreversible conformational collapse has been noted for some polysaccharide simulations with the GLYCAM06j force field. Here, we investigate the cause of this phenomenon through comparative simulations of a range of saccharides with both the GLYCAM06j and the CHARMM36 carbohydrate force fields. We find that conformational collapse in GLYCAM06j occurs for saccharide chains containing the deoxy sugar α-l-rhamnose after relatively long simulation intervals. Further, we explore the mechanism of conformational collapse and show that this phenomenon arises because of the anomalous low energy in GLYCAM06j (as compared to quantum mechanical calculations) of a specific orientation of α-l-Rha to α-l-Rha glycosidic linkages, which are subsequently sustained by intramolecular interactions in the saccharide chain. We identify the lack of partial charges on aliphatic hydrogens in GLYCAM as the source of this anomaly, demonstrating that addition of small partial atomic charges on the aliphatic protons in rhamnose removes the conformational collapse phenomenon. This work reveals the large cumulative impact that small partial charges may have on the dynamic behavior of polysaccharides and indicates that future reparameterization of the GLYCAM06j force field should investigate the addition of partial charges on all aliphatic hydrogens.
Collapse
Affiliation(s)
- Ryan D Lazar
- Department of Computer Science, University of Cape Town, Rondebosch, Cape Town 7701, South Africa
| | - Farideh B Akher
- Department of Computer Science, University of Cape Town, Rondebosch, Cape Town 7701, South Africa
| | - Neil Ravenscroft
- Department of Chemistry, University of Cape Town, Rondebosch, Cape Town 7701, South Africa
| | - Michelle M Kuttel
- Department of Computer Science, University of Cape Town, Rondebosch, Cape Town 7701, South Africa
| |
Collapse
|
17
|
Jitapunkul K, Toochinda P, Lawtrakul L. Molecular Dynamic Simulation Analysis on the Inclusion Complexation of Plumbagin with β-Cyclodextrin Derivatives in Aqueous Solution. Molecules 2021; 26:molecules26226784. [PMID: 34833875 PMCID: PMC8622708 DOI: 10.3390/molecules26226784] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2021] [Revised: 11/07/2021] [Accepted: 11/09/2021] [Indexed: 11/16/2022] Open
Abstract
Stable encapsulation of medically active compounds can lead to longer storage life and facilitate the slow-release mechanism. In this work, the dynamic and molecular interactions between plumbagin molecule with β-cyclodextrin (BCD) and its two derivatives, which are dimethyl-β-cyclodextrin (MBCD), and 2-O-monohydroxypropyl-β-cyclodextrin (HPBCD) were investigated. Molecular dynamics simulations (MD) with GLYCAM-06 and AMBER force fields were used to simulate the inclusion complex systems under storage temperature (4 °C) in an aqueous solution. The simulation results suggested that HPBCD is the best encapsulation agent to produce stable host-guest binding with plumbagin. Moreover, the observation of the plumbagin dynamic inside the binding cavity revealed that it tends to orient the methyl group toward the wider rim of HPBCD. Therefore, HPBCD is a decent candidate for the preservation of plumbagin with a promising longer storage life and presents the opportunity to facilitate the slow-release mechanism.
Collapse
|
18
|
Shadfar Z, Yahiaoui O, Collier TA, Fallon T, Allison JR. Illustration of a computational pipeline for evaluating cyclodextrin host-guest complex formation through conformational capture of bullvalene. J Chem Phys 2021; 154:154105. [PMID: 33887942 DOI: 10.1063/5.0045115] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Cyclodextrins have a diverse range of applications, including as supramolecular hosts, as enzyme active-site analogs, in improving drug solubility and delivery, and in molecular selection. We have investigated their ability to form stable complexes with bullvalenes, unusual organic cage molecules that spontaneously interconvert between numerous degenerate isomers. The shape-shifting nature of substituted bullvalenes raises the potential for dynamic adaptive binding to biological targets. We tested whether β- and γ-cyclodextrins can capture particular bullvalene isomers and whether the preferred binding mode(s) differ between isomers. We first applied our computational host-guest interaction potential energy profiling to determine the best binding mode(s) of unsubstituted bullvalene and each isomer of methylenehydroxybullvalene to β- and γ-cyclodextrin. Subsequent molecular dynamics simulations of the predicted host-guest complexes showed that while unsubstituted bullvalene has a single, albeit ill-defined, binding mode with either cyclodextrin, each isomer of methylenehydroxybullvalene has two possible modes of binding to β-cyclodextrin but only a single, nebulous mode of binding to γ-cyclodextrin. Experimental determination of the binding free energy of each methylenehydroxybullvalene-cyclodextrin complex showed that methylenehydroxybullvalene is more likely to bind to β-cyclodextrin than to γ-cyclodextrin, despite its smaller cavity. Together, our results suggest that β-cyclodextrin, but not γ-cyclodextrin, shows promise for conformational capture of mono-substituted bullvalenes. More broadly, our computational pipeline should prove useful for rapid characterization of cyclodextrin host-guest complexes, avoiding the need for costly synthesis of guest molecules that are unlikely to bind stably, as well as providing detailed atomic-level insight into the nature of complexation.
Collapse
Affiliation(s)
- Zahra Shadfar
- Centre for Theoretical Chemistry and Physics, Institute of Natural and Mathematical Sciences, Massey University, Auckland 0632, New Zealand
| | - Oussama Yahiaoui
- Department of Chemistry, School of Physical Sciences, University of Adelaide, Adelaide, SA 5005, Australia
| | - Thomas A Collier
- Centre for Theoretical Chemistry and Physics, Institute of Natural and Mathematical Sciences, Massey University, Auckland 0632, New Zealand
| | - Thomas Fallon
- Department of Chemistry, School of Physical Sciences, University of Adelaide, Adelaide, SA 5005, Australia
| | - Jane R Allison
- Maurice Wilkins Centre for Molecular Biodiscovery and School of Biological Sciences, University of Auckland, Auckland 1010, New Zealand
| |
Collapse
|
19
|
Arsiccio A, Rospiccio M, Shea JE, Pisano R. Force Field Parameterization for the Description of the Interactions between Hydroxypropyl-β-Cyclodextrin and Proteins. J Phys Chem B 2021; 125:7397-7405. [PMID: 34210121 PMCID: PMC8287564 DOI: 10.1021/acs.jpcb.1c04033] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Cyclodextrins are cyclic oligosaccharides, widely used as drug carriers, solubilizers, and excipients. Among cyclodextrins, the functionalized derivative known as hydroxypropyl-β-cyclodextrin (HPβCD) offers several advantages due to its unique structural features. Its optimal use in pharmaceutical and medical applications would benefit from a molecular-level understanding of its behavior, as can be offered by molecular dynamics simulations. Here, we propose a set of parameters for all-atom simulations of HPβCD, based on the ADD force field for sugars developed in our group, and compare it to the original CHARMM36 description. Using Kirkwood-Buff integrals of binary HPβCD-water mixtures as target experimental data, we show that the ADD-based description results in a considerably improved prediction of HPβCD self-association and interaction with water. We then use the new set of parameters to characterize the behavior of HPβCD toward the different amino acids. We observe pronounced interactions of HPβCD with both polar and nonpolar moieties, with a special preference for the aromatic rings of tyrosine, phenylalanine, and tryptophan. Interestingly, our simulations further highlight a preferential orientation of HPβCD's hydrophobic cavity toward the backbone atoms of amino acids, which, coupled with a favorable interaction of HPβCD with the peptide backbone, suggest a propensity for HPβCD to denature proteins.
Collapse
Affiliation(s)
- Andrea Arsiccio
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| | - Marcello Rospiccio
- Molecular Engineering Laboratory, Department of Applied Science and Technology, Politecnico di Torino, 24 corso Duca degli Abruzzi, Torino 10129, Italy
| | - Joan-Emma Shea
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States.,Department of Physics, University of California, Santa Barbara, California 93106, United States
| | - Roberto Pisano
- Molecular Engineering Laboratory, Department of Applied Science and Technology, Politecnico di Torino, 24 corso Duca degli Abruzzi, Torino 10129, Italy
| |
Collapse
|
20
|
Yan X, Wang Y, Meng T, Yan H. Computational Insights Into the Influence of Substitution Groups on the Inclusion Complexation of β-Cyclodextrin. Front Chem 2021; 9:668400. [PMID: 34095084 PMCID: PMC8176092 DOI: 10.3389/fchem.2021.668400] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Accepted: 05/03/2021] [Indexed: 11/21/2022] Open
Abstract
Cyclodextrins (CDs) and their derivatives have good prospects in soil remediation application due to their ability to enhance the stability and solubility of low water-soluble compounds by inclusion performance. To investigate the effect of different structural properties of cyclodextrin and its derivatives on the inclusion complexation, molecular dynamic (MD) simulations were performed on the inclusion complexes formed by three kinds of CDs with polycyclic aromatic hydrocarbons (PAHs). Based on neutral β-CD, the other two CDs were modified by introducing substitutional groups, including 2-hydroxypropyl and sulfonated butyl (SBE) functional groups in the ring structure, called HP-CD and SBE-CD. MD results show that PAH can merely enter into the cavity of SBE–β-CD from its wide rim. The substitutional groups significantly affect the structure of CDs, which may also cause the flipping of the glucose units. However, the substitutional groups can also enlarge the volume of the hydrophobic cavity, resulting in a tight combination with the guest molecules.
Collapse
Affiliation(s)
- Xianghua Yan
- School of Pharmaceutical Sciences, Liaocheng University, Liaocheng, China.,School of Chemistry and Chemical Engineering, Liaocheng University, Liaocheng, China
| | - Yue Wang
- School of Pharmaceutical Sciences, Liaocheng University, Liaocheng, China
| | - Tong Meng
- School of Pharmaceutical Sciences, Liaocheng University, Liaocheng, China
| | - Hui Yan
- School of Pharmaceutical Sciences, Liaocheng University, Liaocheng, China
| |
Collapse
|
21
|
Rospiccio M, Arsiccio A, Winter G, Pisano R. The Role of Cyclodextrins against Interface-Induced Denaturation in Pharmaceutical Formulations: A Molecular Dynamics Approach. Mol Pharm 2021; 18:2322-2333. [PMID: 33999634 PMCID: PMC8289300 DOI: 10.1021/acs.molpharmaceut.1c00135] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Protein-based pharmaceutical
products are subject to a variety
of environmental stressors, during both production and shelf-life.
In order to preserve their structure, and, therefore, functionality,
it is necessary to use excipients as stabilizing agents. Among the
eligible stabilizers, cyclodextrins (CDs) have recently gained interest
in the scientific community thanks to their properties. Here, a computational
approach is proposed to clarify the role of β-cyclodextrin (βCD)
and 2-hydroxypropyl-β-cyclodextrin (HPβCD) against granulocyte
colony-stimulating (GCSF) factor denaturation at the air–water
and ice–water interfaces, and also in bulk water at 300 or
260 K. Both traditional molecular dynamics (MD) simulations and enhanced
sampling techniques (metadynamics, MetaD) are used to shed light on
the underlying molecular mechanisms. Bulk simulations revealed that
CDs were preferentially included within the surface hydration layer
of GCSF, and even included some peptide residues in their hydrophobic
cavity. HPβCD was able to stabilize the protein against surface-induced
denaturation in proximity of the air–water interface, while
βCD had a destabilizing effect. No remarkable conformational
changes of GCSF, or noticeable effect of the CDs, were instead observed
at the ice surface. GCSF seemed less stable at low temperature (260
K), which may be attributed to cold-denaturation effects. In this
case, CDs did not significantly improve conformational stability.
In general, the conformationally altered regions of GCSF seemed not
to depend on the presence of excipients that only modulated the extent
of destabilization with either a positive or a negative effect.
Collapse
Affiliation(s)
- Marcello Rospiccio
- Molecular Engineering Laboratory, Department of Applied Science and Technology, Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino 10129, Italy
| | - Andrea Arsiccio
- Molecular Engineering Laboratory, Department of Applied Science and Technology, Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino 10129, Italy
| | - Gerhard Winter
- Department of Pharmacy, Ludwig-Maximilians-University, 81377 Munich, Germany
| | - Roberto Pisano
- Molecular Engineering Laboratory, Department of Applied Science and Technology, Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino 10129, Italy
| |
Collapse
|
22
|
Gieroba B, Kalisz G, Sroka-Bartnicka A, Płazińska A, Płaziński W, Starek M, Dąbrowska M. Molecular Structure of Cefuroxime Axetil Complexes with α-, β-, γ-, and 2-Hydroxypropyl-β-Cyclodextrins: Molecular Simulations and Raman Spectroscopic and Imaging Studies. Int J Mol Sci 2021; 22:ijms22105238. [PMID: 34063471 PMCID: PMC8156438 DOI: 10.3390/ijms22105238] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2021] [Revised: 05/11/2021] [Accepted: 05/13/2021] [Indexed: 11/16/2022] Open
Abstract
The formation of cefuroxime axetil+cyclodextrin (CA+CD) complexes increases the aqueous solubility of CA, improves its physico-chemical properties, and facilitates a biomembrane-mediated drug delivery process. In CD-based tablet formulations, it is crucial to investigate the molecular details of complexes in final pharmaceutical preparation. In this study, Raman spectroscopy and mapping were applied for the detection and identification of chemical groups involved in α-, β-, γ-, and 2-hydroxypropyl-β-CD (2-HP- β-CD)+CA complexation process. The experimental studies have been complemented by molecular dynamics-based investigations, providing additional molecular details of CA+CD interactions. It has been demonstrated that CA forms the guest–host type inclusion complexes with all studied CDs; however, the nature of the interactions is slightly different. It seems that both α- and β-CD interact with furanyl and methoxy moieties of CA, γ-CD forms a more diverse pattern of interactions with CA, which are not observed in other CDs, whereas 2HP-β-CD binds CA with the contribution of hydrogen bonding. Apart from supporting this interpretation of the experimental data, molecular dynamics simulations allowed for ordering the CA+CD binding affinities. The obtained results proved that the molecular details of the host–guest complexation can be successfully predicted from the combination of Raman spectroscopy and molecular modeling.
Collapse
Affiliation(s)
- Barbara Gieroba
- Department of Biopharmacy, Medical University of Lublin, ul. Chodzki 4a, 20-093 Lublin, Poland; (G.K.); (A.S.-B.); (A.P.)
- Correspondence: (B.G.); (W.P.)
| | - Grzegorz Kalisz
- Department of Biopharmacy, Medical University of Lublin, ul. Chodzki 4a, 20-093 Lublin, Poland; (G.K.); (A.S.-B.); (A.P.)
| | - Anna Sroka-Bartnicka
- Department of Biopharmacy, Medical University of Lublin, ul. Chodzki 4a, 20-093 Lublin, Poland; (G.K.); (A.S.-B.); (A.P.)
- Department of Genetics and Microbiology, Institute of Microbiology and Biotechnology, Maria Curie-Skłodowska University, ul. Akademicka 19, 20-033 Lublin, Poland
| | - Anita Płazińska
- Department of Biopharmacy, Medical University of Lublin, ul. Chodzki 4a, 20-093 Lublin, Poland; (G.K.); (A.S.-B.); (A.P.)
| | - Wojciech Płaziński
- Jerzy Haber Institute of Catalysis and Surface Chemistry, ul. Niezapominajek 8, 30-239 Krakow, Poland
- Correspondence: (B.G.); (W.P.)
| | - Małgorzata Starek
- Department of Inorganic and Analytical Chemistry, Faculty of Pharmacy, Jagiellonian University Medical College, ul. Medyczna 9, 30-688 Kraków, Poland; (M.S.); (M.D.)
| | - Monika Dąbrowska
- Department of Inorganic and Analytical Chemistry, Faculty of Pharmacy, Jagiellonian University Medical College, ul. Medyczna 9, 30-688 Kraków, Poland; (M.S.); (M.D.)
| |
Collapse
|
23
|
Stasiłowicz A, Rosiak N, Tykarska E, Kozak M, Jenczyk J, Szulc P, Kobus-Cisowska J, Lewandowska K, Płazińska A, Płaziński W, Cielecka-Piontek J. Combinations of Piperine with Hydroxypropyl-β-Cyclodextrin as a Multifunctional System. Int J Mol Sci 2021; 22:4195. [PMID: 33919582 PMCID: PMC8072981 DOI: 10.3390/ijms22084195] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2021] [Revised: 04/09/2021] [Accepted: 04/16/2021] [Indexed: 12/12/2022] Open
Abstract
Piperine is an alkaloid that has extensive pharmacological activity and impacts other active substances bioavailability due to inhibition of CYP450 enzymes, stimulation of amino acid transporters and P-glycoprotein inhibition. Low solubility and the associated low bioavailability of piperine limit its potential. The combination of piperine with 2-hydroxypropyl-β-cyclodextrin (HP-β-CD) causes a significant increase in its solubility and, consequently, an increase in permeability through gastrointestinal tract membranes and the blood-brain barrier. X-ray powder diffraction (XRPD), differential scanning calorimetry (DSC), Fourier-transform infrared spectroscopy (FT-IR), nuclear magnetic resonance (NMR) were used to characterize interactions between piperine and HP-β-CD. The observed physicochemical changes should be combined with the process of piperine and CD system formation. Importantly, with an increase in solubility and permeability of piperine as a result of interaction with CD, it was proven to maintain its biological activity concerning the antioxidant potential (2,2-diphenyl-1-picryl-hydrazyl-hydrate assay), inhibition of enzymes essential for the inflammatory process and for neurodegenerative changes (hyaluronidase, acetylcholinesterase, butyrylcholinesterase).
Collapse
Affiliation(s)
- Anna Stasiłowicz
- Department of Pharmacognosy, Faculty of Pharmacy, Poznan University of Medical Sciences, Swiecickiego 4, 60-781 Poznań, Poland; (A.S.); (N.R.)
| | - Natalia Rosiak
- Department of Pharmacognosy, Faculty of Pharmacy, Poznan University of Medical Sciences, Swiecickiego 4, 60-781 Poznań, Poland; (A.S.); (N.R.)
| | - Ewa Tykarska
- Department of Chemical Technology of Drugs, Poznan University of Medical Sciences, Grunwaldzka 6, 60-780 Poznan, Poland;
| | - Maciej Kozak
- Department of Macromolecular Physics, Faculty of Physics, Adam Mickiewicz University, Uniwersytetu Poznańskiego 2, 61-614 Poznań, Poland;
| | - Jacek Jenczyk
- NanoBioMedical Centre, Adam Mickiewicz University, Wszechnicy Piastowskiej 3, 61-614 Poznań, Poland;
| | - Piotr Szulc
- Department of Agronomy, Poznań University of Life Sciences, Dojazd 11, 60-632 Poznań, Poland;
| | - Joanna Kobus-Cisowska
- Department of Gastronomy Sciences and Functional Foods, Faculty of Food Science and Nutrition, Poznań University of Life Sciences, Wojska Polskiego 31, 60-624 Poznań, Poland;
| | - Kornelia Lewandowska
- Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego 17, 60-179 Poznan, Poland;
| | - Anita Płazińska
- Department of Biopharmacy, Faculty of Pharmacy, Medical University of Lublin, Chodzki 4a, 20-093 Lublin, Poland;
| | - Wojciech Płaziński
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, 30-239 Krakow, Poland;
| | - Judyta Cielecka-Piontek
- Department of Pharmacognosy, Faculty of Pharmacy, Poznan University of Medical Sciences, Swiecickiego 4, 60-781 Poznań, Poland; (A.S.); (N.R.)
| |
Collapse
|
24
|
Erdős M, Frangou M, Vlugt TJH, Moultos OA. Diffusivity of α-, β-, γ-cyclodextrin and the inclusion complex of β-cyclodextrin: Ibuprofen in aqueous solutions; A molecular dynamics simulation study. FLUID PHASE EQUILIBRIA 2021; 528:112842. [PMID: 33024350 PMCID: PMC7529625 DOI: 10.1016/j.fluid.2020.112842] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2020] [Revised: 09/03/2020] [Accepted: 09/18/2020] [Indexed: 05/28/2023]
Abstract
Cyclodextrins (CDs) are widely used in drug delivery, catalysis, food and separation processes. In this work, a comprehensive simulation study on the diffusion of the native α-, β- and γ-CDs in aqueous solutions is carried out using Molecular Dynamics simulations. The effect of the system size on the computed self-diffusivity is investigated and it is found that the required correction can be as much as 75% of the final value. The effect of the water force field is examined and it is shown that the q4md-CD/TIP4P/2005 force field combination predicts the experimentally measured self-diffusion coefficients of CDs very accurately. The self-diffusion coefficients of the three native CDs were also computed in aqueous-NaCl solutions using the Joung and Cheatham (JC) and the Madrid-2019 force fields. It is found that Na+ ions have higher affinity towards the CDs when the JC force field is used and for this reason the predicted diffusivity of CDs is lower compared to simulations using the Madrid-2019 force field. As a model system for drug delivery and waste-water treatment applications, the diffusion of the β-CD:Ibuprofen inclusion complex in water is studied. In agreement with experiments for similar components, it is shown that the inclusion complex and the free β-CD have almost equal self-diffusion coefficients. Our analysis revealed that this is most likely caused by the almost full inclusion of the ibuprofen in the cavity of the β-CD. Our findings show that Molecular Dynamics simulation can be used to provide reasonable diffusivity predictions, and to obtain molecular-level understanding useful for industrial applications of CDs.
Collapse
Affiliation(s)
- Máté Erdős
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - Michalis Frangou
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - Thijs J H Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| | - Othonas A Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, Delft 2628CB, The Netherlands
| |
Collapse
|
25
|
Lewis-Atwell T, Townsend PA, Grayson MN. Comparisons of different force fields in conformational analysis and searching of organic molecules: A review. Tetrahedron 2021. [DOI: 10.1016/j.tet.2020.131865] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
26
|
Tedizolid-Cyclodextrin System as Delayed-Release Drug Delivery with Antibacterial Activity. Int J Mol Sci 2020; 22:ijms22010115. [PMID: 33374358 PMCID: PMC7795824 DOI: 10.3390/ijms22010115] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2020] [Revised: 12/21/2020] [Accepted: 12/22/2020] [Indexed: 12/18/2022] Open
Abstract
Progressive increase in bacterial resistance has caused an urgent need to introduce new antibiotics, one of them being oxazolidinones with their representative tedizolid. Despite the broad spectrum of activity of the parent tedizolid, it is characterized by low water solubility, which limits its use. The combination of the active molecule with a multifunctional excipient, which is cyclodextrins, allows preservation of its pharmacological activity and modification of its physicochemical properties. Therefore, the aim of the study was to change the dissolution rate and permeability through the model membrane of tedizolid by formation of solid dispersions with a cyclodextrin. The research included identification of tedizolid-hydroxypropyl-β-cyclodextrin (tedizolid/HP-β-CD) inclusion complex by thermal method (Differential Scanning Colorimetry), spectroscopic methods (powder X-ray diffraction, Fourier-Transform Infrared spectroscopy), and molecular docking. The second part of the research concerned the physicochemical properties (dissolution and permeability) and the biological properties of the system in terms of its microbiological activity. An increase in the dissolution rate was observed in the presence of cyclodextrin, while maintaining a high permeation coefficient and high microbiological activity. The proposed approach is an opportunity to develop drug delivery systems used in the treatment of resistant bacterial infections, in which, in addition to modifying the physicochemical properties caused by cyclodextrin, we observe a favorable change in the pharmacological potential of the bioactives.
Collapse
|
27
|
Sandilya A, Natarajan U, Priya MH. Molecular View into the Cyclodextrin Cavity: Structure and Hydration. ACS OMEGA 2020; 5:25655-25667. [PMID: 33073091 PMCID: PMC7557249 DOI: 10.1021/acsomega.0c02760] [Citation(s) in RCA: 62] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/11/2020] [Accepted: 08/14/2020] [Indexed: 05/21/2023]
Abstract
We find, through atomistic molecular dynamics simulation of native cyclodextrins (CDs) in water, that although the outer surface of a CD appears like a truncated cone, the inner cavity resembles a conical hourglass because of the inward protrusion of the glycosidic oxygens. Furthermore, the conformations of the constituent α-glucose molecules are found to differ significantly from a free monomeric α-glucose molecule. This is the first computational study that maps the conformational change to the preferential hydrogen bond donating capacity of one of the secondary hydroxyl groups of CD, in consensus with an NMR experiment. We have developed a simple and novel geometry-based technique to identify water molecules occupying the nonspherical CD cavity, and the computed water occupancies are in close agreement with the experimental and density functional theory studies. Our analysis reveals that a water molecule in CD cavity loses out about two hydrogen bonds and remains energetically frustrated but possesses higher orientational degree of freedom compared to bulk water. In the context of CD-drug complexation, these imply a nonclassical, that is, enthalpically driven hydrophobic association of a drug in CD cavity.
Collapse
Affiliation(s)
- Avilasha
A. Sandilya
- Department
of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Upendra Natarajan
- Department
of Chemical Engineering, Indian Institute
of Technology Madras, Chennai 600036, India
| | - M. Hamsa Priya
- Department
of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
- . Phone: +91-44-22574132
| |
Collapse
|
28
|
Triolo A, Lo Celso F, Russina O. Structural Features of β-Cyclodextrin Solvation in the Deep Eutectic Solvent, Reline. J Phys Chem B 2020; 124:2652-2660. [DOI: 10.1021/acs.jpcb.0c00876] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Alessandro Triolo
- Laboratorio Liquidi Ionici, Istituto Struttura della Materia, Consiglio Nazionale delle Ricerche (ISM-CNR), 00133 Rome, Italy
| | - Fabrizio Lo Celso
- Laboratorio Liquidi Ionici, Istituto Struttura della Materia, Consiglio Nazionale delle Ricerche (ISM-CNR), 00133 Rome, Italy
- Dipartimento di Fisica e Chimica ‘Emilio Segrè’, Università degli studi di Palermo, 90128 Palermo, Italy
| | - Olga Russina
- Laboratorio Liquidi Ionici, Istituto Struttura della Materia, Consiglio Nazionale delle Ricerche (ISM-CNR), 00133 Rome, Italy
- Department of Chemistry, University of Rome Sapienza, 00185 Rome, Italy
| |
Collapse
|
29
|
Turupcu A, Blaukopf M, Kosma P, Oostenbrink C. Molecular Conformations of Di-, Tri-, and Tetra- α-(2→8)-Linked Sialic Acid from NMR Spectroscopy and MD Simulations. Int J Mol Sci 2019; 21:ijms21010030. [PMID: 31861593 PMCID: PMC6981865 DOI: 10.3390/ijms21010030] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2019] [Revised: 12/12/2019] [Accepted: 12/14/2019] [Indexed: 12/26/2022] Open
Abstract
By using molecular dynamics simulations with an efficient enhanced sampling technique and in combination with nuclear magnetic resonance (NMR) spectroscopy quantitative structural information on α-2,8-linked sialic acids is presented. We used a bottom-up approach to obtain a set of larger ensembles for tetra- and deca-sialic acid from model dimer and trimer systems that are in agreement with the available J-coupling constants and nuclear Overhauser effects. The molecular dynamic (MD) simulations with enhanced sampling are used to validate the force field used in this study for its further use. This empowered us to couple NMR observables in the MD framework via J-coupling and distance restraining simulations to obtain conformations that are supported by experimental data. We used these conformations in thermodynamic integration and one-step perturbation simulations to calculate the free-energy of suggested helical conformations. This study brings most of the available NMR experiments together and supplies information to resolve the conflict on the structures of poly-α-2,8-linked sialic acid.
Collapse
Affiliation(s)
- Aysegül Turupcu
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences Vienna, 1190 Vienna, Austria;
| | - Markus Blaukopf
- Department of Chemistry, University of Natural Resources and Life Sciences Vienna, 1190 Vienna, Austria; (M.B.); (P.K.)
| | - Paul Kosma
- Department of Chemistry, University of Natural Resources and Life Sciences Vienna, 1190 Vienna, Austria; (M.B.); (P.K.)
| | - Chris Oostenbrink
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences Vienna, 1190 Vienna, Austria;
- Correspondence: ; Tel.: +43-1-47654-89411
| |
Collapse
|
30
|
Wu J, Ma H, Bu X, Ma C, Zhu L, Hao B, Zhao B, Tian Y. SERS determination of the antihypertensive drugs prazosin and losartan by using silver nanoparticles coated with β-cyclodextrin. Mikrochim Acta 2019; 186:801. [DOI: 10.1007/s00604-019-3946-z] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2019] [Accepted: 10/14/2019] [Indexed: 01/06/2023]
|
31
|
Slochower DR, Henriksen NM, Wang LP, Chodera JD, Mobley DL, Gilson MK. Binding Thermodynamics of Host-Guest Systems with SMIRNOFF99Frosst 1.0.5 from the Open Force Field Initiative. J Chem Theory Comput 2019; 15:6225-6242. [PMID: 31603667 PMCID: PMC7328435 DOI: 10.1021/acs.jctc.9b00748] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Designing ligands that bind their target biomolecules with high affinity and specificity is a key step in small-molecule drug discovery, but accurately predicting protein-ligand binding free energies remains challenging. Key sources of errors in the calculations include inadequate sampling of conformational space, ambiguous protonation states, and errors in force fields. Noncovalent complexes between a host molecule with a binding cavity and a druglike guest molecule have emerged as powerful model systems. As model systems, host-guest complexes reduce many of the errors in more complex protein-ligand binding systems, as their small size greatly facilitates conformational sampling, and one can choose systems that avoid ambiguities in protonation states. These features, combined with their ease of experimental characterization, make host-guest systems ideal model systems to test and ultimately optimize force fields in the context of binding thermodynamics calculations. The Open Force Field Initiative aims to create a modern, open software infrastructure for automatically generating and assessing force fields using data sets. The first force field to arise out of this effort, named SMIRNOFF99Frosst, has approximately one tenth the number of parameters, in version 1.0.5, compared to typical general small molecule force fields, such as GAFF. Here, we evaluate the accuracy of this initial force field, using free energy calculations of 43 α and β-cyclodextrin host-guest pairs for which experimental thermodynamic data are available, and compare with matched calculations using two versions of GAFF. For all three force fields, we used TIP3P water and AM1-BCC charges. The calculations are performed using the attach-pull-release (APR) method as implemented in the open source package, pAPRika. For binding free energies, the root-mean-square error of the SMIRNOFF99Frosst calculations relative to experiment is 0.9 [0.7, 1.1] kcal/mol, while the corresponding results for GAFF 1.7 and GAFF 2.1 are 0.9 [0.7, 1.1] kcal/mol and 1.7 [1.5, 1.9] kcal/mol, respectively, with 95% confidence ranges in brackets. These results suggest that SMIRNOFF99Frosst performs competitively with existing small molecule force fields and is a parsimonious starting point for optimization.
Collapse
Affiliation(s)
- David R Slochower
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California, San Diego , La Jolla , California 92093 , United States
| | - Niel M Henriksen
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California, San Diego , La Jolla , California 92093 , United States
| | - Lee-Ping Wang
- Department of Chemistry , University of California , Davis , California 95616 , United States
| | - John D Chodera
- Computational and Systems Biology Program, Sloan Kettering Institute , Memorial Sloan Kettering Cancer Center , New York , New York 10065 , United States
| | - David L Mobley
- Department of Pharmaceutical Sciences and Department of Chemistry , University of California , Irvine , California 92697 , United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences , University of California, San Diego , La Jolla , California 92093 , United States
| |
Collapse
|
32
|
He P, Sarkar S, Gallicchio E, Kurtzman T, Wickstrom L. Role of Displacing Confined Solvent in the Conformational Equilibrium of β-Cyclodextrin. J Phys Chem B 2019; 123:8378-8386. [PMID: 31509409 DOI: 10.1021/acs.jpcb.9b07028] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This study investigates the role of hydration and its relationship to the conformational equilibrium of the host molecule β-cyclodextrin. Molecular dynamics simulations indicate that the unbound β-cyclodextrin exhibits two state behavior in explicit solvent due to the opening and closing of its cavity. In implicit solvent, these transitions are not observed, and there is one dominant conformation of β-cyclodextrin with an open cavity. Based on these observations, we investigate the hypothesis that the expulsion of thermodynamically unfavorable water molecules into the bulk plays an important role in controlling the accessibility of the closed macrostate at room temperature. We compare the results of the molecular mechanics analytical generalized Born plus nonpolar solvation approach to those obtained through grid inhomogeneous solvation theory analysis with explicit solvation to elucidate the thermodynamic forces at play. The work illustrates the use of continuum solvent models to tease out solvation effects related to the inhomogeneity and the molecular nature of water and demonstrates the key role of the thermodynamics of enclosed hydration in driving the conformational equilibrium of molecules in solution.
Collapse
Affiliation(s)
- Peng He
- Center for Biophysics & Computational Biology/ICMS, Department of Chemistry , Temple University , Philadelphia , Pennsylvania 19122 , United States
| | - Sheila Sarkar
- Department of Science , Borough of Manhattan Community College, The City University of New York , New York , New York 10007 , United States
| | - Emilio Gallicchio
- Department of Chemistry , Brooklyn College, The City University of New York , Brooklyn , New York 11210 , United States.,Ph.D. Programs in Chemistry & Biochemistry , The Graduate Center of the City University of New York , 365 Fifth Avenue , New York , New York 10016 , United States
| | - Tom Kurtzman
- Department of Chemistry , Lehman College, The City University of New York , Bronx , New York 10468 , United States.,Ph.D. Programs in Chemistry & Biochemistry , The Graduate Center of the City University of New York , 365 Fifth Avenue , New York , New York 10016 , United States
| | - Lauren Wickstrom
- Department of Science , Borough of Manhattan Community College, The City University of New York , New York , New York 10007 , United States
| |
Collapse
|
33
|
Turupcu A, Diem M, Smith LJ, Oostenbrink C. Structural Aspects of the O-glycosylation Linkage in Glycopeptides via MD Simulations and Comparison with NMR Experiments. Chemphyschem 2019; 20:1527-1537. [PMID: 30920077 PMCID: PMC6563056 DOI: 10.1002/cphc.201900079] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2019] [Revised: 03/27/2019] [Indexed: 11/26/2022]
Abstract
A powerful conformational searching and enhanced sampling simulation method, and unbiased molecular dynamics simulations have been used along with NMR spectroscopic observables to provide a detailed structural view of O-glycosylation. For four model systems, the force-field parameters can accurately predict experimental NMR observables (J couplings and NOE's). This enables us to derive conclusions based on the generated ensembles, in which O-glycosylation affects the peptide backbone conformation by forcing it towards to an extended conformation. An exception is described for β-GalNAc-Thr where the α content is increased and stabilized via hydrogen bonding between the sugar and the peptide backbone, which was not observed in the rest of the studied systems. These observations might offer an explanation for the evolutionary preference of α-linked GalNAc glycosylation instead of a β link.
Collapse
Affiliation(s)
- Aysegül Turupcu
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| | - Matthias Diem
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| | | | - Chris Oostenbrink
- Department of Material Sciences and Process Engineering, Institute of Molecular Modeling and SimulationUniversity of Natural Resources and Life Sciences ViennaViennaAustria
| |
Collapse
|
34
|
Verkhnyatskaya SA, de Vries AH, Douma‐de Vries E, Sneep RJL, Walvoort MTC. Direct and Regioselective Di-α-fucosylation on the Secondary Rim of β-Cyclodextrin. Chemistry 2019; 25:6722-6727. [PMID: 30801812 PMCID: PMC6563713 DOI: 10.1002/chem.201806090] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2018] [Revised: 02/15/2019] [Indexed: 01/01/2023]
Abstract
A straightforward glycosylation method is described to regio- and stereoselectively introduce two α-l-fucose moieties directly to the secondary rim of β-cyclodextrin. Using NMR and MS fragmentation studies, the nonasaccharide structure was determined, which was also visualized using molecular dynamics simulations. The reported glycosylation method proved to be robust on gram-scale, and may be generally applied to directly glycosylate β-cyclodextrins to make well-defined multivalent glycoclusters.
Collapse
Affiliation(s)
- Stella A. Verkhnyatskaya
- Stratingh Institute for ChemistryUniversity of GroningenNijenborgh 79747 AGGroningenThe Netherlands
| | - Alex H. de Vries
- Groningen Biomolecular Sciences and Biotechnology InstituteUniversity of GroningenNijenborgh 79747 AGGroningenThe Netherlands
| | - Elmatine Douma‐de Vries
- Stratingh Institute for ChemistryUniversity of GroningenNijenborgh 79747 AGGroningenThe Netherlands
| | - Renze J. L. Sneep
- Stratingh Institute for ChemistryUniversity of GroningenNijenborgh 79747 AGGroningenThe Netherlands
| | - Marthe T. C. Walvoort
- Stratingh Institute for ChemistryUniversity of GroningenNijenborgh 79747 AGGroningenThe Netherlands
| |
Collapse
|
35
|
Ramos MC, Horta VAC, Horta BAC. Molecular Dynamics Simulations of PAMAM and PPI Dendrimers Using the GROMOS-Compatible 2016H66 Force Field. J Chem Inf Model 2019; 59:1444-1457. [DOI: 10.1021/acs.jcim.8b00911] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Mayk C. Ramos
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, 21941-909, Brazil
| | - Vitor A. C. Horta
- Departamento de Ciência da Computação, Universidade Federal de Juiz de Fora, Juiz de Fora, 36036-900, Brazil
| | - Bruno A. C. Horta
- Instituto de Química, Universidade Federal do Rio de Janeiro, Rio de Janeiro, 21941-909, Brazil
| |
Collapse
|
36
|
Baz J, Gebhardt J, Kraus H, Markthaler D, Hansen N. Insights into Noncovalent Binding Obtained from Molecular Dynamics Simulations. CHEM-ING-TECH 2018. [DOI: 10.1002/cite.201800050] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Jörg Baz
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Julia Gebhardt
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Hamzeh Kraus
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Daniel Markthaler
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| | - Niels Hansen
- University of Stuttgart; Institute of Thermodynamics and Thermal Process Engineering; Pfaffenwaldring 9 70569 Stuttgart Germany
| |
Collapse
|