Multimodal inclusion complexes between barbiturates and 2-hydroxypropyl-beta-cyclodextrin in aqueous solution: isothermal titration microcalorimetry, (13)C NMR spectrometry, and molecular dynamics simulation.
Multiple types (structures) of inclusion complexes between barbiturates and 2-hydroxypropyl-beta-cyclodextrin (HPCD) were evaluated by isothermal titration microcalorimetry and (13)C NMR spectroscopy. The geometries of the inclusion complexes were suggested by molecular dynamics simulation. Barbituric acid (BA), barbital (B), amobarbital (AB), pentobarbital (PB), secobarbital (SB), cyclobarbital (CB), and phenobarbital (PHB) were used as barbiturates with different substituents on the barbituric acid ring and compared for inclusion types in aqueous solution. The association constants (K), stoichiometries, and thermodynamic parameters change in free energy (DeltaG) change in enthalpy (DeltaH), and change in entropy [DeltaS] for each type of complex were determined from the calorimetric data. The inclusion complexation was largely entropy driven because of hydrophobic interactions. The values of K increased in the order BA<B<AB<PB<SB<CB<PHB. Barbiturates, except B and BA, form two types of inclusion complex with a 1:1 stoichiometry in the un-ionized forms. The first type of inclusion complex with high affinity (K(1)) was characterized by small negative values of DeltaH(1) and large positive DeltaS(1), where the substituent R2 of the barbiturate was initially inserted into the cavity of HPCD through hydrophobic interactions. There was a good relationship between DeltaG(1) obtained from the calorimetric data for the first type of inclusion complex and DeltaG(R2) calculated from the changes in (13)C Nuclear Magnetic Resonance (NMR) chemical shifts for the substituent R2 of barbiturates. These types were very stable in aqueous solution at various pHs. The second type of complex, with low affinity (K(2)), was characterized by large negative values of DeltaH(2) and small positive DeltaS(2), reflecting van der Waals' interactions in the un-ionized forms of barbiturates at pH values less than pK(a). The values of K(2) were markedly decreased to <10(3) M(-1) as the barbiturates were ionized over pH 8. Thus, in the second type, the barbituric acid ring contributed to forming the complexes. The geometries were stabilized by hydrogen bond formation between the hetero atoms in the barbituric acid ring and the secondary hydroxyl groups on the rim of the cyclodextrin. The (13)C NMR chemical shifts of C4 and C6 carbons in the barbituric acid ring were moved upfield significantly by the inclusion complexation. On the other hand, B and BA could form only one type of complex, the lid-type supramolecular complex with small association constants.