{"title":"Energetic and Non-covalent Interaction Insights in Host–Guest Chemistry: A Case Study of Metronidazole-β-Cyclodextrin Complex","authors":"Dhouha Boulahraf, Amel Bendjeddou, Tahar Abbaz, Didier Villemin","doi":"10.1007/s10953-026-01563-7","DOIUrl":null,"url":null,"abstract":"<div><p>Metronidazole (MNZ) displays notable biological activity; however, it is encumbered by suboptimal physicochemical and pharmacokinetic properties. To address these limitations, researchers have explored the formation of inclusion complexes with cyclodextrins. The present study aims to elucidate the inclusion mechanism of MNZ and β-cyclodextrin (β-CD) in a 1:1 stoichiometry, spanning from electronic structure to non-covalent interactions. Two orientations of MNZ were evaluated: in orientation <i>A</i>, MNZ enters through the wide rim; in orientation <i>B</i>, through the narrow rim. Semi-empirical and density functional theory (DFT) calculations, including dispersion-correction DFT (DFT–D3), were performed in gas and water phases using PW6B95, PW6B95–D3, and B3LYP functionals in conjunction with the 6-31G(d) basis set. The study systematically compared the performance of these methods, with a particular focus on the impact of dispersion correction on the accuracy and reliability of the results. The findings obtained from this study indicate that orientation <i>A</i> is more favorable than <i>B</i>. In both orientations, MNZ is fully encapsulated in the β-CD cavity. The Gauge-Including Atomic Orbital (GIAO) method was employed to calculate the shifts of the proton nuclear magnetic resonance (<sup>1</sup>H NMR), which aligned well with the experimental data. To further analyze host–guest interactions, advanced computational tools: AIM, NBO, NCI–RDG, and IGM were utilized. These analyses revealed that the stability of the MNZ@β-CD inclusion complex is driven by hydrogen bonding and van der Waals interactions, with dispersion corrections (DFT–D3) playing a critical role in accurately modeling these non-covalent forces. This study provides an exhaustive analysis of the molecular interactions driving the formation of MNZ@β-CD complexes, emphasizing the reliability of the PW6B95 and PW6B95–D3 functionals and the necessity of incorporating dispersion corrections for an accurate description of these systems.</p></div>","PeriodicalId":666,"journal":{"name":"Journal of Solution Chemistry","volume":"55 6","pages":"1246 - 1269"},"PeriodicalIF":1.3000,"publicationDate":"2026-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Solution Chemistry","FirstCategoryId":"92","ListUrlMain":"https://link.springer.com/article/10.1007/s10953-026-01563-7","RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q4","JCRName":"CHEMISTRY, PHYSICAL","Score":null,"Total":0}
引用次数: 0
Abstract
Metronidazole (MNZ) displays notable biological activity; however, it is encumbered by suboptimal physicochemical and pharmacokinetic properties. To address these limitations, researchers have explored the formation of inclusion complexes with cyclodextrins. The present study aims to elucidate the inclusion mechanism of MNZ and β-cyclodextrin (β-CD) in a 1:1 stoichiometry, spanning from electronic structure to non-covalent interactions. Two orientations of MNZ were evaluated: in orientation A, MNZ enters through the wide rim; in orientation B, through the narrow rim. Semi-empirical and density functional theory (DFT) calculations, including dispersion-correction DFT (DFT–D3), were performed in gas and water phases using PW6B95, PW6B95–D3, and B3LYP functionals in conjunction with the 6-31G(d) basis set. The study systematically compared the performance of these methods, with a particular focus on the impact of dispersion correction on the accuracy and reliability of the results. The findings obtained from this study indicate that orientation A is more favorable than B. In both orientations, MNZ is fully encapsulated in the β-CD cavity. The Gauge-Including Atomic Orbital (GIAO) method was employed to calculate the shifts of the proton nuclear magnetic resonance (1H NMR), which aligned well with the experimental data. To further analyze host–guest interactions, advanced computational tools: AIM, NBO, NCI–RDG, and IGM were utilized. These analyses revealed that the stability of the MNZ@β-CD inclusion complex is driven by hydrogen bonding and van der Waals interactions, with dispersion corrections (DFT–D3) playing a critical role in accurately modeling these non-covalent forces. This study provides an exhaustive analysis of the molecular interactions driving the formation of MNZ@β-CD complexes, emphasizing the reliability of the PW6B95 and PW6B95–D3 functionals and the necessity of incorporating dispersion corrections for an accurate description of these systems.
期刊介绍:
Journal of Solution Chemistry offers a forum for research on the physical chemistry of liquid solutions in such fields as physical chemistry, chemical physics, molecular biology, statistical mechanics, biochemistry, and biophysics. The emphasis is on papers in which the solvent plays a dominant rather than incidental role. Featured topics include experimental investigations of the dielectric, spectroscopic, thermodynamic, transport, or relaxation properties of both electrolytes and nonelectrolytes in liquid solutions.