On the structure of anhydrous β-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.

Hydrogen bond networks in gas-phase complex anions

Hydrogen bond networks (HBNs) have piqued the interest of the scientific community due to their crucial roles in nature. However, HBNs that are isolated from complicated backgrounds for unraveling their characteristics are still scarce. Herein, we propose that HBNs exist in complex anions formed between α-cyclodextrin (α-CD) and four benzoic acids (RBAs) in the gas phase. The complex anions are facilely extracted from solutions via the electrospray ionization technique, and subsequently activated through collision for the investigation of their transition dynamics. It is revealed that the generation of deprotonated α-CD and neutral RBAs is the unexpected dominant dissociation pathway for all the four complex anions, and the complex anions formed from more acidic RBAs exhibit higher stabilities. These dissociation results are successfully explained by the cooperative stretching dynamics of the proposed HBNs that are formed involving the intramolecular HBN of α-CD and the intermolecular hydrogen bonds (HBs) between α-CD and RBAs. Furthermore, the rarely observed low barrier HBs (LBHBs) are suggested to be present in the HBNs. It is believed that the present complex anions can serve as a facilely accessible and informative model for studying HBNs in the future.

Hydrogen bond networks and low barrier hydrogen bonds are demonstrated in the complex anions formed between α-cyclodextrin and benzoic acids.

Unravelling the structures of sodiated β-cyclodextrin and its fragments

We present cryogenic infrared spectra of sodiated β-cyclodextrin [β-CD + Na]+, a common cyclic oligosaccharide, and its main dissociation products upon collision-induced dissociation (CID). We characterize the parent ions using high-resolution ion mobility spectrometry and cryogenic infrared action spectroscopy, while the fragments are characterized by their mass and cryogenic infrared spectra. We observe sodium-cationized fragments that differ in mass by 162 u, corresponding to Bn/Zm ions. For the m/z 347 product ion, electronic structure calculations are consistent with formation of the lowest energy 2-ketone B2 ion structure. For the m/z 509 product ion, both the calculated 2-ketone B3 and the Z3 structures show similarities with the experimental spectrum. The theoretical structure most consistent with the spectrum of the m/z 671 ions is a slightly higher energy 2-ketone B4 structure. Overall, the data suggest a consistent formation mechanism for all the observed fragments.

Computational Insights Into the Influence of Substitution Groups on the Inclusion Complexation of β-Cyclodextrin

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.

First Steps to Rationalize Host-Guest Interaction between α-, β-, and γ-Cyclodextrin and Divalent First-Row Transition and Post-transition Metals (Subgroups VIIB, VIIIB, and IIB)

Cyclodextrins (CDs) are cyclic oligosaccharides mainly composed of six, seven, and eight glucose units, so-called α-, β-, and γ-CDs, respectively. They own a very particular molecular structure exhibiting hydrophilic features thanks to primary and secondary rims and delimiting a hydrophobic internal cavity. The latter can encapsulate organic compounds, but the former can form supramolecular complexes by hydrogen-bonding or electrostatic interactions. CDs have been used in catalytic processes to increase mass transfer in aqueous-organic two-phase systems or to prepare catalysts. In the last case, interaction between CDs and metal salts was considered to be a key point in obtaining highly active catalysts. Up to now, no work was reported on the investigation of factors affecting the binding of metal to CD. In the study herein, we present the favorable combination of electrospray ionization coupled to mass spectrometry [ESI-MS(/MS)] and density functional theory molecular modeling [B3LYP/Def2-SV(P)] to delineate some determinants governing the coordination of first-row divalent transition metals (Mn²⁺, Co²⁺, Ni²⁺, Cu²⁺, and Fe²⁺) and one post-transition metal (Zn²⁺) with α-, β-, and γ-CDs. A large set of features concerning the metal itself (ionic radius, electron configuration, and spin state) as well as the complexes formed (the most stable conformer, relative abundance in MS, CE₅₀ value in MS/MS, binding energy, effective coordination number, average bond lengths, binding site localization, bond dissociation energies, and natural bond orbital distribution) were screened. Taking into account all of these properties, various selectivity rankings have been delineated, portraying differential association/dissociation behaviors. Nonetheless, unique 3D topologies for each CD-metal complex were emphasized. The combination of these approaches brings a stone for building a compendium of molecular features to serve as a suitable descriptor or predictor for a better first round rationalization of catalytic activities.

On Complex Formation between 5-Fluorouracil and β-Cyclodextrin in Solution and in the Solid State: IR Markers and Detection of Short-Lived Complexes by Diffusion NMR

In this work, the nuclear magnetic resonance (NMR) and IR spectroscopic markers of the complexation between 5-fluorouracil (5-FU) and β-cyclodextrin (β-CD) in solid state and in aqueous solution are investigated. In the attenuated total reflectance(ATR) spectra of 5-FU/β-CD products obtained by physical mixing, kneading and co-precipitation, we have identified the two most promising marker bands that could be used to detect complex formations: the C=O and C-F stretching bands of 5-FU that experience a blue shift by ca. 8 and 2 cm-1 upon complexation. The aqueous solutions were studied by NMR spectroscopy. As routine NMR spectra did not show any signs of complexation, we have analyzed the diffusion attenuation of spin-echo signals and the dependence of the population factor of slowly diffusing components on the diffusion time (diffusion NMR of pulsed-field gradient (PFG) NMR). The analysis has revealed that, at each moment, ~60% of 5-FU molecules form a complex with β-CD and its lifetime is ca. 13.5 ms. It is likely to be an inclusion complex, judging from the independence of the diffusion coefficient of β-CD on complexation. The obtained results could be important for future attempts of finding better methods of targeted anticancer drug delivery.

Cavity Closure of 2-Hydroxypropyl-β-Cyclodextrin: Replica Exchange Molecular Dynamics Simulations

2-Hydroxypropyl-β-cyclodextrin (HPβCD) has unique properties to enhance the stability and the solubility of low water-soluble compounds by inclusion complexation. An understanding of the structural properties of HPβCD and its derivatives, based on the number of 2-hydroxypropyl (HP) substituents at the α-d-glucopyranose subunits is rather important. In this work, replica exchange molecular dynamics simulations were performed to investigate the conformational changes of single- and double-sided HP-substitution, called 6-HPβCDs and 2,6-HPβCDs, respectively. The results show that the glucose subunits in both 6-HPβCDs and 2,6-HPβCDs have a lower chance of flipping than in βCD. Also, HP groups occasionally block the hydrophobic cavity of HPβCDs, thus hindering drug inclusion. We found that HPβCDs with a high number of HP-substitutions are more likely to be blocked, while HPβCDs with double-sided HP-substitutions have an even higher probability of being blocked. Overall, 6-HPβCDs with three and four HP-substitutions are highlighted as the most suitable structures for guest encapsulation, based on our conformational analyses, such as structural distortion, the radius of gyration, circularity, and cavity self-closure of the HPβCDs.

A molecular dynamics study of conformations of beta-cyclodextrin and its eight derivatives in four different solvents

Understanding the atomic level interactions and the resulting structural characteristics is required for developing beta-cyclodextrin (βCD) derivatives for pharmaceutical and other applications. The effect of four different solvents on the structures of the native βCD and its hydrophilic (methylated βCD; MEβCD and hydroxypropyl βCD; HPβCD) and hydrophobic derivatives (ethylated βCD; ETβCD) was explored using molecular dynamics (MD) simulations and solvation free energy calculations. The native βCD, 2-MEβCD, 6-MEβCD, 2,6-DMβCD, 2,3,6-TMβCD, 6-HPβCD, 2,6-HPβCD and 2,6-ETβCD in non-polar solvents (cyclohexane; CHX and octane; OCT) were stably formed in a symmetric cyclic cavity shape through their intramolecular hydrogen bonds. In contrast, βCDs in polar solvents (methanol; MeOH and water; WAT) exhibited large structural changes and fluctuations leading to significant deformations of their cavities. Hydrogen bonding with polar solvents was found to be one of the major contributors to this behavior: solvent-βCD hydrogen bonding strongly competes with intramolecular bonding leading to significant changes in the structural stability of βCDs. An exception to this is the hydrophobic 2,6-ETβCD which retained its spherical cavity in all solvents. Based on this, it is proposed that the 2,6-ETβCD can act as a sustained release drug carrier.

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.

Approximate quantum chemical methods for modelling carbohydrate conformation and aromatic interactions: β-cyclodextrin and its adsorption on a single-layer graphene sheet

Adsorption of carbohydrates on graphene has the potential to improve graphene dispersibility in water. Here we assess the ability of DFTB-based and NDDO-based quantum chemical methods to model β-cyclodextrin conformations and interactions with graphene.

Quantum chemical study and isothermal titration calorimetry of β-cyclodextrin complexes with mianserin in aqueous solution

The structures, interaction energies and thermodynamics of the complex formation between mianserin (MIA) and β-cyclodextrin (β-CD) are investigated using computational methods and calorimetric measurements.

Metadynamics supports molecular dynamics simulation-based binding affinities of eucalyptol and beta-cyclodextrin inclusion complexes

The development of various molecular dynamics methods enables the detailed investigation of association processes, like host–guest complexes, including their dynamics and, additionally, the release of the guest compound.

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 complexation of pinostrobin with various cyclodextrin derivatives

Pinostrobin (PNS) is one of the important flavonoids and can be abundantly found in the rhizomes of fingerroot (Boesenbergia rotrunda) and galangal (Alpinia galangal and Alpinia officinarum), the herbal basis of Southeast Asian cooking. Similar to other flavonoids, PNS exhibits anti-oxidative, anti-inflammatory and anti-cancer properties. However, this compound has an extremely low water solubility that limits its use in pharmaceutical applications. Beta-cyclodextrin (βCD) and its derivatives, 2,6-dimethyl-βCD (2,6-DMβCD) and the three hydroxypropyl-βCDs (2-HPβCD, 6-HPβCD and 2,6-DHPβCD), have unique properties that enhance the stability and solubility of such low-soluble guest molecules. In the present study, molecular dynamics simulations were applied to investigate the dynamics and stability of PNS inclusion complexes with βCD and its derivatives (2,6-DMβCD, 2,6-DHPβCD, 2-HPβCD and 6-HPβCD). PNS was able to form complexes with βCD and all four of its derivatives by either the chromone (C-PNS) or phenyl (P-PNS) ring dipping toward the cavity. According to the molecular mechanics-generalized Born surface area binding free energy values, the stability of the different PNS/βCD complexes was ranked as 2,6-DHPβCD>2,6-DMβCD>2-HPβCD>6-HPβCD>βCD. These theoretical results were in good agreement with the stability constants that had been determined by the solubility method.

Probing the common alkali metal affinity of native and variously methylated β-cyclodextrins by combining electrospray-tandem mass spectrometry and molecular modeling

In the study herein, we investigated the solution and gas phase affinity of native and variously methylated β-cyclodextrins (CDs) as hosts towards three common alkali metals as guests namely lithium, sodium and potassium. For this purpose, two complementary approaches have been employed: electrospray-tandem mass spectrometry (ESI-MS/MS) with two energetic regimes: Collision Induced Dissociation (CID) and Higher Collision Dissociation (HCD), respectively, and DFT molecular modeling. These approaches have been achieved by taking into account the interaction of either one or two alkali metals with the host molecules. The results showed a good agreement between experimental and theoretical data. It was demonstrated that increasing the methylation degree strengthened the gas phase affinity towards all studied alkali metals. Furthermore, it was established that the cation selectivity was Na(+) > Li(+) > K(+) and Li(+) > Na(+) > K(+) for the solution and gas phase, respectively.

Infrared multiple photon dissociation spectroscopy and density functional theory (DFT) studies of protonated permethylated β-cyclodextrin-water non-covalent complexes

We present infrared multiple photon dissociation (IRMPD) spectroscopy and quantum chemical calculation results for the protonated permethylated β-cyclodextrin (CD)-water non-covalent complex, the simplest β-CD non-covalent complex, in the gas-phase. The IRMPD spectrum in the region 2700-3750 cm(-1) consisted of three strong peaks at 3096, 3315, and 3490 cm(-1). These spectral features in the experimental IRMPD spectrum were compared with a large set of infrared absorption spectra predicted using density functional theory (DFT) calculations for the protonated β-CD-water complex. Complex III (see ), in which the water molecule (at the primary rim) and the proton (at the secondary rim) were separated, was found to suitably reflect the main spectral characteristics found in the experimental IRMPD spectrum. The absence of the homodromic hydrogen bond ring, due to replacement of hydroxyl groups with methoxy groups in permethylated β-CD, rendered the primary rim open compared with the unmodified β-CD 'one-gate-closed' lowest energy conformer. This study demonstrates that IRMPD studies combined with DFT theoretical calculations can be a good method for studying molecular interactions of large host-guest pairs.

Binding mode and free energy prediction of fisetin/β-cyclodextrin inclusion complexes

In the present study, our aim is to investigate the preferential binding mode and encapsulation of the flavonoid fisetin in the nano-pore of β-cyclodextrin (β-CD) at the molecular level using various theoretical approaches: molecular docking, molecular dynamics (MD) simulations and binding free energy calculations. The molecular docking suggested four possible fisetin orientations in the cavity through its chromone or phenyl ring with two different geometries of fisetin due to the rotatable bond between the two rings. From the multiple MD results, the phenyl ring of fisetin favours its inclusion into the β-CD cavity, whilst less binding or even unbinding preference was observed in the complexes where the larger chromone ring is located in the cavity. All MM- and QM-PBSA/GBSA free energy predictions supported the more stable fisetin/β-CD complex of the bound phenyl ring. Van der Waals interaction is the key force in forming the complexes. In addition, the quantum mechanics calculations with M06-2X/6-31G(d,p) clearly showed that both solvation effect and BSSE correction cannot be neglected for the energy determination of the chosen system.

Enhanced stability of a naringenin/2,6-dimethyl β-cyclodextrin inclusion complex: molecular dynamics and free energy calculations based on MM- and QM-PBSA/GBSA

The structure, dynamic behavior and binding affinity of the inclusion complexes between naringenin and the two cyclodextrins (CDs), β-CD and its 2,6-dimethyl derivative (DM-β-CD), were theoretically studied by multiple molecular dynamics simulations and free energy calculations. Naringenin most likely prefers to bind with CDs through the phenyl ring. Although a lower hydrogen bond formation of naringenin with the 3-hydroxyl group of DM-β-CD (relative to β-CD) was observed, the higher cavity could encapsulate almost the whole naringenin molecule. In contrast for the naringenin/β-CD complex, the phenyl ring feasibly passed through the primary rim resulting in the chromone ring binding inside instead. MM-PBSA/GBSA and QM-PBSA/GBSA binding free energies strongly suggested a greater stability of the naringenin/DM-β-CD inclusion complex. Van der Waals force played an important role as the key guest-host interaction for the complexation between naringenin and each cyclodextrin.

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.

Charge sensitivity approach to mutual polarization of reactants: molecular mechanics perspective

Charge sensitivity analysis (CSA) in force-field atoms resolution was applied to describe the mutual polarization of reactants as well as charge-transfer (CT) effects. An inclusion complex of β-cyclodextrin with salicylic acid was used as a model system. Three CSA models were taken into account and verified on a Born–Oppenheimer molecular dynamics (BOMD) trajectory. The models differed in terms of the equilibrium conditions imposed on the system. It was demonstrated that mutual polarization is an important source of stabilization, in contrast to the results obtained from static charge calculations. The energy lowering induced by CT was small and comparable to the CT stabilization that occurs in hydrogen-bonded systems. All models correctly described the main topological features of the BOMD energy surface. CSA in force-field atoms resolution qualitatively reproduced the charge reorganization accompanying hydrogen-bond formation. It was shown that CSA parameters are very sensitive to the bond formation process, which suggests that they could be applied in reactive force fields as detectors of newly formed chemical bonds.

Figure

Intramolecular hydrogen bond energy and cooperative interactions in α-, β-, and γ-cyclodextrin conformers

Accurate estimation of individual intramolecular hydrogen bond (H-bond) energies is an intricate task for multiply H-bonded systems. In such cases, the hydrogen bond strengths could be highly influenced by the cooperative interactions, for example, those between hydroxyl groups in sugars. In this work, we use the recently proposed molecular tailoring approach-based quantification (Deshmukh, Gadre, and Bartolotti, J Phys Chem A 2006, 110, 12519) to the extended systems of cyclodextrins (CDs). Further, the structure and stability of different conformers of α-, β-, and γ-CDs are explained based on the energetics and cooperative contribution to the strength of these H-bonds. The estimated O-H···O H-bond energies in the various CD conformers are found to vary widely from 1.1 to 8.3 kcal mol(-1). The calculated energy contributions to cooperativity toward the H-bond strengths fall in the range of 0.25-2.75 kcal mol(-1).

Theoretical studies on hydrogen bonding, NMR chemical shifts and electron density topography in alpha, beta and gamma-cyclodextrin conformers

Hydrogen-bonded interactions in alpha-, beta-, and gamma-CD conformers are investigated from the molecular electron density topography and chemical shift in the nuclear magnetic resonance (NMR) spectra calculated by using the Gauge Invariant Atomic Orbital (GIAO) method within the framework of density functional theory. For the lowest-energy CD conformers in the gas phase, the O3-H...O2' hydrogen-bonding interactions are present. Calculated 1H NMR chemical shifts (delta H) correlate well with the hydrogen-bond distance as well as electron density at the bond critical point in the molecular electron density (MED) topography. The conformers of beta- and gamma-CD comprised of relatively strong secondary hydroxyl interactions are stabilized by solvation from polar solvents.