Complexation Thermodynamics of Modified Cyclodextrins: Extended Cavities and Distorted Structures

ITC analysis of polydisperse systems: Unravelling the impact of sample heterogeneity

Binding interactions often involve heterogeneous samples displaying a distribution of binding sites that vary in affinity and binding enthalpy. Examples include biological samples like proteins and chemically produced samples like modified cyclodextrins. Experimental studies often ignore sample heterogeneity and treat the system as an interaction of two homogeneous species, i.e. a chemically well-defined ligand binding to one type of site. The present study explores, by simulations and experiments, the impact of heterogeneity in isothermal titration calorimetry (ITC) setups where one of the binding components is heterogeneous. It is found that the standard single-site model, based on the assumption of two homogeneous binding components, provides excellent fits to simulated ITC data when the binding free energy is normally distributed and all sites have similar binding enthalpies. In such cases, heterogeneity can easily go undetected but leads to underestimated binding constants. Heterogeneity in the binding enthalpy is a bigger problem and may result in enthalpograms of increased complexity that are likely to be misinterpreted as two-site binding or other complex binding models. Finally, it is shown that heterogeneity can account for previously observed experimental anomalies. All simulations are accessible in Google Colab for readers to experiment with the simulation parameters.

Complexation of Lignin Dimers with β-Cyclodextrin and Binding Stability Analysis by ESI-MS, Isothermal Titration Calorimetry, and Molecular Dynamics Simulations

Lignin derived from lignocellulosic biomass is the largest source of renewable bioaromatics present on earth and requires environmentally sustainable separation strategies to selectively obtain high-value degradation products. Applications of supramolecular interactions have the potential to isolate lignin compounds from biomass degradation fractions by the formation of variable inclusion complexes with cyclodextrins (CDs). CDs are commonly used as selective adsorbents for many applications and can capture guest molecules in their internal hydrophobic cavity. The strength of supramolecular interactions between CDs and lignin model compounds that represent potential lignocellulosic biomass degradation products can be characterized by assessing the thermodynamics of binding stability. Consequently, the inclusion interactions of β-CD and lignin model compounds G-(β-O-4')-G, G-(β-O-4')-truncG (guaiacylglycerol-β-guaiacyl ether), and G-(β-β')-G (pinoresinol) were investigated empirically by electrospray ionization mass spectrometry and isothermal titration calorimetry, complemented by molecular dynamics (MD) simulations. Empirical results indicate that there are substantial differences in binding stability dependent on the linkage type. The lignin model β-β' dimer showed more potential bound states including 1:1, 2:1, and 1:2 (guest:host) complexation and, based on binding stability determinations, was consistently the most energetically favorable guest. Empirical results are supported by MD simulations that reveal that the capture of G-(β-β')-G by β-CD is promising with a 66% probability of being bound for G-(β-O-4')-truncG compared to 88% for G-(β-β')-G (unbiased distance trajectory and explicit counting of bound states). These outcomes indicate CDs as a promising material to assist in separations of lignin oligomers from heterogeneous mixtures for the development of environmentally sustainable isolations of lignin compounds from biomass fractions.

Thermodynamic Studies of Interactions between Sertraline Hydrochloride and Randomly Methylated β-Cyclodextrin Molecules Supported by Circular Dichroism Spectroscopy and Molecular Docking Results

The interaction between sertraline hydrochloride (SRT) and randomly methylated β-cyclodextrin (RMβCD) molecules have been investigated at 298.15 K under atmospheric pressure. The method used-Isothermal Titration Calorimetry (ITC) enabled to determine values of the thermodynamic functions like the enthalpy (ΔH), the entropy (ΔS) and the Gibbs free energy (ΔG) of binding for the examined system. Moreover, the stoichiometry coefficient of binding (n) and binding/association constant (K) value have been calculated from the experimental results. The obtained outcome was compared with the data from the literature for other non-ionic βCD derivatives interacting with SRT and the enthalpy-entropy compensation were observed and interpreted. Furthermore, the connection of RMβCD with SRT was characterized by circular dichroism spectroscopy (CD) and complexes of βCD derivatives with SRT were characterized through the computational studies with the use of molecular docking (MD).

Cyclodextrin complexation studies as the first step for repurposing of chlorpromazine

The antipsychotic drug chlorpromazine (CPZ) has potential for the treatment of acute myeloid leukemia, if central nervous system side-effects resulting from its passage through the blood-brain barrier can be prevented. A robust drug delivery system for repurposed CPZ would be drug-in-cyclodextrin-in-liposome that would redirect the drug away from the brain while avoiding premature release in the circulation. As a first step, CPZ complexation with cyclodextrin (CD) has been studied. The stoichiometry, binding constant, enthalpy, and entropy of complex formation between CPZ and a panel of CDs was investigated by isothermal titration calorimetry (ITC). All the tested CDs were able to include CPZ, in the form of 1:1, 1:2 or a mixture of 1:1 and 1:2 complexes. In particular, a substituted γ-CD, sugammadex (the octasodium salt of octakis(6-deoxy-6-S-(2-carboxyethyl)-6-thio)cyclomaltooctaose), formed exclusively 1:2 complexes with an extremely high association constant of 6.37 × 10⁹ M^-2. Complexes were further characterized by heat capacity changes, one- and two-dimensional (ROESY) nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics simulations. Finally, protection of CPZ against photodegradation by CDs was assessed. This was accelerated rather than reduced by complexation with CD. Altogether these results provide a molecular basis for the use of CD in delayed release formulations for CPZ.

Complexation Kinetics of Cyclodextrins with Bile Salt Anions: Energy Barriers for the Threading of Ionic Groups

Binding constants for thousands of cyclodextrin complexes have been reported in the literature, but much less is known about the kinetics of these host-guest complexes. In the present study, inclusion complexes of bile salts with β-cyclodextrin, γ-cyclodextrin, and a methylated β-cyclodextrin were studied by nuclear magnetic resonance (NMR) lineshape analysis to explore the structural factors that govern the complexation kinetics. For complexes with β-cyclodextrin, the association rate constants ranged from 2 × 10⁶ to 2 × 10⁷ M^-1 s^-1 while the dissociation rate constants ranged from 12 to 6000 s^-1 at 25 °C. The kinetics were thus significantly slower than for any other β-cyclodextrin complex reported in the literature, due to the large energy barrier for threading the ionic sidechains of the bile salt anions. Bile salts with taurine and glycine sidechains had identical binding affinities, but the kinetics differed by a factor of 10. Introduction of a single hydroxyl group at the binding site of the bile salts reduced the lifetimes and binding constants of the complexes by more than 50 times. The strong temperature dependence of the rate constants revealed that the large activation energies were mainly enthalpic with a small contribution from entropy. The larger γ-cyclodextrin was threaded by the nonionic end of the bile salts, and the kinetics were too fast to be accurately determined. The study demonstrates that ionic groups on guest molecules constitute significant energy barriers for the threading and dethreading of β-cyclodextrin hosts.

Identifying absolute configurations of PCB atropisomers by comparison of their experimental specific rotations with their DFT calculated values

Nineteen enantiomer pairs of polychlorinated biphenyls (PCBs) with three or four chloro substituents about the central carbon–carbon bond form a stable subclass of compounds whose biological effects vary with their chirality. Optical rotations for this group of PCBs were determined from density functional calculations employing extended atomic orbital gauge invariant basis sets. A comparison of these results with the experimental ones found from the literature for 10 of the pairs enabled the identification of their absolute configurations as analytes in gas chromatography studies.

Extracavity Effect in Cyclodextrin/Surfactant Complexation

Cyclodextrin (CD) complexation is a convenient method to sequester surfactants in a controllable way, for example, during membrane-protein reconstitution. Interestingly, the equilibrium stability of CD/surfactant inclusion complexes increases with the length of the nonpolar surfactant chain even beyond the point where all hydrophobic contacts within the canonical CD cavity are saturated. To rationalize this observation, we have dissected the inclusion complexation equilibria of a structurally well-defined CD, that is, heptakis(2,6-di- O-methyl)-β-CD (DIMEB), and a homologous series of surfactants, namely, n-alkyl- N, N-dimethyl-3-ammonio-1-propanesulfonates (SB3- x) with chain lengths ranging from x = 8 to 14. Thermodynamic parameters obtained by isothermal titration calorimetry and structural insights derived from nuclear magnetic resonance spectroscopy and molecular dynamics simulations revealed that, upon inclusion, long-chain surfactants with x = ≥10 extend beyond the canonical CD cavity. This enables the formation of hydrophobic contacts between long surfactant chains and the extracavity parts of DIMEB, which make additional favorable contributions to the stability of the inclusion complex. These results explain the finding that the stability of CD/surfactant inclusion complexes monotonously increases with the surfactant chain length even for long chains that completely fill the canonical CD cavity.

Validation and Comparison of Force Fields for Native Cyclodextrins in Aqueous Solution

Molecular dynamics simulations of native α-, β-, and γ-cyclodextrin in aqueous solution have been conducted with the goal to investigate the performance of the CHARMM36 force field, the AMBER-compatible q4md-CD force field, and five variants of the GROMOS force field. The properties analyzed are structural parameters derived from X-ray diffraction and NMR experiments as well as hydrogen bonds and hydration patterns, including hydration free enthalpies. Recent revisions of the torsional-angle parameters for carbohydrate systems within the GROMOS family of force fields lead to a significant improvement of the agreement between simulated and experimental NMR data. Therefore, we recommend using the variant 53A6_GLYC instead of 53A6 and 56A6_{CARBO_R} or 2016H66 instead of 56A6_CARBO to simulate cyclodextrins in solution. The CHARMM36 and q4md-CD force fields show a similar performance as the three recommended GROMOS parameter sets. A significant difference is the more flexible nature of the cyclodextrins modeled with the CHARMM36 and q4md-CD force fields compared to the three recommended GROMOS parameter sets.

Free-energy patterns in inclusion complexes: the relevance of non-included moieties in the stability constants

A MD/PMF-based procedure is designed for quantification of the interaction and respective components, guiding complex formation in water between β-CD and several naphthalene derivatives, highlighting the relevance of substituents.

Thermodynamic investigation of the interaction between cyclodextrins and preservatives - Application and verification in a mathematical model to determine the needed preservative surplus in aqueous cyclodextrin formulations

Preservatives are inactivated when added to conserve aqueous cyclodextrin (CD) formulations due to complex formation between CDs and the preservative. To maintain the desired conservation effect the preservative needs to be added in apparent surplus to account for this inactivation. The purpose of the present work was to establish a mathematical model, which defines this surplus based upon knowledge of stability constants and the minimal concentration of preservation to inhibit bacterial growth. The stability constants of benzoic acid, methyl- and propyl-paraben with different frequently used βCDs were determined by isothermal titration calorimetry. Based upon this knowledge mathematical models were constructed to account for the equilibrium systems and to calculate the required concentration of the preservations, which was evaluated experimentally based upon the USP/Ph. Eur./JP monograph. The mathematical calculations were able to predict the needed concentration of preservation in the presence of CDs; it clearly demonstrated the usefulness of including all underlying chemical equilibria in a mathematical model, such that the formulation design can be based on quantitative arguments.