Origin of Stereoselectivity of the Photoinduced Asymmetric Phase

It is well-known that the theoretical chemistry can offer alternative opportunity to study the mechanism,(19) especially the untouchable transition st...
0 downloads 0 Views 3MB Size
Article pubs.acs.org/joc

Origin of Stereoselectivity of the Photoinduced Asymmetric PhaseTransfer-Catalyzed Perfluoroalkylation of β‑Ketoesters Chen Yang,†,‡ Wei Zhang,† Yi-He Li,‡ Xiao-Song Xue,† Xin Li,*,† and Jin-Pei Cheng†,‡ †

State Key Laboratory of Elemento-organic Chemistry, College of Chemistry, Nankai University, Collaborative Innovation Center of Chemical Science and Engineering, Tianjin 300071, P.R. China ‡ Department of Chemistry and Biology, College of Science, National University of Defense Technology, Changsha 410073, China S Supporting Information *

ABSTRACT: The sources of asymmetric induction in the photoinduced phase-transfer-catalytic perfluoroalkylation of βketoesters were interpreted by density functional theory calculations. The calculations indicated that multiple hydrogen-bonding interactions mode rather than π−π stacking interactions model better described the transition structure of the reaction. The quantitative estimation of these noncovalent interactions within the key transition states identified that the O−H···O and specific C−H···O hydrogen bonds play key roles for the stereocontrol.



INTRODUCTION Though it is in infancy, the asymmetric catalysis driven by visible light1 has exhibited thriving vitality past a decade.2 Among the known enantioselective photocatalysis, the key chiral intermediates are generated mostly through covalent catalysis with photogenerated active species, such as chiral enamines,3 iminium ions,4 Breslow intermediates,5 metal complexes,6 and enolates.7 Recently, more challenging noncovalent strategy to control stereoselectivity in visible-light photochemistry has emerged. 8−16 As a result, several asymmetric transformations were achieved, such as αhydroxylation,8 aza-Pinacol cyclization,9 intramolecular11 and intermolecular12 cyclizations, oxidative C−H functionalization,13 redox neutral C−C bond forming reaction,14 and antiMarkovnikov hydroetherification.15 Besides the excellent examples mentioned above, recently, Melchiorre and co-workers exploited an attractive externalphotocatalyst-free strategy to achieve the phase-transfercatalyzed (PTC) enantioselective perfluoroalkylation of βketoesters,16 which relied on the formation of photoactive electron donor−acceptor (EDA) complexes17 between the chiral enolates and electron-deficient perfluoroalkyl iodides (Scheme 1). On basis of conducted diagnostic experiments on mechanism, they proposed a radical chain propagation pathway to interpret the asymmetric photochemical perfluoroalkylation. The transformation was triggered by irradiation of the chiral EDA complex, leading to the formation of the perfluoroalkyl radical RF· via the reductive cleavage of the perfluoroalkyl iodide RfI. The selective attack of the chiral enolate by the electrophilic RF· was believed to be the stereocontrol step. The resulting ketyl intermediate reacted with another RfI, yielding the product, the PTC catalyst and another RF·. The PTC catalyst and RF· entered into cycles to continuously provide products. © 2017 American Chemical Society

Although the PTC approaches, especially the cinchonaderived PTC catalysts, have found lots of applications on asymmetric transformations,18 we’re astonished that a full picture of the origin of stereoinduction remains obscure by screening literatures. It is well-known that the theoretical chemistry can offer alternative opportunity to study the mechanism,19 especially the untouchable transition states, besides experimental procedures. Recently, overcoming the large and flexible disadvantage of cinchona-derived20 PTC catalysts, some calculations were performed to throw new insights on PTC reactions.21−27 Lipkowitz and O’Donnell conducted a precedent theoretical study on the interaction between the N-benzylcinchoninium cation and Schiff basederived enolates with molecular dynamics and semiempirical quantum theory in 1991.21 They found the Z-enolates selective binding on the front side of the N-benzylcinchoninium may account for the observed enantioselectivity. The Palomo group proposed a dual activation model to fit the observed enantioselectivity in the N-benzylquininium ion catalyzed azaHenry reaction, where the catalyst 9-OH group interacts with the NO2 group of the nitro-compound and the catalyst +N− CH site interacts with the imine.22 Andrus and Ess disclosed the oxy-anion-ammonium interactions in transition states rather than π-face interactions in the ground states are responsible for the asymmetric allylation of Schiff base catalyzed by a 9-allyated cinchona derived PTC catalyst.23 Combination of density functional theory calculations and kinetic model, Pliego et al. pointed out both the hydrogen bond and electrostatic interaction involving the leaving chloride anion in the favored transition state are the key factors for the enantioselective inducement in the cinchoninium ion catalyzed alkylation of Received: May 9, 2017 Published: August 22, 2017 9321

DOI: 10.1021/acs.joc.7b01130 J. Org. Chem. 2017, 82, 9321−9327

Article

The Journal of Organic Chemistry Scheme 1. Photoinduced Enantioselectvie Pefluoroalkylation of β-Ketoesters

indanone.24a,b González and co-workers believed the hydrogenbond interaction dominated the favorable transition state but the π−π stacking interaction between substrate and the quinoline moiety of the catalyst destabilized the unfavorable transition state for it constrained the relative orientations of the substrates and catalyst, which led to the enantiopure epoxidation products of nitroalkenes.25 Capobianco and Palombi stressed the strength difference of hydrogen bond network established between 9-OH and 6′-OH of the quininium ion and the intermediate facilitated the desired product in an asymmetric tandem hemiaminal-heterocyclization-aza-Mannich reaction of 2-formylbenzonitriles.26 These theoretically enriched the understandings of the asymmetric induction mechanism of cinchona-derived phasetransfer catalysts, however, to our knowledge, phase-transfer catalytic examples involved in photochemical radical progress have not been reported by far. The researches on mechanisms and origins of stereoselectivity28,29 will undoubtedly accelerate the rational development of photochemistry. Herein, we reported the mechanism calculations on the photo-organocatalytic perfluoroalkylation of β-ketoesters.



Their relative free energies were calculated at M06-2X-D3/ def2TZVPP-SMD(PhCl)// M06-2X/6-31G(d)-SMD(PhCl) level.



RESULTS AND DISCUSSION Since the Melchiorre group proposed a radical chain propagation of the perfluoroalkylation reaction initiated by EDA complexes,16 we first evaluated the existence of the EDA complex by calculations. Simulating the catalyst with N-methyl quinuclidinium ion, eight possible ground EDA complexes were located (Figures 1 and S1). Through a TD-DFT calculation40 at M06-2X/6-31G(d)+Lanl2dz(I) level, we found that the EDA complexes through a C−I bond paralleled the plane of the enolate interaction (EDA1−2 and EDAS1−2) rather than halogen-bonding interaction EDA3−4 and EDAS3−4, between the CF3I and the enolate, are able to absorb visible light (wavelength in green). These results agreed well with the observed optical absorption spectra.16 The maximum absorption wavelength calculated at TD-CAM-B3LYP/6-311++G(d,p)+SDD(I) also gave the similar results (wavelength in blue). The calculated structures of EDA1−2 are very similar to the EDA complex suggested by Melchiorre and Maseras in their previous computational study.28d Meanwhile, the halogen-bonding complex seemed to be irradiated just by ultraviolet light, which is in line with the Chen’s recent study.41 The orbital analysis of EDA2 complex indicated the charge-transfer transition from HOMO (in red and green) to LUMO (in purple and yellow) should be responsible for the visible light absorption (λ = 434.7 nm). We also tried to obtain the structure of the singlet state of EDA complex, however, it failed because of the difficulty to keep the excited state in the same state during the optimization. Next, we came to consider the transition states (TSs) in the stereocontrolling step. Balancing the computational cost and precision, we adopted N-3,4,5-trifluorobenzyl cinchoninium ion as the catalyst and the reaction of methyl 1-oxo-2,3-dihydro1H-indene-2-carboxylate and CF3I as the model (Scheme 2). The existence of torsional freedom enabled cinchoninium cation to possess many conformations. Recently, the study from Pliego et al. pointed out that the N-4-trifluoromethylphenyl cinchoninium cation adopted the anti-open conformation in the alkylation reaction of indanone after they examined many

COMPUTATIONAL METHODS

All calculations were performed with Gaussian09 packages at T = 298.15 K and p = 1 atm.30 Unless noted, all energetics are reported in kcal/mol, and the bond lengths are reported in angstroms (Å). Structures were generated using CYLview,31 Multiwfn,32 and VMD.33 Geometry optimizations were conducted using M06-2X functional34 with 6-31G(d) basis set for C, H, O, N, and F atoms, and Lanl2dz35 for I atom at gas phase. The frequency calculations were conducted at the same level of theory to confirm the nature of stationary points and obtain the thermal corrections. The vertical excitations of EDA complexes were calculated at TD-M06-2X/6-31G(d)+Lanl2dz(I) and TD-CAM-B3LYP/6-311++G(d,p)+SDD(I) level without vibronic corrections.36 Single-point energies of the transition states were calculated at the M06-2X-D3/6-311++G(d,p) and M06-2X-D3/ def2TZVPP level, respectively.37,38 The solvent effect of chlorobenzene was evaluated with the SMD method.39 The thermal corrections derived from the vibrational frequency calculations were then added to the electronic energies calculated at M06-2X-D3/6-311++G(d,p) or M06-2X-D3/def2TZVPP level to obtain Gibbs free energies. Several transition states were also reoptimized in solution phase at M06-2X/6-31G(d) level based on the structures obtained in gas phase. 9322

DOI: 10.1021/acs.joc.7b01130 J. Org. Chem. 2017, 82, 9321−9327

Article

The Journal of Organic Chemistry

calculations. Other conformations of the PTC catalyst, including different conformations of the catalyst’s 3,4,5trifluorobenzyl and vinyl group, however, were not considered in the following calculations. As described in Scheme 2, three possible transition-state models were proposed with respect to the stereocontrolling step. In model a, the TSs were stabilized by hydrogen bond interactions established between 9-OH and +N−CH of the quininium ion and the enolate, and π−π stacking interaction between enolate and the phenyl group of the catalyst. In the model b, the TSs were stabilized by hydrogen bond interactions established between 9-OH and +N−CH of the quininium ion and the enolate, and π−π stacking interaction between enolate and the quinoline moiety of the catalyst. In model c, the TSs were stabilized only by hydrogen bond network between the two oxygens of the β-ketoester and the catalyst. According to these TS models, we calculated corresponding transition states (Figure 2). The relative free energies were obtained at the M06-2X-D3/6-311++G(d,p)-SMD(PhCl)// M06-2X/6-31G(d) level (values in green). The ts1−2 belonged to model a, the ts3−4 belonged to model b, and the ts5−8 belonged to model c. Given the product is S enantiomer, ts1−4 failed to predict the right isomer. The proposed π−π stacking interaction8a,b seemed to have little influence on the enantioselectivity induction in the current reaction.24 Among the TSs corresponded to model c, the ts7-S and ts8-R structures are more stable than ts5-S and ts6-R. Because keto carbonyl contains more negative charge than ester carbonyl, the hydrogen-bonding interactions between OH group of catalyst and keto carbonyl of enolate is stronger than that between OH group of catalyst and ester carbonyl of enolate. The ts7-S, leading to the S product, is the most stable transition state among all TSs. The ts8-R is more unfavorable than ts7-S by 1.3 kcal/mol. The 1.3 kcal/mol difference in energy indicated an 81% ee value of S enantiomer, which is in good agreement with the observed 86% ee.42 Even if we analyze all calculated TSs with the Boltzmann equation, the predicted ee would be 71%, which is not in disagreement with the experimental result. Thus, the origin of the enantioselectivity can be attributed to the difference between ts7-S and ts8-R. Additional calculations were performed to confirm the importance of model c in the enantioselective inducement. In order to attenuate the basis set superposition error when energy calculations were conducted with 6-311++G(d,p), the triple-ζ

Figure 1. EDA complexes and maximum absorption wavelength.

structures for an anion-catalyst complex.24a,b The Houk group also found that the most stable comformer of a bis-quaternized cinchoninium cation adopted the anti-open conformation.24c Besides, the X-ray structure of the optimal cinchoninium PTC catalyst with a 2′-position modification of the quinoline ring reported in Melchiorre’s experiments is anti-open conformation, which can serve as a direct reference.16 Considering its resemblance to the N-4-trifluoromethylphenyl cinchoninium cation and bis-quaternized cinchoninium cation and the experimental X-ray data, the most stable structure of N-3,4,5trifluorobenzyl cinchoninium cation should be the anti-open conformation and was used as the initial structure in our Scheme 2. Model Reaction and Putative Transition-State Models

9323

DOI: 10.1021/acs.joc.7b01130 J. Org. Chem. 2017, 82, 9321−9327

Article

The Journal of Organic Chemistry

experiments (see the Supporting Information of ref 16).16 If the proposed model c was indeed in charge of the stereoselectivity and different cinchona-derived PTC catalysts shared almost the same transition states, the observed enantioselectivites with different PTC catalysts could be reproduced via transition state calculations. For example, the 3,5-ditertbutylbenzyl cinchonidine-derived PTC catalyst provided product with only 1% ee in their original reports. We posited that the pseudoenantiomeric 3,5-dimethylbenzyl cinchonine-derived PTC catalyst should also lead to a low level of stereocontrol (Scheme 3). With the Scheme 3. Further Examination of the Hydrogen-Bonding Network Model c with a Cinchonine-Derived PTC Catalyst

established method (M06-2X-D3/6-311++G(d,p)-SMD(CH2Cl2)//M06-2X/6-31G(d)), ts9-S and ts10-R were located (Figure 3). Here again, it should be noted that C6F13I was

Figure 2. Stereodetermining transition state structures and their relative free energies (units: kcal/mol).

basis set def2-TZVPP was first applied.43 The resulting relative free energies were presented with blue color in Figure 2. This more precise method gave similar relative energies. The predicted ee value is 83% with the Boltzmann equation, also indicating model c is the preferred activation model. Performing the geometry optimization in solution phase with implicit solvent model did not change the geometry of transistion states much as well (Figure S2). At the M06-2X-D3/def2TZVPPSMD(PhCl)//M06-2X/6-31G(d)-SMD(PhCl) level, the energy difference between ts7-S1 and ts8-R1 reached up to 2.1 kcal/mol. And the energies of TSs belonging to model a and b are higher at least by 2.5 kcal/mol than that of ts7-S1, meaning the model c is prominent in control of enantioselectivity. Besides, another two model-c-type TSs (syn-ts7 and syn-ts8) with syn-open catalyst conformation were calculated (Figure S3). It turned out energies of transitions states catalyzed by synopen catalyst are much higher. For example, syn-ts7 is 6.9 kcal/ mol above ts7-S in energy. In the end, we noticed the substitution pattern at the benzyl group of the PTC catalyst had a profound influence on the enantioselectivity in Melchiorre’s

Figure 3. Transition state structures of the cinchonine-derived PTC catalyst and their relative free energies (units: kcal/mol).

replaced by CF3I in our calculations. Indeed, the difference of the two TSs in energy is only 0.02 kcal/mol, indicating the enantioselectivity of product is very low (3% ee) when the reaction was conducted in CH2Cl2. This result lent much support to the hydrogen-bonding network model in control of enantioselectivity. The influence of electronic modification of the benzyl group of the PTC catalyst on the enantioselectivity could be explained as follows. The introduction of withdrawing groups into the aromatic ring of the catalyst would strengthen the related hydrogen-bonding interaction between the enolate and the catalyst. For example, the distance of hydrogen bond (C22−H42···O77) in ts7-S is shorter than that (C22−H42··· O74) in ts9-S (2.38 Å versus 2.61 Å), meaning it is stronger in ts7-S. The importance of such a hydrogen bond has been proposed in some cinchona alkaloid-derived PTC catalyzed transformations24a,27b and also in bifunctional thiourea9324

DOI: 10.1021/acs.joc.7b01130 J. Org. Chem. 2017, 82, 9321−9327

Article

The Journal of Organic Chemistry catalyzed reactions.44 Meanwhile, the introduction of withdrawing groups can increase the strength of electrostatic interaction between the positive nitrogen center of the catalyst and the negative enolate, which contributed to the stabilization of the transition states.23,24a The distance (DN45−O77) in ts7-S catalyzed by 3,4,5-trifluorobenzyl cinchoninium is 3.50 Å, while it is 3.53 Å (D N45−O74 ) in ts9-S catalyzed by 3,5ditertbutylbenzyl cinchoninium. Forming a tight complex for the introduction of withdrawing group to the aromatic ring of catalyst would ensure discrimination of one of the enantiotopic faces in the TSs, indicating a better performance of stereocontrolling. In the following, we tried to figure out which hydrogen bond is dominant. Treating ts7-S and ts8-R with Yang’s noncovalent analytic method,45 we can clearly see the existence of hydrogen bond network in these two transition states (Figure 4). Taking

estimated the energies of the hydrogen bonds within key transition states using Espinosa equation. As we show in Table 1, the O−H···O hydrogen bond contributed most to stabilize Table 1. Estimated Strength of Hydrogen Bonds in the Transition States TS

type

ts7-S

O47−H44···O77 C16−H40···O77 C17−H41···O77 C22−H42···O77 C17−H41···O74 C19−H39···O74 C7−H35···O74 O47′-H44′···O77′ C16′-H40′···O77′ C17′-H41′···O77′ C22′-H42′···O77′ C17′-H41′···O74′ C19′-H39′···O74′ C7′-H35′···O74′

ts8-R

length 1.69 2.27 2.07 2.38 2.21 2.15 2.34 1.74 2.37 2.04 2.31 2.19 2.28 2.38

Å Å Å Å Å Å Å Å Å Å Å Å Å Å

energy (kcal/mol) −14.0 −2.8 −5.0 −3.8 −2.2 −3.7 −2.2 −12.1 −2.3 −5.3 −4.1 −2.5 −2.7 −2.1

the corresponding transition states among all the hydrogen bonds. The energy of this type hydrogen bond is −14.0 kcal/ mol in ts7-S and −12.1 kcal/mol in ts8-R, respectively, which is in line with the distance (O47−H44···O77 1.69 Å vs O47′H44′···O77′ 1.74 Å). And the difference of O−H···O hydrogen bond in ts7-S and ts8-R reached 1.9 kcal/mol, reflecting the importance of O−H···O hydrogen bond in control of enantioselectivity. The C19−H39···O74 hydrogen bond in ts7-S (2.15 Å) is 1.0 kcal/mol stronger than that in ts8-R (2.28 Å). The C16−H40···O77 hydrogen bond in ts7-S (2.27 Å) is 0.5 kcal/mol stronger than that in ts8-R (2.37 Å). These two interactions are other factors in control of enantioselectivity. The energy difference of remaining hydrogen bonds is 0.1−0.3 kcal/mol as the tiny difference of distance. Taking all the energy difference of hydrogen bonds together, these interactions make ts7-S more stable than ts8-R by 2.6 kcal/ mol. Though this value (2.6 kcal/mol) derived from hydrogenbonding interactions is somewhat bigger than the difference (1.3 kcal/mol) between ts7-S and ts8-R, the quantitative estimation of hydrogen-bonding strength is a very powerful tool to interpret the origin of enantioselectivity. Compared to the Re face attack, the attack of chiral enolate by trifluoromethyl radical from the Si face is more favorable resulting in multiple hydrogen-bonding interactions, leading to the S product through ts7-S.

Figure 4. Noncovalent interaction analysis of transition states ts7-S and ts8-R.

ts7-S as an example, there are O47−H44···O77, C16−H40··· O77, C17−H41···O77, C22−H42···O77, C17−H41···O74, C19−H39···O74, and C7−H35···O74 hydrogen bonds. The existence of multiple hydrogen bonds makes it still not easy to pick up which hydrogen bond contributes more for the enantioselectivity with the noncovalent analysis, while the quantitative description of the strength of the hydrogen bond will lend much help to give a clear picture of the selectivity. In 1998, Espinosa et al. disclosed the energy of a given hydrogen bond (Eint) is proportionate to the local electron potential energy density V(rcp) of the bond critical point with a coefficient of 0.5.46 With this relationship (Eint = 0.5* V(rcp)), Yu and Cheng have analyzed the strength of hydrogen bonds in the transition states of the asymmetric olefin isomerization47 catalyzed by a cinchona alkaloid derivative in 2013 and successfully revealed the importance of multiple C− H···O hydrogen bonds in the selectivity-control.48 Here we



CONCLUSIONS The stereoselective origin of the photoinduced, asymmetric phase-transfer catalytic perfluoroalkylation of β-ketoesters was interpreted by density functional theory calculations. By means of TD-DFT calculations at M06-2X/6-31G(d)+Lanl2dz(I) level, EDA complexes through a C−I bond paralleled the plane of the enolate interaction rather than halogen-bonding interaction between the CF3I and the enolate are proposed to be visible light photoactive. The examination of three putative transition state models suggests that the simple multiple hydrogen-bonding mode rather than the other two mixed π−π stacking interaction and hydrogen-bonding modes better describes the transition structures of the perfluoroalkylation reaction. The reason why introduction of withdrawing groups 9325

DOI: 10.1021/acs.joc.7b01130 J. Org. Chem. 2017, 82, 9321−9327

Article

The Journal of Organic Chemistry

14642. (g) Yang, Q.; Zhang, L.; Ye, C.; Luo, S.; Wu, L.-Z.; Tung, C.-H. Angew. Chem., Int. Ed. 2017, 56, 3694. (4) (a) Murphy, J. J.; Bastida, D.; Paria, S.; Fagnoni, M.; Melchiorre, P. Nature 2016, 532, 218. (b) Dell'Amico, L.; Fernández-Alvarez, V. M.; Maseras, F.; Melchiorre, P. Angew. Chem., Int. Ed. 2017, 56, 3304. (5) DiRocco, D. A.; Rovis, T. J. Am. Chem. Soc. 2012, 134, 8094. (6) For selected examples, see: (a) Huo, H.; Shen, X.; Wang, C.; Zhang, L.; Rose, P.; Chen, L.-A.; Harms, K.; Marsch, M.; Hilt, G.; Meggers, E. Nature 2014, 515, 100. (b) Du, J.; Skubi, K. L.; Schultz, D. M.; Yoon, T. P. Science 2014, 344, 392. (c) Kainz, Q. M.; Matier, C. D.; Bartoszewicz, A.; Zultanski, S. L.; Peters, J. C.; Fu, G. C. Science 2016, 351, 681. (7) Wei, G.; Zhang, C.; Bureš, F.; Ye, X.; Tan, C.-H.; Jiang, Z. ACS Catal. 2016, 6, 3708. (8) (a) Lian, M.; Li, Z.; Cai, Y.; Meng, Q.; Gao, Z. Chem. - Asian J. 2012, 7, 2019. (b) Wang, Y.; Zheng, Z.; Lian, M.; Yin, H.; Zhao, J.; Meng, Q.; Gao, Z. Green Chem. 2016, 18, 5493. (c) Wang, Y.; Yin, H.; Tang, X.; Wu, Y.; Meng, Q.; Gao, Z. J. Org. Chem. 2016, 81, 7042. (9) Rono, L. J.; Yayla, H. G.; Wang, D. Y.; Armstrong, M. F.; Knowles, R. R. J. Am. Chem. Soc. 2013, 135, 17735. (10) (a) Ding, W.; Zhou, Q.-Q.; Xuan, J.; Li, T.-R.; Lu, L.-Q.; Xiao, W.-J. Tetrahedron Lett. 2014, 55, 4648. (b) Hepburn, H. B.; Melchiorre, P. Chem. Commun. 2016, 52, 3520. (11) Alonso, R.; Bach, T. Angew. Chem., Int. Ed. 2014, 53, 4368. (12) Tröster, A.; Alonso, R.; Bauer, A.; Bach, T. J. Am. Chem. Soc. 2016, 138, 7808. (13) Bergonzini, G.; Schindler, C. S.; Wallentin, C.-J.; Jacobsen, E. N.; Stephenson, C. R. J. Chem. Sci. 2014, 5, 112. (14) (a) Uraguchi, D.; Kinoshita, N.; Kizu, T.; Ooi, T. J. Am. Chem. Soc. 2015, 137, 13768. (b) Kizu, T.; Uraguchi, D.; Ooi, T. J. Org. Chem. 2016, 81, 6953. (15) Yang, Z.; Li, H.; Li, S.; Zhang, M.-T.; Luo, S. Org. Chem. Front. 2017, 4, 1037. (16) Woźniak, Ł.; Murphy, J. J.; Melchiorre, P. J. Am. Chem. Soc. 2015, 137, 5678. (17) For a review on photoactive electron donor−acceptor complexes, see: (a) Lima, C. G. S.; de M Lima, T.; Duarte, M.; Jurberg, I. D.; Paixão, M. W. ACS Catal. 2016, 6, 1389. For selected examples of electron donor−acceptor complexes in photochemitry, see: (b) Arceo, E.; Jurberg, I. D.; Á lvarez-Fernández, A.; Melchiorre, P. Nat. Chem. 2013, 5, 750. (c) Beatty, J. W.; Douglas, J. J.; Miller, R.; McAtee, R. C.; Cole, K. P.; Stephenson, C. R. J. Chem. 2016, 1, 456. (d) Quint, V.; Morlet-Savary, F.; Lohier, J. F.; Lalevee, J.; Gaumont, A. C.; Lakhdar, S. J. Am. Chem. Soc. 2016, 138, 7436. (18) (a) Hashimoto, T.; Maruoka, K. Chem. Rev. 2007, 107, 5656. (b) Shirakawa, S.; Maruoka, K. Angew. Chem., Int. Ed. 2013, 52, 4312. (c) Kaneko, S.; Kumatabara, Y.; Shirakawa, S. Org. Biomol. Chem. 2016, 14, 5367. (d) Liu, S.; Kumatabara, Y.; Shirakawa, S. Green Chem. 2016, 18, 331. (19) For selected reviews, see: (a) Cheong, P. H.-Y.; Legault, C. Y.; Um, J. M.; Ç elebi-Ö lçüm, N.; Houk, K. N. Chem. Rev. 2011, 111, 5042. (b) Houk, K. N.; Cheong, P. H.-Y. Nature 2008, 455, 309. (c) Lam, Y.-h.; Grayson, M. N.; Holland, M. C.; Simon, A.; Houk, K. N. Acc. Chem. Res. 2016, 49, 750. (d) Sperger, T.; Sanhueza, I. A.; Schoenebeck, F. Acc. Chem. Res. 2016, 49, 1311. (20) Tanriver, G.; Dedeoglu, B.; Catak, S.; Aviyente, V. Acc. Chem. Res. 2016, 49, 1250. (21) Lipkowitz, K. B.; Cavanaugh, M. W.; Baker, B.; O'Donnell, M. J. J. Org. Chem. 1991, 56, 5181. (22) Gomez-Bengoa, E.; Linden, A.; López, R.; Múgica-Mendiola, I.; Oiarbide, M.; Palomo, C. J. Am. Chem. Soc. 2008, 130, 7955. (23) Cook, T. C.; Andrus, M. B.; Ess, D. H. Org. Lett. 2012, 14, 5836. (24) (a) de Freitas Martins, E.; Pliego, J. R. ACS Catal. 2013, 3, 613. (b) Martins, E. F.; Pliego, J. R., Jr J. Mol. Catal. A: Chem. 2016, 417, 192. (c) He, C. Q.; Simon, A.; Lam, Y.-h.; Brunskill, A. P. J.; Yasuda, N.; Tan, J.; Hyde, A. M.; Sherer, E. C.; Houk, K. N. J. Org. Chem. 2017, 82, 8645. ́ (25) Vidal-Albalat, A.; Swiderek, K.; Izquierdo, J.; Rodríguez, S.; Moliner, V.; González, F. V. Chem. Commun. 2016, 52, 10060.

into the aromatic ring of the PTC catalyst would gain better enantioselectivity was attributed to a better molecular recognization as the strength of the related hydrogen-bonding and electrostatic interactions between the enolate and the PTC catalyst increased. In the case of N-3,4,5-trifluorobenzyl cinchoninium cation, we believed that multiple hydrogenbonding interactions between the 9-OH moiety and +N−CH site of the cinchoninium ion and enolate are the key factors to obtain high enantioselectivity. These noncovalent interactions within the key transition states were estimated quantitatively and clarified that the O−H···O and specific C−H···O interactions governed the stereocontrolling step. These findings enriched the understanding of Cinchona alkaloids based phasetransfer catalysis in photochemistry and may lend help to design novel phase-transfer catalysts and related reactions.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.joc.7b01130. EDA complexes with TD-DFT descriptions, Gibbs free energies of calculated transition states, stereo-determining transition state structures, transition state structures, and Cartesian coordinates of all the optimized structures (PDF)



AUTHOR INFORMATION

Corresponding Author

*[email protected]. ORCID

Xiao-Song Xue: 0000-0003-4541-8702 Xin Li: 0000-0001-6020-9170 Jin-Pei Cheng: 0000-0001-8822-1577 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We are grateful to the National Natural Science Foundation (Grant Nos. 21390400, 21421062) and the State Key Laboratory on Elemento-organic Chemistry for financial support.



REFERENCES

(1) For selected reviews on photochemistry, see: (a) Prier, C. K.; Rankic, D. A.; MacMillan, D. W. C. Chem. Rev. 2013, 113, 5322. (b) Romero, N. A.; Nicewicz, D. A. Chem. Rev. 2016, 116, 10075. (c) Schultz, D. M.; Yoon, T. P. Science 2014, 343, 1239176. (2) For selected reviews on asymmeric photochemistry, see: (a) Brimioulle, R.; Lenhart, D.; Maturi, M. M.; Bach, T. Angew. Chem., Int. Ed. 2015, 54, 3872. (b) Wang, C.; Lu, Z. Org. Chem. Front. 2015, 2, 179. (c) Wang, D.; Zhang, L.; Luo, S. Huaxue Xuebao 2017, 75, 22. (d) Miyabe, H.; Kawashima, A.; Yoshioka, E.; Kohtani, S. Chem. - Eur. J. 2017, 23, 6225. (e) Zhang, L.; Meggers, E. Acc. Chem. Res. 2017, 50, 320. (3) (a) Nicewicz, D. A.; MacMillan, D. W. C. Science 2008, 322, 77. (b) Cecere, G.; König, C. M.; Alleva, J. L.; MacMillan, D. W. C. J. Am. Chem. Soc. 2013, 135, 11521. (c) Nagib, D. A.; Scott, M. E.; MacMillan, D. W. C. J. Am. Chem. Soc. 2009, 131, 10875. (d) Shih, H.W.; Vander Wal, M. N.; Grange, R. L.; MacMillan, D. W. C. J. Am. Chem. Soc. 2010, 132, 13600. (e) Welin, E. R.; Warkentin, A. A.; Conrad, J. C.; MacMillan, D. W. C. Angew. Chem., Int. Ed. 2015, 54, 9668. (f) Zhu, Y.; Zhang, L.; Luo, S. J. Am. Chem. Soc. 2014, 136, 9326

DOI: 10.1021/acs.joc.7b01130 J. Org. Chem. 2017, 82, 9321−9327

Article

The Journal of Organic Chemistry

Piquemal, J.-P.; Beratan, D. N.; Yang, W. J. Chem. Theory Comput. 2011, 7, 625. (46) Espinosa, E.; Molins, E.; Lecomte, C. Chem. Phys. Lett. 1998, 285, 170. (47) Wu, Y.; Singh, R. P.; Deng, L. J. Am. Chem. Soc. 2011, 133, 12458. (48) Xue, X.-S.; Li, X.; Yu, A.; Yang, C.; Song, C.; Cheng, J.-P. J. Am. Chem. Soc. 2013, 135, 7462.

(26) Capobianco, A.; Di Mola, A.; Intintoli, V.; Massa, A.; Capaccio, V.; Roiser, L.; Waser, M.; Palombi, L. RSC Adv. 2016, 6, 31861. (27) (a) Petrova, G. P.; Li, H.-B.; Maruoka, K.; Morokuma, K. J. Phys. Chem. B 2014, 118, 5154. (b) Kamachi, T.; Yoshizawa, K. Org. Lett. 2014, 16, 472. (c) Kamachi, T.; Yoshizawa, K. J. Chem. Inf. Model. 2016, 56, 347. (d) Cannizzaro, C. E.; Houk, K. N. J. Am. Chem. Soc. 2002, 124, 7163. (28) (a) Zhang, Q.; Zhang, Z.-Q.; Fu, Y.; Yu, H.-Z. ACS Catal. 2016, 6, 798. (b) Le Vaillant, F.; Wodrich, M. D.; Waser, J. Chem. Sci. 2017, 8, 1790. (c) Gutierrez, O.; Tellis, J. C.; Primer, D. N.; Molander, G. A.; Kozlowski, M. C. J. Am. Chem. Soc. 2015, 137, 4896. (d) FernandezAlvarez, V. M.; Nappi, M.; Melchiorre, P.; Maseras, F. Org. Lett. 2015, 17, 2676. (e) Demissie, T. B.; Ruud, K.; Hansen, J. H. Organometallics 2015, 34, 4218. (29) (a) Bahamonde, A.; Melchiorre, P. J. Am. Chem. Soc. 2016, 138, 8019. (b) Johnson, M. W.; Hannoun, K. I.; Tan, Y.; Fu, G.; Peters, J. C. Chem. Sci. 2016, 7, 4091. (c) Cismesia, M. A.; Yoon, T. P. Chem. Sci. 2015, 6, 5426. (d) Tsudaka, T.; Kotani, H.; Ohkubo, K.; Nakagawa, T.; Tkachenko, N. V.; Lemmetyinen, H.; Fukuzumi, S. Chem. - Eur. J. 2017, 23, 1306. (e) Pan, X.; Fang, C.; Fantin, M.; Malhotra, N.; So, W. Y.; Peteanu, L. A.; Isse, A. A.; Gennaro, A.; Liu, P.; Matyjaszewski, K. J. Am. Chem. Soc. 2016, 138, 2411. (f) Bahamonde, A.; Murphy, J. J.; Savarese, M.; Brémond, É.; Cavalli, A.; Melchiorre, P. J. Am. Chem. Soc. 2017, 139, 4559. (30) 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.; Keith, T.; 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, O.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09, revision D.01; Gaussian, Inc., 2013. (31) Legault, C. Y. CYLview, 1.0b; Université de Sherbrooke: Sherbrooke (Québec) Canada, 2009. (32) Lu, T.; Chen, F. J. Comput. Chem. 2012, 33, 580. (33) Humphrey, W.; Dalke, A.; Schulten, K. J. Mol. Graphics 1996, 14, 33. (34) (a) Zhao, Y.; Truhlar, D. G. Acc. Chem. Res. 2008, 41, 157. (b) Zhao, Y.; Truhlar, D. Theor. Chem. Acc. 2008, 120, 215. (35) Hay, P. J.; Wadt, W. R. J. Chem. Phys. 1985, 82, 299. (36) Yanai, T.; Tew, D. P.; Handy, N. C. Chem. Phys. Lett. 2004, 393, 51. (37) Grimme, S. WIREs: Comput. Mol. Sci. 2011, 1, 211. (38) Weigend, F.; Ahlrichs, R. Phys. Chem. Chem. Phys. 2005, 7, 3297. (39) Marenich, A. V.; Cramer, C. J.; Truhlar, D. G. J. Phys. Chem. B 2009, 113, 6378. (40) Adamo, C.; Jacquemin, D. Chem. Soc. Rev. 2013, 42, 845. (41) Wang, Y.; Wang, J.; Li, G.-X.; He, G.; Chen, G. Org. Lett. 2017, 19, 1442. (42) The experimental value (86% ee) was quoted from entry 7 in Table 1 from Melchiorre’s original publication, see ref 16. This value is obtained from the reaction of methyl 1-oxo-2,3-dihydro-1H-indene-2carboxylate and C6F13I. It should be noted that we use CF3I to replace C6F13I in our calculations. (43) Armstrong, A.; Boto, R. A.; Dingwall, P.; Contreras-Garcia, J.; Harvey, M. J.; Mason, N. J.; Rzepa, H. S. Chem. Sci. 2014, 5, 2057. (44) Lippert, K. M.; Hof, K.; Gerbig, D.; Ley, D.; Hausmann, H.; Guenther, S.; Schreiner, P. R. Eur. J. Org. Chem. 2012, 2012, 5919. (45) (a) Johnson, E. R.; Keinan, S.; Mori-Sanchez, P.; ContrerasGarcia, J.; Cohen, A. J.; Yang, W. J. Am. Chem. Soc. 2010, 132, 6498. (b) Contreras-Garcia, J.; Johnson, E. R.; Keinan, S.; Chaudret, R.; 9327

DOI: 10.1021/acs.joc.7b01130 J. Org. Chem. 2017, 82, 9321−9327