The β-cyclodextrin/benzene complex and its hydrogen bonds – a theoretical study using molecular dynamics, quantum mechanics and COSMO-RS
Four highly ordered hydrogen-bonded models of β-cyclodextrin (β-CD) and its inclusion complex with benzene were investigated by three different theoretical methods: classical quantum mechanics (QM) on AM1 and on the BP/TZVP-DISP3 level of approximation, and thirdly by classical molecular dynamics si...
Main Authors: | , |
---|---|
Format: | Article |
Language: | English |
Published: |
Beilstein-Institut
2013-01-01
|
Series: | Beilstein Journal of Organic Chemistry |
Subjects: | |
Online Access: | https://doi.org/10.3762/bjoc.9.15 |
_version_ | 1818655414038822912 |
---|---|
author | Jutta Erika Helga Köhler Nicole Grczelschak-Mick |
author_facet | Jutta Erika Helga Köhler Nicole Grczelschak-Mick |
author_sort | Jutta Erika Helga Köhler |
collection | DOAJ |
description | Four highly ordered hydrogen-bonded models of β-cyclodextrin (β-CD) and its inclusion complex with benzene were investigated by three different theoretical methods: classical quantum mechanics (QM) on AM1 and on the BP/TZVP-DISP3 level of approximation, and thirdly by classical molecular dynamics simulations (MD) at different temperatures (120 K and 273 to 300 K). The hydrogen bonds at the larger O2/O3 rim of empty β-CDs prefer the right-hand orientation, e.g., O3-H…O2-H in the same glucose unit and bifurcated towards …O4 and O3 of the next glucose unit on the right side. On AM1 level the complex energy was −2.75 kcal mol−1 when the benzene molecule was located parallel inside the β-CD cavity and −2.46 kcal mol−1 when it was positioned vertically. The AM1 HOMO/LUMO gap of the empty β-CD with about 12 eV is lowered to about 10 eV in the complex, in agreement with data from the literature. AM1 IR spectra displayed a splitting of the O–H frequencies of cyclodextrin upon complex formation. At the BP/TZVP-DISP3 level the parallel and vertical positions from the starting structures converged to a structure where benzene assumes a more oblique position (−20.16 kcal mol−1 and −20.22 kcal mol−1, resp.) as was reported in the literature. The character of the COSMO-RS σ-surface of β-CD was much more hydrophobic on its O6 rim than on its O2/O3 side when all hydrogen bonds were arranged in a concerted mode.This static QM picture of the β-CD/benzene complex at 0 K was extended by MD simulations. At 120 K benzene was mobile but always stayed inside the cavity of β-CD. The trajectories at 273, 280, 290 and 300 K certainly no longer displayed the highly ordered hydrogen bonds of β-CD and benzene occupied many different positions inside the cavity, before it left the β-CD finally at its O2/O3 side. |
first_indexed | 2024-12-17T03:09:18Z |
format | Article |
id | doaj.art-f7174d2442084851962c32b2bc6009f7 |
institution | Directory Open Access Journal |
issn | 1860-5397 |
language | English |
last_indexed | 2024-12-17T03:09:18Z |
publishDate | 2013-01-01 |
publisher | Beilstein-Institut |
record_format | Article |
series | Beilstein Journal of Organic Chemistry |
spelling | doaj.art-f7174d2442084851962c32b2bc6009f72022-12-21T22:05:53ZengBeilstein-InstitutBeilstein Journal of Organic Chemistry1860-53972013-01-019111813410.3762/bjoc.9.151860-5397-9-15The β-cyclodextrin/benzene complex and its hydrogen bonds – a theoretical study using molecular dynamics, quantum mechanics and COSMO-RSJutta Erika Helga Köhler0Nicole Grczelschak-Mick1Wacker Chemie AG, Consortium für elektrochemische Industrie, Zielstattstrasse 20, D-81379 München, GermanyWacker Chemie AG, Consortium für elektrochemische Industrie, Zielstattstrasse 20, D-81379 München, GermanyFour highly ordered hydrogen-bonded models of β-cyclodextrin (β-CD) and its inclusion complex with benzene were investigated by three different theoretical methods: classical quantum mechanics (QM) on AM1 and on the BP/TZVP-DISP3 level of approximation, and thirdly by classical molecular dynamics simulations (MD) at different temperatures (120 K and 273 to 300 K). The hydrogen bonds at the larger O2/O3 rim of empty β-CDs prefer the right-hand orientation, e.g., O3-H…O2-H in the same glucose unit and bifurcated towards …O4 and O3 of the next glucose unit on the right side. On AM1 level the complex energy was −2.75 kcal mol−1 when the benzene molecule was located parallel inside the β-CD cavity and −2.46 kcal mol−1 when it was positioned vertically. The AM1 HOMO/LUMO gap of the empty β-CD with about 12 eV is lowered to about 10 eV in the complex, in agreement with data from the literature. AM1 IR spectra displayed a splitting of the O–H frequencies of cyclodextrin upon complex formation. At the BP/TZVP-DISP3 level the parallel and vertical positions from the starting structures converged to a structure where benzene assumes a more oblique position (−20.16 kcal mol−1 and −20.22 kcal mol−1, resp.) as was reported in the literature. The character of the COSMO-RS σ-surface of β-CD was much more hydrophobic on its O6 rim than on its O2/O3 side when all hydrogen bonds were arranged in a concerted mode.This static QM picture of the β-CD/benzene complex at 0 K was extended by MD simulations. At 120 K benzene was mobile but always stayed inside the cavity of β-CD. The trajectories at 273, 280, 290 and 300 K certainly no longer displayed the highly ordered hydrogen bonds of β-CD and benzene occupied many different positions inside the cavity, before it left the β-CD finally at its O2/O3 side.https://doi.org/10.3762/bjoc.9.15AM1benzeneCOSMO-RScyclodextrinhydrogen bondsinclusion complexmolecular dynamicsquantum mechanics |
spellingShingle | Jutta Erika Helga Köhler Nicole Grczelschak-Mick The β-cyclodextrin/benzene complex and its hydrogen bonds – a theoretical study using molecular dynamics, quantum mechanics and COSMO-RS Beilstein Journal of Organic Chemistry AM1 benzene COSMO-RS cyclodextrin hydrogen bonds inclusion complex molecular dynamics quantum mechanics |
title | The β-cyclodextrin/benzene complex and its hydrogen bonds – a theoretical study using molecular dynamics, quantum mechanics and COSMO-RS |
title_full | The β-cyclodextrin/benzene complex and its hydrogen bonds – a theoretical study using molecular dynamics, quantum mechanics and COSMO-RS |
title_fullStr | The β-cyclodextrin/benzene complex and its hydrogen bonds – a theoretical study using molecular dynamics, quantum mechanics and COSMO-RS |
title_full_unstemmed | The β-cyclodextrin/benzene complex and its hydrogen bonds – a theoretical study using molecular dynamics, quantum mechanics and COSMO-RS |
title_short | The β-cyclodextrin/benzene complex and its hydrogen bonds – a theoretical study using molecular dynamics, quantum mechanics and COSMO-RS |
title_sort | β cyclodextrin benzene complex and its hydrogen bonds a theoretical study using molecular dynamics quantum mechanics and cosmo rs |
topic | AM1 benzene COSMO-RS cyclodextrin hydrogen bonds inclusion complex molecular dynamics quantum mechanics |
url | https://doi.org/10.3762/bjoc.9.15 |
work_keys_str_mv | AT juttaerikahelgakohler thebcyclodextrinbenzenecomplexanditshydrogenbondsatheoreticalstudyusingmoleculardynamicsquantummechanicsandcosmors AT nicolegrczelschakmick thebcyclodextrinbenzenecomplexanditshydrogenbondsatheoreticalstudyusingmoleculardynamicsquantummechanicsandcosmors AT juttaerikahelgakohler bcyclodextrinbenzenecomplexanditshydrogenbondsatheoreticalstudyusingmoleculardynamicsquantummechanicsandcosmors AT nicolegrczelschakmick bcyclodextrinbenzenecomplexanditshydrogenbondsatheoreticalstudyusingmoleculardynamicsquantummechanicsandcosmors |