Xuexiang
Ma
,
Aili
Feng
,
Chengbu
Liu
and
Dongju
Zhang
*
Key Lab of Colloid and Interface Chemistry, Ministry of Education, Institute of Theoretical Chemistry, School of Chemistry and Chemical Engineering, Shandong University, Jinan 250100, P. R. China. E-mail: zhangdj@sdu.edu.cn; Tel: +86 0531 88365833
First published on 12th November 2021
Density functional theory calculations were performed on a prototype of three-component reactions involving aryl iodides, 2,6-substituted aryl bromides, and acrylates to understand the construction of axially chiral biaryls through the palladium/chiral norbornene cooperative catalysis (Nat. Catal.2020, 3, 727–733). The reaction is established to occur in sequence through oxidative addition of aryl iodide to Pd(0), norbornene insertion, C–H activation, oxidative addition of aryl bromide to Pd(II), C–C reductive elimination, norbornene extrusion, alkene insertion, and β-H elimination. The C–H activation is characterized as the rate-determining step, and the oxidative addition of aryl bromide on Pd(II) control the enantioselectivity of the reaction. The construction of more favorable (S)-C–C axial chirality is attributed to the smaller steric repulsion between the ortho-methyl group of aryl bromide and the norbornene moiety. The theoretical results rationalize the reaction sequence of aryl iodide and aryl bromide on both Pd(0) and Pd(II) and clarify the substantial role of the ortho-ester group in the aryl bromide.
The palladium/norbornene (Pd/NBE) cooperative catalysis pioneered by Catellani17 and later modified by Lautens,18 Yu,19,20 Dong,21,22 and others,23–25 is expected to be a promising approach to the construction of biaryl motifs. However, the asymmetric version of the Catellani reaction has barely been reported, until the first enantioselective construction of axially chiral biaryls was achieved by using a chiral phosphine ligand.26 Inspired by the recent progress, Zhou et al.27 successfully achieved axially chiral biaryls by using palladium and chiral norbornene (NBE*) catalyzed three-component cross-coupling reactions. As indicated in Scheme 1 by a representative example, this new strategy transforms aryl iodide 1, 2,6-substituted aryl bromide 2 and acrylate 3 (terminating reagent) into the axially chiral biaryl P in good yield (63%) with excellent enantioselectivity (97% ee). The optimal reaction conditions used Pd(OAc)2 as the catalyst, (1S,4R)-2-ethyl ester-substituted NBE* as the mediator, and K2CO3 as the base at 105 °C in acetonitrile (MeCN).
![]() | ||
Scheme 1 Pd/NBE* cooperative three-component cascade reaction for the construction of axially chiral biaryls. |
To understand the reactivity, Zhou and co-workers27 have tentatively proposed a plausible mechanism. As shown in Scheme 2, the reaction is assumed to start from the oxidative addition of aryl iodide 1 on Pd(0) species, generating the aryl Pd(II) intermediate A. Subsequently, the chiral norbornene (NBE*) inserts into Pd–C bond to give intermediate B, which then evolves into the five-membered aryl–norbornyl–palladacycle (ANP) intermediate Cvia C–H bond activation. Oxidative addition of 2,6-substituted aryl bromide 2 on Pd(II) center of C leads to Pd(IV) species D, from which C–C reductive elimination occurs to form axially chiral Pd(II) species E. Then ipso-carbopalladation intermediate F can be obtained via extrusion of NBE* from E (β-carbon elimination). Finally, F reacts with the terminating reagent acrylate 3 to furnish the final axially chiral biaryl P with regeneration of Pd(0) catalyst. The proposed mechanism provides a fundamental understanding of the construction of axial chirality. However, the mechanism details keep unexplored, and the observed excellent enantioselectivity still has not been rationalized well. In addition, it is noted that the proposed mechanism involves two oxidation addition steps: the first is for the substrate (aryl iodide) on Pd(0) and the second is for the arylation reagent (2,6-substituted aryl bromide) on Pd(II). Why is it such a reaction sequence? In other words, with simultaneous presence of two aryl halides, why the first oxidation addition gives preference to aryl iodide and the second to aryl bromide. Herein, we report a density functional theory (DFT) based computational study, from which we hope to clarify these essential questions.
![]() | ||
Scheme 2 Plausible mechanism for the Pd/NBE* cooperative three-component cascade process to construct axial chirality biaryls. |
According to Scheme 2, the next step is the C–H bond activation to form the ANP intermediate. Based on previous studies by groups of Lan41 and Hong42 on the similar Pd/NBE cooperative Catellani reactions, this process is expected to occur with the assistance of a base. Here K2CO3 (2.5 equiv., Scheme 1) was used as proton abstraction reagent in the C–H bond activation. Upon coordination of K2CO3, intermediate 8 is brought to a tetracoordinated Pd(II) complex 9 with an energy release of 22.7 kcal mol−1. The subsequent C(sp2)–H activation proceeds according to a concerted metalation–deprotonation (CMD) mechanism43,44 to form the ANP (10) viaTS3 with a barrier of 27.7 kcal mol−1.
Other three potential pathways for the C–H bond activation are also shown in Fig. 1 for comparison. The red describes the pathway with the removal of KI prior the C–H activation to give species 9′, which then carries out the C–H bond activation viaTS3′. However, TS3′ is indicated to be energetically less favorable than TS3 by 10.4 kcal mol−1. This can be attributed to the large energy requirement for cleaving the strong Pd–I bond to release of KI.
The blue and brown pathways correspond to the C–H bond activation processes by using acetate anion as the proton abstraction reagent, as done in several recent works of Pd/NBE cooperative Catellani reactions.41,42 Identified transition states in these two pathways are denoted as TS3′′ and TS3′′′, which are counterparts of TS3 and TS3′, respectively. Our calculations indicate that these two pathways are also energetically less favorable than the black one. This fact, together with the low concentration of the acetate anion from a catalytic amount of Pd(OAc)2, suggests that the acetate anion is not appropriate for assisting the proton abstraction process. Thus, we believe that the C–H metalation occurs with assistance of K2CO3 which was used in the experiment in a stoichiometric amount (Scheme 1). The situation is similar to previous theoretical reports by Bi45 and Liu.46
After clarifying the formation mechanism of ANP intermediate, we now consider the question mentioned above: why the aryl iodide (as substrate) rather than the aryl bromide (as arylation reagent) carries out the initial oxidation addition on the Pd(0) center. To clarify the question, we carried out calculations for the oxidative addition of aryl bromide and the results are compared with that of the aryl iodide in Fig. 2a. It is found that the oxidation addition of aryl bromide viaTS1′ needs to overcome a barrier of 16.1 kcal mol−1, which is larger by 3.3 kcal mol−1 than that for the aryl iodide. Thus the Pd(0) species selectively reacts with aryl iodide instead of with aryl bromide, as confirmed by previous experimental and theoretical studies.47,48 This can be attributed to the larger bond energy of C–Br bond than C–I bond (78.4 vs. 67.3 kcal mol−1 from the present calculations) and the strong electron-withdrawing capacity of the –CO2Me group in the aryl bromide, which remarkably reduces the negative charge of the carbon atom of the C–Br bond (−0.08 vs. −0.23 e in the C–I bond) and thus reduces its nucleophilic attack ability on the positively charged palladation center.
![]() | ||
Fig. 2 (a) Aryl iodide vs. aryl bromide oxidative addition, and (b) NBE*vs. alkene insertion into the Pd(II)–C(aryl) bond. Calculated NBO charges are given in e. |
It should be noted that instead of NBE*, the alkene (3) can also insert into the Pd–C bond of aryl–palladium iodine intermediate 6, resulting in the coupling product of 1 with 3. This pathway would result in the direct coupling product of 1 with 3. The calculated results are shown in Fig. 2b, and the reaction occurs viaTS2′. This pathway is located to be 1.6 kcal mol−1 less favorable than the NBE* insertion pathway viaTS2. Inserting NBE* or 3 into the Pd–C bond involves the charge transfer from the highest occupied molecular orbital (HOMO) of NBE* or 3 to the lower unoccupied molecular orbital (LUMO) of 6. As shown in Fig. 3, the calculated HOMO–LUMO gap is 4.24 eV for the NBE* insertion and 4.74 eV for the alkene insertion, implying that the NBE* insertion is more favorable than the alkene insertion. This is in agreement with the experimental obsevation27 where the direct Heck coupling product between 1 and 3 is in very low yield.
![]() | ||
Fig. 4 Calculated Gibbs energy profiles for the formation of (S)-C–C chiral axis (the right profile) in comparison with that of (R)-C–C chiral axis (the left profile). |
In contrast, the calculated result for the construction of (R)-C–C axial chirality is shown in the left profile of Fig. 4. The barriers of the C–Br bond oxidative oxidation viaTS6 and the C–C reductive elimination viaTS7 are 22.2 and 30.5 kcal mol−1 (from 12 to TS6/TS7). Thus, the construction of (R)-C-C axial chirality involves much higher barrier than that of the (S)-C–C axial chirality, 30.5 vs. 20.2 kcal mol−1, exclusively favoring the latter. The optimized geometries show that the enantioselectivity mainly originates from the different steric repulsions in transition states involved in these two branches. As shown in Fig. 5, in TS7, there is a distinct steric repulsion between the ortho-methyl group of aryl bromide and the bridged methylene of the norbornene (H–H distances: 1.90 and 2.13 Å), while in TS4 and TS5, no remarkable steric hindrance is observed (the shortest H–H distance is 2.16 Å in TS5). Therefore, the (S)-C–C axial chirality pathway is more favored than the (R)-C–C axial chirality pathway. We also investigated the other two oxidation addition pathways (TS4′ and TS6′), which involve the same reaction scheme but with opposite Br, O coordination position compared to TS4 and TS6. It is found that TS4′ and TS6′ are energetically less favorable by ∼8 kcal mol−1 than TS4 and TS6, respectively (Fig. S2†).
![]() | ||
Fig. 5 Optimized transition state geometries involved in the construction of C–C chiral axis. The selected bond distances are given in Å. |
Furthermore, the other three potential pathways from ANP according to Scheme 3 were also calculated: the C–C reductive elimination (Fig. 6a), the oxidative addition of aryl iodide to the Pd(II) center (Fig. 6b), and the terminating reagent insertion into one of two Pd–C bonds in ANP (Fig. S3 and S4†). The C–C reductive elimination viaTS8 is found to involve a barrier as high as 41.9 kcal mol−1, explaining the fact that no the C–C reductive elimination by-product was observed in the experiment. This may be attributed to the large strain of the 3,4-bridged ring and the steric repulsion between the palladium center and –CO2Et group of the norbornene moiety in TS8. The oxidative addition of aryl iodide to the Pd(II) center viaTS9 requires to overcome a barrier of 21.2 kcal mol−1, which is higher by 7.7 kcal mol−1 in comparison with the oxidative addition of aryl bromide viaTS4, implying that this process is also inaccessible. This situation is different from the oxidative addition of aryl iodide and aryl bromide on the Pd(0) center discussed above, where the aryl iodide is more prone to oxidative addition than aryl bromide due to the weaker C–I bond. Furthermore, the possibility for the terminating reagent insertion into Pd–C bond must also be ruled out due to high barriers involved (Fig. S3 and S4†). These results indicate that from ANP, the energetically most favorable reaction is the oxidation addition of aryl bromide on Pd(II).
![]() | ||
Fig. 6 Calculated Gibbs energies (in kcal mol−1) of the direct reductive elimination (a) and oxidative addition (b) pathways from ANP. |
The different reactivity of these two aryl halides towards the Pd(II) center can be understood by comparing the optimized geometries of TS4 and TS9. In TS4, the ortho-ester group of reacting aryl bromide coordinates to Pd(II) center to form a fused five–five-membered palladacycle structure, while such a fused bis-palladacycle structure is absent in TS9. Therefore, we assume that the bis-palladacycle structure plays a substantial role for stabilizing the oxidative addition transition state of aryl bromide (TS4). To test this hypothesis, we replaced the ortho-ester group in TS4 with a hydrogen atom and added an acetonitrile (solvent) molecule which enables the Pd(II) to keep its quasi-planar tetracoordinated geometry. Upon these substitutions, the optimized oxidative addition transition state is indicated as TS4′′, a mono-palladacycle structure and is compared with the bis-palladacyclic TS4 in Fig. 7. It is found that the barrier of the oxidative addition increases to 24.6 kcal mol−1 (from ANP to TS4′′), confirming the intrinsic role of the bis-palladacycle for stabilizing the oxidative addition transition state. On the other hand, the calculated barriers crossing TS4′′ and TS9, 24.6 and 21.2 kcal mol−1 respectively, support the generally accepted view that the oxidative addition of aryl bromide on a metal center is more difficult than aryl iodine. It is noted that in TS4, there is another fused 5/6-ring, which presents a quasi-planar motif and thus remarkably increase the rigidity and stability of the structure. The result is consistent with the previous study by our group on the palladium-catalyzed carboxylate-assisted C(sp3)–H arylation, wherein the unique 5/6-ring transition state was found to account for the facile C–C reductive elimination.49
![]() | ||
Fig. 7 Optimized geometries of TS4 and TS4′′ with selected bond distances (red values, in Å). All hydrogen atoms have been omitted for clarity. |
From these discussion, it is clear that the ortho-ester group in the aryl bromide is of crucial importance, which makes the oxidative addition of the aryl bromide on the Pd(II) more facile than that of the aryl iodine. So the Pd(II) species can selectively react with the aryl bromide instead of the aryl iodide. This substrate characteristics of aryl bromine have been also shown in several reports.12,41,42
Subsequently, 3 replaces the NBE* in 17, giving the olefin-Pd complex 18, from which the CC inserts into the Pd–C(sp2) bond. Both the 1,2- and 2,1-insertion modes have been taken into account, as shown in Scheme 4, where each mode contains two different pathways depending on –CO2tBu group in 3 pointing to or away from the chiral axis. It is found that two 2,1-insertion modes involve almost isoenergetic transition states (−55.9 kcal mol−1 for TS11 and −56.0 kcal mol−1 for TS14), and are energetically more favorable than two corresponding 1,2-insertion modes (−52.2 kcal mol−1 for TS11′ and −46.7 kcal mol−1 for TS14′). This can be attributed to both the larger π electron density of C2 than C1 due to the electron-withdrawing property of –CO2tBu group, favoring the nucleophilic attack of C2 to Pd(II), and the smaller steric hindrance between the phenyl moiety and the –CO2tBu group in the 2,1-insertion mode than that in the 1,2-insertion.
![]() | ||
Scheme 4 Four different modes of tert-butyl acrylate insertion into Pd–C(sp2) bond with calculated relative Gibbs energies (in kcal mol−1). |
In Fig. 8, only two favorable 2,1-insertion modes viaTS11 and TS14 are given, and the calculated barriers are 5.3 and 4.9 kcal mol−1, respectively. The resultant intermediates are 19 and 22. Finally, the catalytic cycle is closed through the β-H elimination process. From 19, β-H1 elimination viaTS12 (−69.3 kcal mol−1) affording the E-configuration product P is calculated to be more favorable by 4.1 kcal mol−1 than β-H2 elimination viaTS13 (−65.2 kcal mol−1) affording the Z-configuration product P′. The fact can be attributed to the greater steric hindrance between the phenyl moiety and the ester group –CO2tBu in TS13. Similarly, from 22, transition states resulting P and P′ are denoted as TS15 and TS16, which are less stable by 2.3 and 6.8 kcal mol−1 than TS12 and TS13, respectively. Thus 22 is expected to be less possible for providing the desired axially chiral product. From Fig. 8, it is clear that the β-H elimination viaTS12, offering the (S)-type, E-configuration axially chiral biaryl product, is the most favorable pathway.
To show a complete mechanistic picture of the Pd/NBE* cooperative three-component cascade reaction, we combine the main results in Fig. 1, 4, and 8 into Fig. 9. The construction of axial chirality is achieved via eight steps: Ar–I oxidative addition to Pd(0), norbornene insertion, C–H activation, Ar–Br oxidative addition to Pd(II), C–C reductive elimination, norbornene extrusion, alkene insertion, and β-H elimination. The C–H activation is identified as the rate-determining step of the reaction with a barrier of 27.7 kcal mol−1, which matches the experimental conditions (105 °C). The oxidative addition of Ar–Br on Pd(II) followed by the C–C reductive elimination favors the construction of the (S)-C–C axial chirality. This situation is different from the recent report by Lan et al.41 for a similar three-component Catellani reaction, where the rate-determining step of the reaction was identified as the oxidative addition of aryl bromide.
![]() | ||
Fig. 9 Complete mechanistic picture of the Pd/NBE* cooperative three-component cascade reaction for construction of the (S)-C–C axially chiral biaryl. |
Finally, we try to rationalize the experimentally observed enantiomeric excess (ee) value of the product, 97%. From our calculations (Fig. 4), the barriers for constructing the (S)- and (R)-C–C chiral axes are 20.2 and 30.5 kcal mol−1, respectively. The large difference (10.3 kcal mol−1) implies the reaction exclusively favors the (S)-C–C axially chiral biaryl. In this case, how is a small amount of R-type product produced? We conjecture that the R-type product may originate from the S–R isomerization. The calculated barrier of the C–C chiral axis rotation is 34.1 kcal mol−1 (Fig. 10). The present reaction was carried out under 105 °C, which provides a potential possibility of the S–R isomerization.
![]() | ||
Fig. 10 Calculated Gibbs energies (in kcal mol−1) for the S–R isomerization by rotating the chiral axis. |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1cy01863a |
This journal is © The Royal Society of Chemistry 2022 |