Mechanistic Exploration of the Competition Relationship between

2 days ago - The introduction of C=O, C=C, C=S, or C=N bond has emerged as effective strategy for carbocycle synthesis. A computational mechanistic st...
0 downloads 9 Views 1MB Size
Subscriber access provided by Queen Mary, University of London

Mechanistic Exploration of the Competition Relationship between Ketone and C=C, C=N or C=S Bond in the Rh(III)-Catalyzed Carbocyclization Reactions Yang-Yang Xing, Jian-Biao Liu, Chuan-Zhi Sun, Fang Huang, and De-Zhan Chen J. Org. Chem., Just Accepted Manuscript • DOI: 10.1021/acs.joc.8b00292 • Publication Date (Web): 23 Mar 2018 Downloaded from http://pubs.acs.org on March 24, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Organic Chemistry

Mechanistic Exploration of the Competition Relationship between Ketone and C=C, C=N or C=S Bond in the Rh(III)-Catalyzed Carbocyclization Reactions Yang-Yang Xing, Jian-Biao Liu,* Chuan-Zhi Sun, Fang Huang and De-Zhan Chen* College of Chemistry, Chemical Engineering and Materials Science, Collaborative Innovation Center of Functionalized Probes for Chemical Imaging in Universities of Shandong, Shandong Normal University, Jinan 250014, P. R. China Abstract The introduction of C=O, C=C, C=S, or C=N bond has emerged as effective strategy for carbocycle synthesis. A computational mechanistic study of Rh(III)-catalyzed coupling of alkynes with enaminones, sulfoxonium ylides or α-carbonyl-nitrones was carried out. Our results uncover the roles of dual directing groups in the three substrates and confirm that the ketone acts as the role of the directing group while the C=C, C=N, or C=S bond serves as the cyclization site. By comparing the coordination of ketone versus C=C, C=N, or C=S bond, as well as the chemoselectivity concerning the six- versus five-membered formation, a competition relationship is revealed within the dual directing groups. Furthermore, after the alkyne insertion, instead of the originally proposed direct reductive elimination mechanism, the ketone enolization is found to be essential prior to the reductive elimination. The following C(sp2)−C(sp2) reductive elimination is more favorable than C(sp3)−C(sp2) formation, which can be explained by the aromaticity difference in the corresponding transition states. The substituent effect on controlling the selectivity was also discussed. 1.Introduction Transition-metal-catalyzed C−H functionalization has served as a powerful tool for the construction of cyclic scaffolds, which enhances the efficiency and scope for the synthesis of diverse heterocyclic and carbocyclic compounds.1 To date, owing to the widely-used N-, S- and P-containing directing groups (DGs), the cyclic products are mainly heterocycles, in which the moiety or whole of the DGs is retained in the target frameworks.2 However, due to the competing heteroatom cyclization reaction, synthesis of carbocyclic compounds is more challenging, which is noteworthy in view of the fact that carbocycles are carbon skeleton motifs present in various kinds of organics, pharmaceutical molecules, and natural compounds. Therefore, the development of efficient method to synthesize carbocycles has attracted considerable attention in recent years. Utilization of a weaker coordinating moiety is usually required to promote the carbocycles formation, and the most widely-applied synthetic strategy is the introduction of ketone group that exhibits high electrophilic reactivity.3 In the pioneering work by Glorius, a method for the comprehensive and selective preparation of diversely functionalized indenol and fulvene derivatives using ketone as DG was achieved.3a In another key example, In 2011, Cheng et al. developed Rh(III)-catalyzed carbocyclization of aryl ketones and alkynes (eq 1, Scheme 1).3b Similar reactions involving Rh(III) or Ru(II)-catalyzed C−H functionalization by annulation with alkenes or alkynes are also reported.3 Since ketone group acts as both directing and cyclization 1

ACS Paragon Plus Environment

The Journal of Organic Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

sites in these reactions, only the five-membered carbocyclic ring-containing targets are obtained. Cheng et al. further reported Pd(OAc)2-catalyzed cyclization of allylarenes with alkynes, in which the weaker alkene moiety is used as DG (eq 2).4 However, the use of alkene group brings about several known drawbacks such as requiring strong acidic condition. Recently, Zhu et al. proposed an enaminone-directed C–H functionalization method, and Rh(III)-catalyzed synthesis of naphthalenes via coupling of enaminones with alkynes was realized (eq 3).5 The novel reaction offers an important method for facile synthesis of six-membered carbocycles. The obtained naphthalene derivatives contain both a nucleophilic aldehyde unit and an electrophilic hydroxyl unit in the same molecule, which makes them interesting synthons in further transformations. Similar reactions for the Cp*Rh(III)-catalyzed coupling of sulfoxonium ylides (eq 4) or α-carbonyl-nitrones (eq 5) with alkynes were reported by Li et al..6 The carbocyclization strategies in eqs 3-5 have similar catalytic cycles. Directing and cyclization sites are separated in the functional groups, and ketone acts as the role of the directing site while the C=C, C=N, or C=S bond serves as the cyclization site. Because of the different electronic structures of C=C, C=N and C=S bonds, various naphthalene derivatives were obtained in these reactions. The unit of the aldehyde is incorporated into the target naphthalenes (eq 3), while the dimethylsulfoxide or nitrone group dissociates completely from the final products (eqs 4 and 5). Although plausible catalytic cycles were proposed respectively for the novel transformations, however, the detailed mechanisms especially the cooperative pathway of the directing and cyclization sites have not been explored at the outset of this work. In addition, the actual sequence of the dimethylsulfoxide or nitrone group dissociation step and migratory insertion remains unclear. Herein, we performed density functional theory (DFT) calculations on the three representative rhodium(III)-catalyzed carbocyclization reactions to understand the exact cooperative mechanism of the dual directing groups and the underlying factors resulting in the product difference in the three similar catalytic cycles.

2

ACS Paragon Plus Environment

Page 2 of 18

Page 3 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Organic Chemistry

Scheme 1. Representative carbocyclization reactions employing various DGs.

2. Computational Methods Geometry optimizations and frequency calculations were performed at B3LYP/BS1 level7 in dichloroethane (DCE) using the SMD solvation model.8 BS1 designates a mixed basis set of ECP28MDF_AVDZ for Rh9 and 6-31G(d,p) for other atoms.10 The results of frequency calculations were further examined to confirm each structure is a local minimum (no imaginary frequency) or a transition state (only one imaginary frequency). IRC calculations were performed to verify the right connections among a transition state and its forward and reverse minima. Thermodynamic corrections at 353.15 K and 1 atm for all structures in DCE were obtained by harmonic frequency calculations at the B3LYP/BS1 level. To obtain solvation-corrected free energies, we performed single-point energy calculations at 11 M06 /ECP28MDF_AVDZ/6-311++G(d,p) level in DCE with SMD on the B3LYP/BS1-optimized geometries. All calculations were performed using the Gaussian 09 program package.12 Natural bond order (NBO) calculations were performed by GenNBO 5.0 program13 using the wavefunction got from B3LYP/BS1 level in DCE with the SMD model on the selected systems. The 3D structures were prepared using CYLView 14 and VMD.15

3. Results and Discussion

3

ACS Paragon Plus Environment

The Journal of Organic Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 18

OH O O

OH

HOAc HNMe2

1a or 1b + 2a or 1c X=OAc or OPiv

O HX

R

R

[Cp*RhX]+

+

+

Ph

Ph

Ph

alkyne insertion

Ph

Ph

+

β−Η

elimination

HOAc HNMe2

P1 NMe2 R= C

Ph Ph hydrolysis deamination O R

+

[Rh]

H-[Rh]

Ph Ph

R

O

[Rh]

C− H activation

HOAc H 2O O

+ [Rh]

+

[Rh]

reductive elimination

[Rh] reductive Ph elimination

+

hydrolysis deamination +

H-[Rh]

[Rh]

β−Η

Ph

Ph

elimination

Ph

O

HOAc H 2O NMe2

Ph

isomerization O R= S

DMSO or tBuNO

;

O N

DMSO or tBuNO tBu

O O

1a

NMe2 C

O

O S

1b

O

O N

P2

+ tBu Ph

1c

[Rh] reductive elimination Ph

protonolysis O [Rh] +

[Cp*RhX]+ HX

Ph

[Rh]=Cp*Rh Ph

Scheme 2. Proposed catalytic cycle for the Rh(III)-catalyzed coupling of enaminones, sulfoxonium ylides or α-carbonyl-nitrones with alkynes. As shown in Scheme 2, the three reactions have similar mechanisms, which are all initiated by the C−H activation to form the five-membered rhodacyclic intermediates. Coordination of the substrate 2a to the metal center results in the π-complexes. The insertion intermediates undergo the O-bound to C-bound tautomerization to generate the seven-membered rings, after which the reactions bifurcate onto two pathways. In one route, the release of dimethylamino, dimethylsulfoxide or nitrone group occurs first, then the C−C reductive elimination happens. Alternatively, the reductive elimination might occur first. In the following sections, we will discuss the processes in detail. 3.1 C− −H activation and alkyne insertion

4

ACS Paragon Plus Environment

Page 5 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Organic Chemistry

Figure 1. Highest occupied molecular orbital (HOMO) of substrates 1a and 1b, and HOMO-1 of substrate 1c. Numbers in red are NPA charges. In addition to the dative bond between the metal center with the oxygen atom of the ketone moiety, the substrate 1a can also coordinate with rhodium through the π electrons of the C=C moiety (Figure 1). The frontier MOs of substrates 1b and 1c possess similar characters. As revealed by the MO shapes, for the initial C−H activation and alkyne insertion steps in the reaction, either the ketone or C=C group of substrate 1a can play the role of coordinating. The NPA charges of O and C atoms in 1a are −0.61 and −0.46, respectively, suggesting that O atom is more nucleophilic. Similar result was observed in 1c while the C atom in 1b becomes more electronegative compared to O atom. For clarity, we only listed the results of C−H activation and alkyne insertion steps of substrate 1a in this subsection, and the results of substrate 1b or 1c are given in Figures S2-S3. It’s noteworthy that 1a could rotate around the C−C bond to generate a thermodynamically less stable rotamer 1a′′ via a barrier of 12.6 kcal/mol (Figure S4). Herein, we chose the more stable substrate 1a as the starting point to compute the overall catalytic cycle.

Figure 2. Free energy profiles for C−H activation and alkyne insertion of substrate 1a. Energies are relative to [Cp*RhOAc]+ + 1a and are mass balanced (similarly hereinafter). Active species [Cp*RhOAc]+ interacts firstly with 1a to form the O-coordinated enolate 2 or its C-coordinated isomer 2′′.16 The results in Figure 2 show two competitive pathways, that is, 2 to 3 to 4 to 5 versus 2′′ to 3′′ to 4′′ to 5′′. In the black path, the C−H activation of the enolate 2 proceeds via the concerted metalation-deprotonation (CMD) mechanism, with an energy barrier of 20.4 kcal/mol (TS1 relative to 2). The resulting five-membered rhodacycle 3 subsequently undergoes alkyne insertion to afford 5 via TS2, with an energy barrier of 14.5 kcal/mol (TS2 relative to 4) . In the red path, 2′′ undergoes C−H activation via TS1′′ to generate the corresponding rhodacycle 3′′, and the energy barrier for this step is calculated to be 22.6 kcal/mol (TS1′′ relative to 2′′). The 5

ACS Paragon Plus Environment

The Journal of Organic Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

subsequent alkyne insertion proceeds via TS2′′ with an energy barrier of 16.6 kcal/mol (TS2′′ relative to 4′′). Thus, in the black path of 2 to 5, the rate determining step is C−H activation with a barrier of 20.4 kcal/mol, and the rate determining step is also C−H activation with a barrier of 22.6 kcal/mol in the red path of 2′′ to 5′′. As a result, the black path is more favorable than the red one, and the possibility of C-coordination could be ruled out at this stage. On the basis of the calculated free energy profiles of the whole catalytic cycle, C−H activation is the rate-determining step (vide infra), and this is consistent with the large experimental kinetic isotope effects (KIE) (kH/kD = 5.3). Several examples of the similar C−H activation and olefin/alkyne insertion processes have been reported in recent years.17 As shown in Figure 2, the intermediate 4 is lower in energy by 11.5 kcal/mol than 4′′. To obtain further insight into the energy difference between intermediates 4 and 4′′, we performed NBO analysis. It is known that the π-back donation exist in the metal-alkyne complexes. The sum of the stabilization energies E(2) arising from the donation from πC≡C to Rh-5s and the π-back donation are nearly identical in intermediates 4 and 4′′ (Figure 3). In 4, the interactions between the lone pair of O and the vacant Rh-5s orbital were also observed and the total stabilization energy is 37.5 kcal/mol. However, in 4′′, the stabilization energy arising from the donation from πC=C to Rh-5s is relatively small (27.6 kcal/mol). These results reveal that the coordination of lone pair of electrons on the O atom to Rh brings about more stable complex, which makes the ketone rather than the alkene group the preferred directing group. Similar results were also observed for substrates 1b and 1c (see Figures S2-S3).

Figure 3. Key orbital interactions in 4 and 4′′. Alternatively, the reaction could also proceed via alkyne insertion into Rh−O bond to deliver the desired complex 5 (see Figure S5 for details). However, this path is kinetically unfavorable, with a forbidden high barrier. Additionally, the C−H activation happening at the site of alkene carbon atom adjacent to the dimethylamino group was also excluded in view of the high stationary points along this path (see Figure S6 for details). 6

ACS Paragon Plus Environment

Page 6 of 18

Page 7 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Organic Chemistry

3.2 Reductive elimination

Figure 4. Free energy profiles for the reductive elimination of 6. For the optimized structures of TS3 and TS3′′, hydrogen atoms are omitted for clarity. Intermediate 5 undergoes isomerization to afford a seven-membered rhodacycle 6. The ensuing reductive elimination leads to the formation of intermediate 9, an RhI species. As shown in Figure 4, the direct C(sp3)-C(sp2) reductive elimination via the transition state TS3′′ is found to be unfavorable (see the red path in Figure 4). This result differs markedly from the previous theoretical results, which indicates the pivotal influence of dimethylamino and diphenyl alkenyl groups of 6 on the process of reductive elimination.16 A previous theoretical study by Xu has demonstrated that ketone enolization is necessary before the reductive elimination in the Rh(III)-catalyzed cascade oxidative annulation of benzoylacetonitrile with alkynes.18 However, the underlying reason has not been deeply explored yet. As shown in Figure 4, ketone enolization converts 6 to a less stable complex 7, and the ensuing C(sp2)-C(sp2) reductive elimination is relatively facile, 19 with lower in free energy transition state TS3. Subsequently, complex 8 undergoes ketone enolization again to reach the intermediate 9. The transition state TS3′′ is higher in free energy than TS3 by 14.7 kcal/mol. A comparison of structures of the two transition states indicates two factors may be responsible for the remarkable energy difference. Firstly, TS3 involves a ligand π-Rh coordination, which is widely present in π-conjugative structure of the six-membered ring18,20 In TS3′′, the lack of π-system on the ligand makes the C(sp3)-C(sp2) reductive elimination unassisted. Secondly, the dihedral angle φC1-C2-C3-C4 in TS3′′ is 50.0°, which is significantly larger than the value in TS3 (31.1°), implying a stronger distortion within the formation of six-membered ring. Nucleus-independent chemical shifts (NICS) has been a popular method for the evaluation of magnetic aromaticity. To gain deep insights into the activation energy differences between TS3 and TS3′′, the spatial magnetic properties of the two transition states have been calculated by the gauge-invariant atomic orbitals (GIAO) perturbation method employing the NICS concept.21 The through-space NMR shieldings (TSNMRS) visualized as the iso-chemical-shielding surfaces (ICSSs) of TS3 and TS3′′ are given in Figure 5. The obvious π-electron delocalization due to the conjugation of C1=C2−C3=C4−C5=C6 bond in TS3 strengthens the anisotropy effect. As shown in 7

ACS Paragon Plus Environment

The Journal of Organic Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5, the six-membered ring in TS3 obviously exhibits the tendency to generate paratropic ring currents and therefore aromaticity. In contrast, the classical anisotropy of the corresponding six-membered ring in TS3′′ displays deshielding. In summary, our results show that the aromaticity of TS3 is stronger than that of TS3′′, which makes the C(sp2)−C(sp2) bond formation more favorable.

Figure 5. TSNMRS of TS3 and TS3′′ visualized as ICSSs (blue, red, orange, yellow, green represent 2 ppm, 4 ppm, 8 ppm, 16 ppm, 32 ppm shielding respectively; white represents −1 ppm deshielding). The free energy profiles for the reductive elimination of 6b and 6c are displayed in Figure 6. The C−C bond formations of the two intermediates are similar with that of intermediate 6a. The ketone enolization step is again necessary prior to the direct C(sp3)-C(sp2) reductive elimination. It is worth noting that, different from 6a and 6c, the reaction for 6b undergoes C−S bond cleavage, simultaneously forming new C−C bond. The result indicates that the DMSO dissociation pathway is more favorable, which can be attributed to the lower BDE of C=S bond in 1b compared to that of the C=C bond in 1a and C=N bond in 1c. The bond strength for the three kinds of double bonds is well revealed by the calculated Wiberg bond index (WBI): 1.6, 1.4 and 1.1 for 1a, 1b and 1c respectively.

8

ACS Paragon Plus Environment

Page 8 of 18

Page 9 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Organic Chemistry

Figure 6. Free energy profiles for the reductive elimination of 6b and 6c. Energies are relative to [Cp*RhOAc]+ + 1b and [Cp*RhOPiv]+ + 1c, respectively. 3.3 Hydroxyl/aldehyde functionalization Upon formation of Rh(I) complex 9, the ensuing β-H elimination occurs facilely via TS4, only crossing a barrier of 7.0 kcal/mol. The [Cp*RhH]+ released in this step reacts with the external oxidant AgOAc to afford [Cp*Rh]2+, which allows further occupation of anion OAc- to regenerate the catalytically active species [Cp*RhOAc]+. In the presence of H2O, the resultant intermediate 10 undergoes stepwise functionalization, including hydrolysis and deamination to generate aldehyde and hydroxyl units (see the red path in Figure 7). In TS5, the concerted hydrogen transfer to oxygen atom and attack of carbon atom occur, and the process spans an energy barrier of 17.8 kcal/mol (TS5 relative to 11). Subsequently, the HOAc-assisted intramolecular hydrogen transfer occurs with a negligible barrier, releasing the by-product HNMe2. The naphthalene derivative product contains both a nucleophilic aldehyde unit and an electrophilic hydroxyl unit in the same molecule after this stepwise functionalization. The energy profiles for the hydroxyl functionalization of 9b and 9c are shown in Figures S7-S8. Experimentally, in the absence of H2O, the target product P1 was formed in 25% yield. Our results reveal that HOAc molecule could also transfer its hydrogen atom to oxygen atom of carbonyl in intermediate 10 (see the blue path in Figure 7). This process is found to be facile, with an energy barrier of 7.8 kcal /mol (TS5′′ relative to 11′′). However, the following deamination with the assistance of HOAc is kinetically less favorable, with a higher energy barrier of 23.9 kcal/mol (TS6′′ relative to 13′′), indicating the important role of the surrounding water molecules in the hydroxyl/aldehyde functionalization process. Our theoretical results agree well with the 9

ACS Paragon Plus Environment

The Journal of Organic Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 18

experimental fact that water molecules can improve the product yield.

OH

H NMe 2 + [Rh] H Ph

Ph

9 1.5

OH O

O

TS4 5.5

NMe2

Ph

Ph

Ph

Ph

H [Rh]

+

Ph C

H H O O NMe2

10 10.1

H2O

HOAc

11 6.9

Ph

11' 8.0

NMe 2

Ph

Ph Ph

TS5' 0.2 C O

12' 1.9

O

OH

O

NHMe 2 C

Ph OH O

Ph

N H Me2 Ph O

Ph

TS6 0.2

13 0.2

HOAc OH

O

13' 4.9

HOAc

12 1.6

TS6' 28.8

O

N H O Me2 Ph

OH

TS5 10.9

O

C O H O O NMe2

OH

Ph

H

C O O

O N H Me2 Ph

H O NMe 2

O O

O

OH

O

O H N H Me Ph 2

P1+(MeCO)2O+HNMe2 0.7 O O

P1+HOAc+HNMe2 10.9

Ph

NMe 2 Ph

Ph Ph

Ph

hydrolysis

deamination

acidification

deamination

elimination

Figure 7. Free energy profiles for the hydroxyl/aldehyde functionalization of 9. Another pathway starting from rhodacycle 6, namely hydrolysis prior to reductive elimination, was also investigated (see Figure S9 for details). Accordingly, there are two possible routes after the formation of hydrate 15: (i) deamination followed by C(sp3)-C(sp2) reductive elimination, then β-H elimination to afford [Cp*RhH]+; (ii) direct C(sp3)-C(sp2) reductive elimination, followed by deamination and β-H elimination. The results in Figure S9 reveal that both routes are kinetically less favorable, with high barriers. Therefore, the possibility of hydrolysis prior to cyclization can be ruled out. At this stage, intermediates 6b/6c could alternatively undergo elimination of the DMSO/tBuNO before the reductive elimination (see Figure 8). For 6b, release of DMSO is facile, with a barrier of 12.9 kcal/mol. The following C(sp2)-C(sp2) reductive elimination also easily occurs, with a low activation energy (9.2 kcal/mol). However, 6c undergoes tBuNO elimination via TS7c with an energy barrier of 29.0 kcal/mol, generating a highly unstable intermediate 15b. This step is thermodynamically unfavorable with ∆G0 = 21.8 kcal/mol (15b relative to 6c). The differences between the results of 6b/6c are probably due to the intrinsically higher stability of C=N than C=S bond.

10

ACS Paragon Plus Environment

Page 11 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Organic Chemistry

Figure 8. Elimination of DMSO/tBuNO prior to reductive elimination in 6b/6c, with energies in kcal/mol. 3.4 The formation of six-membered versus five-membered ring To understand the competition between ketone and C=C, C=N, C=S bond, we studied the free energy profiles for the formation of six-membered and five-membered intermediates from the alkyne insertion intermediates 5, 5b and 5c (see Figure 9). The formation of six-membered rings has been discussed above. Alternatively, intermediate 5, 5b or 5c may undergo direct C(sp2)-C(sp2) reductive elimination via TS13, TS7b, or TS8c to generate the corresponding five-membered ring 20, 16b or 16c. Comparing the two routes from 5, the formation of six-membered ring via TS3 is 9.4 kcal/mol more accessible than the five-membered formation pathway via TS13 (7 → TS3, 8.4 kcal/mol vs 5 → TS13, 17.8 kcal/mol). In complex 5b or 5c, the formation of six-membered carbocycle is also kinetically favorable (6b → TS5b, 12.9 kcal/mol vs 5b → TS7b, 24.1 kcal/mol; 7c→TS3c, 3.9 kcal/mol vs 5c→TS8c, 12.8 kcal/mol). Therefore, the possibility for the generation of five-membered products from 5, 5b and 5c could be excluded. It should be noted that the five-membered product was formed in the Rh(III)-catalyzed carbocyclization of aryl ketones and alkynes (see eq 1 in Scheme 1), in which the external oxidant Cu(OAc)2 helps release the Rh catalyst via transmetalation and the resulting copper(I) five-membered ring undergoes protonation to generate the final product. We have also considered the possibility of 20, 16b and 16c evolving through ring expansion to 9 ,9b and 9c, respectively. However, all the three paths are kinetically unfavorable, with high energy barriers.

11

ACS Paragon Plus Environment

The Journal of Organic Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

five-membered formation

six-membered formation O

H

O

NMe2 + [Rh]

O +

[Rh] Ph

Ph Ph

9b

9

Ph

Ph

Page 12 of 18

O + N tBu [Rh] Ph

O [Rh] NMe2 +

O

O

O

Ph Ph

tBu Ph Ph TS8c

Ph

TS8b

TS14 44.0 TS9c 29.4

TS3c 19.5

9c 9.8

O

N +

G (kcal/mol)

9 1.5

[Rh]

S +

Ph TS14

9c

[Rh]

7c 15.6

TS3 17.4

8c 3.5 TS6b 3.1

7 9.0

15b 6.1

TS5b 1.5

8 9.9

TS7b 14.2

6c 8.3

TS13 9.0

6 5.5 6b 14.4

16b 4.9 20 3.4 16c 6.1

TS8c 5.5

5c 5 7.3 8.8 5b

TS8b 28.7

9.9 NMe2

NMe2

tBu

S

9b 85.5

O

+

O

[Rh] Ph Ph

TS13

[Rh] Ph Ph

20

Ph

TS7b

N

S

+ OO [Rh] Ph

+

+ OO [Rh] Ph

Ph

16b

O

O +

tBu

[Rh] Ph Ph

TS8c

N

O O + [Rh] Ph

Ph

16c

Figure 9. Free energy profiles of the competing pathways for the formation of six-membered and five-membered rings from intermediates 5, 5b and 5c. 3.5 The substituent effect on the reactivity and regioselectivity of cyclization of enaminones. As seen in Scheme 3, the cyclization of meta-substituted enaminones could principally occur at either para- or ortho-site. Experimental observations indicated that the methyl at the meta position of substrate 1d generated a single para-coupled product P3 with a yield of 75%, while the possible ortho-coupled product P3′′ was not observed. Interestingly, a mixture of para- and ortho-coupled products P4 and P4′′ were obtained in the cyclization of meta-fluoro-substituted substrate 1e with a high regioisomeric ratio (P4 : P4′′ = 1 : 0.15). Figure 10 shows the free energy profiles of the full catalytic cycle for para- and ortho-coupled with both meta-methyl- and meta-fluoro-substituted enaminones (see Figures S10−S11 for complete results). According to the energetic span model introduced by Kozuch and Shaik,22 the energetic span (δE) that serves as the apparent activation energy of catalytic cycle is defined by the energy difference between turnover-frequency(TOF)-determining transition state (TDTS), the TOF-determining intermediate (TDI), and the reaction driving force. For meta-methyl-substituted enaminones, the TDTSs are the reductive elimination transition states TS3d and TS3d′′, respectively, and the TDIs are the alkyne insertion species 5d and 5d′′, respectively. The apparent activation energy for para-coupled is 27.2 kcal/mol (TS3d − 5d), and the value is 30.1 kcal/mol (TS3d′′ − 5d′′) for ortho-coupled. By comparing the energetic spans of the two free energy surfaces (∆∆G‡ = 2.9 kcal/mol), the result was in good agreement with experiment, in which only product P3 was formed. The final apparent activation energies for para- and ortho-coupled with meta-fluoro-substituted enaminones are 27.2 kcal/mol and 28.0 kcal/mol, respectively, indicating the possibility of cyclization at the ortho position. The calculated activation energy also agrees qualitatively well with the experimental 12

ACS Paragon Plus Environment

Page 13 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Organic Chemistry

value ( P4 : P4′′=1 : 0.15).

ortho H

OH

OH O

O

O

Me

Me

NMe2 para 1d

ortho H F

Ph

Ph Ph

H

NMe2 H

OH O

OH O

O

para

Ph Ph P3' not observed Me

Ph P3 75%

[Rh]

F Ph

Ph Ph

Ph [Rh]

Ph P4 64%

1e

F

Ph P4'

(P4:P4' = 1 : 0.15)

Scheme 3. Substituent effect on reactivity and regioselectivity of cyclization of enaminones to naphthalenes.

Figure 10. Free energy profiles of the full catalytic cycle for the cyclization of meta-methy-substituted and meta-fluoro-substituted enaminones.

Figure 11. Noncovalent interactions of meta-methy-substituted and meta-fluoro-substituted enaminones. To gain deep insights into the difference between the two substituents on regioselectivity, the 13

ACS Paragon Plus Environment

The Journal of Organic Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

noncovalent interactions in TDTSs are analyzed by the reduced density gradient (RDG) method (see Figure 11). In the case of the bulky methyl group, weaker repulsion between the phenyl ring and hydrogen atom is observed in TS3d. The corresponding transition state TS3d′′ suffers a stronger repulsion, which accounts for the energy difference between TS3d and TS3d′′. In contrast, the less sterically hindered fluoro group could not completely prevent the formation of ortho-coupled regioisomer due to the nearly identical weaker repulsion existence in both TS3e and TS3e′′. 5. Conclusions In the current study, we have performed a DFT mechanistic study on the Rh(III)-catalyzed coupling of alkynes with enaminones, sulfoxonium ylides or α-carbonyl-nitrones via dual directing groups. All the three reactions begin with C−H activation and the resulting cyclometalated intermediates undergo alkyne insertion to generate the seven-membered rings. It should be noted that the ketone instead of the double bond acts as the directing group in these steps. The O-bound to C-bound tautomerization is crucial prior to the subsequent reductive elimination. C(sp2)-C(sp2) reductive elimination is more favorable than the originally proposed direct reductive elimination mechanism in the three reactions. For substrate 1a or 1c, C(sp2)-C(sp2) reductive elimination prefers taking place before the elimination of dimethylamino or nitrone group, whereas substrate 1b prefers undergoing the release of dimethylsulfoxide prior to the reductive elimination. The mechanistic difference can be explained by the fact that the C=C and C=N bonds are more stable than C=S bond. The C−H activation is the rate-determining step in all the three catalytic cycles. Our computational results agree well with experimental kinetic isotope effects (KIE) study. By comparing the coordination of ketone versus C=C, C=N, or C=S bond, as well as the chemoselectivity concerning the six- versus five-membered formation, the roles of dual directing groups in the three substrates were uncovered. The present results confirm that the ketone acts as the role of the directing group while the C=C, C=N, or C=S bond serves as the cyclization site. The competition relationship of the dual directing groups is expected to be useful for substrates design in the synthesis of diverse carbocycles. Associated Content Supporting Information Additional computational results, energies, and Cartesian coordinates of all optimized structures. This material is available free of charge via the Internet at http://pubs.acs.org. Author Information Corresponding Authors *E-mail for J.-B. L.: [email protected] *E-mail for D.-Z. C.: [email protected] Notes The authors declare no competing financial interest. Acknowledgments 14

ACS Paragon Plus Environment

Page 14 of 18

Page 15 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Organic Chemistry

This work was financially supported by the National Natural Science Foundation of China (NSFC No. 21375082) of China. Reference (1) (a) Segawa, Y.; Maekawa, T.; Itami, K. Angew. Chem., Int. Ed. 2015, 54, 66−81. (b) Huang, H.-W.; Ji, X.-C.; Wu, W.-Q.; Jiang, H.-F. Chem. Soc. Rev. 2015, 44, 1155−1171. (c) Huang, Z.-X.; Lim, H. N.; Mo, F.-Y.; Young, M. C.; Dong, G.-B. Chem. Soc. Rev. 2015, 44, 7764−7786. (d) Fukui, Y.; Liu, P.; Liu, Q.; He, Z.-T.; Wu, N.-Y.; Tian, P.; Lin, G.-Q. J. Am. Chem. Soc. 2014, 136, 15607−15614. (e) Dateer, R. B.; Chang, S. J. Am. Chem. Soc. 2015, 137, 4908−4911. (f) Mazuela, J.; Banerjee, D.; Backvall, J. J. Am. Chem. Soc. 2015, 137, 9559−9562. (g) Kawaguchi, Y.; Yasuda, S.; Kaneko, A.; Oura, Y.; Mukai, C. Angew. Chem., Int. Ed. 2014, 53, 7608−7612. (h) McNally, A.; Haffemayer, B.; Collins, B. S. L.; Gaunt, M. J. Nature 2014, 510, 129−133. (i) Huang, Z.-X.; Lim, H. N.; Mo, F.-Y.; Young, M. C.; Dong, G.-B. Chem. Soc. Rev. 2015, 44, 7764−7786. (j) Li, B.-J.; Shi, Z.-J. Chem. Soc. Rev. 2012, 41, 5588−5598. (k) Balcells, D.; Clot, E.; Eisenstein, O. Chem. Rev. 2010, 110, 749−823. (l) Kuhl, N.; Hopkinson, M. N.; WencelDelord, J.; Glorius, F. Angew. Chem., Int. Ed. 2012, 51, 10236−10254.(m) Engle, K. M.; Mei, T. S.; Wasa, M.; Yu, J.-Q. Acc. Chem. Res. 2012,45, 788−802. (2) (a) Jayakumar, J.; Parthasarathy, K.; Cheng, C. Angew. Chem., Int. Ed. 2012, 51, 197−200. (b) Lian, Y.-J.; Huber, T.; Hesp, K. D.; Bergman, R. G.; Ellman, J. A. Angew. Chem., Int. Ed. 2013, 52, 629−633. (c) Qi, Z.-S.; Li, X.-W. Angew. Chem., Int. Ed. 2013, 52, 8995−9000. (d) Burns, D. J.; Lam, H. W. Angew. Chem., Int. Ed. 2014, 126, 10089−10093. (e) Lian, Y.-J.; Bergman, R. G.; Lavis, L. D.; Ellman, J. A. J. Am. Chem. Soc. 2013, 135, 7122−7125. (f) Shi, Z.-Z.; Koester, D. C.; Boultadakis-Arapinis, M.; Glorius, F. J. Am. Chem. Soc. 2013, 135, 12204−12207. (g) Liu, B.-Q.; Song, C.; Sun, C.; Zhou, S.-G.; Zhu, J. J. Am. Chem. Soc. 2013, 135, 16625−16631. (h) Ueura, K.; Satoh, T.; Miura, M. Org. Lett. 2007, 9, 1407−1409. (i) Unoh, Y.; Hashimoto, Y.; Takeda, D.; Hirano, K.; Satoh, T.; Miura, M. Org. Lett. 2013, 15, 3258−3261. (3) (a) Muralirajan, K.; Parthasarathy, K.; Cheng, C. Angew. Chem., Int. Ed. 2011, 50, 4169−4172. (b) Patureau, F. W.; Besset, T.; Kuhl, N.; Glorius, F. J. Am. Chem. Soc. 2011, 133, 2154−2156. (c) Yu, S.-J.; Liu, S.; Lan, Y.; Wan, B.-S.; Li, X.-W. J. Am. Chem. Soc. 2015, 137, 1623−1631. (d) Mehta, V. P.; García-López, J.; Greaney, M. F. Angew. Chem., Int. Ed. 2014, 126, 1555−1559. (e) Shi, X.-Y.; Li, C.-J. Org. Lett. 2013, 15, 1476−1479. (4) Gandeepan, P.; Cheng, C. Org. Lett. 2013, 15, 2084−2087. (5) Zhou, S.-G.; Wang, J.-H.; Wang, L.-L.; Song, C.; Chen, K.-H.; Zhu, J. Angew. Chem., Int. Ed. 2016, 55, 9384−9388. (6) (a) Wang, Q.; Xu, Y.-W.; Yang, X.-F.; Li, Y.-Y.; Li, X.-W. Chem. Commun. 2017, 53, 9640−9643. (b) Xu, Y.-W.; Yang, X.-F.; Zhou, X.-K.; Kong, L.-H.; Li, X.-W. Org. Lett. 2017, 19, 4307−4310. (7) (a) Becke, A. D. J. Chem. Phys. 1993, 98, 5648–5652. (b) Lee, C.; Yang. W.; Parr, R. G. Phys. Rev. B: Condens. Matter 1988, 37, 785–789. (8) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. J. Phys. Chem. B 2009, 113, 6378−6396. (9) (a) Figgen, D.; Peterson, K. A.; Dolg, M.; Stoll, H. J. Chem. Phys. 2009, 130, 164108. (b) Peterson, K. A.; Figgen, D.; Dolg, M.; Stoll, H. J. Chem. Phys. 2007, 126, 124101.

15

ACS Paragon Plus Environment

The Journal of Organic Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(10)The influence of basis sets on the results was considered, and the corresponding results were given in Figure S1. (11) (a) Zhao, Y.; Truhlar, D. G. J. Chem. Phys. 2006, 125, 194101. (b) Zhao, Y.; Truhlar, D. G. Theor. Chem. Acc. 2008, 120, 215−241. (12)Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. J. A.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford, CT, 2009. (13)Glendening, E. D.; Badenhoop, J. K.; Reed, A. E.; Carpenter, J. E.; Bohmann, J. A.; Morales, C. M.; Weinhold, F. NBO 5.0 ; Theoretical Chemistry Institute, University of Wisconsin, Madison, 2001. (14)Legault, C. Y. CYLview, 1.0b, Université de Sherbrooke, Canada, 2009, http://www.cylview.org. (15)Humphrey, W.; Dalke, A.; Schulten, K. J. Mol. Graphics 1996, 14, 33−38. (16)Yu, S.-J.; Liu, S.; Lan, Y.; Wan, B.-S.; Li, X.-W. J. Am. Chem. Soc. 2015, 137, 1623−1631. (17) (a) Ajitha, M. J.; Huang, K.-W. Organometallics 2016, 35, 450−455. (b) Zhang, Z.-C.; Yang, S.-W.; Li, J.; Liao, X.-J. J. Org. Chem. 2016, 81, 9639−9646. (c) Xing, Y.-Y.; Liu, J.-B.; Tian, Y.-Y.; Sun, C.-Z.; Huang, F.; Chen, D.-Z. J. Phys. Chem. A 2016, 120, 9151−9158. (18) Fu, X.-N.; Shang, Z.-F.; Xu, X.-F. J. Org. Chem. 2016, 81, 8378−8385. (19) (a) Mann, G.; Baranano, D.; Hartwig, J. F.; Rheingold, A. L.; Guzei, I. A. J. Am. Chem. Soc. 1998, 120, 9205−9219. (b) Cohen, R.; Milstein, D.; Martin, J. M. L. Organometallics 2004, 23, 2336−2342. (c) Culkin, D. A.; Hartwig, J. F. Organometallics 2004, 23, 3398−3416. (d) Pérez-Rodríguez, M.; Braga, A. A. C.; Garcia-Melchor, M.; Pérez-Temprano, M. H.; Casares, J. A.; Ujaque, G.; de Lera, A. R.; Álvarez, R.; Maseras, F.; Espinet, P. J. Am. Chem. Soc. 2009, 131, 3650−3657. (e) Racowski, J. M.; Dick, A. R.; Sanford, M. S. J. Am. Chem. Soc. 2009, 131, 10974−10983. (f) Pérez-Rodríguez, M.; Braga, A. A. C.; de Lera, A. R.; Maseras, F.; Álvarez, R.; Espinet, P. Organometallics 2010, 29, 4983−4991. (g) Ji, C.-L.; Hong, X. J. Am. Chem. Soc. 2017, 139, 15522−15529. (20) (a) Yu, Z.-X.; Cheong, P. H.; Liu, P.; Legault, C. Y.; Wender, P. A.; Houk, K. N. J. Am. Chem. Soc. 2008, 130, 2378−2379. (b) Zuidema, E.; van Leeuwen, P. W. N. M.; Bo, C. Organometallics 2005, 24, 3703−3710. (21) (a) Klod, S.; Kleinpeter, E. J. Chem. Soc., Perkin Trans. 2 2001, 1893−1898. (b) Kleinpeter, E.; Koch, A. J. Phys. Chem. A 2012, 116, 5674−5680. (c) 16

ACS Paragon Plus Environment

Page 16 of 18

Page 17 of 18 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Organic Chemistry

Fallah-Bagher-Shaidaei, H.; Wannere, C. S.; Corminboeuf, C.; Puchta, R.; Schleyer, P. V. R. Org. Lett. 2006, 8, 863−866. (22) (a) Kozuch, S.; Shaik, S. Acc. Chem. Res. 2011, 44, 101−110. (b) Kozuch, S. WIREs Comput. Mol. Sci. 2012, 2, 795−815.

17

ACS Paragon Plus Environment

The Journal of Organic Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table of Contents Graphic

18

ACS Paragon Plus Environment

Page 18 of 18