Valentine P.
Ananikov
* and
Evgeniy G.
Gordeev
Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, Leninsky Pr. 47, Moscow, 119991, Russia. E-mail: val@ioc.ac.ru
First published on 15th September 2011
A new efficient approach was developed for the synthesis of aromatic and heteroaromatic compounds based on [4 + 2] cycloaddition of unsubstituted and heteroatom-substituted alkyne and enyne units. The developed approach provides a practical Green chemical route to several types of important bicyclic products (indane, cyclopentapyridines, indole, isoindole, indolizine, isophosphindole, benzofuran, benzothiophene, benzoselenophene and corresponding dihydro derivatives) starting from simple linear compounds. The mechanism of the reactions was revealed by theoretical calculations using different methods, including CCSD(T) and MP4(SDTQ) for energy calculations and B3LYP, M052X, B3PW91, BLYP and MP2 levels for evaluation of molecular structures.
The classical approach for the synthesis of (hetero)aromatic compounds employs a substitution pathway, which is always accompanied by production of waste (as with any substitution reaction). Cycloaddition methodology is an attractive Green alternative with 100% atom efficiency of the addition reaction, which avoids waste.1 The cycloaddition strategy has an important practical advantage, because it is a single step route to substituted products where the position of the substituent can be easily controlled by the structure of the initial reagent. In fact, a challenging problem of selective preparation of cyclic (as well as bicyclic and even polycyclic) products with all the substituents in place can be solved using a cycloaddition strategy.2,3
[4 + 2] Cycloaddition of enyne and alkyne units is a fascinating direct route to construct six-membered aromatic rings. In a simple model case the intermolecular cycloaddition reaction of vinylacetylene and acetylene leads to the formation of benzene (Scheme 1). Experimental and theoretical studies were carried out to investigate the nature of the reaction and have suggested intermediate involvement of strained cyclic alleneIII (1,2,4-cyclohexatriene).4 In spite of great mechanistic interest, synthetic application of this reaction was limited due to rather harsh experimental conditions required to carry out the intermolecular transformation.
![]() | ||
Scheme 1 Intermolecular [4 + 2] cycloaddition reaction between acetylene and vinylacetylene (but-1-en-3-yne). |
Excellent practical utilization of this cycloaddition strategy was demonstrated by Danheiser and co-workers via introducing yne-yne link and facilitating intramolecular cycloaddition (Scheme 2).5 Further studies in the field resulted in the development of an efficient and flexible modular approach for the synthesis of important classes of aromatic and heteroaromatic compounds.6 The mechanistic studies of both intermolecular (Scheme 1) and intramolecular (Scheme 2) reactions were carried out and have revealed preferable factors for the cycloaddition step involving an intramolecular pathway.7,8 Smaller activation barriers due to lowered activation entropy and favorable intrinsic structural arrangements of reagent I afforded formation of cycloaddition products in good to high yields under simple experimental conditions (Scheme 2).5,6
![]() | ||
Scheme 2 Intramolecular [4 + 2] cycloaddition reaction governed by yne-yne link in non-1-ene-3,8-diyne (I). |
Inspired by Schreiner's recent concept on the description of reaction space of addition reactions9 and in light of our previous studies7,8 we have carried out a systematic theoretical investigation of the interconversions of alkyne and enyne key-units. Surprisingly, we have found six new intramolecular [4 + 2] cycloaddition reactions governed by an ene-yne link (Scheme 3). A detailed study of the reaction mechanism was carried out using density functional theory (B3LYP, BLYP, B3PW91 and M052X), Møller–Plesset methods (MP2, MP4(SDTQ)) and coupled cluster CCSD(T) calculations. Reliable synthetic pathways to aromatic and heteroaromatic products were designed based on the theoretical mechanistic study.
![]() | ||
Scheme 3 Proposed in the present study an intramolecular [4 + 2] cycloaddition reaction governed by ene-yne link in non-3-ene-1,8-diyne and corresponding heteroatom derivatives |
![]() | ||
Scheme 4 The studied mechanism of the intramolecular [4 + 2] cycloaddition reaction of E- (1) and Z- (11) isomers of non-3-ene-1,8-diyne. |
Entry | Level | 1 | 2-TS | 3 | 4-TS | 5 | 6-TS | 7-TS | 8 | 9-TS | 10 | 11 | 12-TS | 13 | 14-TS | Av. Dev.e |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
a For structures see Scheme 4, for optimized molecules see Fig. 2. b Energy evaluation at indicated theory level using 6-311+G(d) basis set, geometry optimization at B3LYP/6-311+G(d) level. c Energy evaluation and geometry optimization at indicated theory level using 6-311+G(d) basis set. d Energy evaluation and geometry optimization at B3LYP/6-311++G(d,p) level. e Average deviation of corresponding ΔE values with those obtained at CCSD(T) level; calculated for points 1–10.31 | ||||||||||||||||
1 | B3LYP (ΔE)b | 0.0 | 31.4 | −16.3 | 33.9 | 4.2 | 17.2 | 8.7 | −11.0 | −7.6 | −100.0 | −1.3 | 35.6 | −13.8 | 2.3 | 2.9 |
2 | CCSD(T) (ΔE)b | 0.0 | 29.9 | −23.9 | 28.0 | −4.7 | 15.2 | 9.3 | −10.6 | −8.1 | −99.9 | −1.5 | 32.6 | −19.9 | 1.8 | 0.0 |
3 | MP2 (ΔE)c | 0.0 | 23.0 | −20.5 | 22.2 | −3.5 | 16.8 | 6.1 | −7.4 | −7.4 | −104.1 | −1.8 | 29.6 | −17.1 | 2.1 | 2.8 |
4 | MP4(SDTQ) (ΔE)b | 0.0 | 25.9 | −22.6 | 24.7 | −4.2 | 18.6 | 7.5 | −9.8 | −8.6 | −100.8 | −1.5 | 30.6 | −18.8 | 0.7 | 1.6 |
5 | BLYP (ΔE)c | 0.0 | 27.6 | −13.0 | 33.7 | 6.2 | 15.4 | 12.3 | −7.4 | −2.6 | −93.0 | −1.2 | 29.6 | −11.4 | 5.9 | 4.9 |
6 | B3PW91 (ΔE)c | 0.0 | 27.5 | −25.8 | 21.0 | −6.1 | 6.8 | −2.0 | −19.6 | −17.7 | −109.7 | −1.3 | 30.9 | −23.2 | −6.6 | 5.7 |
7 | M052X (ΔE)c | 0.0 | 28.1 | −24.1 | 28.0 | −2.8 | 13.0 | 2.6 | −15.4 | −14.0 | −108.3 | −1.9 | 36.0 | −20.3 | −1.7 | 3.1 |
8 | B3LYP (ΔE)d | 0.0 | 32.5 | −14.7 | 34.2 | 5.6 | 18.1 | 9.8 | −9.0 | −6.4 | −98.7 | −1.3 | 37.1 | −12.1 | 4.2 | 3.9 |
1 | 2-TS | 3 | 4-TS | 5 | 6-TS | 7-TS | 8 | 9-TS | 10 | 11 | 12-TS | 13 | 14-TS | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
a Imaginary frequencies for the transition states: 476.6 i cm−1 (2-TS), 1542.0 i cm−1 (4-TS), 625.1 i cm−1 (6-TS), 1001.8 i cm−1 (7-TS), 783.8 i cm−1 (9-TS), 205.3 i cm−1 (12-TS), 2340.0 i cm−1 (14-TS). | ||||||||||||||
C1–C2 | 1.207 | 1.252 | 1.331 | 1.408 | 1.494 | 1.471 | 1.409 | 1.423 | 1.444 | 1.396 | 1.206 | 1.284 | 1.341 | 1.423 |
C2–C3 | 1.421 | 1.386 | 1.329 | 1.359 | 1.364 | 1.435 | 1.422 | 1.497 | 1.459 | 1.397 | 1.422 | 1.387 | 1.327 | 1.390 |
C3–C4 | 1.341 | 1.370 | 1.510 | 1.483 | 1.479 | 1.408 | 1.475 | 1.485 | 1.447 | 1.392 | 1.343 | 1.373 | 1.528 | 1.488 |
C4–C5 | 1.499 | 1.494 | 1.538 | 1.520 | 1.507 | 1.517 | 1.516 | 1.507 | 1.509 | 1.514 | 1.500 | 1.496 | 1.545 | 1.547 |
C5–C6 | 1.540 | 1.541 | 1.541 | 1.546 | 1.551 | 1.547 | 1.547 | 1.551 | 1.551 | 1.548 | 1.540 | 1.546 | 1.549 | 1.543 |
C6–C7 | 1.544 | 1.541 | 1.549 | 1.552 | 1.552 | 1.547 | 1.551 | 1.552 | 1.551 | 1.548 | 1.543 | 1.563 | 1.565 | 1.550 |
C7–C8 | 1.461 | 1.472 | 1.514 | 1.509 | 1.511 | 1.511 | 1.512 | 1.514 | 1.513 | 1.514 | 1.461 | 1.449 | 1.512 | 1.506 |
C8–C9 | 1.203 | 1.242 | 1.343 | 1.401 | 1.486 | 1.424 | 1.385 | 1.426 | 1.419 | 1.392 | 1.203 | 1.264 | 1.352 | 1.361 |
C4–C8 | 3.220 | 2.447 | 1.558 | 1.452 | 1.378 | 1.410 | 1.414 | 1.359 | 1.369 | 1.400 | 3.954 | 2.841 | 1.519 | 1.480 |
C1–C9 | 5.129 | 2.054 | 1.489 | 1.417 | 1.363 | 1.384 | 1.410 | 1.385 | 1.378 | 1.397 | 4.886 | 1.673 | 1.466 | 1.431 |
C1–H | 1.064 | 1.073 | 1.087 | 1.085 | 1.085 | 1.090 | 1.090 | 1.089 | 1.088 | 1.086 | 1.064 | 1.081 | 1.088 | 1.093 |
C3–H | 1.089 | 1.088 | 1.087 | 1.090 | 1.096 | 1.093 | 1.090 | 1.103 | 1.090 | 1.087 | 1.088 | 1.089 | 1.087 | 1.094 |
1.315 | 1.107 | 1.203 | ||||||||||||
C4–H | 1.088 | 1.086 | 1.102 | 1.463 | 2.434 | 2.243 | 1.286 | — | — | — | 1.089 | 1.091 | 1.103 | 1.117 |
C2–H | — | — | — | 1.391 | 1.096 | 1.137 | 2.268 | — | 1.457 | 1.086 | — | — | — | — |
C9–H | 1.064 | 1.073 | 1.087 | 1.086 | 1.088 | 1.088 | 1.089 | 1.091 | 1.089 | 1.087 | 1.064 | 1.085 | 1.087 | 1.088 |
C1–C2–C3 | 178.3 | 138.4 | 130.8 | 123.0 | 112.0 | 118.1 | 113.2 | 113.3 | 111.8 | 120.4 | 177.6 | 134.0 | 129.3 | 111.3 |
C2–C3–C4 | 124.4 | 116.1 | 111.6 | 94.7 | 106.5 | 115.3 | 122.8 | 121.2 | 123.6 | 119.2 | 125.6 | 124.1 | 115.2 | 125.0 |
C9–C8–C7 | 178.2 | 163.0 | 129.7 | 130.7 | 124.9 | 125.6 | 131.1 | 128.7 | 128.9 | 129.2 | 178.9 | 171.9 | 127.5 | 130.8 |
![]() | ||
Fig. 2 B3LYP/6-311+G(d) optimized molecular structures of initial compounds, transition states, intermediates and products of the studied reaction. Displacement vectors corresponding to imaginary frequency are shown for each transition state (see Scheme 4 for structures and Table S2† for optimized geometry parameters). |
Calculations at the B3LYP level will be considered first using a mechanistically reliable ΔG energy surface (Sections 2.1), followed by discussion of the CCSD(T) energy surface, calculations at various theory levels (MP2, MP4(SDTQ), BLYP, B3PW91 and M052X) as well as computations applying a larger basis set (Section 2.2).
Comparison with available experimental and theoretical data on similar reactions will be provided next (Section 2.3). Finally, the reaction leading to different heterocyclic products will be revealed using the studied [4 + 2] cycloaddition strategy (Sections 2.4 and 2.5).
Starting from 3 one of the possible pathways involves a 1,3-hydrogen shift (4-TS) leading to formation of structure 5 with a “Möbius benzene”7,8 arrangement of the six-membered ring. The isomerization of 5 through the transition state 6-TS led to the final product 10. Another pathway includes a 1,2-hydrogen shift (7-TS) with intermediate formation of cyclic carbene8 and affords the same final product 10 after the second 1,2-H shift through the 9-TS.
The calculated energy surface indicated the activation barrier of ΔG≠ = 34.5 kcal mol−1 for the cycloaddition step, and formation of allene3 resulted in an energy gain of ΔG = −9.7 kcal mol−1 (Fig. 1). Comparison of the two possible pathways of further transformation of 3via the 4-TS, 6-TS or via the 7-TS, 9-TS led to selection of the latter as more energetically favored. The activation barrier of 1,2-H shift ΔG≠ = 23.1 kcal mol−1 (3 → 7-TS) was calculated as much lower compared to the barrier of 1,3-H shift ΔG≠ = 47.2 kcal mol−1 (3 → 4-TS). Cyclic carbene8 was more stable by ΔG = −6.6 kcal mol−1 relative to initial compound 1 and less stable by ΔG = 3.1 kcal mol−1 relative to cyclic allene3. The last 1,2-H shift through the 8 → 9-TS step was characterized with a very small activation barrier of ΔG≠ = 3.8 kcal mol−1. Formation of the final product 10 was calculated to be highly exothermic ΔG = −91.4 kcal mol−1, thus providing the necessary driving force for the transformation.
Starting from Z-non-3-ene-1,8-diyne 11 as the initial compound the transition state of the [4 + 2] cycloaddition reaction 12-TS was located. The calculated energy barrier was somewhat higher, ΔG≠ = 40.5 kcal mol−1, and formation of cyclic allene13 was less exothermic, ΔG = −6.0 kcal mol−1, compared to the reaction involving E-geometry of the double bond (ΔG≠ = 34.5 kcal mol−1, ΔG = −9.7 kcal mol−1). Rearrangement of the hydrogen atoms at C3 and C4 centers from syn- to anti- orientation took place via the inversion-like transition state 14-TS, which connected 13 with more stable isomer 3. The calculated 13 → 14-TSactivation barrier was relatively small, ΔG≠ = 15.1 kcal mol−1. The further transformation involving 1,2-H shifts led to the final product 10 as discussed above.
Optimized structures of 2-TS and 12-TS and IRC calculations confirmed concerted nature of the [4 + 2] cycloaddition step: both carbon-carbon bonds C1–C9 and C4–C8 were formed upon movement along the reaction coordinate toward allene intermediates 3 and 13. However, molecular geometries of the studied cycloaddition transition states significantly differ (Table 2 and Fig. 2). A more synchronous transition state 2-TS was located for the cycloaddition involving E-isomer 1, C1–C9 = 2.054 Å and C4–C8 = 2.447 Å, whereas a highly asynchronous transition state 12-TS was found for the cycloaddition of Z-isomer 11, C1–C9 = 1.673 Å and C4–C8 = 2.841 Å. Clearly, the degree of asynchronicity with the Z-geometry of the double bond was much larger compared to the same reaction starting with E-geometry: Δ = 0.393 Å and Δ = 1.168 Å for 2-TS and 12-TS, respectively.10
Formation of a six-membered ring resulted in elongation of the triple bonds which was more pronounced for the Z-isomer, C1–C2 = 1.252 Å (2-TS), 1.284 Å (12-TS) and C8–C9 = 1.242 Å (2-TS), 1.264 Å (12-TS). In contrast, the CC unit was more conservative and the length of the double bond changed only slightly, C3–C4 = 1.370 Å (2-TS), 1.373 Å (12-TS). Thus, a less favorable Z-geometry arrangement of the conservative C
C unit was compensated by more flexible terminal triple bonds. In agreement with calculated geometry trends, the earlier transition state 2-TS was characterized by a smaller activation barrier compared to more late transition state 12-TS (Scheme 4).
Further transformation of cyclic allene3via an energetically favored pathway involved two transition states (7-TS and 9-TS) with geometry corresponding to a 1,2-H shift process (Table 2). In the final product 10 a benzene ring was formed with C–C bonds lengths 1.39–1.40 Å typical for an aromatic system. Selected geometry parameters of all optimized stationary points are shown in Table 2 and will not be discussed in detail here. It is interesting to point out that related allowed 1,5-shifts for the out-of-plane π-system were reported by Alabugin and co-workers.11 Such shifts accompanied with aromatization are known to be relatively fast.11
In terms of average deviation between the calculated energy values, very good agreement was obtained between CCSD(T)//B3LYP and MP4(SDTQ)//B3LYP surfaces (1.6 kcal mol−1). The M052X level performed well (3.1 kcal mol−1), while worse agreement was observed with BLYP (4.9 kcal mol−1) and B3PW91 (5.7 kcal mol−1) functionals (Table 1). Re-optimization of the energy surface with a larger basis set led to only a minor change of the calculated energy surface (Entries 1, 8; Table 1).
Thus, the mechanistic conclusions were found to be valid at all theory levels used in the present study. The studied [4 + 2] cycloaddition reaction takes place with the activation barrier of 26–32 kcal mol−1 and the formation of the allene intermediate is an exothermic process. Further transformation of the allene3via the 1,2-H shifts is a more favorable pathway and formation of the final product 10 is a highly exothermic and irreversible reaction.
From experiment it is known that cycloaddition reactions corresponding to Scheme 2 took place at 180–250 °C and resulted in product formation in good to high yields.5 Several other examples of this cycloaddition strategy were demonstrated under similar conditions and provided a useful practical route to substituted bicyclic derivatives.4,6 The proposed mechanism of the reaction involving the cyclic allene and carbene intermediates was confirmed by trapping experiments and fragmentation reactions.4–6 For convenient practical application microwave-accelerated experimental procedures were reported for heating reaction mixtures of cycloaddition reactions (100–250 °C).3
The reaction explored in the present study (Schemes 3, 4) was characterized by the calculated activation energy for the cycloaddition step and stability of the allene key-intermediate very similar to the experimentally proven transformations (Schemes 1, 2). Thus, it is natural to expect the cycloaddition approach designed in the present study and governed by the ene-yne link to take place as a thermal reaction similar to the process involving an yne-yne link. As a further practical extension of new transformations described in the present study, known reagents and metal catalysts may be added to carry out the reactions under more mild conditions.2,3
![]() | ||
Scheme 5 Cycloaddition reactions leading to heterocyclic products with calculated ΔE values at CCSD(T)/6-311+G(d)//B3LYP/6-311+G(d) level (without parenthesis) and ΔG values at B3LYP/6-311+G(d) level (in parenthesis); in kcal mol−1 (see Supporting Information† Table S4 and Scheme S1 for another representation of the energy surface).13,14 |
Transition states 16-TS, 20-TS, and 24-TS were located for the [4 + 2] cycloaddition of amines with the activation barrier at CCSD(T) level of 29.1–30.3 kcal mol−1 (Scheme 5a–5c). In spite of similar activation energies, stability of cyclic allene intermediates was different: −18.5 kcal mol−1 (17), −25.7 kcal mol−1 (21), and −33.0 kcal mol−1 (25). Indoline 18 and isoindoline 22 were calculated as the final heteroaromatic products, which can be formed after a series of H-shifts.12 Comparison with the 1 → 2-TS → 3cycloaddition step discussed above, the activation barrier of 29.9 kcal mol−1 and stability of cyclic allene of −23.9 kcal mol−1 were calculated in the absence of nitrogen atoms (Entry 2, Table 1).
The [4 + 2] cycloaddition of nitriles26 and 30 was characterized with high activation barriers of 43.4 and 39.0 kcal mol−1, respectively (Scheme 5d–5e). Formation of allene intermediates was found to be endothermic by 14.3 and 14.0 kcal mol−1, respectively. The reaction leading to pyridine products was calculated to be exothermic by −58.9 and −59.5 kcal mol−1 with a smaller energy gain when compared to formation of indoline (18), isoindoline (22) and indane (10) products −98.5–−105.5 kcal mol−1 (Scheme 5, 4 and Table 1).
For the synthesis of the pyridine system 33 the reaction involving imine34 should be a more useful alternative than transformation of the nitrile30. A smaller activation energy of 32.0 kcal mol−1 (35-TS) and better stability of the intermediate allene −23.3 kcal mol−1 (36) were found to be more favorable for this transformation. Larger energy gains in the formation of product 33 (cf.Scheme 5e and 5f) reflects a contribution from relative stability of the reagents 30 and 34 (see Figure S1†).
Cycloaddition of imine37 was found to be unique among the studied reactions (Scheme 5g). The most stable cyclic allene39 (−40.2 kcal mol−1) and the most stable N-heterocylcic carbene intermediate 40 (−56.1 kcal mol−1) were found in the calculations. Possible rearrangement of the carbene may involve H-shifts from five-membered rings leading to dihydroindolizine derivates 41 and 42.
Comparing calculated energy surfaces at CCSD(T) and B3LYP levels we may notice the tendency outlined already for the reaction leading to product 10 (Table 1). Calculated activation energies for the cycloaddition step were in a good agreement, while cyclic allene intermediates were found to be more stable at the CCSD(T) level. Both theory levels confirm the exothermic nature of the studied cycloaddition reactions on the final product formation stage. Optimized geometries of all stationary points for the transformations involving heteroatom-substituted molecules as well as reactions of cis-34 and cis-37 are provided in the Supporting Information (Scheme S2, Figures S4–S6, Table S5†).
It is interesting to point out, that compounds 36 and 39 significantly differ in relative stability: ΔE = E(39) − E(36) = −24.2, −16.9 and −23.0 kcal mol−1 at the B3LYP/6-311+G(d), CCSD(T)/6-311+G(d)//B3LYP/6-311+G(d) and MP2(full)/6-311++G(d,p) levels of theory, respectively. Thus, compound 39 with a tertiary nitrogen atom is much more thermodynamically preferred, than molecule 36 with two carbon substituents at the nitrogen atom. Compound 39 has more formal “aromatic” attributes compared to 36: a six-membered heterocycle is planar (36 is not planar) and interatomic distances of six-membered heterocycles are close to one another (in 36 there is a marked difference). We have applied Bader quantum chemical theory of atoms in molecules15 to reveal the difference between structures 36 and 39 and have found that molecule 39 possess more similar electronic features to the aromatic pyridine molecule, rather than molecule 36, in total agreement with observed thermodynamic stability.16 This electronic structure also suggests that the reaction is somewhat mechanistically different from the other cases discussed in the present study. It may be also considered as a pseudopericyclic process, which starts with nucleophilic attack of the nitrogen at the yne terminus. A similar consideration has been reported for radicals “masquerading as electrophiles” and this mechanism would make the process a new and quite interesting type of cycloaromatization reaction.17,18
It should be pointed out that indoline (18), isoindoline (22), and dihydroindolizine (41, 42) can be easily converted to indole, isoindole and indolizine derivatives using well-established experimental methods.19 Thus, providing an easy route to bicyclic heteroaromatic products.
The calculated energy surfaces for the reactions of phosphorus (43), oxygen (47), sulfur (51) and selenium (55) analogs indicated activation energy 28–31 kcal mol−1 and relative stability of the allene intermediate of about −25 kcal mol−1 at CCSD(T) level (Scheme 6). All transformations were characterized with a large energy gain of about −103 kcal mol−1. Higher energy barriers and slightly lower stability of cyclic allenes and products were found on the ΔG energy surface calculated at the B3LYP level. The energy parameters well agree with those calculated for the nitrogen derivatives (Scheme 5b).
![]() | ||
Scheme 6 Cycloaddition reactions leading to heterocyclic products calculated ΔE values at CCSD(T)/6-311+G(d)//B3LYP/6-311+G(d) level (without parenthesis) and ΔG values at B3LYP/6-311+G(d) level (in parenthesis); in kcal mol−1 (see Supporting Information† Table S6 for energy surface, and Figure S7 and Table S7 for optimized structures).20,14 |
Indeed, the calculations have shown that the studied cycloaddition reactions should be suitable for preparation of a series of heterocyclic compounds, such as 2,3-dihydro or 1,3-dihydro derivatives of isophosphindole (46), 2-benzofuran (50), 2-benzothiophene (54), and 2-benzoselenophene (58). Removal of two hydrogen atoms and formation of heteroaromatic rings can be carried out using known experimental procedures.19
The mechanism of the reaction of non-3-ene-1,8-diyne (1) (Scheme 3) was proposed to involve cycloaddition and a series of H-shifts with the final formation of the indane core (10). Both E- and Z- derivatives can be involved in the reaction and furnish the same final product. According to the calculations, the E-isomer was predicted to be more reactive, although the difference in reactivity was small and the E-/Z- mixture may be utilized in the reaction as a starting material. The reliability of the calculations was confirmed at various theory levels.
Using nitrogen-substituted reagents all possible cycloaddition reactions were investigated to develop a new strategy for the synthesis of heteroaromatic products (Scheme 4). The calculations revealed the following order of relative reactivity in the studied cycloaddition reaction: amines > imines > nitriles. Thus, the reaction of amines with the final formation of indoline (18) and isoindoline (22) cores should be possible in more mild conditions at lower temperature compared to the synthesis of pyridine systems (29, 33) from imines and even less reactive nitriles.
Cycloaddition reactions involving imines are of much interest and can result in the formation of pyridine (33) or dihydroindolizine (41, 42) cores depending on the position of nitrogen atom in the starting reagent. It is worthwhile mentioning that synthesis of compounds 18, 22, 41 and 42 may be followed by hydrogen abstraction with well-known experimental procedures to prepare indole, isoindole and indolizine products.
Important to note, according to the theoretical predictions formation of stable cyclic allene (39) and carbene (40) intermediates can be expected in the reaction of imine37. In a certain case such intermediate species can be directly observed by experimental methods or analyzed using known trapping techniques.
A series of new reactions were developed in the present study based on [4 + 2] cycloaddition of alkyne and enyne key units as well as various substituted units. The designed reactions utilize Green chemistry potential of cycloaddition approach and provide a convenient synthetic route to important classes of aromatic and heteroaromatic products. We encourage implementing the developed reactions and investigating their scope under experimental conditions.
Geometry optimization of the initial compounds, transition states and products of the reactions was carried out using the B3LYP hybrid density functional method.22 For the studied structures normal coordinate analysis was performed to characterize the nature of the stationary points and to calculate thermodynamic properties (298.15 K and 1 atm). Transition states were confirmed with IRC (Intrinsic Reaction Coordinate) calculations using the standard method.23
Geometry optimization and energy calculations were also carried out with MP2 (Entry 3, Table 1),24BLYP (Entry 5, Table 1),22B3PW91 (Entry 6, Table 1),22,25 and M052X (Entry 7, Table 1)26 methods using 6-311+G(d) basis set.
Single point energy calculations were carried out at CCSD(T)/6-311+G(d) level27 were performed for all studied reactions (Entry 2, Table 1; Sections 2.3–2.5) and single point calculations at MP4(SDTQ)/6-311+G(d) level28 were executed for the cycloaddition of non-3-ene-1,8-diyne (Entry 4, Table 1). In both cases B3LYP/6-311+G(d) optimized geometries were used for the single point calculations.
In the previous studies it was established B3LYP, MP4(SDTQ) and CCSD(T) levels of theory reasonably well describe the energy and geometry parameters of the cycloaddition reactions of interest,8 the utility of M052X for studying cycloaddition reactions was reported recently.29
All calculations were performed without any symmetry constraints using the Gaussian09 program.30
Footnote |
† Electronic supplementary information (ESI) available: Tables S1–S7, Figures S1–S7, and Schemes S1–S2. See DOI: 10.1039/c1sc00380a |
This journal is © The Royal Society of Chemistry 2011 |