Francesca
Corbellini
a,
Fijs W. B.
van Leeuwen
a,
Hans
Beijleveld
a,
Huub
Kooijman
b,
Anthony L.
Spek
b,
Willem
Verboom
a,
Mercedes
Crego-Calama
*a and
David N.
Reinhoudt
*a
aLaboratory of Supramolecular Chemistry and Technology, MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE, Enschede, The Netherlands. E-mail: m.cregocalama@utwente.nl; Fax: +31 53 489 4645; Tel: +31 53 4892984
bDepartment of Crystal and Structural Chemistry, Utrecht University, H. R. Kruytgebouw, N-827 Padualaan 8, 3584 CH, Utrecht, The Netherlands. E-mail: h.kooijman@chem.uu.nl; Fax: +31 30 2533940; Tel: +31 30 2532533
First published on 8th December 2004
The formation and characterization of molecular capsules resulting from the self-assembly between oppositely charged (thia)calix[4]arenes in polar solvents like MeOH and MeOH–H2O are reported. The multiple ionic interactions allow the self-assembly of the complementary (thia)calix[4]arenes 1 and 2a–d into 1∶1 complexes as revealed by 1H NMR and mass spectrometry (ESI-MS). Isothermal titration calorimetry (ITC) was used to determine the association constants, which, depending on the ionic groups involved in the complexes, vary between 103 and 106 M−1. An X-ray structure of the assembly 1•2d was also obtained. Unlike in solution, in the solid state 1•2d forms a 1∶1 three-dimensional network in which 2d adopts a 1,2-alternate conformation.
Ionic interactions are often employed as an important attractive force in both biological9 and artificial molecular recognition.10 Nevertheless, these strong interactions have been only marginally used for building molecular capsules. Surprisingly, before our initial studies only a few examples of stable cage-like complexes using multiple ionic interactions have been reported in the literature and only very recently a few other publications have followed these first examples.11
Previously, we proved the effectiveness of multiple ionic interactions as an alternative to hydrogen bonds and metal–ligand coordination for the synthesis of noncovalent supramolecular capsules in aqueous solutions. Calix[4]arenes functionalized at their upper rim with oppositely charged groups form molecular capsules in polar organic solvent where guest encapsulation can be observed.12 In parallel with our study, the synthesis of molecular capsules based on ionic interactions between oppositely charged calix[4]- and calix[6]arenes has been published by the group of Schrader.13
Here, a more extensive study is presented on the use of ionic interactions for the synthesis of supramolecular capsules in polar solvents. We report that strong association between the oppositely charged calix[4]arenes 1 and 2a–c, functionalized at the upper rim with amidinium and either sulfonate, carboxylic acid and phosphonate moieties, respectively, results in the formation of 1∶1 molecular capsules (Chart 1). Capsule formation between 1 and the larger and less preorganized thiacalix[4]arene 2d has also been investigated. The four sulfur atoms in the thiacalix[4]arene 2d provide an enlarged skeleton compared to the calix[4]arene scaffold.14 As suggested by molecular modelling studies, the resulting capsule 1•2d is suitable for the encapsulation of larger guest molecules.
![]() | ||
Chart 1 |
Crystal data for 1•2d: C60H92N8O8·C24H12O16S8·8C2H6OS·0.776H2O, Mr
= 2505.39, colorless, block-shaped crystal (0.1 × 0.2 × 0.3 mm), triclinic, space group P
(no. 2) with a
= 18.061(2), b
= 20.378(2), c
= 25.193(4)
Å, α
= 69.128(15)°, β
= 69.681(15)°, γ
= 69.840(18)o, U
= 7862.2(17)
Å3, Z
= 2, Dc
= 1.058 g cm−3, F(000)
= 2656, μ(Mo Kα)
= 0.279 mm−1, 155874 reflections measured, 28453 independent, Rint
= 0.0831.†
![]() | ||
Scheme 1 Reagents and conditions: (a) Et2AlCl, 1-propylamine, fluorobenzene, 5 days, 80 °C, yield 65%; (b) t-BuLi, CO2, 30 min, rt, HCl 6M, yield 72%; (c) p-HCOH, HCl, CH3COOH, H3PO4, 1,4-dioxane, 3 days, reflux, 45%; (d) (EtO)3P, 24 h reflux; (e) TMSBr, 12 h rt, then MeOH, yield 54%. |
The synthesis of compound 6 was carried out in an alternative way to the ones reported in the literature,18,22 following a similar route as described for the synthesis of the analogous 5,11,17,23-tetrakis(chloromethyl)-25,26,27,28-tetrakis(2-ethoxyethoxy)-calix[4]arene (1,3-alternate).21
Reaction of the chloromethylated calix[4]arene 6 with triethylphosphite, under Arbuzov conditions, and subsequent hydrolysis with TMSBr afforded calix[4]arene 2c in 54% yield (Scheme 1).
Compounds 1 and 2a–c give sharp 1H NMR spectra in CD3OD, showing that they do not aggregate in solution. Moreover, the presence of one pair of doublets for the bridge methylene protons clearly reveals the C4v symmetry common to tetrasubstituted cone calix[4]arenes. Unlike the analogous tetrasulfonatecalix[4]arene 2a, which is fixed in the cone conformation by the alkylation of the phenolic oxygens, tetrasulfonatethiacalix[4]arene 2d can exist in solution in four different conformations. However, according to density functional theory calculations of the structure and conformational equilibrium of thiacalix[4]arene, the most stable conformer is the cone, which is stabilized by a cyclic array of hydrogen bonds at the lower rim.28 The presence of only a single peak (∼8 ppm) for the aromatic protons in the 1H NMR spectrum of the thiacalix[4]arene 2d in CD3OD indicates indeed the cone conformation.
Interestingly, the proton signals of the amidinium propyl chains of compound 1, Hα, Hβ and Hγ, shifted upfield upon assembly formation (Fig. 1 and Table 1). In analogy to what was observed before,12,29 this upfield shift is attributed to the inclusion of the propyl side chain of calix[4]arene 1 in the cavity of the formed capsule 1•2.
![]() | ||
Fig. 1 1H NMR spectra (CD3OD, 298 K) of (a) 2a, (b) 1 and (c) capsule 1•2a obtained by precipitation in H2O upon mixing solutions of 1 and 2a. (*) indicates solvent signals. (For proton labelling see Chart 1.) |
For assembly 1•2a only small changes were observed for all the other signals (Δδ < 0.1 ppm), while for the assembly 1•2b downfield changes were observed also for the methylene bridge protons (Δδ = 0.12 ppm). Downfield shifts accompanied by broadening were also detected for the first CH2 group of the ethoxyethyl chains at the lower rim of 2b (Δδ = 0.15 ppm). Assembly 1•2c shows only minor changes in the resonances of the propyl side chain protons (Table 1). Most likely the hindrance caused by the bulkier phosphonate groups precludes the access of the propyl chain to the cavity of the calix[4]arene 1. Nevertheless, small changes are detectable for the signals of the aromatic protons and the methylene bridge hydrogens of 2c (Δδ ∼ 0.05 ppm). Analogously to what is observed for the assemblies based on calix[4]arenes, upfield shifts in the resonances of the protons of the amidinium chains of 1 (Table 1) and of the aromatic protons of both building blocks [Δδ (1) = 0.23 ppm, Δδ (2d) = 0.11 ppm] were observed for assembly 1•2d. This is a good indication that a molecular capsule is formed where one of the alkyl chains of 1 is pointing inside the internal cavity.
Crystallization of assembly 1•2d in DMSO–MeOH provided single crystals suitable for X-ray diffraction studies. In the solid state 1 and 2d form an infinite three-dimensional network. The unit cell contains, besides 1 and 2d in a 1∶1 ratio, eight molecules of DMSO, one molecule of water and a number of completely disordered solvent molecules, which are located in a region around the crystallographic inversion centre. Thiacalix[4]arene 2d adopts a 1,2-alternate conformation, thus presenting two adjacent sulfonate moieties on each side of the plane. Calix[4]arene 1 adopts instead a pinched-cone conformation in which the two opposite aromatic rings are bent inwards and the other two aromatic rings are forced to bend outwards. A very complicated structure is formed in which hydrogen bonds are formed either between the two building blocks or between one of the building blocks and solvent molecules. Fig. 2(a) depicts one type of interaction where two opposite amidinium groups of 1 form hydrogen bonds to one of the sulfonate groups of the thiacalix[4]arene 2d. Fig. 2(b) shows hydrogen bonds between two adjacent amidinium groups of 1 and two sulfonate moieties of 2d. Interestingly, the 1,2-alternate conformation in 2d is stabilized by bifurcated hydrogen bonds, involving an intramolecular hydrogen bond between the hydroxyl group and a sulfur atom of 2d as well as an intermolecular bond donated to the oxygen of a DMSO molecule (not shown).
![]() | ||
Fig. 2 Portions of the X-ray structure of the 3D network showing the hydrogen-bonding interactions between calixarene 1 and thiacalix[4]arene 2d. See text for details of views (a) and (b). |
This unexpected result prompted us to further investigate the assembly of 1•2d in solution. According to the X-ray structure, compound 2d in 1•2d adopts a 1,2-alternate conformation. If the same conformation exists in solution, two signals for the aromatic protons of 2d would be expected, while, for symmetry reasons, one signal would be expected for the cone conformation. As mentioned before, in CD3OD solution at 25 °C, the 1H NMR spectrum of 1•2d shows a single resonance for the aromatic protons of 2d, which accounts for a symmetrical conformer (cone). However, the singlet could also be an averaged signal of different conformers that are rapidly interconverting. Temperature-dependent 1H NMR spectra of 1•2d were recorded in CD3OD. At −50 °C a singlet for the aromatic protons was still observed while a further decrease in the temperature resulted in the precipitation of the thiacalix[4]arene 2d. This result indicates that if different conformers are present at 25 °C, at −50 °C the flexibility of 2d is still too high for the individual conformers to be revealed by NMR spectroscopy, which is not very likely, thus confirming the cone conformation.
To study the effect of the concentration on the self-assembly of 1 and 2d, 1H NMR studies in DMSO-d6, in which a better solubility of the complex is observed, were performed. The 1H NMR spectrum of a 20 mM solution of 1•2d in DMSO-d6 showed no shifts in the propyl chains of the amidinium groups and only small upfield changes for the –NH signals of 1 (Fig. 3).
![]() | ||
Fig. 3 1H NMR spectra (DMSO-d6, 298 K) of (a) 2d, (b) 1 and (c) 1•2d (c = 20 mM). (*) indicates solvent signals. (For proton labelling see Chart 1.) |
At a much lower concentration (1 mM), a DMSO-d6 solution of 1•2d showed a different spectrum, which is analogous to the one observed in CD3OD (spectra not shown). The signals of the protons of the propyl side chain of 1 were all upfield-shifted (Table 1). Furthermore, significant upfield shifts were also observed for the –NH signals of 1 (upfield shifts, Δδ = 0.62, 0.64, 0.30 ppm). These results show that the concentration somehow affects the geometry of the assembly. The structure obtained in the solid state could exist in concentrated solutions of the two building blocks, while at lower concentration (1 mM) a well-defined complex is most probably formed. This assumption is corroborated by analysis of the assembly 1•2d by 2D NMR where NOE connectivities between the aromatic protons of 2d and the protons of the propyl amidinium chain of 1 were found. Perhaps in the solid state the molecules organize by optimizing both space and intramolecular interactions (Fig. 2), while in solution the flexible thiacalix[4]arene 2d assembles to give a well-defined 1•2d molecular capsule.
Equimolar solutions of the different assemblies were analyzed by electrospray mass spectrometry (ESI-MS) and in all cases the presence of a doubly charged peak for the capsule was observed. Peaks for the isolated building blocks were also detected. The spectrum of an equimolar solution of 1 and 2a in MeOH showed the doubly charged signals of the capsule at m/z 1064.8 corresponding to [(1•2a) + 2Na]2+. A peak at m/z 717.5 for [(1•2a) + 3Na]3+ was also present. The spectrum of an equimolar solution of 1 and 2b in MeOH showed the doubly charged signals of the capsule at m/z 969.7 and m/z 987.8 corresponding to [(1•2b) + 2H]2+ and [(1•2b) + H + Na]2+, respectively, together with a peak at m/z 1049.5 for [1 − 4HCl + H]+ and two peaks at 911.2 and 933.8 for [2b + Na]+ and [(2b − H + Na) + H]+, respectively (Fig. 4).
Analogously, the measured ESI-MS spectrum of an equimolar solution of 1 and 2c in MeOH showed the [(1•2c) + Na + H]2+ and [(1•2c) + 2Na]2+ ions at m/z 1079.5 and m/z 1090.5 together with signals at m/z 1049 and 911.3 for [1 − 4HCl + H]+ and [2c + Na]+, respectively. For assembly 1•2d analyses at different concentrations and using different mass spectrometric techniques were performed. A 0.5 mM solution of 1•2d in MeOH was analyzed by ESI-MS. The spectrum showed a peak at m/z 955.4 corresponding to [(1•2d) + 2Na]2+. The FAB-MS spectrum showed the presence of a peak at m/z 1865.9 corresponding to [(1•2d) + H]+. The same results were obtained from the ESI-MS and FAB-MS of both a 1 mM and 20 mM solution of 1•2d in DMSO. No evidence for high m/z aggregates was observed. In agreement with the results obtained using ESI-MS and FAB-MS, the MALDI-TOF spectrum of a 20 mM solution of 1•2d in DMSO showed only the presence of the expected 1∶1 complex.
Isothermal titration calorimetry (ITC) was used to study the thermodynamics of the association between 1 and 2a–d. The titrations were carried out by adding aliquots of calix[4]arene tetraamidinium 1 to a solution of calix[4]arenes 2a–d at 298 K in MeOH–H2O (xwater = 0.4) in the presence of either tetrabutylammonium perchlorate (Bu4NClO4, I = 0.01 M) or borate buffer (Na2B4O7, I = 0.01 M and I = 0.03 M).
The titration curves showed an inflection point after addition of 1.0 equiv. of 1, which confirms the formation of complexes with a 1∶1 stoichiometry (Fig. 5). The titration data were fitted to a 1∶1 binding model using a nonlinear least-squares fitting procedure. The corresponding association constants and thermodynamic parameters are listed in Table 2.
![]() | ||
Fig. 5 Calorimetric titration of 2a (0.1 mM) with 1 (1.0 mM) in MeOH–H2O (xwater = 0.4) at 298 K in the presence of Bu4NClO4, I = 0.01 M (left), and Na2B4O7, I = 0.01 M, pH = 9.2 (right). (a) and (c) show the heat evolution per injection of 1; (b) and (d) are the resulting binding curve and best fit curve. |
Assembly | K a/M−1 | ΔH°/kJ mol−1 | ΔS°/J K−1 mol−1 |
---|---|---|---|
a Background electrolyte: 1.0 × 10−2 M Bu4NClO4 (I = 0.01 M). b Background electrolyte: 0.33 × 10−2 M Na2B4O7 (I = 0.01 M, pH = 9.2). c Background electrolyte: 1.0 × 10−2 M Na2B4O7 (I = 0.03 M, pH = 9.2). d Could not be determined because of the small heat effect. | |||
1•2aa | (8.5 ± 1.4) × 106 | 14.1 ± 0.1 | 180 ± 2 |
1•2ab | (1.2 ± 0.1) × 106 | −69.9 ± 0.2 | −118 ± 1 |
1•2ac | (1.1 ± 0.1) × 105 | −58.2 ± 0.4 | −98 ± 2 |
1•2ba | (8.6 ± 1.9) × 106 | 17.1 ± 0.2 | 190 ± 2 |
1•2bb | (1.5 ± 0.1) × 106 | −57.9 ± 0.1 | −76 ± 1 |
1•2bc | (1.8 ± 0.1) × 105 | −49.7 ± 0.4 | −66 ± 2 |
1•2ca | —d | —d | —d |
1•2cb | (4.5 ± 0.2) × 104 | −55.3 ± 0.1 | −96 ± 1 |
1•2cc | (6.7 ± 0.1) × 103 | −40.9 ± 0.3 | −63 ± 1 |
1•2da | (1.0 ± 0.1) × 107 | 25.9 ± 0.1 | 221 ± 1 |
In the presence of 1.0 × 10−2 M (I = 0.01 M) of Bu4NClO4 as background electrolyte, association constants in the order of 106–107 M−1 for 1•2a and 1•2b were determined. The association process is strongly entropy driven, most likely due to the release of highly ordered solvent molecules into the bulk solvent. The positive values obtained for the enthalpies are ascribed to the energy needed to desolvate the charged groups, which overrides the negative enthalpic contribution due to the formation of the ionic interactions of the assemblies.30 The binding curve obtained for assembly 1•2c displayed only a very small endothermic effect that precluded accurate curve fitting. This was attributed to possible intramolecular interactions between the phosphonate groups of 2c. These interactions could compete with the assembly of the molecular capsule. The formation of capsule 1•2d shows the most favorable entropic contribution and the least favorable enthalpic term. Since the interaction between sulfonate and amidinium groups should be similar for calix[4]arene 2a or thiacalix[4]arene 2d with 1, the more positive values observed for both ΔH° and ΔS° must be the result of a higher degree of solvation of the larger thiacalix[4]arene 2d.
ITC measurements for assemblies 1•2a–c were also performed in the presence of borate buffer (Na2B4O7). The titrations were carried out at two different concentrations of borate buffer (c = 1.0 × 10−2 M and c = 0.33 × 10−2 M), leading to ionic strengths of I = 0.03 M and I = 0.01 M, respectively. Compared to 1•2a and 1•2b, assembly 1•2c shows lower association constants (Ka = 4.5 × 104 M−1 or Ka = 6.7 × 103 M−1, Table 2). This is most probably the result of a lower degree of preorganization of 2c due to the freedom of rotation around the C–P bond. Probably the formation of the capsule causes a reduction in the torsional entropy that overrides the gain in entropy, caused by the desolvation of the charged groups as reflected in the negative values for ΔS°.31
The enthalpograms for titrations performed in the presence of the two different background electrolytes (Bu4NClO4vs. borate buffer) but at the same ionic strength (I = 0.01 M) show opposite signs of the heat effects, which indicates a change in the thermodynamics of association. Going from Bu4NClO4 to borate buffer the overall process occurring in solution changes from endothermic and entropy-driven to exothermic and enthalphy-driven with the overall curve showing a negative contribution both from ΔH° and ΔS° (Fig. 5). It is well-known that the heat effect varies depending on the particular electrolyte employed. Most probably the different solvation modes of the charged moieties are of importance for the differences in the thermodynamics of binding.
More interesting is the variation of the heat effect depending on the specific electrolyte concentration in solution. The association constant decreases on increasing the ionic strength. For example, for assembly 1•2aKa decreases from 1.2 × 106 M−1 at I = 0.01 M to Ka = 1.1 × 105 M−1 at I = 0.03 M. Increasing the ionic strength in general weakens the electrostatic interactions. When a higher concentration of salt is present in solution, the Coulombic attractive forces between the two building blocks are considerably shielded and the binding is weakened. This is reflected in a less negative value of ΔH°. Although the observed Ka decreases upon increasing the salt concentration, a change in the mode of binding is not expected as no change in the stoichiometry of the interaction was observed.
Footnote |
† CCDC reference number 253844. See http://www.rsc.org/suppdata/nj/b4/b412409j/ for crystallographic data in .cif or other electronic format. |
This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2005 |