Susanne
Löffler
a,
Jens
Lübben
b,
Axel
Wuttke
c,
Ricardo A.
Mata
c,
Michael
John
b,
Birger
Dittrich
*d and
Guido H.
Clever
*a
aFaculty for Chemistry and Chemical Biology, TU Dortmund University, Otto-Hahn-Str. 6, 44227 Dortmund, Germany. E-mail: guido.clever@tu-dortmund.de
bInstitute for Inorganic Chemistry, Georg-August-University Göttingen, Tammannstr. 4, 37077 Göttingen, Germany
cInstitute for Physical Chemistry, Georg-August-University Göttingen, Tammannstr. 6, 37077 Göttingen, Germany
dInstitute for Inorganic Chemistry and Structural Chemistry, Heinrich-Heine University Düsseldorf, Universitätsstr. 1, 40225 Düsseldorf, Germany
First published on 4th April 2016
Substituent control in self-assembled host systems allows for a fine-tuning of structure, dynamics and guest preference. Flat banana-shaped ligands L1 assemble with Pd(II) cations into the interpenetrated coordination cage dimer [3BF4@Pd4L18], capable of sequential guest uptake. In contrast, the introduction of bulky adamantyl groups in ligand L2 prevents dimerization and results in the clean formation of monomeric cage species [Pd2L24]. Owing to steric crowding, the adamantyl substituent is considerably bent sideways with respect to the ligand backbone, and is rapidly flipping between both faces of the free ligand giving rise to two energetically degenerate conformers. Surprisingly, the flipping is preserved in the cage, albeit at a lower rate due to entropic reasons. Despite the very dense packing within the self-assembled structure, the cage is able to encapsulate a series of bis-anionic guests in an induced-fit fashion. Electronic structure calculations revealed a substantial contribution from dispersion interactions between the guest and the surrounding adamantyl groups that stabilize the host–guest complex. Guest exchange kinetics were quantified and the influence that encapsulated guests imparted on the ligand flipping dynamics was examined by a series of 2D NMR experiments. Four synchrotron X-ray structures of the cage and its host–guest complexes are presented, allowing for unprecedented insight into the host–guest interactions of a sterically overcrowded host and its guest-induced distortion.
In the field of metallo-supramolecular self-assembly of discrete rings and cages,22 several examples of dynamic functionality fused to rigid architectures in exohedral positions have been developed by Yang23 and Stang.23,24 Other self-assembled structures were designed containing flexible functionalities as integral parts of their framework. Noteworthy is the tetrahedral cage consisting of six rotaxane edges that has been recently reported by Nitschke, Sanders and Schalley.25 A further degree of complexity is achieved in cages with inbuilt switching mechanisms, controlled by external stimuli such as light. Fujita and coworkers have reported a spherical cage in which the photoisomerization of endohedrally arranged azobenzenes modifies the hydrophobicity of the cavity environment.26 In contrast, Zhou have reported a coordination cage containing azobenzene substituents in exohedral positions, in which the photoswitching affects supramolecular aggregation processes.27 We have recently introduced self-assembled coordination cages based on dithienylethene (DTE) photoswitches in which the inherent flexibility of the ligand backbones can be switched by light with drastic effects on guest binding and assembly size.28 Dynamic host–guest systems have not only been examined with respect to guest affinity and uptake/release kinetics, but also regarding the motility of guests inside the cavity in terms of tumbling or flipping motions.
For example, Raymond and coworkers have studied the dynamics of various benzyl phosphonium guest derivatives when encapsulated in self-assembled tetrahedral cages.29 The group of Schmittel has investigated the tumbling of a DABCO molecule inside a supramolecular cavity flanked with two Zn-porphyrin complexes.7 Similar effects have also been observed in purely organic supramolecular container molecules.30 The processes observed in these non-covalent host–guest systems share similarities with molecular gyroscopes in which a central, dynamic element is covalently bound to a surrounding cage or ring structure.31–33 Steric interactions between the guest and surrounding host structure, as well as the associated solvent molecules, have been identified as important factors that modulate the dynamics of such systems by what may be called molecular friction.31
Several examples of self-assembled cages have been reported, in which steric bulk was systematically incorporated into the framework in order to modulate the guest binding properties. Severin and coworkers have recently reported a family of voluminous ditopic ligands featuring an iron-clathrochelate core that assemble into sterically congested [Re4L4] rings and [Pd6L12] cages.34,35 Yoshizawa has employed a range of metal-mediated and organic capsules using large anthracene panels and studied their aggregation and guest binding properties.36 We have installed steric bulk into the backbone of bis-monodentate ligands to control the dimerization and guest binding selectivity of interpenetrated double cages.37 Furthermore, the Fujita group has shown that also remotely installed steric bulk can influence the uptake of guests inside self-assembled cages.38
The direct consequence of introducing steric bulk into a host system is the reduction of the size of the internal cavity and the width of the portals that allow for entry and release of guest molecules, thereby restricting the size of the encapsulated guest and slowing down their exchange kinetics.39 However, the common belief that steric crowding is always repulsive in nature (and thus generally detrimental to guest binding) seems to be wrong. Even in the absence of functional groups that supposedly dominate intermolecular interactions such as hydrogen bond donors or acceptors, attractive interactions between closely associated molecules are always promoted by London dispersion interactions.40 Though individually weak, dispersive forces can add up substantially in larger systems such as supramolecular assemblies.41 Although dispersion interactions have been known since the 1930's, this insight has only recently become actively discussed again by the broader community, as evidenced by coining the definition of “dispersion energy donors” (DEDs) and the emphasis given to the explicit treatment of dispersion in modern methods of electronic structure calculation.42 This motivated us to synthesize self-assembled host systems equipped with bulky but otherwise unfunctionalized groups in order to study the interplay between molecular crowding, self-assembly and host–guest chemistry.
Fig. 2 1H NMR spectra of the cage and some host–guest complexes. (a) Ligand L2, (b) cage [Pd2L24], (c) [G1@Pd2L24], and (d) [G2@Pd2L24] (400 MHz, 298 K, CD3CN). The signals of the encapsulated guest molecules are highlighted in red. Hi, Hj and Hk represent outside (exo) and Hi′, Hj′ and Hk′ inside (endo) pointing hydrogen atoms of the adamantyl substituent with respect to the acridone backbone (compare Fig. 4a). Empty circle: [K(18-crown-6)]+. |
The high-resolution ESI mass spectrum of the self-assembly product clearly supports the formation of the cage [Pd2L24] by exhibiting a series of species [Pd2L24 + nBF4](4−n)+ (n = 0–2) containing a variable number of BF4− counter anions (Fig. 3a).
Fig. 3 ESI mass spectra of the cage and some host–guest complexes. (a) [Pd2L24], (b) [G1@Pd2L24], and (c) [G2@Pd2L24]. |
1H–1H NOESY NMR experiments on ligand L2 show correlations between the exo (Hi, Hj and Hk) and the endo (Hi′, Hj′ and Hk′) protons of the adamantyl group (Fig. SI 16†). The corresponding cross-peaks intensify with increasing temperature, hence indicating a dynamic movement of the adamantyl group within the ligand structure. As rotation of the adamantyl group can be excluded, due to the presence of the double bond, we propose that a flipping motion occurs in the ligand that allows the endo and exo faces of the adamantyl substituent to swap positions (Fig. 2 and 4a). The rate constant for the flipping was determined from the 2D NMR experiment to be k ≈ 150 s−1 at 298 K (for comparison, all rate constants were determined at the same temperature of 298 K). A comparable case of dynamic exchange between two degenerate conformations was previously described for a structurally related dicyanomethylene-substituted acridone, although the electronic situation in this push–pull system differs considerably from electron-rich ligand L2.47
The formation of the monomeric cage [Pd2L24] was expected to have a tremendous effect on the motility of the adamantyl groups. To our surprise, the flipping dynamics of the adamantyl groups can still be observed in the self-assembled cage, albeit at much lower rates. VT (variable temperature) 1H–1H NOESY NMR experiments were performed for cage [Pd2L24] in order to obtain the flipping rates at different temperatures, from which the activation parameters were calculated by the Eyring plot method (ESI†). The flipping process inside the cage is slowed down significantly with respect to the free ligand, with a rate constant of k ≈ 0.32 s−1 (Fig. 4b, SI 17 and SI 18†). This is mainly caused by an unfavourable entropic contribution (ΔS‡ligand = 11.4 J mol−1 K−1vs. ΔS‡cage = −164.5 J mol−1 K−1). The slower flipping can be explained by the high degree of steric crowding inside the supramolecular self-assembly. Interestingly, from an enthalpic point of view, the flipping inside the cage is favoured over flipping in the free ligand (ΔH‡ligand = 64.0 kJ mol−1vs. ΔH‡cage = 26.9 kJ mol−1). A computational study on the free ligand supported the proposed flipping motion: at ambient temperatures (298 K) the activation barrier was calculated to be ΔG‡298,calc = 56.6 kJ mol−1, comparing well with the experimental value of ΔG298,exp = 60.6 kJ mol−1 extracted from the NMR data. The calculated structures of the conformational minimum and the transition state of the flipping motion are depicted in Fig. 4a (for further details see the Computational section).
To explain the mechanism of the adamantyl flipping motion within the densely packed coordination cage, four possible mechanisms were proposed and investigated: (a) full decomplexation of one ligand, subsequent flipping of all adamantyl groups in the tentative, sterically relaxed [Pd2L23] fragment, followed by recomplexation; (b) one-side decomplexation of one ligand, then flipping of all four adamantyl groups and recomplexation; (c) concerted flipping of all adamantyl groups passing each other simultaneously in the centre of the cage and (d) sequential flipping of each adamantyl group of the intact cage. A series of experiments was conducted in order to deduce the most plausible mechanism. A 1H EXSY NMR measurement of a 3:1 mixture of free ligand L2 and the [Pd2L24] cage showed no exchange between those two species on the NMR time scale (Fig. SI 25†). Therefore, a mechanism following a complete dissociation of one of the ligands (suggestion a) could be excluded. A subsequent experiment showed that the addition of up to five equivalents of pyridine as a competitive ligand to the [Pd2L24] cage does not lead to any significant changes in the rate constant of the flipping motion (Fig. SI 26†). If one-side dissociation of a ligand would be involved in the cage's flipping mechanism, the addition of a competing ligand should increase the flipping rate, since additional donors can be used to accelerate ligand exchange of square-planar metal centres that proceed through an associative mechanism.48 As the pyridine–platinum bond is kinetically much more robust than the palladium analogue,49 also a platinum-based cage was prepared by heating ligand L2 with [Pt2Cl2(CH3CN)2] and AgClO4 for 3 days at 80 °C. Again, 1H–1H NOESY NMR experiments on cage [Pt2L24] show flipping of the adamantyl groups inside the cage, with an exchange rate (k ≈ 0.13 s−1) similar in magnitude to that of the palladium cage (Fig. SI 19†).
Both of the latter two experiments rule out that a partial dissociation pathway is involved in the flipping mechanism (suggestion b). Furthermore, the concerted flipping of all adamantyl groups (suggestion c) is very unlikely for steric reasons (there is not enough space in the interior of the cage to harbour the four adamantyl groups). Therefore, we hypothesize that the flipping is a sequential movement (suggestion d).
Fig. 5 X-ray crystal structures of (a) [Pd2L24], (b) [G1@Pd2L24], (c) [G4@Pd2L24] and (d) [G5@Pd2L24]. Color scheme: C, grey; N, blue; O, red; S, yellow; F, green; B, brown; Pd, tan; Fe, orange. For clarity, solvent molecules and free (non encapsulated) anions are omitted. The ESI and CCDC numbers 1053080–83 contain details of the crystallographic data. |
A case worth discussing in more detail is the encapsulation of the 1,1′-ferrocene bis(sulfonate) guest G1. The 1H NMR spectrum of the host–guest complex (Fig. 2c) revealed that the two signals, assigned to the cyclopentadienyl protons in α and β position relative to the sulfonate substituent, are further split twofold upon encapsulation by the host.50 This splitting pattern arises from the desymmetrization of the C2h-symmetric ferrocene derivative (anti conformation) once transferred into the C4h-symmetric coordination cage. The number of the cage's NMR signals, however, does not change, indicating that the guest is rapidly rotating inside the cavity. In accordance with our previous studies, cyclic voltammetry showed that the Fe(II)/Fe(III) redox-potential is anodically shifted upon encapsulation of the ferrocene inside the cationic cage (Fig. SI 27†). This can be explained by the effect of the tetracationic cage disfavouring the removal of an electron from the encapsulated ferrocene guest.45
Due to the substantial size of the guests, these are in close contact with the surrounding adamantyl residues in the equatorial region of the cavity as can be observed in the X-ray structures (Fig. 5) and also in the NOESY NMR spectra recorded for the host–guest complexes. Fig. SI 22 in the ESI† highlights close contacts between the cage's adamantyl residues and the protons of guest G2. After one equivalent of the guests G1–G5 has been added to the cage, the signals of the empty [Pd2L24] cage disappeared and only signals of the host–guest systems [G1–5@Pd2L24] remained in the 1H NMR spectra. Hence, guests G1–G5 are strongly bound inside the cage (with association constants estimated to be larger than 104 to 105 L mol−1) and the guest exchange with the solvent exterior is slow on the NMR time scale. In contrast, addition of guests G6–G9 to cage [Pd2L24] resulted in a gradual shifting of the cage's 1H NMR signals, indicating a fast exchange of these guests. Since more than one equivalent of guests G6–G9 was required to saturate the cages, these guests bind more weakly to the cage than those previously discussed.51 The different binding behaviour is a result of differences in the guest structures. Guest molecules G1–G5 are bis-anionic and have a good size match for encapsulation,46 while the other guests are either too small (G8 and G9), too large (G7) or lack a second sulfonate group (G6) to be encapsulated with high affinity inside the coordination cage. An NMR exchange spectroscopy (EXSY) experiment of a 50:50 sample of empty cage [Pd2L24] and host–guest complex [G1@Pd2L24] or [G2@Pd2L24] was employed to determine the exchange rates of the guests as k ≈ 7 s−1 for [G1@Pd2L24] and k ≈ 4 s−1 for [G2@Pd2L24] (Fig. 4 and ESI Fig. SI 23 and SI 24†). In comparison, the previously reported monomeric cage composed of a dibenzocycloheptatriene-based ligand L3, which is of similar length and backbone structure but does not have a bulky adamantyl group, shows an exchange rate of k ≈ 90 s−1 with the bis(sulfonate) guest G1 (Fig. SI 15†).52 This led us to conclude that the bulky adamantyl groups of ligand L2 congest the portals of cage [Pd2L24] and slow down the rate of guest exchange. As we discuss below, we hypothesize that they also contribute to the thermodynamic stabilization of the guest molecules inside the cavity.
The double bond between the backbone and the adamantyl substituent itself is not significantly bent, with all six carbon atoms constituting the tetrasubstituted double bond system occupying a common plane. All adamantyl groups are bowed sideways with radial symmetry and occupy the four portals of the cage. This results in a highly packed structure with short inter-ligand distances, measuring only 3.11 Å from proton Ha to proton Hj′ of the neighbouring ligand. The cavity provides a rather spacious environment near the Pd(pyridine)4-planes (where the anions are found), while the bulky adamantyl groups reach into the centre of the cage leaving a narrow channel measuring about 6.3 Å in diameter (see the ESI† for further views of the X-ray structures). In addition, single crystals of the host–guest systems [G1@Pd2L24], [G4@Pd2L24] and [G5@Pd2L24] were obtained via slow vapour diffusion of methanol (for [G1@Pd2L24]), benzene (for [G4@Pd2L24]) and diethyl ether (for [G5@Pd2L24]) into acetonitrile solutions of these complexes.
In all three cases, the relatively large guest molecules were encapsulated inside the cages, with each guest orienting its anionic substituents in close proximity to the Pd(II) cations. This is most pronounced for guests G4 and G5 but less so for shorter ferrocene-based guest G1 that is only able to closely approach one of the palladium atoms (Fig. 5). Interestingly, guest encapsulation forces the adamantyl groups of the ligands to bend even further towards the outside of the cage (Fig. 5b–d). According to DFT computational results, this guest-induced overstretching of the ligand should not cost an energetic penalty of more than 4 kJ mol−1.
Table 1 summarizes a selection of distances extracted from the X-ray structures of the cage and its host–guest complexes. An interesting trend can be observed in this data: the more bulky guests push the adamantyl residues further away from the cavity, which in turn leads to a gradual shortening of the cage along the Pd–Pd axis. In this sense, the cage behaves similar to a flexible container that slightly shrinks from top to the bottom when the cargo is so bulky that it bulges out the vessel to the sides.
[2BF4@Pd2L24] | [G4@Pd2L24] | [G5@Pd2L24] | [G1@Pd2L24] | |
---|---|---|---|---|
a Calculated from the DFT (B3LYP/6-31G*)53 structure of the anions. b As defined by the middle of the line connecting the two palladium atoms. c The longest distance each is highlighted in bold font. | ||||
Mol. volume of the guest [Å3]a | 109.7 | 226.2 | 239.0 | 256.1 |
Distance [Å] between adamantyl sp2-C atom and center of cageb | 5.4 | 6.1 | 6.9 | 6.8 |
5.3 | 5.5 | 5.5 | 5.8 | |
5.2 | 5.0 | 5.4 | 5.7 | |
5.2 | 5.0 | 5.1 | 5.6 | |
Average distance [Å] | 5.3 | 5.4 | 5.7 | 6.0 |
Pd–Pd distance [Å] | 16.2 | 16.2 | 15.9 | 16.0 |
Electronic structure calculations were then carried out in order to obtain further insight into the binding of the guest molecules, in particular the role of dispersion forces. The parent cage containing only BF4− counter anions as well as two host–guest systems were considered, namely [G4@Pd2L24] and [G5@Pd2L24]. First, the contribution of dispersion interactions to the cage stability was computed. Calculations at the SCS-LMP2/aug′-cc-pVTZ level of theory55 gave a dispersion contribution of 14.3 kJ mol−1 for the [Pd2L24] system. The dispersion contribution in the [G4@Pd2L24] system was calculated to be 9.3 kJ mol−1, somewhat smaller given the larger distances between the four ligands that result from the guest-induced expansion of the cage structure. This is not a particularly large effect and although it will counterbalance the steric repulsion of the bulky adamantyl moieties, it should not be a major factor in the cage formation itself.
Interestingly, dispersion forces seem to play a significant role in binding of the guest molecules inside the cages interior. The overall interaction energies between the adamantyl residues and the guests (excluding the cation–anion interactions) for the systems [G4@Pd2L24] and [G5@Pd2L24] were computed as 52.8 kJ mol−1 and 67.2 kJ mol−1. The dispersion energy contributions obtained are 86.3 kJ mol−1 and 64.0 kJ mol−1 for [G4@Pd2L24] and [G5@Pd2L24], respectively. These correspond to a strong contribution in binding, even exceeding the total interaction in the case of [G4@Pd2L24]. In [G5@Pd2L24], the adamantyl units are further apart, reducing the relative weight of the dispersion forces. Consequently, we anticipate the main contributions for the adamantyls' interaction with the guests to be Pauli repulsion and dispersion, the two of opposite sign. These results clearly show that dispersion forces can easily add up to large values in such a supramolecular construct. To better visualize this effect, we have plotted the Dispersion Interaction Densities (DIDs), calculated between the guest and the four adamantane units (Fig. 6c and d; the definition of DID is provided in the ESI†). The images show that the π-systems of the guests interact strongly with the surface of the pocket which is in accordance with the results of the Hirshfeld analysis. In particular, some hot spots are identifiable where the adamantyl groups are in close contact with the encapsulated guests.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 1053080–1053083. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c6sc00985a |
This journal is © The Royal Society of Chemistry 2016 |