Stereocontrol through Synergistic Catalysis in the Enantioselective α

Jan 8, 2018 - We speculate that the computational model presented in this paper may provide a guideline to a rational design of catalyst combinations ...
0 downloads 0 Views 2MB Size
Article Cite This: J. Org. Chem. 2018, 83, 1304−1311

pubs.acs.org/joc

Stereocontrol through Synergistic Catalysis in the Enantioselective α‑Alkenylation of Aldehyde: A Computational Study Mahendra Patil* UM-DAE Centre for Excellence in Basic Sciences, Health Centre, University of Mumbai, Vidyanagari Campus, Kalina, Santacruz (East), Mumbai 400098, India S Supporting Information *

ABSTRACT: We describe the computational study of an interesting class of reactions involving synergistic action of chiral amine and Cu catalyst. The stereoselectivity-determining step of the enantioselective α-alkenylation of aldehyde has been investigated using the density functional theory methods to gain insight into the origin of the product selectivity. We found that the catalytic activation of reactants in the form of enamine and alkenyl Cu(III) intermediates significantly reduces the activation barrier of the addition step through an improved interaction between these two intermediates at the transition state (referred as enamine and Cu catalyst fragments of the transition state in the text). The transition-state stabilization through interaction between catalytic fragments, as demonstrated by the interaction/distortion model, clearly outperforms destabilization incurred due to the distortion of catalytic fragments and hence is recognized as a major factor contributing to the high stereoselectivity of reaction. Furthermore, the metal−enamine interaction described through the Cu···C7 distance is identified as a vital noncovalent interaction at the transition state. Our calculations show that the catalytic (covalent) activations and metal−enamine interaction can operate in tandem to amplify the net interaction between two catalytic fragments. The cooperative nature of these interactions is also reflected in the trend of interaction energies, which show a large variation with a subtle change in the metal−enamine interaction. Our computational model verified for the different catalytic combinations of chiral amine and Cu catalysts successfully rationalizes the experimentally observed enantioselectivity.



Asymmetric α-alkenylation of aldehydes, an elaboration of olefin framework adjacent to a carbonyl group, represents one of the most desirable, however very challenging processes for the synthesis of optically active α-substituted carbonyl compounds (Scheme 1). So far, several novel strategies have

INTRODUCTION In recent years, synergistic catalysis involving simultaneous activation of reacting partners using two discrete catalysts has spurred newer avenues to control the stereochemical outcome of chemical reactions.1 A principal feature of synergistic catalysis, which underscores its differences compared to the other modes of catalysis, is the nature of activation offered by two independent catalysts to the different reactants, allowing bimolecular reactions via HOMO(nucleophile) as well as LUMO(electrophile) activation in the same chemical environment. Owing to unprecedented reactivity achieved by the synergistic combination of catalysts, this catalytic protocol has been successfully applied for the development of new asymmetric transformations.2 While the application of this multicatalytic approach in chemical synthesis has witnessed explosive growth in recent times, the efforts directed toward elucidating its mechanism using computational methods are rather very limited.3 Herein, we investigate the selectivitydetermining step of the α-alkenylation of aldehydes which involves the synergistic combination of chiral amine and Cu catalyst. The mechanistic details of synergistic functioning of chiral amine and Cu catalyst would be crucial, in particular to understand how two catalysts work together to achieve the high enantioselectivity in the reaction. © 2018 American Chemical Society

Scheme 1. Enantioselective α-Alkenylation of Aldehydes with Alkenyl Boronic Acids5 or Alkenyl Iodonium Salts6

been designed to add an alkenyl group at the α-position of carbonyl group.4 Very recently, MacMillan and co-workers elegantly demonstrated the scope of dual catalysis for the αfunctionalization of aldehydes by combining the copper and chiral amine catalysts in a single process. In these reports, Macmillan and his group have described the enantioselective αalkenylation of aldehydes using suitable coupling partners such Received: November 7, 2017 Published: January 8, 2018 1304

DOI: 10.1021/acs.joc.7b02822 J. Org. Chem. 2018, 83, 1304−1311

Article

The Journal of Organic Chemistry

effective core potential and the associated double-ζ basis set for Br.15 Solvent effects were incorporated in the calculations using the SMD solvation model16 and ethyl ethanoate (dielectric constant: 5.98) as solvent. We have also included empirical dispersion corrections in calculations.17 This approach is designated as “M062X” in the text. Since the discussion in this paper involves the subtle changes in the structural features and associated energy variations of the diastereomeric transition states, we avoid including energetics obtained for the different basis set (6-311+G**) through single point calculations in the text. The relevant information is provided in the Supporting Information. Geometry optimizations were carried out in the solvent phase without any constraints. The optimized stationary points were characterized as local minima or transition structures by harmonic force constant analysis. In addition, intrinsic reaction coordinate (IRC) calculations were performed to verify the connection between reactant and product via transition structures.18 The lowest energy conformers of reactants, intermediates, and transition structures are reported in the text. Free energies of all stationary points were obtained from calculated thermal and entropic corrections at 298 K using unscaled vibrational frequencies. Relative free energies of all stationary points such as reactants, intermediates, and transition states were also computed at the following DFT levels: (1) L1: SMD(ethylethanoate)/M06-2X/6-311+G**//SMD(ethylethanoate)/M06-2X/6-311G** (2) L2: SMD(ethylethanoate)/B3LYP/6-311G**//SMD(ethylethanoate)/M06-2X/6-311G** (3) L3: SMD(ethylethanoate)/wB97XD/6-311G**//SMD(ethylethanoate)/M06-2X/6-311G** The LANL2DZ basis set for Br was used in calculations whenever it was applicable. The empirical dispersion corrections were included throughout the calculations. All these results are documented in the Supporting Information. Natural charges were computed using the natural bond order (NBO) analysis.19 The topological analysis of transition states is performed using Multiwfn program.20

as alkenyl boronic acids5 or alkenyl iodonium salts.6 Mechanistically, copper catalyst reacts with the alkenyl boronic acid or iodonium salt to generate alkenyl Cu(III) intermediate,7 whereas chiral amine reacts with the aldehyde and forms an enamine intermediate as shown in Scheme 2. The coupling of Scheme 2. Putative Intermediates in Enantioselective αAlkenylation of Aldehydesa

a

The catalysts investigated in this study are shown in the inset.

intermediates, i.e., enamine and alkenyl Cu(III) intermediate generated in two discrete catalytic cycles, is most likely to occur in the steric-free site of the enamine intermediate, avoiding hindrance provided by bulky groups on the chiral amine.8 In such cases, one can speculate that the addition of enamine to the alkenyl Cu(III) intermediate may experience somewhat similar local environment at the competing transition states, resulting into low enantioselectivity. On the contrary, very high enantioselectivity was reported for these reactions. Interestingly enough, changes in the achiral copper catalyst (Cu(OAc)2 vs CuBr) and the chirality of chiral amine catalyst in these reactions provided products with opposite selectivity.5,6 This experimental evidence indicates an active role for the chiral amine as well as achiral copper catalysts, apart from functioning as a HOMO or LUMO activator, in causing differentiation at the diastereomeric transition states. The precise role of copper and amine catalysts in the stereoselectivity-determining step of these reactions has not yet been delineated. In this paper, we present computational insights on the chiral amine−copper synergistic catalysis in the enantioselective α-alkenylation of aldehyde reactions. It is worth mentioning that MacMillan and his group have also employed the synergistic catalysis protocol with considerable success in new asymmetric bond-forming transformations such as the enantioselective α-trifluoromethylation,9 α-arylation,10 and α-oxidation11 of aldehyde. Hence, we anticipate that this study may help us to understand the role of catalysts in a wide range of asymmetric reactions involving chiral amine−transition-metal synergistic catalysis.





RESULTS AND DISCUSSION The α-alkenylation of aldehydes catalyzed by combination of chiral amine and Cu catalyst may proceed through the pathway shown in Scheme 2. This pathway involves four key steps: formation of enamine intermediate using organocatalysis, formation of Cu(III) intermediate using Cu catalysis, addition of enamine intermediate to Cu(III) intermediate leading to iminium ion intermediate, and subsequent hydrolysis to form a final product.21 In general, the step involving addition of enamine intermediate to the electrophile (C−C or C− heteroatom bond formation) is considered as the stereoselectivity-determining step of enamine-mediated reactions. However, Blackmond and co-workers in recent studies have proposed that the relative stability and reactivity of downstream diastereomeric intermediates formed in the catalytic cycle can dictate the enantioselectivity of chiral amine catalyzed reaction of aldehyde with electrophile rather than the diastereomeric transition states of the addition step in which a new stereogenic center is created by the attack of enamine to electrophile.22 In another interesting study, Pápai and co-workers have suggested that the stereoselectivity and rate of enamine mediated reaction can be determined in two different steps of the catalytic cycle.23 They have also shown that the addition of electrophile to the enamine intermediate (C−C bond forming step) controls the stereoselectivity of reaction, whereas subsequent protonation step determines the rate of enamine mediated reaction. It is worth noting that the catalytic system investigated in this study is completely different from the catalytic system that has been used by Blackmond et al. and Pápai et al. in their studies. Contrary to the diarylprolinol ether catalyst investigated in the

COMPUTATIONAL DETAILS

All electronic structure calculations were performed using the Gaussian 09 quantum chemical program.12 Density functional theory (DFT) was applied using the M06-2X functional13 in combination with the 6311G** basis set for C, H, O, N, and Cu14 and the LANL2DZ 1305

DOI: 10.1021/acs.joc.7b02822 J. Org. Chem. 2018, 83, 1304−1311

Article

The Journal of Organic Chemistry

Table 1. Relative Free Energiesa (kcal/mol) of Reaction Complexes (Rc) and Transition States (Ts) Computed at the M062X Level catalyst combination

Rc-(R)

Ts-(R)

ΔG⧧

Rc-(S)

Ts-(S)

ΔG⧧

% eeb

%eec

I: 5A + 6a II: 5B + 6a III: 5C + 6a IV: 5D + 6a V: 5E + 6a VI: 5E + 6b

−7.5 −5.8 −10.2 −8.1 −1.7 −2.3

−4.3 −3.2 −8.4 −6.5 0.8 1.5

3.2 2.6 1.8 1.6 2.5 3.8

−7.9 −5.8 −5.2 −7.2 −5.8 −7.1

−4.4 −2.6 −2.4 −3.9 −2.5 −5.0

3.5 3.2 2.8 3.3 3.3 2.1

8 (R) 47 (R) >99 (R) 98 (R) 99 (S) >99 (S)

d 93 (R) 86 (R) 93 (R) d 91 - 97e (S)

a The free energies computed relative to separated intermediates (5 and 6a or 6b). bEnantiomeric excess computed using difference in free energies of the diastereomeric transition states. cExperimentally observed enantiomeric excess. dExperimental enantiomeric excess are not available. eRange of enantioselctivity (expresess by %ee) reported for the reaction of aldehydes with different alkenyl iodonium salts (see ref 6).

recent mechanistic studies,22,23 the present study explores the role of imidazolidinone based chiral amine catalyst bearing two bulky substituents on the C-2 and C-5 position of the ring. Such catalytic framework can restrict the formation of downstream intermediates which were proposed by Blackmond et al. to explain the stereochemical outcome of α-chlorination of aldehyde-catalyzed pyrrolidine-based chiral amine catalyst.22 Moreover, the rate-limiting protonation step discussed by Pápai et al.23 is not relevant in the catalytic cycle proposed for αalkenylation of aldehydes reaction in Scheme 2. Considering all these facts, the step involving the addition of enamine intermediate (5) to the alkenyl Cu(III) intermediate (6) is expected to be the stereoselectivity-determining step of reaction. Hereafter, the stereoselectivity-determining step is referred as the “addition step”. The addition step is investigated using DFT methods. The five imidazolidinone-based catalysts with different α-substituents are employed as chiral amine catalysts in this study (Scheme 2). These catalysts are denoted by the letters A, B, C, and D which have −methyl (−CH3), −phenyl (−C6H5), −phenylmethyl (−CH2C6H5), and −naphthylen-1-ylmethyl (−CH2C12H11) groups at the C-5 of imidazolidinone ring, respectively. The tert-butyl substituent at the C-2 of imidazolidinone ring is common to these catalysts. On the other hand, the catalyst E has −phenyl (−C6H5) and tert-butyl groups on the C-2 and C-5 positions of the ring, respectively. The configuration of chiral centers (C-2 and C-5) of catalysts A, B, C, and D are (2S,5S), whereas catalyst E has the (2R,5R) configuration on chiral centers. Addition of enamine intermediates (5) formed by the reaction of chiral amines (A−E) with the propanal (1), and the alkenyl Cu(III) intermediate (6a) derived from the reaction of Cu(OAc)2 and boronic acid (2) were considered for the computational investigations (I−V, Table 1). The addition step with alkenyl Cu(III) intermediate (6b) derived from the reaction of CuBr and iodonium salt (3) is also investigated. The energetics of addition of 5E to 6b which is relevant to the experimental studies is reported in the text (VI, Table 1), and information about addition of other enamines (5A−D) to 6b is provided in the Supporting Information. To simplify our work, triflate ion of 6b was replaced by the acetate ion in calculations. Enamine intermediate (5) has two major isomers arising from the orientation of enamine double bond relative to the bulky substituents (tert-butyl and R1 substituents) on the imidazolidinone ring. For the sake of convenience, syn and anti notations (as shown in Scheme 3) are used to depict the orientation of the enamine double bond relative to the R1 group. A detailed conformational analysis of syn-(5) and anti(5) was performed to identify the conformational possibilities

Scheme 3. Two Faces of Enamine Intermediate Available for the Addition to Electrophile

pertaining to different orientations of R1-group. In general, anti-(5) is found to be slightly unfavorable as compared to syn(5) presumably because of the proximity of enamine double bond to the bulky tert-butyl group at the 2-position of the ring. Since bulky groups placed at the same side of the imidazolidinone ring can effectively hinder one face of enamine olefinic bond, the other face of enamine bond is relatively exposed for the addition step. As a result, for the enamine derived from catalysts A−D, the re-face of syn-(5) and si-face of anti-(5) is preferred for the C−C bond formation with the alkenyl Cu(III) intermediate (6a or 6b). The addition of syn(5) to 6a through the re face will lead to the (R)-enantiomer, whereas the si face addition will result in the (S)-enantiomer of 4. In the case of enamine derived from catalyst E, the si face of syn-(5E) and re-face of anti-(5E) will be relatively exposed in the addition step. Consequently, the addition of 6a (or 6b) to the si-face of syn-(5E) and re-face of anti-(5E) will lead to the (R) and (S)-enantiomers of 4, respectively. The corresponding transition states with different catalyst combinations were located at the M062X level. The free energies of reaction complexes (Rc) and transition states (Ts) for each catalyst computed relative to separated intermediates (5 and 6a or 6b) are provided in Table 1.24 The channel leading to (R)- and (S)enantiomers of 4 is designated by (R) and (S), respectively. Comparison of relative free energies of TSs involving different catalysts shows a distinct preference for the addition of syn-(5) to 6a. While the predicted product selectivity is in agreement with experimentally observed selectivity, calculated enantiomeric excess (ee) varies significantly with the catalyst combinations. Our calculations predict a very low enantioselectivity for the catalyst combination I(5A + 6a) whereas for II(5B + 6a), moderate enantioselectivity is estimated. The best results (high enantioselectivities) are obtained for the catalyst combinations III (5C + 6a), IV(5D + 6a), V(5E + 6a), and VI(5E + 6b). The representative catalyst combinations (IV and VI), for which predicted enantioselectivities are found to be in 1306

DOI: 10.1021/acs.joc.7b02822 J. Org. Chem. 2018, 83, 1304−1311

Article

The Journal of Organic Chemistry Table 2. Distortion/Interaction Analysis for the Addition Step with Different Catalyst Combinations (I−VI)a Ts-(R)

Ts-(S)

catalyst combination

ΔE⧧int

ΔE⧧dis(5)

ΔE⧧dis(6)

ΔE⧧dis

ΔE⧧

ΔE⧧int

ΔE⧧dis(5)

ΔE⧧dis(6)

ΔE⧧dis

ΔE⧧

I: 5A + 6a II: 5B + 6a III: 5C + 6a IV: 5D + 6a V: 5E + 6a VI: 5E + 6b

−69.3 −65.2 −74.4 −76.9 −76.3 −77.0

15.2 14.0 16.6 20.0 19.3 23.0

31.9 30.2 33.6 31.0 40.4 34.0

47.1 44.2 50.2 51.0 59.7 57.0

−22.2 −21.0 −24.2 −25.9 −16.6 −20.0

−69.5 −69.7 −67.3 −72.2 −65.4 −82.2

16.1 18.4 15.5 17.0 14.1 17.7

30.9 30.2 30.8 30.8 30.3 39.9

47.0 48.6 46.3 47.8 44.4 57.6

−22.5 −21.1 −21.0 −24.4 −21.0 −24.6

ΔE⧧int: interaction energy. ΔE⧧dis(5): distortion energy for enamine fragment. ΔE⧧dis(6): distortion Cu catalyst fragment. ΔE⧧dis: total distortion energy. ΔE⧧: activation energy. All energies are given in kcal/mol. a

Figure 1. Optimized transition states for the addition step with catalyst combinations IV and VI. See Scheme 3 for the atom numbering. The values in parentheses indicate geometric distortion of the enamine moiety in the transition state as compared to the parent enamine intermediate. Select bond distances are given in Å and dihedral angles are in degrees. Atom colors: H = white, C = silver, N = blue, O = red and Br = brick red.

corresponds to enantiomeric excess (ee) of 98% in favor of Rproduct, which reasonably agrees with the experimentally observed %ee (93%).25 On the other hand, for the catalyst combination VI, >99% ee in favor of S-product is estimated. The change in trend of enantioselectivity from (R)- to (S)product with a variation in the catalyst combination IV to VI is in agreement with experimental results.5,6 The vital question here is how these catalyst combinations are able to tune asymmetric environment for the addition of 5 to 6a (or 6b) in

close agreement with experimental enantioselectivities, discussed succinctly in the following section. In the case of catalyst combination IV, relative to separated intermediates, free energies of reaction complexes Rc-(R)-IV and Rc-(S)-IV are −8.1 and −7.2 kcal/mol, respectively. The transition states Ts-(R)-IV and Ts-(S)-IV are found to be very close to the reacting complexes with activation barriers of 1.6 and 3.3 kcal/mol, respectively, rendering an early transition state for the addition step. The difference in free energies 1307

DOI: 10.1021/acs.joc.7b02822 J. Org. Chem. 2018, 83, 1304−1311

Article

The Journal of Organic Chemistry

dihedral angle ω3: N−C6−C7−C(Me) which exhibits distortion at the C7 during C7−C8 bond formation as shown in Figure 1. In Ts-(R)-IV and Ts-(S)-IV, the ω1 and ω2 show a modest deviation while, ω3 shows a considerable deviation compared to the corresponding angles of the enamine intermediate (5) geometry. The deviation in these dihedral angles is more in Ts-(R)-IV than in Ts-(S)-IV which perfectly correlates with the trend of distortion energies predicted for these transition states. Similar to the distortion energies, the interaction energies between two catalytic fragments at the transition state also vary with the different catalyst combinations. It is also clear from Table 2 that the stabilization offered by the interaction between enamine fragment and Cu catalyst fragment is much higher than the destabilization incurred due to distortion of these fragments. In order to identify the source of high interaction energies, we have investigated catalyst-free model reactions involving α-alkenylation of aldehydes (Scheme 4).29 In the first

diastereomeric transition states even though, such additions prefer a less hindered site of the chiral amine catalyst. We next sought to identify the steric as well as electronic factors which govern the stereoselectivity of the reaction by examining lowest energy transition states of the addition step. In general, various factors play crucial role in determining the stereochemical outcome of amino catalysis.26 These factors include (i) steric interactions offered by α-substituents of chiral amine catalyst, (ii) hydrogen-bond interactions between catalyst and substrate, and (iii) geometric distortion of catalyst and substrates at the transition states. In the current reaction, bulky α-substituents of the chiral amine catalysts, which are positioned at the same side of the imidazolidinone ring, provide effective shielding to one of the faces of enamine double bond. Such disposition of bulky groups on the catalyst reduces the number of possible approaches of alkenyl Cu(III) intermediate to the enamine 5 and thus limits the effective interactions of reacting partners. Furthermore, the interaction and the distortion in catalytic fragment namely, enamine fragment and Cu catalyst fragment, can also contribute to the energy differentiation in the diastereomeric transition states. The nature and energy costs associated with the distortion of catalytic fragments at the transition states were determined using the distortion/interaction analysis.27 The distortion energies Edis associated with distortion of enamine fragment Edis(5) and Cu catalyst fragment (Edis(6)) at the transition state compared to the ground state geometries of 5 and 6a, respectively, were calculated. This analysis also provides interaction energies Eint between enamine fragment and Cu catalyst fragment of the transition states. The energy components for each catalyst combinations (I−VI) and corresponding values are summarized in Table 2. The difference of ΔEact of diastereomeric transition states in Table 2 show a reasonable correlation with the trend of selectivity (Table 1).28 It is apparent from the data shown in the Table 2 that the distortion energies of alkenyl Cu fragment are higher than the distortion energies of enamine fragment. However, the distortion energies of alkenyl Cu fragment remain qualitatively same in the transition states, except for the catalyst combination VI which is discussed separately (vide infra). Thus, the difference in the total distortion energies of diastereomeric transition states mainly arises from the differential distortion of enamine fragment at the transition states. The optimized geometries of TSs of representative catalyst combination IV (Ts-(R)-IV and Ts-(S)-IV) were considered for a detailed analysis (Figure 1). It is noticed that the coordination of acetate ions to Cu at these transition states shows major deviation when it is compared with the ground state geometry of alkenyl Cu(III) intermediate (6a). The distortion in the coordination of acetate ions to Cu at the transition states presumably occurs to minimize the steric clash with the imdazolidinone ring. However, distortion energies of Cu catalytic fragment in Ts-(R)-IV and Ts-(S)-IV are comparable. Next, we compare geometries of Ts-(R)-IV and Ts-(S)-IV with the geometry of enamine intermediate (5). The C−C bond formation and enamine-iminium ion conversion at the transition state is accompanied by the change in conformation of imidazolidinone ring. To account the distortion in enamine fragments, dihedral angles ω1: C2−N− C6−C7 and ω2: C5−N−C6−C7 which describe changes in the conformation of imidazolidinone ring relative to enamine double bond C6−C7 were analyzed. We have also noted the

Scheme 4. Transition State of Catalyst-Free Model Reactionsa

a

Select bond distances are shown in Å. Similar to catalyzed reactions, the carbon atoms involved in the C−C bond formation are designated as C7 and C8.

model reaction (M1), addition of enolate of propanal to 3phenylpropene is considered. The free energy of activation for this reaction is found to be as high as 27.7 kcal/mol, while the energy of interaction between two reactants at the transition state is 11.7 kcal/mol. In another model reaction, we have considered addition of enol form of propanal to the alkenyl Cu intermediate. In this case, free energy of activation and interaction energy at the transition state are 14.0 and 33.8 kcal/mol, respectively. These results clearly demonstrate that high interaction energies stemming from the augmented interaction between catalytic fragments can be attributed to the catalytic (HOMO and LUMO) activation of reactants achieved in the synergistic catalysis. Such dual catalytic activation of reactants, ultimately, helps to reduce the activation barrier of addition step. It is worthwhile to note that free energies of activation predicted for the addition step (C−C bond formations) in present work are significantly lower as compared to the free energies of activation reported for the C− C bond-forming step in various enamine mediated reactions.30 However, the effect of HOMO-raising and LUMO-lowering activation achieved through synergistic catalysis is expected to remain same in the diastereomeric transition states, and likely to contribute equally to the interaction energies. Interestingly, interaction energies computed for the diastereomeric transition states of various catalytic combinations are not in the same range. For instance, the interaction energy in Ts-(R)-IV is 4 kcal/mol higher than that of Ts-(S)-IV which suggests, that the 1308

DOI: 10.1021/acs.joc.7b02822 J. Org. Chem. 2018, 83, 1304−1311

Article

The Journal of Organic Chemistry

computed for the transition states of various catalyst combinations (Table 2) exhibit interesting variations. A plot of distortion vs interaction energies of transition states of catalyst combinations (I−VI) display a reasonable correlation with R2 = 0.782 (Figure 2a).31 This trend in the distortion and

stabilizing noncovalent interactions between two catalytic fragments may also contribute to the net interaction energy. The comparison of transition-state geometries reveals the possibility of two stabilizing noncovalent interactions at the transition state. The first one includes the interaction between Cu and the enamine moiety (referred as metal−enamine interaction). In Ts-(R)-IV, Cu is located at 2.12 and 2.48 Å from C7 and C6, respectively (Figure 1). Natural population analysis of the transitions state shows Cu bears a partial positive charge (1.24 e), while C7 and C6 of enamine bond have −0.44 e and 0.21 e, respectively. This charge distribution suggests that Cu may interact with the C7 of enamine moiety which is also involved in the C−C bond formation process at the transition state. The corresponding frontier molecular orbital showing Cu···C7 orbital interaction is shown in Figure 1. Along with the Cu···C7 interaction, a hydrogen-bonding interaction between the oxygen of acetate and hydrogen on the C2 or C5 of imidazolidinone ring is also identified in the transition state. The (acetate) O···H distance at the Ts-(R)-IV is 2.14 Å. These two interactions Cu···C7 and (acetate) O···H(Cα) are relatively weak in Ts-(S)-IV, which is evident from the slightly longer Cu···C7 and (acetate) O···H distances in Ts-(S)-IV (2.14 and 2.24 Å, respectively). We have also analyzed topological features of electron density distribution at the transition state to quantify the extent of non covalent interactions using multiwfn program. The electron densities of bond critical points (BCP) for the Cu···C7 and (acetate) O···H interactions in Ts-(R)-IV are 0.073 and 0.016 au, respectively. On the other hand, in Ts-(S)-IV, the electron density values of BCP for these interactions are found to be lower by 0.004−0.005 au as compared to Ts-(R)-IV. It is clear from topological analysis of transition states that Cu···C7 interaction is more prominent than the (acetate) O···H (Cα) hydrogen-bonding interaction and could be a major contributing factor that leads to the differences in interaction energies at the diastereomeric transition states. To obtain further support to our hypothesis, we have also examined transition-state geometries (Figure 1) of the catalyst combination VI in which (S)-product selectivity is predicted. In this catalyst combination, one of the acetate ions of 6a is replaced by bromide ion to model the alkenyl Cu(III) intermediate (6b) derived from the reaction of iodonium salt and CuBr. Even with this catalyst combination, the Cu···C7 distance in the lower energy Ts-(S)VI is shorter by 0.03 Å compared to Ts-(R)-VI (2.10 vs 2.13 Å). Accordingly, interaction energy of the catalytic fragments in the preferred Ts-(S)-VI is found to be 5 kcal/mol higher than that of Ts-(R)-VI. The distortion energies of Cu catalyst fragment, however, follow a different trend in comparison to other catalyst combinations. The distortion energies of Cu catalyst fragment in Ts-(R)-VI and Ts-(S)-VI, and net distortion energies of these transition states are higher than that of other transition states. This is mainly because of the large displacement of bromide ion in the transition states relative to the original position in the ground-state geometry of 6b. Interestingly, the net distortion energies in Ts-(R)-VI and Ts-(S)-VI qualitatively remained the same, and therefore, we reasoned that the difference in interaction energies which primarily arises from the Cu···C7 interaction might be responsible for the energy difference in these transition states. After having analyzed stabilizing and destabilizing factors at the transition states, we further discussed the salient features of computational model of synergistic catalysis presented in this study. It is noticed that the interaction and distortion energies

Figure 2. Correlation plot of (a) distortion energy vs interaction energy (b) Cu···C7 bond distance vs interaction energy. Energies are given in kcal/mol and bond distances in Å.

interaction energies conveys that catalytic fragments may experience a sizable distortion to maximize the interaction between these fragments at the transition states. However, sufficiently strong stabilizing interaction between two catalytic fragments undertakes destabilization caused by the distortion in catalytic fragments, resulting into low activation energy for the transition state. Our computational model also suggests that a subtle change in the metal−enamine interaction described through the Cu···C7 distance may determine the variation in the interaction energies of diastereomeric transition states. Consistent with this hypothesis, a linear relationship between the Cu···C7 distance (at the transition state) and interaction energies computed for the transition states of catalyst combinations (I−VI) is obtained (Figure 2b).32 The plot shows that the deviation of 0.06 Å in Cu···C7 distance may cause approximately 15 kcal/mol variation in the interaction energy. Such a large variation in the interaction energy is possible only if the metal-enamine interaction operate in tandem with catalytic (covalent) activation. Careful examination of the metal-enamine interaction allows envisaging a threecentered transition state for the addition step wherein, the metal activates both reacting centers, C7(nucleophilic) and C8(electrophilic). Thus, the enamine fragment may experience additional nucleophilic activation through metal−enamine interaction. The cooperative nature of these interactions may facilitate C7···C8 bond formation, rendering additional stabilization to the transition state. A further support for the 1309

DOI: 10.1021/acs.joc.7b02822 J. Org. Chem. 2018, 83, 1304−1311

The Journal of Organic Chemistry



coaction of catalytic activations and metal-enamine interaction comes from the comparison of Cu···C7 distance in the transition state of catalyzed vs catalyst-free model (M2) reaction. In catalyst-free transition state of model reaction (M2), the Cu···C7 distance is increased by about 0.14−0.19 Å relative to that in the transition state of catalyzed reactions (Scheme 4). The weakening of Cu···C7 interaction without an enamine activation of aldehyde manifests the cooperation between the catalytic activations and metal-enamine interaction at the transition state. Besides, a topology analysis of Cu···C7 interaction, which recognized it as a relatively strong noncovalent interaction than the hydrogen bonding interaction, justifies the large fluctuation in interaction energies with a delicate change in the Cu···C7 distance. In summary, the differential metal−enamine interaction at the diastereomeric transition states within a synergistic catalytic framework primarily dictates the enantioselectivity of reaction.

Article

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.joc.7b02822. Optimized Cartesian coordinates and total energies of all computed stationary points (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Mahendra Patil: 0000-0001-9631-4407 Notes

The author declares no competing financial interest.





ACKNOWLEDGMENTS I thank Department of Science and Technology (DST), India (SB/FT/CS-086/2012), for the partial funding.

CONCLUSIONS Computational investigations of the stereoselectivity-determining step of the enantioselective α-alkenylation of aldehyde reveal interesting aspects of the synergistic catalysis in the asymmetric reaction. We have examined different combinations of chiral amines (A−E) and Cu catalysts in the computational model to provide a comprehensive view on the origin of the high enantioselectivity in the reaction. In general, bulky substituents on the chiral amine, as demonstrated in the various activation modes of aminocatalysis such as iminium, enamine and SOMO activation catalysis, play a crucial role in the asymmetric induction, providing a differential steric environment for the addition of nucleophile/electrophile.33 However, in the present reaction, the addition of the enamine intermediate derived from the chiral amines (A−E) to alkenyl Cu(III) intermediate (electrophile) prefers a face opposite to bulky groups and thus, minimizes the contribution of steric factor in the chiral induction. In the absence of steric discrimination at the competing transition states, our calculations show that the interaction between catalytic fragments (enamine and alkenyl Cu catalyst fragments) in the synergistic catalytic framework can be decisive for the stereoinduction in reaction. The high interaction energies, in the range of 65 to 80 kcal/mol, are predicted using distortion/ interaction analysis. Furthermore, at the transition state, the Cu is seen to activate the enamine fragment through the Cu···C7 interaction. The typical Cu···C7 distance at the transition state is found to be in the range of 2.10 to 2.16 Å. This interaction allows Cu catalyst to engage in nucleophilic (at C7) and electrophilic (at C8) activation simultaneously, via a threecentered transition state. Owing to cooperative nature of these interactions, a subtle change in the Cu···C7 interaction leads to a substantial variation in the interaction energies of diastereomeric transition states. Taken together, our computational results demonstrate that the basis for enantioselectivity in this reaction may be ascribed to the different degrees of interaction among catalytic fragments at the competing transition states, and the interaction energies are mainly determined by the cooperation of catalytic (covalent) activations and the metal−enamine interaction at the transition state. We speculate that the computational model presented in this paper may provide a guideline to a rational design of catalyst combinations for other asymmetric reactions that rely on the synergistic action of organo- and metal catalysis.



REFERENCES

(1) (a) Afewerki, S.; Córdova, A. Chem. Rev. 2016, 116, 13512. (b) Deng, Y.; Kumar, S.; Wang, H. Chem. Commun. 2014, 50, 4272. (c) Hopkinson, M. N.; Sahoo, B.; Li, J.-L.; Glorius, F. Chem. - Eur. J. 2014, 20, 3874. (d) Denard, C. A.; Hartwig, J. F.; Zhao, H. ACS Catal. 2013, 3, 2856. (e) Du, Z.; Shao, Z. Chem. Soc. Rev. 2013, 42, 1337. (f) Allen, A. E.; MacMillan, D. W. C. Chem. Sci. 2012, 3, 633. (g) Rueping, M.; Koenigs, R. M.; Atodiresei, I. Chem. - Eur. J. 2010, 16, 9350. (h) Zhong, C.; Shi, X. Eur. J. Org. Chem. 2010, 2010, 2999. (i) Shao, Z.; Zhang, H. Chem. Soc. Rev. 2009, 38, 2745. (2) (a) Meazza, M.; Rios, R. Synthesis 2016, 48, 960. (b) Wang, M. H.; Scheidt, K. A. Angew. Chem., Int. Ed. 2016, 55, 14912. (c) Inamdar, S. M.; Shinde, V. K.; Patil, N. T. Org. Biomol. Chem. 2015, 13, 8116. (d) Deng, Y.; Kumar, S.; Wang, H. Chem. Commun. 2014, 50, 4272. (e) Cohen, D. T.; Scheidt, K. A. Chem. Sci. 2012, 3, 53. (f) Ramachary, D. B.; Jain, S. Org. Biomol. Chem. 2011, 9, 1277. (g) Ambrosini, L. M.; Lambert, T. H. ChemCatChem 2010, 2, 1373. (h) Wasilke, J.-C.; Obrey, S. J.; Baker, R. T.; Bazan, C. G. Chem. Rev. 2005, 105, 1001. (3) (a) Lam, Y.-h.; Grayson, M. N.; Holland, M. C.; Simon, A.; Houk, K. N. Acc. Chem. Res. 2016, 49, 750. (b) Jindal, G.; Kisan, H. K.; Sunoj, R. B. ACS Catal. 2015, 5, 480. (4) (a) Li, Y.; Ibsen, L.; Jørgensen, K. A. Org. Lett. 2017, 19, 1200. (b) Guo, J.; Lin, L.; Liu, Y.; Liu, X.; Feng, X. Org. Lett. 2016, 18, 5540. (c) Xiao, J. Org. Lett. 2012, 14, 1716. (d) Lou, S.; Fu, G. C. J. Am. Chem. Soc. 2010, 132, 5010. (e) Taylor, A. M.; Altman, R. A.; Buchwald, S. L. J. Am. Chem. Soc. 2009, 131, 9900. (f) Dai, X.; Strotman, N. A.; Fu, G. C. J. Am. Chem. Soc. 2008, 130, 3302. (g) Kim, H.; MacMillan, D. W. C. J. Am. Chem. Soc. 2008, 130, 398. (h) Poulsen, T. B.; Bernardi, L.; Bell, M.; Jørgensen, K. A. Angew. Chem., Int. Ed. 2006, 45, 6551. (i) Chieffi, A.; Kamikawa, K.; Ahman, J.; Fox, J. M.; Buchwald, S. L. Org. Lett. 2001, 3, 1897. (5) Stevens, J. M.; MacMillan, D. W. C. J. Am. Chem. Soc. 2013, 135, 11756. (6) Skucas, E.; MacMillan, D. W. C. J. Am. Chem. Soc. 2012, 134, 9090. (7) (a) King, A. E.; Brunold, T. C.; Stahl, S. S. J. Am. Chem. Soc. 2009, 131, 5044. (b) King, A. E.; Ryland, B. L.; Brunold, T. C.; Stahl, S. S. Organometallics 2012, 31, 7948. (8) (a) Um, J. M.; Gutierrez, O.; Schoenebeck, F.; Houk, K. N.; MacMillan, D.W. C. J. Am. Chem. Soc. 2010, 132, 6001. (b) Gordillo, R.; Houk, K. N. J. Am. Chem. Soc. 2006, 128, 3543. (9) Allen, A. E.; MacMillan, D. W. C. J. Am. Chem. Soc. 2010, 132, 4986. (10) Allen, A. E.; MacMillan, D. W. C. J. Am. Chem. Soc. 2011, 133, 4260. 1310

DOI: 10.1021/acs.joc.7b02822 J. Org. Chem. 2018, 83, 1304−1311

Article

The Journal of Organic Chemistry

2014, 43, 4953. (d) Ess, D. H.; Houk, K. N. J. Am. Chem. Soc. 2008, 130, 10187. (28) The distortion/interaction analysis has been frequently used to elucidate the origin of stereoselectivity in asymmetric reactions. See: (a) Duarte, F.; Paton, R. S. J. Am. Chem. Soc. 2017, 139, 8886. (b) Champagne, P. A.; Houk, K. N. J. Am. Chem. Soc. 2016, 138, 12356. (c) Ess, D. H.; Houk, K. N. J. Am. Chem. Soc. 2008, 130, 10187. (d) Wolf, L. M.; Thiel, W. J. Org. Chem. 2014, 79, 12136. (e) Green, A. G.; Liu, P.; Merlic, C. A.; Houk, K. N. J. Am. Chem. Soc. 2014, 136, 4575. (29) Optimized geometries of transition states of model reactions are shown in the Supporting Information. See Figure S1. (30) (a) Simon, A.; Lam, Y.-h.; Houk, K. N. J. Org. Chem. 2017, 82, 8186. (b) Dieckmann, A.; Breugst, M.; Houk, K. N. J. Am. Chem. Soc. 2013, 135, 3237. (c) Duarte, F. J. S.; Santos, A. G. J. Org. Chem. 2012, 77, 3252. (d) Bertelsen, S.; Marigo, M.; Brandes, S.; Dinér, P.; Jørgensen, K. A. J. Am. Chem. Soc. 2006, 128, 12973. (e) Clemente, F. R.; Houk, K. N. Angew. Chem., Int. Ed. 2004, 43, 5766. (f) Bahmanyar, S.; Houk, K. N. J. Am. Chem. Soc. 2001, 123, 11273. (31) A plot of distortion vs interaction energies of transition states of all catalyst combinations (I−X) is shown in the Supporting Information. See Figure S2. (32) A plot of the Cu···C7 distance (at the transition state) vs interaction energies of transition states of all catalyst combinations (I− X) is shown in Supporting Information. See Figure S2. (33) (a) Bertelsen, S.; Jørgensen, K. A. Chem. Soc. Rev. 2009, 38, 2178. (b) MacMillan, D. W. C. Nature 2008, 455, 304. (c) Gaunt, M. J.; Johansson, C. C. C.; McNally, A.; Vo, N. T. Drug Discovery Today 2007, 12, 8. (d) Dalko, P. I.; Moisan, L. Angew. Chem., Int. Ed. 2004, 43, 5138. (e) Dalko, P. I.; Moisan, L. Angew. Chem., Int. Ed. 2001, 40, 3726.

(11) Simonovich, S. P.; Van Humbeck, J. F.; MacMillan, D. W. C. Chem.Sci. 2012, 3, 58. (12) Gaussian 09, Revision D.01, 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. A., Jr.; 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, Inc., Wallingford, CT, 2013. (13) Zhao, Y.; Truhlar, D. G. Theor. Chem. Acc. 2008, 120, 215. (14) (a) McLean, A. D.; Chandler, G. S. J. Chem. Phys. 1980, 72, 5639. (b) Krishnan, R.; Binkley, J. S.; Seeger, R.; Pople, J. A. J. Chem. Phys. 1980, 72, 650. (c) Blaudeau, J.-P.; McGrath, M. P.; Curtiss, L. A.; Radom, L. J. Chem. Phys. 1997, 107, 5016. (d) Wachters, A. J. H. J. Chem. Phys. 1970, 52, 1033. (e) Hay, P. J. J. Chem. Phys. 1977, 66, 4377. (f) Raghavachari, K.; Trucks, G. W. J. Chem. Phys. 1989, 91, 1062. (g) Binning, R. C., Jr.; Curtiss, L. A. J. Comput. Chem. 1990, 11, 1206. (h) McGrath, M. P.; Radom, L. J. Chem. Phys. 1991, 94, 511. (i) Curtiss, L. A.; McGrath, M. P.; Blaudeau, J.-P.; Davis, N. E.; Binning, R. C., Jr.; Radom, L. J. Chem. Phys. 1995, 103, 6104. (15) (a) Hay, P. J.; Wadt, W. R. J. Chem. Phys. 1985, 82, 270. (b) Wadt, W. R.; Hay, P. J. J. Chem. Phys. 1985, 82, 284. (c) Hay, P. J.; Wadt, W. R. J. Chem. Phys. 1985, 82, 299. (16) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. J. Phys. Chem. B 2009, 113, 6378. (17) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. J. Chem. Phys. 2010, 132, 154104. (18) (a) Gonzalez, C.; Schlegel, H. B. J. Chem. Phys. 1989, 90, 2154. (b) Gonzalez, C.; Schlegel, H. B. J. Phys. Chem. 1990, 94, 5523. (19) Glendening, E. D.; Reed, A. E.; Carpenter, J. E.; Weinhold, F. NBO Version 3.1. (20) Lu, T.; Chen, F. J. Comput. Chem. 2012, 33, 580. (21) The steps involved in the formation of enamine intermediate and alkenyl Cu(III) intermediate are also investigated. The relevant information is provided in the Supporting Information. See Figure S3 and Scheme S1. (22) (a) Burés, J.; Armstrong, A.; Blackmond, D. G. J. Am. Chem. Soc. 2011, 133, 8822. (b) Burés, J.; Armstrong, A.; Blackmond, D. G. J. Am. Chem. Soc. 2012, 134, 6741. (23) Földes, T.; Madarász, Á .; Révész, Á ; Dobi, Z.; Varga, S.; Hamza, A.; Nagy, P. R.; Pihko, P. M.; Pápai, I. J. Am. Chem. Soc. 2017, 139, 17052. (24) All possible conformers of syn-(5) as well as anti-(5) were considered to identify the lowest energy diastereomeric transition states of the addition step. (25) The activation barriers for diastereomeric transition states computed using different DFT functionals are provided in the Supporting Information; see Tables S2−S4. (26) (a) Halskov, K. S.; Donslund, B. S.; Paz, B. M.; Jørgensen, K. A. Acc. Chem. Res. 2016, 49, 974. (b) Walden, D. M.; Ogba, O. M.; Johnston, R. C.; Cheong, P. H.-Y. Acc. Chem. Res. 2016, 49, 1279. (c) Sunoj, R. Acc. Chem. Res. 2016, 49, 1019. (d) Wheeler, S. E.; Seguin, T. J.; Guan, Y.; Doney, A. C. Acc. Chem. Res. 2016, 49, 1061. (e) Cheong, P. H.-Y.; Legault, C. Y.; Um, J. M.; Nihan, Ç .-Ö .; Houk, K. N. Chem. Rev. 2011, 111, 5042. (27) (a) Bickelhaupt, F. M.; Houk, K. N. Angew. Chem., Int. Ed. 2017, 56, 10070. (b) Wolters, L. P.; Bickelhaupt, F. M. WIREs Comput. Mol. Sci. 2015, 5, 324. (c) Fernández, I.; Bickelhaupt, F. M. Chem. Soc. Rev. 1311

DOI: 10.1021/acs.joc.7b02822 J. Org. Chem. 2018, 83, 1304−1311