Sandra
Kaabel
ab,
Jasper
Adamson
c,
Filip
Topić
b,
Anniina
Kiesilä
b,
Elina
Kalenius
b,
Mario
Öeren
a,
Mart
Reimund
a,
Elena
Prigorchenko
a,
Aivar
Lõokene
a,
Hans J.
Reich
d,
Kari
Rissanen
*b and
Riina
Aav
*a
aDepartment of Chemistry, Tallinn University of Technology, Akadeemia tee 15, 12618 Tallinn, Estonia. E-mail: riina.aav@ttu.ee
bUniversity of Jyvaskyla, Department of Chemistry, Nanoscience Center, P.O. Box. 35, FI-40014 Jyvaskyla, Finland. E-mail: kari.t.rissanen@jyu.fi
cNational Institute of Chemical Physics and Biophysics, Akadeemia tee 23, 12618 Tallinn, Estonia
dDepartment of Chemistry, University of Wisconsin, Madison, WI 53706, USA
First published on 30th November 2016
A novel eight-membered macrocycle of the hemicucurbit[n]uril family, chiral (all-R)-cyclohexanohemicucurbit[8]uril (cycHC[8])‡ binds anions in a purely protic solvent with remarkable selectivity. The cycHC[8] portals open and close to fully encapsulate anions in a 1:1 ratio, resembling a molecular Pac-Man™. Comprehensive gas, solution and solid phase studies prove that the binding is governed by the size, shape and charge distribution of the bound anion. Gas phase studies show an order of SbF6− ≈ PF6− > ReO4− > ClO4− > SCN− > BF4− > HSO4− > CF3SO3− for anion complexation strength. An extensive crystallographic study reveals the preferred orientations of the anions within the octahedral cavity of cycHC[8] and highlights the importance of the size- and shape-matching between the anion and the receptor cavity. The solution studies show the strongest binding of the ideally fitting SbF6− anion, with an association constant of 2.5 × 105 M−1 in pure methanol. The symmetric, receptor cavity-matching charge distribution of the anions results in drastically stronger binding than in the case of anions with asymmetric charge distribution. Isothermal titration calorimetry (ITC) reveals the complexation to be exothermic and enthalpy-driven. The DFT calculations and VT-NMR studies confirmed that the complexation proceeds through a pre-complex formation while the exchange of methanol solvent with the anion is the rate-limiting step. The octameric cycHC[8] offers a unique example of template-controlled design of an electroneutral host for binding large anions in a competitive polar solvent.
Exemplary studies on anion binding by cyclic hexapeptides show that introducing a confined cavity to the structure of the receptor significantly increases the anion binding ability, as abundant 2:1 host–guest complexes were observed with halides and SO42− where the anion was enclosed in the cavity formed by two cyclopeptide units.10,11 This led to the design of sandwich-like bis(cyclopeptide) receptors that reached association constants of up to 106 M−1 for binding SO42− in a water–methanol mixture.12,13 Employing a hydrophobic pocket for anion binding in water has been illustrated by the Gibb's octa-acid cavitand,14 where the binding of partially hydrated anions yielded association constants up to 103 M−1 at pH 11.5. Likewise, the size-dependent complexation of dodecaborate dianions within the hydrophobic cavity of γ-cyclodextrin was studied by Nau and co-workers, with the strongest association in the case of B12Br122− (Ka = 9.6 × 105 M−1).15
Cucurbituril family members have been explored as ion receptors due to their well-defined hydrophobic cavity.16,17 Hemicucurbiturils,18 cyclohexanohemicucurbit[6]urils,19 bambus[6]urils20 and biotin[6]urils21–23 have been shown to bind anions within their electron-deficient (i.e. partially positively charged) hydrophobic cavities. Presently, bambus[6]urils hold the record for the strongest anion binding (Ka = 5.5 × 107 M−1) by a neutral host in exclusively protic solvent.24 Based on the crystal structures of six-membered hemicucurbiturils, bambus[6]uril25–28 and biotin[6]uril21 complexes the central cavities readily accommodate halide anions. Larger anions have been shown to prefer the formation of 1:2 host:guest complexes with double-cone-shaped bambus[6]urils,29,30 with the anions bound away from the centre of the cavity.
We have previously demonstrated the synthesis of the first 8-membered hemicucurbituril, (all-R)-cyclohexanohemicucurbit[8]uril (cycHC[8]), by an approach using anion-templating.31 Given that larger anions such as PF6− and CF3CO2− were observed to act as effective templates, we decided to investigate the binding of other larger inorganic anions by cycHC[8]. Such anions are for example used in ionic liquids32 (BF4−, PF6−, SbF6−, CF3SO3−) and as oxidizing agents33 (ClO4−, IO4−). On the other hand, they are considered as environmental pollutants.33,34 Binding of these anions in protic media is important from a biological and environmental point of view.
Fig. 1 Molecular structure of (all-R)-cyclohexanohemicucurbit[8]uril, cycHC[8] (left), and the X-ray structure of an inclusion complex with SbF6− (right). |
The scope of anionic guests that form complexes with cycHC[8] was determined by mass spectrometry (see ESI†). Only 1:1 complexes were observed in ESI-MS spectra and abundant complexes were observed with SbF6− ≈ PF6− > ReO4− > ClO4− > SCN− > BF4− > HSO4− > CF3SO3−, ranked by decreasing affinity for cycHC[8]. The anions H2PO4−, AcO−, Br−, Cl−, I−, F−, NO3− and NO2− showed only weak complexation, while no complexes were formed with AuBr4−, Br3− and CN−. The order of affinity was ascertained through competition experiments, performed on three-component mixtures of the host with two competing anions in 1:1:1 molar ratio (ESI Fig. S2 and S3†). Halide anions (13 to 35 Å3 in volume36), while readily forming complexes with 6-membered hemicucurbiturils,37–41 appear to have very low affinity towards cycHC[8], presumably due to a mismatch in size with the cavity of cycHC[8]. As expected, the affinity towards more heavily solvated anions was found to be lower than for weakly solvated ones. On the other hand, tetrahedral and octahedral anions falling into the volume range of 50 to 80 Å3 form stronger complexes with cycHC[8], presumably due to a better size fit. This is also in good agreement with Rebek's rule suggesting a packing coefficient (PC)42 of 0.55 ± 0.09 for optimal fit in the cycHC[8] cavity with a volume of 123 Å3.31 The MS/MS collision-induced dissociation (CID) experiments (Fig. S4†) on isolated complexes revealed the most efficient dissociation (lowest kinetic stability) for the host–guest complexes with highest PCs (>0.55). This results from the sensitive interplay of the attractive and repulsive forces between the anion and the cavity walls, and the lack of void space. More importantly, it also indicates the full encapsulation of the anions in the complexes in the gas phase. The same phenomenon has been studied in detail with cucurbiturils and azoalkanes,43 but is reported here for the first time with anionic complexes.
The crystal structures of the host–guest complexes, obtained by single crystal X-ray diffraction, demonstrate unambiguously the 1:1 stoichiometry of the anion inclusion complexes in the solid state (Fig. 2). Single crystals of the complexes were obtained from solutions of cycHC[8] in methanol with the guest added as a tetrabutylammonium (TBA) or tetrabutyl-phosphonium (TBP) salt. The TBA or TBP cations and a number of solvent molecules fill the space between the capsule-like moieties, affording isostructural crystals regardless of the guest anion used. The guest anions are situated at the center of the cycHC[8] cavity, in a manner depending on their size and shape. The octahedral anions SbF6− and PF6− are locked in a fixed position showing no disorder. Four fluorine atoms of SbF6−/PF6− lying on the equatorial plane of the macrocycle point to the four corners of its square-shaped belt, while the two remaining fluorines point to the opposite portals (Fig. 2a and b). The Hirshfeld surface44 plotted for the encapsulated octahedral anion SbF6− indicates that the host–guest C–H⋯F interactions are responsible for the fixed orientation of these guests (ESI Fig. S6†). The shortest C–H⋯F distances are found between the four equatorial fluorine atoms and the axial 2ax protons in each corner of the macrocyclic cavity (Fig. 3B and Table S3†). In contrast to the octahedral anions, the tetrahedral anions BF4−, ClO4−, ReO4− and IO4− have more space inside the cavity (Fig. 2c–f) and show disorder in the crystal. As the host cavity is symmetric and therefore offers a number of equal interaction sites, several orientations of the tetrahedral anions are equally favored. The Hirshfeld surface of encapsulated IO4− (Fig. S7†), together with the close contact analysis of other tetrahedral anions (Tables S5–S8†) reveals that these anions can only form a limited number of interactions with the host cavity wall in a given orientation. Smaller anions like BF4− and ClO4− can form even fewer host–guest interactions simultaneously. The large CF3SO3− anion has two orientations (Fig. 2g and Table S9†). Given that cycHC[8] preserves its conformation almost fully regardless of the guest anion encapsulated, its cavity can be considered as an octahedrally shaped void, which encapsulates guests in a manner depending on their shape, volume and complementarity of the interactions.
Next, we examined anion complexation with cycHC[8] in solution using 1H NMR spectroscopy. The complex with SbF6− showed the largest complexation-induced chemical shift changes (downfield shifts of 0.52, 0.23 and 0.64 ppm for 2ax, 4ax and 6ax, respectively). The signals for 2ax, 4ax and 6ax in the complex with SbF6− also showed significant signal broadening at room temperature, indicating a slow guest exchange rate for SbF6−. The guest exchange slows down at low temperature resulting in the signals of SbF6−@cycHC[8] and excess free cycHC[8] being separate (Fig. 3C(d)).
The Job plot analysis confirmed 1:1 stoichiometry of binding for all studied guests in methanol (ESI†), as also observed by crystallography in the solid state and by mass spectrometry in the gas phase. Association constants were determined by NMR titrations (Table 1), simultaneously following three cycHC[8] proton signals 1ax, 2ax and 3eq (Fig. 3A). The association constant for SbF6− was, due to the broadening of the 2ax signal, determined only from the 1ax and 3eq proton signals. The range of association constants varied over five orders of magnitude, strongly dependent on the size and shape of the guest. The poor water solubility of cycHC[8] prevented measurements in pure water, but no significant decrease in binding strength was observed when pure methanol was changed to a 1:1 methanol/water mixture (Table 1, rows 2–3).
Anion | Cation | V anion (Å3) | PC | K a (M−1) |
---|---|---|---|---|
a Anion volume is based on optimized anion geometries (BP86-D/def2-TZVPD) and calculated using a triangulated sphere model (based on CSD default atomic radii) through the Olex2 program package.45 PC is defined as the ratio between the Vanion to Vcavity(host).42Vcavity(cycHC[8]) = 123.0 Å3 is measured from the crystal structure of cycHC[8].31 b 1H NMR in 1:1 MeOD/D2O. | ||||
SbF6− | Na+ | 81.8 | 0.67 | (2.5 ± 0.7) × 105 |
PF6− | Bu4N+ | 70.6 | 0.57 | (2.8 ± 0.4) × 104 |
PF6−b | Bu4N+ | 70.6 | 0.57 | (2.6 ± 0.2) × 104 |
PF6− | Na+ | 70.6 | 0.57 | (2.0 ± 0.2) × 104 |
ReO4− | Bu4N+ | 64.8 | 0.53 | (4.7 ± 0.4) × 103 |
IO4− | Na+ | 64.3 | 0.52 | (1.8 ± 0.2) × 103 |
ClO4− | Bu4N+ | 54.7 | 0.45 | (4.7 ± 0.2) × 102 |
BF4− | Bu4N+ | 51.6 | 0.42 | (4.8 ± 0.4) × 10 |
CF3SO3− | Bu4N+ | 82.3 | 0.67 | (3.9 ± 0.5) × 10 |
CF3CO2− | Bu4N+ | 68.7 | 0.56 | <10 |
For tetrahedral and octahedral anions, the affinity to the host grows exponentially with the increasing size of the guest, ranging from 48 M−1 for the smallest tested anion BF4− to 250000 M−1 for the largest tested octahedral anion SbF6− (Table 1 and Fig. 4). The anion size dependency for anion binding has also been discussed by Sindelar and co-workers for a bambus[6]uril host in chloroform,24 with the highest selectivity towards ClO4−. Given the double-cone shape of the dodecabenzylbambus[6]uril host and the restricted diameter of the central cavity of 6-membered hemicucurbiturils, it seems that anions larger than perchlorate are bound away from the center of the bambus[6]uril macrocycle, inside the cone-shaped pockets formed by the extending substituents.25,30 With cycHC[8], a single binding site is suggested by the crystal structures, with the anion in each case fully encapsulated in the center of the cavity. The selectivity of this 8-membered host is therefore determined by the parameters of its central cavity. Based on the crystal structures, the correlation between the binding strength and the size of the anion can arguably be ascribed to the number of host–guest interactions an anion can form simultaneously. Thus smaller anions BF4− and ClO4−, which are able to form only one or two C–H⋯anion interactions in a given orientation, are bound with considerably lower affinities.
Fig. 4 The association constant dependency on the anion size and the electrostatic surface potential of the studied anions. Surface potential calculated using Gaussian 09,46 visualized using GaussView5,47 red to blue surface color range spans from −0.2 to 0.2. Pale red dashed lines represent the 68 ± 11 Å3 anion volume range (PC = 0.55 ± 0.09). |
Surprisingly, however, the association constant for the binding of roughly octahedral CF3SO3− is dramatically lower compared to the similarly sized octahedral guest SbF6−, regardless of the several interactions between the host and the encapsulated CF3SO3− apparent in the crystal structure (Table S9†). Given the inherent symmetry of cycHC[8], one can argue that the binding might be stronger with anions having the charge equally distributed over the surface. To assess whether the low binding affinity of CF3SO3− is caused by the asymmetric charge distribution, a control experiment was conducted with CF3CO2−, similar in volume to PF6−, but with a charge distribution resembling that in CF3SO3−. The fact that CF3CO2−, which effectively templates the synthesis of cycHC[8] (in acetonitrile) and has been shown by diffusion NMR to bind to cycHC[8] in chloroform,31 does not bind to cycHC[8] in methanol (Ka < 10), indicates that the binding of anions to cycHC[8] in methanol is sensitive to charge distribution around the anion.
According to the range of optimal packing coefficients, 0.55 ± 0.09, the association constants within the tested group of tetrahedral and octahedral anions with roughly spherical charge distribution should follow a statistical distribution around the optimal guest fit at PC = 0.55. For cycHC[8], with a cavity volume 123.0 Å3, the optimal guest volume should therefore centre in the range of 68 ± 11 Å3 (Fig. 4). However, the strongest association was in fact observed for SbF6− (Vanion = 81.8 Å3), suggesting a shift from the optimal PC to higher values which might be due to the stabilizing effect of host–guest interactions with tighter-fitting guests like SbF6− or the conformational flexibility of the macrocycle itself.
Next, isothermal titration calorimetry (ITC) was used to determine thermodynamic parameters for the complexes of cycHC[8] with SbF6− and PF6− in methanol. The raw thermogram and the binding isotherm for SbF6− are presented in Fig. 5 with the calculated parameters given in Table 2. The Ka values obtained by ITC were comparable with those obtained by NMR spectroscopy. Similarly to NMR, ITC showed higher affinity of cycHC[8] for SbF6− than PF6−, with binding being exothermal, enthalpy driven and entropically disfavored for both. Comparing the two anions revealed that the binding of SbF6− was accompanied by a greater change in both enthalpy and entropy. This is in line with tighter binding of the larger SbF6− anion, giving rise to a stronger host–guest interaction (enthalpic term) but also a greater loss in degrees of freedom (entropic term). Generally, the thermodynamic profile of cycHC[8] binding resembles the behavior of other hemicucurbiturils.37
Fig. 5 Raw thermogram (above) for SbF6− titration with cycHC[8] and the binding isotherm (below) from the integrated thermogram fit using the one-site model. |
Guest | ITC titration | ΔG° | VT-NMR | |||||||
---|---|---|---|---|---|---|---|---|---|---|
Activation parameters for complex formation | Association rate constantak × 103 s−1 | |||||||||
K a (M−1) | ΔH° | TΔS° | ITC titration | NMR titration | DFT calc. | ΔH‡ | TΔS‡ | ΔG‡ | ||
a Determined at 291 K. | ||||||||||
NaSbF6 | (1.02 ± 0.03) × 105 | −56.2 ± 0.3 | −27.7 | −28.5 | −30.8 | — | 38.6 | −13.4 | 51.9 | 2.6 |
NaPF6 | (1.29 ± 0.04) × 104 | −43.8 ± 0.2 | −20.4 | −23.4 | −24.5 | −22 | 42.1 | −5.5 | 47.5 | 17.5 |
The complexation kinetics of normal cucurbiturils has been thoroughly48–56 studied and the exchange of neutral guests has been shown to proceed through a single step,55,56 while the complexation of cationic guests proceeds through a number of intermediates.50–54 To the best of our knowledge, the kinetics of anion binding by hemicucuribiturils has not been previously studied. A particularly interesting aspect of the cycHC[8] behavior is the conformational dynamics of guest encapsulation. Given the bulkiness of the encapsulated guests (Fig. 1), the conformation of the host clearly has to change considerably for the guest to pass through the narrow portals of cycHC[8]. The flexibility of the host is seemingly crucial to the encapsulation and, likewise, to the ejection of the guests, and the reaction pathway of host–guest complex formation was therefore studied computationally. Density functional theory (DFT) calculations were used to model the complexation of cycHC[8] with PF6− utilizing COSMO solvation model for methanol.57,58
By positioning up to four methanol molecules inside the cavity of cycHC[8] we found that, at temperatures above 100 K, a single molecule of methanol is preferably accommodated within the cavity, close to its center. As the association constants derived from NMR titration with sodium and tetrabutylammonium PF6− salts were very similar, cation influence was assumed to be negligible and was not studied. Modelling the exchange of methanol in MeOH@cycHC[8] with PF6− showed the initial formation of a pre-complex with PF6− at the portal of cycHC[8] (MeOH@cycHC[8] + PF6−; Fig. 6). Next, through a transition state involving the dissociation of a hydrogen bond between methanol and cycHC[8], the inclusion complex PF6−@cycHC[8] is formed (the complexation reaction pathway is visualized in a video, ESI†). The DFT-derived Gibbs free energy difference indicates that the methanol-solvated cycHC[8] is 22 kJ mol−1 higher in energy compared to its inclusion complex with PF6−, which is in good agreement with the experimental ΔG values calculated from equilibrium constants obtained by NMR and ITC analyses (Table 2).
Fig. 6 Minimum energy geometries of cycHC[8] complexes with PF6− in MeOH and the associated relative Gibbs free energy values. |
Additional experimental insight into the kinetics of the complexation and the reaction pathway was gained by variable temperature NMR (VT-NMR) studies of SbF6− and PF6−, using a 2:1 host-to-guest ratio (Fig. 7). The complexation reaction order was determined by dilution experiments near the coalescence temperature (241 K and 253 K for SbF6− and PF6−, respectively).
Reaction rates remained constant upon dilution, indicating that the complexation process follows first order kinetics, characteristic for a unimolecular reaction. This suggests that the overall complexation reaction occurs via a low-energy pre-complex. Moreover, the complexation reaction rate constants for both SbF6− and PF6− (Table 2) were determined, allowing us to derive the activation parameters (Table 2) using the Eyring equation (for details, see ESI†). As expected, the complexation rate constant was an order of magnitude higher for PF6− than for the larger SbF6−, also in good agreement with our initial solution studies by NMR. The entropy of activation of the complexation is negative for both anions, as expected for the host and guest forming one host–guest complex. Other entropic contributions that play a role in the host–guest formation reaction are the entropic penalty of orientation of the guest within the macrocycle and the desolvation of the cavity of the host and the collapse of the methanol cluster around the chaotropic guests.
Both the computational and kinetic studies propose the existence of pre-complexes, although the computations suggest them to be higher in energy than the precursors, whereas the VT-NMR results require this pre-complex to be more stable than the starting components. This reflects a complicated complexation reaction pathway that includes several steps. We propose that either the anion pre-complex formation or the reorganization of the methanol solvation shell around anions and the macrocycle play an important role, which merits further investigation.
This work, besides being the first comprehensive anion-binding study on an eight-membered hemicucurbituril-type macrocycle, also demonstrates the unique anion binding properties of a macrocyclic host easily accessible through a simple templated synthetic protocol.31 Moreover, it proposes a pathway to and encourages the preparation of new hemicucurbiturils for anion binding and transport, catalysis in confined space and other supramolecular applications.
Footnotes |
† Electronic supplementary information (ESI) available: MS, NMR, dynamic NMR and computational details and a DFT-based video of complexation. CCDC 1514736–1514741, 1521388. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c6sc05058a |
‡ The name cyclohexylhemicucurbituril, previously used for these macrocycles, is changed in accordance with the IUPAC nomenclature for fused cycles, as the cyclohexane substituents are fused with the parent hemicucurbituril. |
This journal is © The Royal Society of Chemistry 2017 |