Dorian S.
Sneddon
,
Paul V.
Kevorkian
and
Thomas R.
Hoye
*
Department of Chemistry, University of Minnesota, 207 Pleasant St. SE, Minneapolis, MN 55455, USA. E-mail: hoye@umn.edu
First published on 7th January 2025
The hexadehydro-Diels–Alder (HDDA) reaction is a cycloisomerization between a conjugated diyne and a tethered diynophile that generates ortho-benzyne derivatives. Considerable fundamental understanding of aryne reactivity has resulted from this body of research. The multi-yne cycloisomerization substrate is typically pre-formed and the (rate-limiting) closure of this diyne/diynophile pair to produce the isomeric benzyne generally requires thermal input, often requiring reaction temperatures of >100 °C and times of 16–48 h to achieve near-full conversion. We report here that diynoic acids can be dimerized and that the resulting substrate, having a 3-atom anhydride linker (i.e., OCOC
O), then undergoes HDDA cyclization within minutes at or below room temperature. This allows for the novel in situ assembly and cyclization of HDDA benzyne precursors in an operationally simple protocol. Experimental kinetic data along with DFT computations are used to identify the source of this surprisingly huge rate acceleration afforded by the anhydride linker: >107 faster than the analogous multi-yne having, instead, a CH2OCH2 ether linker.
It was nearly 30 years later that Diels and Alder published their seminal work producing adducts such as 5, in what we regard today as the classical Diels–Alder reaction (Fig. 1b).7 A more recent variant of a net cycloaddition reaction to produce six-membered carbocycles is the hexadehydro-Diels–Alder (HDDA) reaction (Fig. 1c).8 For example, a poly-yne substrate joined by a three-atom tether such as the ether linker in 6 undergoes net intramolecular [4 + 2] cycloaddition to form highly reactive, yet selective, ortho-benzyne derivatives such as 7 under purely thermal conditions.9 Subsequent experiments done by the Ward group10 in which phenylpropiolic acid (1) was converted to the anhydride 2 using DCC were shown to produce the known TDDA adduct 4 even at 0 °C.11 Considering that most TDDA and HDDA cyclizations proceed by stepwise mechanisms,12 we reasoned that an anhydride linked multi-yne substrate could potentially produce an HDDA benzyne also at ambient temperature. To date, there are few examples of HDDA substrates that cyclize at or below room temperature.8,13 One final reaction that influenced our thinking and design of the project we describe here is shown in Fig. 1d. Nicolaou and coworkers described the formation of symmetrical anhydrides under very mild conditions via transient, mixed carboxylic–sulfonic anhydrides.14
By incubating 10a with MsCl and pyridine in a non-interacting solvent (DCM instead of the THF used in Fig. 2), we were able to establish reaction conditions that allowed for generalized benzyne trapping. The first example employing these conditions used furan as an external trap (Fig. 3a). This led to formation of the oxanorbornadiene derivative 14a in excellent yield. We observed that trifluoroacetic anhydride and methanesulfonic anhydride also were effective dehydrating agents, although we did not further explore this. As shown in Fig. 3b, this anhydride-HDDA process tolerates the presence of electron-rich and -poor aryl (10b–d), heteroaryl (10e), alkyl (10f), and silyl (10g–h) substituents on the terminus of the diyne carboxylic acid substrate.
We also used several other trapping agents (Fig. 3c): benzyl azide provided the triazole 15; thioanisole gave the sulfides 17; ethyl diazoacetate led to the pyrazole 19; mesityl nitrile oxide gave the isoxazole 20; a munchnone derivative gave rise to the isoindoles 21; and the simple dienes N-Boc-pyrrole, cyclopentadiene, and anthracene afforded the [4 + 2] cycloadducts 22, 23, and 24 respectively. We note that the tetrayne-derived products that contain an arylethynyl substituent on the phthalic anhydride scaffold show blue fluorescence, which represents further opportunity for investigation.
Not surprisingly, the phthalic anhydride moiety in the products of these one-pot, in situ, anhydride assembly, cycloisomerization, and trapping cascades could be efficiently derivatized with primary amines to afford the phthalimide derivatives 25a–f (Fig. 4).16 Phthalic anhydrides are substrates for a considerable array of additional types of transformation.16 A recent report of anti-bacterial activities of similar phthalimide derivatives against plant pathogens is an example of contemporary interest in this class of compound.17
Recognizing that the applicability of this strategy would be enhanced if we could selectively make and cyclize mixed (i.e., unsymmetric) anhydrides, we explored conditions that might allow this. Using the MsCl-mediated conditions for anhydride formation from a mixture of two different diynoic acids, one would expect, of course, a relatively non-selective formation of products from the three possible homo- and hetero-anhydrides. This was confirmed in a preliminary experiment (1H NMR analysis of the crude product mixture18). Thus, we needed to find more-controlled, milder conditions for selectively forming a mixed anhydride that would persist sufficiently long for its HDDA cycloisomerization to take place. That is, conditions were needed that would not induce anhydride exchange of the unsymmetric anhydride, resulting in formation of unwanted symmetric anhydrides. We explored the use of an acid chloride/acid pair of substrates to generate the mixed anhydride. First, a dichloromethane solution of the preformed acid chloride 26 (structure in Fig. 5a) and the diyne acid 10f was treated with a slight excess of pyridine. Again18 we observed a more or less statistical mixture of the crossed and symmetrical products 27f, 28f, 14b, and 14f. This implied that pyridine and/or pyridinium chloride was catalyzing an unwanted anhydride metathesis event. We also examined the use of a silver(I) carboxylate salt to couple with the acid chloride 26.19 However, anhydride formation was thwarted by decarboxylation to give a terminal alkyne. We deemed other nucleophilic amine bases (e.g., DBU, DMAP, or Et3N) to be poor choices because of the susceptibility of ynoyl chlorides or anhydrides to undergo either 1,2- or 1,4-addition to the ynoyl moiety by species containing a basic nitrogen atom.
Finally, we chose to explore the use of Proton-sponge® [1,8-bis(dimethylamino)naphthalene], a less nucleophilic, yet more basic amine than pyridine.20 This choice proved very rewarding (Fig. 5). Now the products derived from the unsymmetrical mixed anhydrides, arising from reaction of the acid chloride 26 with the acid 10f or 10h, were produced in good yield (Fig. 5b). Likewise, the triynes from 26 + 10i–k were efficiently formed as shown by product formation. These results indicate that scrambling of the initial mixed anhydride was slower than its HDDA cyclization using Proton-sponge®. In the case of the mixed tetraynes, two different benzynes were formed, accounting for the formation of product pairs 27f/28f and 27h/28h.
The triynes derived from the propiolic acid derivatives 10i–k also cyclized, now each to the furan adduct 27i–k from trapping of the only possible benzyne intermediate. Notice that formation of the mesityl-containing product 27i arises from an HDDA process even though a TDDA adduct could have been formed from the diradical intermediate common to both pathways. The presence of non-hydrogen substituents on C2 and C6 of the mesityl group significantly slows the TDDA mode of cyclization.21
A particularly striking result from the series of reactions involving of triyne anhydrides was the formation of 27k from propiolic acid. To the best of our knowledge, this is the first example of a triyne containing a terminal alkyne to undergo an HDDA cyclization at room temperature. The additional remote alkyne in a tetrayne enhances the rate of a HDDA cycloisomerization compared with an analogous triyne substrate.22–24 This is because the additional radical stabilizing energy (RSE) of the alkyne lowers the barrier in the rate-limiting, benzyne-forming event.25 Therefore, we wondered if the ability to make mixed anhydrides would provide an opportunity to now observe the formation of the triyne anhydride by in situ1H NMR analysis. Indeed, reaction of 10k and 26 in CDCl3 containing furan and Proton-sponge® allowed us to observe an acyclic anhydride (see 29a, Fig. 6a) for the first time. This species appeared and competitively disappeared as a transient intermediate enroute to the furan-trapped product 27k.
The ability to observe the mixed anhydride in the relatively less reactive triyne anhydride 29a (from 10k and 26, Fig. 6a) meant that it should be possible to assess the magnitude of the rate enhancement afforded by the anhydride linkage. A stack plot of spectra over time shows changes in the methoxy resonances of the PMP groups in the acid chloride, mixed anhydride, and final product 26, 29a, and 27k, respectively (Fig. 6b). Under the indicated conditions, the rate of formation of the anhydride, a bimolecular process, was observed to compete with that of its cyclization to the benzyne, which is not observed, of course, because of its extremely fast reaction with furan. The changes in the relative amounts of the three PMP-containing species over time were simulated with COPASI,26 which provided the data represented in Fig. 6c. From this the extracted half-life for the HDDA cyclization of 29a to the benzyne was 9 minutes.
We next asked how the rate of the anhydride linked substrate 29a compared with that of the analogous ester- and ether-containing substrates 29b/c and 29d (Fig. 7). Because these latter three all cycloisomerize much more slowly than the anhydride analog, each was easily synthesized and purified. Each was then first converted preparatively into the HDDA benzyne-trapped products 30b/c (from reaction with furan) and 30d (from reaction with cyclooctane27).
![]() | ||
Fig. 7 Large rate enhancement (red) for the HDDA cyclization of the triyne 29a having an anhydride tether compared to the rates of substates 29b–d having ester or ether tethers. |
With NMR data for authentic samples of products 30b–d in hand, we proceeded to measure the half-life for the rate-limiting cycloisomerization of each of the triynes 29b–d in CDCl3 (see ESI† for the detailed protocol). Using a convenient temperature for the reaction rate of these much less reactive substrates (relative to the anhydride 29a), we used NMR spectroscopy to monitor the disappearance of each. The half -lives at 120 °C for 29b and 29c were very similar to one another (10.5 and 9.5 h, respectively) as well as to that seen for an analog of 29b in an earlier study.13,23 Substrate 29d having the ether linkage proved to be even less reactive. (All of the measured half-lives for 29a–d are also listed in column two of the tabulated data in Fig. 8b).
Clearly, the anhydride linker present in substrate 29a and, by extension, all of the other anhydrides involved in the reactions in Fig. 2, 3, and 5 are imparting a dramatic rate-acceleration on the HDDA cycloisomerization. Why? To address this question, we carried out a revealing distortion/interaction (aka, activation strain) analysis (Fig. 8a).28 Generally speaking and in the well-stated words of Bickelhaupt and Houk, “the activation strain or distortion/interaction model is a tool to analyze activation barriers that determine reaction rates. For bimolecular reactions, the activation energies are the sum of the energies to distort the reactants into geometries they have in transition states plus the interaction energies between the two distorted molecules. The energy required to distort the molecules is called the activation strain or distortion energy. This energy is the principal contributor to the activation barrier. The transition state occurs when this activation strain is overcome by the stabilizing interaction energy.”28
The specific approach used here involved first computing the energetic barrier for the initial C–C bond formation in each of 29a–d leading to the corresponding diradical (see ESI†) enroute to each benzyne. These computed activation energies (ΔE‡) are given in column 3 of Fig. 8b. They correlate well with the approximate relative rates of reaction of each of these four substrates (extrapolated from their reactions at different temperatures29). The tethering atoms in each computed structure of 31a–d‡ were then excised and replaced by a hydrogen atom to give rise to the four different structures of 34‡.30 In each of these, the internuclear distance l in 31‡ is maintained. The two “halves” of 31‡ were then separated to infinite distance, giving 32dist and 33dist, the distorted but non-interacting, free alkyne components. To identify the distortion energy developed in the transition structures (TSs) 31‡, the tether was also excised from 32 to produce, following optimization, the undistorted alkyne components 32 and 33. The difference in (the single point) energies between 32distvs.32 and 33distvs.33 comprise the total distortion energy (ΔE‡dist) Likewise, the difference in energy between 34‡ and 32dist + 33dist gives the interaction energy for each of the four analogs.
Notably, the values of both the ΔE‡dist and the ΔE‡int [as well as, therefore, their sum ΔE‡total] were nearly identical for all four sets of structures (columns 4–6, Fig. 8b). Moreover, the interaction energy in 34‡ for each of the four reactions was destabilizing. This is an uncommon, but not unprecedented,31 phenomenon, including an example seen for a computed, stepwise HDDA reaction.32 The implication is that at the distances and geometries between the alkyne moieties in 34‡ (and 31‡), the alkynes have more destabilizing, repulsive interaction than stabilizing, bonding interaction. The energetic contribution imposed/provided by the tether ΔE‡teth is, therefore, the difference in activation energies between the intact tethered substrates 29a–d proceeding to 31a–d‡ (ΔE‡) vs. the untethered analogs 32 and 33 proceeding to 34a–d‡ ΔE‡total.
These are significantly different across the series of tethers (column 7, Fig. 8b). Relative to the ether in 29d, the anhydride in 29a provides ca. 10 kcal mol−1 greater stabilization to its transition structure.
At a more fundamental level, the question remains: why is the anhydride so different? Two factors may be contributing. First, nearly all acyclic anhydrides, including the HDDA substrates here, are composed (as judged by computation as well as microwave and rotational–vibrational spectroscopies) of a family of non-coplanar conformations33 (the exception being formic anhydride34). In the phthalic anhydride products produced here, the enforced planarity of the anhydride could lead to a greater degree of stabilization (resonance and/or strain?)35 relative to the ester-to-lactone and ether-to-ether cyclizations. To probe that possibility, we examined the simple set of homodesmotic reactions (note the same number and type of functional groups on each side of the reaction equation) shown in Fig. 9. The enthalpy change between the (all-anti conformer of the) ring-opened vs. the ring-closed forms for each of the three functional groups (i.e., ether, vs. ester, vs. anhydride) was computed with DFT. The magnitudes of the differences in the ΔH's across this series of cyclizations correlate well with the differences in the rates of the HDDA cyclizations. This reinforces the conclusions from the distortion/interaction analysis. That is, that the primary differentiating factor in the large differences in reactivity for the series 29a–d is the enthalpy change associated with each tether as it undergoes the cyclization, the anhydride being uniquely favorable.
Finally, a 2023 report from the Chen and Li groups at Guizhou University describes a robust method for preparing symmetrical N,N-bis-(3-aryl)propiolylimides that cyclize at ambient temperature to give naphthalimide derivatives 38 (Fig. 10a).17 Precursor amines 35 were acylated with two equivalents each of an arylpropiolic acid 36 and N,N′-diisopropylcarbodiimide (DIC). The intermediate imides 37 spontaneously underwent a TDDA reaction (cf.Fig. 1a) to produce products 38 in very good to excellent yields. Over four dozen examples were reported. All reactions were performed at ambient temperature for 24 h.
![]() | ||
Fig. 10 Tetradehydro-Diels–Alder (TDDA) reactions of propiolic acid chloride-derived anhydride- and imide-linked TDDA substrates 2 and 37. The imide cyclizes faster than the anhydride. |
A reviewer of our manuscript asked whether there was any other type of linker that would accelerate a cyclization by virtue of the planarization phenomenon seen with the anhydride. Because the imide-templated reaction from Chen and Li represents such a linker and especially because it bears considerable similarity to the anhydride work disclosed here, we were curious about the actual extent to which an imide would accelerate this type of cycloisomerization. Toward that end, we prepared, in situ, the symmetrical TDDA substrates, phenylpropiolic anhydride (2) and N-(p-tolyl)-N,N-bis-(3-phenylpropiolyl)imide (41) (Fig. 10b). We used 1H NMR spectroscopy to monitor reactions in CDCl3 between phenylpropiolic acid chloride (39) with either 3-phenylpropiolic acid (1) in the presence of Proton-sponge® or 3-phenyl-N-(p-tolyl)propiolamide (40) in the presence DMAP/Et3N.36 The rate of the TDDA cycloisomerization of the imide 41 to produce the naphthalimide derivative 42 was significantly faster than that of the anhydride 2 to give 4 (see ESI†). Further studies to more fully evaluate and quantify the extent of this even greater rate acceleration by an imide are planned. We note in closing that this strategy for imide formation is readily amenable to the preparation of unsymmetrical N,N-bis-(propiolyl)imides.
We also deduced the reason for the rate enhancement afforded by the anhydride linker. By synthesizing a series of triyne substrates that varied only in their linker structure (ether, ester, and anhydride), we found that the rate of HDDA cycloisomerization increased, dramatically, in the stated order (estimated to be >107 compared to the analogous ether-linked substrate29). To rationalize this trend, we performed a distortion/interaction analysis that revealed that: (i) the distortion/interaction within the diyne–diynophile reacting atoms was nearly identical for all three classes of substrate, and (ii) all four substrates have a destabilizing (i.e., positive) interaction energy at the internuclear distance present in the computed transition structure of each.
A series of simple homodesmotic calculations were performed to assess the enthalpy difference between ring-opened vs. -closed forms for THF, γ-butyrolactone, and succinic anhydride. This showed that it is the energetic contribution associated with the formation of the fused, cyclic five-membered anhydride in the benzyne intermediate that is the principal contributor to the substantially higher rate of the HDDA cycloisomerization. Thus, the differences in stabilization of the transition structures across the series of four substrate classes arise nearly entirely from changes within the reorganization of the tethering atoms. This understanding explains why that, in this setting, the anhydride tether is a kinetically privileged linker. Finally, we have observed that an N-aryl imide analog of the anhydride linker induces an even greater rate acceleration.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 2353613. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4sc07232d |
This journal is © The Royal Society of Chemistry 2025 |