Manuel
Weh‡
a,
Asja A.
Kroeger‡
bcd,
Olga
Anhalt
e,
Amir
Karton
*bf and
Frank
Würthner
*ae
aInstitut für Organische Chemie, Universität Würzburg, Am Hubland, 97074 Würzburg, Germany. E-mail: wuerthner@uni-wuerzburg.de
bSchool of Molecular Sciences, The University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia
cResearch School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
dInstitute for Nanoscale Science & Technology, Flinders University, Adelaide, South Australia 5042, Australia
eCenter for Nanosystems Chemistry, Bavarian Polymer Institute, Universität Würzburg, Theodor-Boveri-Weg, 97074 Würzburg, Germany
fSchool of Science and Technology, University of New England, Armidale, NSW 2351, Australia. E-mail: amir.karton@une.edu.au
First published on 27th November 2023
Corannulene is known to undergo a fast bowl-to-bowl inversion at r.t. via a planar transition structure (TS). Herein we present the catalysis of this process within a perylene bisimide (PBI) cyclophane composed of chirally twisted, non-planar chromophores, linked by para-xylylene spacers. Variable temperature NMR studies reveal that the bowl-to-bowl inversion is significantly accelerated within the cyclophane template despite the structural non-complementarity between the binding site of the host and the TS of the guest. The observed acceleration corresponds to a decrease in the bowl-to-bowl inversion barrier of 11.6 kJ mol−1 compared to the uncatalyzed process. Comparative binding studies for corannulene (20 π-electrons) and other planar polycyclic aromatic hydrocarbons (PAHs) with 14 to 24 π-electrons were applied to rationalize this barrier reduction. They revealed high binding constants that reach, in tetrachloromethane as a solvent, the picomolar range for the largest guest coronene. Computational models corroborate these experimental results and suggest that both TS stabilization and ground state destabilization contribute to the observed catalytic effect. Hereby, we find a “mutual induced fit” between host and guest in the TS complex, such that mutual geometric adaptation of the energetically favored planar TS and curved π-systems of the host results in an unprecedented non-planar TS of corannulene. Concomitant partial planarization of the PBI units optimizes noncovalent TS stabilization by π–π stacking interactions. This observation of a “mutual induced fit” in the TS of a host–guest complex was further validated experimentally by single crystal X-ray analysis of a host–guest complex with coronene as a qualitative transition state analogue.
Aside from the synthetic fine-tuning of the bowl depth, some groups have investigated the impact of supramolecular host–guest encapsulation on corannulene's inversion process. In this context, Nitschke and co-workers recently showed how a curved host and second guest cooperatively inhibit the dynamic motion of complexed corannulene.13 In contrast, Siegel and co-workers showed that supramolecular complexation in a suitable box-shaped cyclophane could be used to stabilize the planar geometry of corannulene's inversion TS, which accelerates the rate of bowl-inversion,14,15i.e. supramolecular catalysis.16,17 Thus, in 2014 the authors used the corannulene inversion catalyzed by ExBox4+, a tetracationic cyclophane consisting of 1,4-phenylene-extended bipyridinium moieties bridged end-to-end by para-xylylene units, to provide a conceptually simple supramolecular example of induced fit catalysis: catalysis originates from the reorganization of the host from a strained conformation in the reactant complex to an energetically favorable geometry in the complex with the complementary planar transition structure of the guest.14 It has further been shown computationally that adsorption onto 2D materials such as graphene may catalyze the bowl-inversion in a similar fashion.18–20 These results show that the kinetics of the inversion process can be controlled both by synthetic modifications and by supramolecular complexation, and that the geometry of corannulene in the ground state can be adjusted by suitable metal coordination or substitution. However, there is still no example of a non-planar transition structure of corannulene's bowl-to-bowl inversion.
We recently reported the synthesis of a PBI cyclophane composed of chiral chromophores,21 whose chemical structure is depicted in Fig. 1. This macrocycle proved to be an excellent host molecule for the encapsulation of small carbohelicenes due to an optimal shape complementarity between the helically twisted PBI units and the guests. In-depth studies showed that this host is able to lower the enantiomerization barrier of [5]helicene by acting as a template catalyst.21,22 This example of an intricate process proceeding via a non-planar TS has drawn our attention to the possible influence of this cyclophane template, composed of rather rigid non-planar PBIs, on the planar TS of corannulene's bowl-to-bowl inversion. Based on the seminal work of Juríček et al. on induced fit catalysis,14 the frequently discussed concepts of anion-π23,24 and cation-π25 catalysis, as well as on theoretical work on π–π catalysis26 exerted by PBI cyclophanes with planar chromophores on the stereoinversions of polycyclic aromatic hydrocarbons (PAHs) with planar TSs,27,28 we were interested in the impact of this host's comparably rigid, non-planar and therefore non-complementary geometry on the dynamic bowl-to-bowl inversion of corannulene.
For the larger guests, i.e. perylene, benzo[ghi]perylene and coronene, fluorescence titration studies were performed. While only a fluorescence turn-off is revealed upon the addition of perylene (20 π-electrons, Fig. S4†), an exciplex band emerges for benzo[ghi]perylene complexation (22 π-electrons, Fig. S5†). This band is further increased in intensity for coronene (24 π-electrons, Fig. 2b and S6†) as the largest guest, supporting a significant charge transfer character for these complexes. The corresponding complexation energies fit the expected linear increase with a larger size of the guests' aromatic surfaces (Fig. 2d). For coronene, the binding constant in chloroform reaches a remarkably high value of Ka = 2.8 × 108 M−1. Upon decreasing the solvent polarity, the binding constants can be further increased for this host if chloroform is replaced by tetrachloromethane as a solvent. For anthracene and pyrene, the association constants could be increased by three orders of magnitude (Fig. S8 and S9†). If a linear free energy relationship is assumed for all PAH guests for the transition from chloroform to tetrachloromethane, the Gibbs free energy can be extrapolated for the larger guests (Fig. 2d). This yields a complexation Gibbs free energy of 66.9 kJ mol−1 for coronene in tetrachloromethane, corresponding to an impressive binding constant of around 1012 M−1 in this organic solvent. Table 1 summarizes the binding constants for the guest complexations by 1-PP and the corresponding Gibbs free energies in chloroform and tetrachloromethane.
Guest | K a (CHCl3) [M−1] | −ΔG0 (CHCl3) [kJ mol−1] | K a (CCl4) [M−1] | −ΔG0 (CCl4) [kJ mol−1] |
---|---|---|---|---|
a Numbers are estimated based on linear extrapolation of the Gibbs free energy in Fig. 2d and calculation of the binding constants according to . | ||||
Anthracene | 3.4 × 104 | 25.6 | 3.3 × 107 | 42.5 |
Pyrene | 1.3 × 105 | 28.9 | 2.0 × 108 | 46.9 |
Triphenylene | 6.1 × 105 | 32.7 | 1.7 × 109a | 52.1a |
Perylene | 4.0 × 107 | 42.9 | 1.2 × 1010a | 57.0a |
Benzo[ghi]perylene | 1.5 × 108 | 46.2 | 9.5 × 1010a | 62.0a |
Coronene | 2.8 × 108 | 47.7 | 7.0 × 1011a | 66.9a |
Corannulene | 2.4 × 103 | 19.1 | 1.4 × 106 | 34.7 |
Thus, despite the apparent structural mismatch between the helically chiral cavity of the host and the planar guests, we find extraordinarily high binding affinities for the complexes with planar PAHs. In contrast, curved corannulene (20 π-electrons) shows a much weaker binding by 1-PP. The optical changes in the corresponding UV/vis titration study (Fig. 2c and S7†) are qualitatively comparable to the one of triphenylene (Fig. S3†), which has 18 π-electrons. However, with a binding constant of Ka = 2.4 × 103 M−1 in chloroform at 295 K, the complexation energy is by 13.6 kJ mol−1 lower, despite the larger aromatic surface of corannulene. Indeed, based on the linear relationships shown in Fig. 2 we can calculate Gibbs free energies for a planar 20 π-electron PAH guest of 42.9 kJ mol−1 in chloroform and 57.0 kJ mol−1 in tetrachloromethane that are significantly larger than those observed for the bowl-shaped corannulene (Table 1).
Fig. 3 Structures of free corannulene and corannulene⊂1-PP optimized at the B3LYP-D3(BJ)/def2-SVP level of theory along with a zoom-in on the guest to illustrate its compression upon complexation. |
The studies on the guest binding by 1-PP revealed that high binding energies can be achieved for the encapsulation of flat PAHs despite the twisted nature of the PBI units in the cyclophane. To obtain further insights into the structures of these complexes, we grew a co-crystal from a mixture of coronene, the most strongly bound guest, and 1-PP by slow evaporation of n-hexane into a chlorobenzene solution. The resulting crystal could be analyzed by single crystal X-ray analysis (Fig. 4a). Crystallographic analysis showed that the complex crystallizes in the triclinic system (for details see the ESI†). The crystal structure not only clearly confirms the intercalation of coronene between the twisted PBIs, but also reveals a distorted geometry of flat coronene within the complex. Furthermore, the twist of the PBI's naphthalene subunits shows a flattening of the perylene core with a small twist angle of 12.8° (average from both chromophores of the two complexes in the unit cell). Hence, the loss of planarity of coronene is accompanied by a reduction in the PBI core twists. The intermolecular contacts between the twisted perylene units and the distorted coronene guest were further qualitatively investigated with the help of NCI (noncovalent interactions) analysis46 of the co-crystal (Fig. 4b).
The green regions in the NCI plot can be attributed to van der Waals interactions between the host and the guest and illustrate how the induced structural adaption in this complex optimizes noncovalent interactions between the aromatic planes of the guest and the perylene bisimides. The additional stabilizing CH–π interactions between coronene and the para-xylylene spacer units further become apparent. When inspecting the NCI plot from a top-down perspective (Fig. 4b, right), the intermolecular binding iso-surface also reveals that almost the full aromatic surface of coronene contributes to the host–guest binding. This explains the continuous increase of the complexation Gibbs free energy with substrate size even beyond perylene, which could otherwise be presumed as the perfect guest for this host.
To study the dynamic bowl-to-bowl inversion process, the alkyl substituent of ethylcorannulene proved a good diagnostic group in order to quantify the inversion of corannulene as the methylene protons become heterotopic at low temperatures, when the bowl-to-bowl motion is slow on the NMR time scale.14 Accordingly, variable temperature NMR studies were performed in deuterated dichloromethane (DCM-d2), which due to its low melting point is more suitable than chloroform or tetrachloromethane to cover the whole dynamic process. We note that the additional ethyl substituent had no impact on the binding affinity, which was tested by host–guest titration and subsequent data evaluation (Fig. S11 and S12†). We further note in this context that the UV/vis absorption spectrum of the complex in tetrachloromethane is in good agreement with the one in DCM (Fig. S14†), which indicates that no structural changes in the host–guest assembly are induced by the less polar solvent but only weak solvatochromism is observed.
The binding affinity for corannulene is weaker in DCM than in chloroform. This was also confirmed for the most strongly bound guest coronene (Fig. S13†), indicating that this is likely to be a general solvent effect. As only one set of signals results for the host protons even at temperatures below 170 K for a host-ethylcorannulene mixture with 2 eq. of 1-PP (Fig. S19†), guest exchange is fast compared to the bowl-to-bowl inversion over the whole temperature range needed for the studies on this process. Hence, the bowl-to-bowl motion is the rate-determining step which allows direct conclusions on the guest inversion from the dynamic NMR data not only in the absence but also in the presence of the host template and we expect the solvent effect to not affect the studies on the kinetics of interest.
At room temperature, a quartet for the methylene protons of the free substrate is apparent, which becomes more and more broadened upon cooling down until the coalescence at 233 K is reached. At lower temperatures these protons split into two signals. For free ethylcorannulene (Fig. 5a and S15†), an inversion barrier of 45.8 kJ mol−1 could be derived from the coalescence temperature and the signal splitting according to eqn (S1†), which is in very good agreement with the literature value for corannulene.8 A similar temperature series was carried out in the presence of 0.5 eq. of 1-PP (Fig. 5b and S16†), revealing this barrier to be decreased by around 25% to 34.2 kJ mol−1 with a coalescence temperature of 176 K for the methylene protons of the guest. This decrease of the barrier towards inversion was further supported by VT experiments with more (Fig. S17†) or less (Fig. S18†) than 0.5 equivalents of the host, which had only minor impact on the catalytic effect: while two equivalents of 1-PP lowered the barrier to 33.7 kJ mol−1, supporting the fact that the guest exchange is fast compared to the bowl-to-bowl inversion (vide supra), we determined a slightly higher barrier of 36.8 kJ mol−1 with 0.1 equivalents of host. These changes might be attributable to the slow diffusion of the reactants at very low temperatures close to the freezing point of dichloromethane, when the solvent gets more viscous.47
Fig. 5 Excerpt from variable temperature 1H NMR of ethylcorannulene in the (a) absence and in the (b) presence of 0.5 equivalents of 1-PP (c (1-PP) ≈ 1 mM, exact temperatures are given in the ESI†). |
In order to rationalize the observed acceleration of the bowl-to-bowl inversion within the cyclophane host, we performed computational studies of the uncatalyzed and host–guest catalyzed processes at the SMD-(CCl4)-PW6B95-D3(BJ)/def2-TZVPP//B3LYP-D3(BJ)/def2-SVP level of theory.37–45Fig. 6 shows the resulting schematic potential energy profiles for the two processes, as well as a breakdown of the contributions to relative noncovalent stabilizations and destabilizations of transition structures obtained from second generation absolutely localized molecular orbital (ALMO) energy decomposition analyses (EDA).48
Fig. 6 (a) Schematic potential energy surface (PES) of the uncatalyzed and the catalyzed bowl-to-bowl inversion of corannulene in the absence and the presence of 1-PP in CCl4 with Gibbs free energies given relative to the free host and guest. Geometries and energies were obtained at the SMD-(CCl4)-PW6B95-D3(BJ)/def2-TZVPP//B3LYP-D3(BJ)/def2-SVP level of theory. (b) Calculated transition structure for the bowl-to-bowl inversion of corannulene and corannulene⊂1-PP. The 1,7-bridging motif in the cyclophane is cut off for clarity. (c) Breakdown of relative noncovalent stabilizations and destabilizations of TSs (ΔΔE = ΔETS − ΔERC) from ALMO-EDA-II analysis48 at the ωB97M-V/def2-SVP level of theory into contributions from Pauli repulsions (ΔΔEPAULI), electrostatic interactions (ΔΔEELEC), dispersion interactions (ΔΔEDISP), polarization (ΔΔEPOL), and charge transfer (ΔΔECT). |
For the bowl-inversion of free corannulene, we obtain a Gibbs free energy barrier of 47.7 kJ mol−1, which is in good agreement with the experimentally obtained value of 45.8 kJ mol−1 (Fig. 6a). Further, in qualitative agreement with the experimental data we find complexation of corannulene with 1-PP to lead to a notable reduction in the bowl-to-bowl inversion barrier. We note that in comparison to the experimentally determined bowl-inversion barrier upon complexation of 34.2 kJ mol−1, we computationally obtain a lower barrier of 21.4 kJ mol−1, indicating that computationally we predict a larger catalytic effect. This discrepancy between theory and experiment is attributed to the inaccuracy of the DFT methods and implicit solvent model employed and to aspects of the experimental system partially captured by the computational simulations. As such, a recent study estimating the catalytic effect of ExCage6+ on the bowl-inversion of indenocorannulene reported that due to issues of low substrate concentration and low binding affinity, the observed rate constant of enantiomerization was by approximately two orders of magnitude smaller than the real rate constant of the inversion occurring within the host, i.e. the process that can be modelled computationally.15 Due to temperature requirements of the variable temperature NMR studies, we here use DCM, a competitive solvent expected to give lower binding affinities, close to its freezing point when determining the inversion barrier experimentally. It is possible that this may have a similar effect on the observed rate and resulting experimental barrier as reported in ref. 15. Most notably, however, we find the TS in the optimized complex with 1-PP to show the same unusual distortion from planarity that we observed experimentally for the complex between 1-PP and the TS analogue coronene (Fig. 4). We equally find that this geometric change in the guest goes along with a flattening of the PBI moieties of the host (Fig. 6b and Table S2†). Energy decomposition analyses (Fig. 6c and Table S3†) show that compared to the reactant complex, noncovalent interactions between the guest and the perylene planes are indeed maximized in the TS complex, in which host and guest geometries mutually adapt to optimize π–π stacking interactions. Such “mutual induced fit transition structure stabilization” is a well conceivable explanation for the observed catalytic effect of this host for the corannulene bowl-inversion, despite the initial non-complementarity between the non-planar PBI chromophores and the planar guest TS. In order to gain quantitative insights into the energetic consequences of this geometric adaptation upon complexation, we obtain an estimate of the contribution of ground state destabilization to the catalytic effect based on our optimized geometries. Table 2 provides the differences in electronic energies between the geometries of the free optimized corannulene equilibrium and transition structures and PBI cyclophane, and the corresponding geometries of the catalyst and substrate resulting from optimizations of the host–guest complexes. The difference between the sums of the total destabilizations to host and guest in the ground and transition state complexes provides an estimate of the contribution of ground state destabilization to catalysis.
As apparent from Table 2, our calculations indicate that apart from transition structure stabilization, ground state destabilization contributes to the observed catalytic effect. While overall distortions upon complexation only result in a total destabilization of 2.5 kJ mol−1 in the TS complex, we obtain a total destabilization of 10.2 kJ mol−1 for the reactant complex. Thus, we arrive at a contribution of ground state destabilization to catalysis of 7.7 kJ mol−1. It is instructive to note that this contribution is larger than the contribution of ground state destabilization of 2.1 kJ mol−1 reported for the ExBox4+ catalyzed corannulene bowl-inversion.14 This is presumably a consequence of the large, rigid PBI π-surface of this host compared to the smaller, more flexible 1,4-phenylene-extended bipyridinium moieties of ExBox4+. Hereby, the energy difference between the freely optimized corannulene bowl and its geometry in the complex with 1-PP is, at 7.7 kJ mol−1, the main contribution to destabilization. Inspection of the corresponding structures reveals that this is mainly the result of a reduction in bowl-depth upon complexation, with the distance between the centroids of the central 5-membered ring to the centroid of the 10 atoms forming the corannulene rim decreasing from 0.901 Å to 0.706 Å (Fig. S21†). By comparison, our calculations reveal that despite the unusual deviation from planarity of the bowl-inversion TS in the complex, this geometric adaptation carries almost negligible energetic penalty (0.5 kJ mol−1), explaining why catalysis is observed despite the need for mutual induced fit binding in the transition structure complex.
Footnotes |
† Electronic supplementary information (ESI) available. CCDC 2300009. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d3sc05341e |
‡ M. W. and A. A. K. contributed equally. |
This journal is © The Royal Society of Chemistry 2024 |