Ab Initio Calculations on Low-Energy Conformers of α-Cyclodextrin

Systematic Study of Different Types of Interactions in α-, β- and γ-Cyclodextrin: Quantum Chemical Investigation

In this work, comprehensive ab initio quantum chemical calculations using the DFT level of theory were performed to characterize the stabilization interactions (H-bonding and hyperconjugation effects) of two stable symmetrical conformations of α-, β-, and γ-cyclodextrins (CDs). For this purpose, we analyzed the electron density using "Atom in molecules" (AIM), "Natural Bond Orbital" (NBO), and energy decomposition method (CECA) in 3D and in Hilbert space. We also calculated the H-bond lengths and OH vibrational frequencies. In every investigated CD, the quantum chemical descriptors characterizing the strength of the interactions between the H-bonds of the primary OH (or hydroxymethyl) and secondary OH groups are examined by comparing the same quantity calculated for ethylene glycol, α-d-glucose (α-d-Glcp) and a water cluster as reference systems. By using these external standards, we can characterize more quantitatively the properties of these bonds (e.g., strength). We have demonstrated that bond critical points (BCP) of intra-unit H-bonds are absent in cyclodextrins, similar to α-d-Glcp and ethylene glycol. In contrast, the CECA analysis showed the existence of an exchange (bond-like) interaction between the interacting O…H atoms. Consequently, the exchange interaction refers to a chemical bond, namely the H-bond between two atoms, unlike BCP, which is not suitable for its detection.

Exploration of the Conformational Scenario for α-, β-, and γ-Cyclodextrins in Dry and Wet Conditions, from Monomers to Crystal Structures: A Quantum-Mechanical Study

Cyclodextrins (CDs) constitute a class of cyclic oligosaccharides that are well recognized and largely applied in the drug delivery field, thanks to their biocompatibility, low cost, and the possibility to be derivatized in order to tune and optimize the complexation/release of the specific drug. The conformational flexibility of these systems is one of their key properties and requires a cost-effective methodology to be studied by combining the accuracy of results with the possibility of exploring a large set of conformations. In the present paper, we have explored the conformational potential energy surface of the monomers and dimers of α-, β-, and γ-cyclodextrins (i.e., 6, 7, and 8 monomeric units, respectively) by means of fast but accurate semiempirical methods, which are then refined by state-of-the-art DFT functionals. Moreover, the crystal structure is considered for a more suitable comparison with the IR spectrum experimentally recorded. Calculations are carried out in the gas phase and in water environments, applying both implicit and explicit treatments. We show that the conformation of the studied molecules changes from the gas phase to the water, even if treated implicitly, thus modifying their complexation capability.

Exploring halide anion affinities to native cyclodextrins by mass spectrometry and molecular modelling

The binding affinities of cyclodextrins complexation with chlorine (Cl^-), bromine (Br^-) and iodine (I^-), were measured by mass spectrometric titrimetry, and the fitting of the binding constants was based on the concentration measurement of the cyclodextrin equilibrium. The binding constants (lg K_a) for α-, β- or γ-cyclodextrin with Cl^- were 3.99, 4.03 and 4.11, respectively. The gas-phase binding affinity of halide anions for native cyclodextrins was probed by collision-induced dissociation. In collision-induced dissociation, the centre-of-mass frame energy results revealed that in the gas phase, for the same type of cyclodextrin, the stability of the complexes decreased in order: Cl > Br > I, and for the same halide anion, the binding stability of the complex with α-, β- or γ-cyclodextrin decreased in the order: γ-cyclodextrin >β-cyclodextrin > α-cyclodextrin. The density functional theory calculations showed that halide anion binding on the primary face had a lower energy than the secondary face and hydrogen bonding was the main driving force for complex formation. The higher stability of the γ-cyclodextrin complex with the Cl anion can be attributed to the higher charge density of the Cl anion and better flexibility of γ-cyclodextrin.

Structures and stabilities of naturally occurring cyclodextrins: a theoretical study of symmetrical conformers

A molecular modeling study of symmetrical conformers of α-, β-, and γ-cyclodextrins in the gas and aqueous phases was carried out using the M06-2X density functional method, with SMD employed as an implicit solvation model. Eight symmetrical conformers were found for each cyclodextrin. Values of geometrical parameters obtained from the modeling study were found to agree well with those obtained from X-ray diffraction structures. A vibrational analysis using harmonic frequencies was performed to determine thermodynamic quantities. The GIAO method was applied to determine proton and carbon-13 NMR chemical shifts, which were then compared with corresponding chemical shifts reported in the literature. Hydrogen-bonding patterns were analyzed using geometrical descriptors, and quantum chemical topology was explored by QTAIM analysis. The results of this study indicated that four of the eight conformers studied for each cyclodextrin are the most populated in aqueous solution. These results provide the foundations for future studies of host-guest complexes involving these cyclodextrins. Graphical abstract δΔG_solvation: variation of free Gibss energy of solvation.

Formation of β-cyclodextrin complexes in an anhydrous environment

The formation of inclusion complexes of β-cyclodextrin was studied at the melting temperature of guest compounds by differential scanning calorimetry. The complexes of long-chain n-alkanes, polyaromatics, and organic acids were investigated by calorimetry and IR spectroscopy. The complexation ratio of β-cyclodextrin was compared with results obtained in an aqueous environment. The stability and structure of inclusion complexes with various stoichiometries were estimated by quantum chemistry and molecular dynamics calculations. Comparison of experimental and theoretical results confirmed the possible formation of multiple inclusion complexes with guest molecules capable of forming hydrogen bonds. This finding gives new insight into the mechanism of formation of host-guest complexes and shows that hydrophobic interactions play a secondary role in this case. Graphical abstract The formation of complexes of β-cyclodextrin with selected n-alkanes, polyaromatics, and organic acids in an anhydrous environment is studied by differential scanning calorimetry, IR spectroscopy, and molecular modeling. The results obtained confirm the possible formation of multiple inclusion complexes with guest molecules capable of forming hydrogen bonds and give a new perspective on the mechanism of formation of host-guest complexes.

Inclusion complex thermodynamics: The β-cyclodextrin and sertraline complex example

Thermodynamic properties for β-cyclodextrin-Sertraline inclusion process was calculated at the density functional theory (DFT) level using the PBE0 functional with 6-31G(d,p), 6-31++G(d,p) and 6-311++G(2df,p) basis sets. Electron correlation was evaluated through Møller-Plesset second-order perturbation theory (MP2). The standard statistical thermodynamic approach was used to assess the entropic contribution to the Gibbs free energy value. According to our results, inclusion of hydration waters to describe the reactants and products in the complex formation reaction model is crucial in order to reproduce the experimental data and seems very coherent with basics thermodynamics yielding good agreement with experiment.

Inclusion complexes of cypermethrin and permethrin with monochlorotriazinyl-beta-cyclodextrin: a combined spectroscopy, TG/DSC and DFT study

The suitable size hydrophobic cavity and monochlorotriazinyl group as a reactive anchor make MCT-β-CD to be widely used in fabric finishing. In this paper, the inclusion complexes of monochlorotriazinyl-beta-cyclodextrin (MCT-β-CD) with cypermethrin (CYPERM) and permethrin (PERM) are synthesized and analyzed by TG/DSC, FT-IR and Raman spectroscopy. TG/DSC reveals that the decomposed temperatures of inclusion complexes are lower by 25-30 °C than that of physical mixtures. DFT calculations in conjunction with FT-IR and Raman spectral analyses are used to study the structures of MCT-β-CD and their inclusion complexes. Four isomers of trisubstituted MCT-β-CD are designed and DFT calculations reveal that 1,3,5-trisubstituted MCT-β-CD has the lowest energy and can be considered as main component of MCT-β-CD. The ground-state geometries, vibrational wavenumbers, IR and Raman intensities of MCT-β-CD and their inclusion complexes were calculated at B3LYP/6-31G (d) level of theory. Upon examining the optimized geometry of inclusion complex, we find that the CYPERM and PERM are inserted into the toroid of MCT-β-CD from the larger opening. The band at 1646 cm(-1) in IR and at 1668 cm(-1) in Raman spectrum reveals that monochloroazinyl group of MCT-β-CD exists in ketone form but not in anion form. The noticeable IR and Raman shift of phenyl reveals that these two benzene rings of CYPERM and PERM stays inside the cavity of MCT-β-CD and has weak interaction with MCT-β-CD. This spectroscopy conclusion is consistent with theoretical predicted structure.

Spin trapping of hydroperoxyl radical by a cyclic nitrone conjugated to β-cyclodextrin: a computational study

Spin trapping of hydroperoxyl radical (HOO^.) by the amide-linked conjugate of 5-carbamoyl-5-methyl-1-pyrroline N-oxide (AMPO) to β-cyclodextrin (β-CD) was studied computationally using a two-layered ONIOM method. From a conformational perspective, the "internal" conformation of 5R-β-CD-AMPO is more favored than the "external" conformation in which the nitrone is located outside of the cavity of the β-CD. When the HOO^. addition product is formed, the most stable isomer has the nitroxyl (N₁-O₁) moiety pointing inside the cavity of the β-CD. Thus, this "internal" conformation might protect the N₁-O₁ moiety of the resulting spin adduct from access by reducing agents, thereby improving the lifetime of the radical adduct. The computed energetic barrier for HOO^. addition to the 5R-β-CD-AMPO is 8.7 kcal/mol, which is marginally smaller than spin trapping by the non-conjugated AMPO (that is, without the β-CD). To optimize the reactivity of the β-CD-AMPO conjugate, the effect of a spacer unit between the AMPO segment and the β-CD moiety with varying methylene units, (CH₂) _n (n = 1, 2, 3), on the energetics of HOO^. addition was evaluated. The structure with only one methylene spacer (n = 1) appears to be optimal as determined by the smaller activation barrier (6.2 kcal/mol) for HOO^. addition to the nitrone moiety. Compared with very time-consuming quantum mechanical methods, the ONIOM method appears to offer significant advantages for evaluation of the best β-CD-AMPO conjugate for trapping of such reactive oxygen species and providing for the rational design of novel nitrones as spin traps.

First principles investigation of noncovalent complexation: a [2.2.2]-cryptand ion-binding selectivity study

A first principles methodology, aimed at understanding the roles of molecular conformation and energetics in host-guest binding interactions, is developed and tested on a system that pushes the practical limits of ab initio methods. The binding behavior between the [2.2.2]-cryptand host (4,7,13,16,21,24-hexaoxa-1,10-diaza-bicyclo[8.8.8]hexacosane) and alkali metal cations (Li(+), Na(+), and K(+)) in gas, water, methanol, and acetonitrile is characterized. Hartree-Fock and density functional theory methods are used in concert with crystallographic information to identify gas phase, energy-minimized conformations. Gas phase free energies of binding, with vibrational contributions, are compared to solution-state binding constants through relative binding selectivity analysis. Calculated relative binding free energies qualitatively correlated with solution state experiments only after gas phase metal desolvation is considered. The B3LYP exchange-correlation functional improves theoretical correlations with experimental relative binding free energies. The relevance of gas phase calculations towards understanding binding in condensed phases is discussed. Natural bond orbital methods highlights previously unidentified intramolecular and intermolecular M(+)(222) chemistries, such as an intramolecular n'(O)-->sigma*(CH) hydrogen bond.