Highly enantioselective construction of CF3-bearing all-carbon quaternary stereocenters: Chiral spiro-fused bisoxazoline ligands with 1, 1′-binaphthyl sidearm for asymmetric Michael-type Friedel-Crafts reaction

Robert Li-Yuan Bao Lei Shi Kang Fu

Citation:  Robert Li-Yuan Bao, Lei Shi, Kang Fu. Highly enantioselective construction of CF3-bearing all-carbon quaternary stereocenters: Chiral spiro-fused bisoxazoline ligands with 1, 1′-binaphthyl sidearm for asymmetric Michael-type Friedel-Crafts reaction[J]. Chinese Chemical Letters, 2022, 33(5): 2415-2419. doi: 10.1016/j.cclet.2021.10.062 shu

Highly enantioselective construction of CF3-bearing all-carbon quaternary stereocenters: Chiral spiro-fused bisoxazoline ligands with 1, 1′-binaphthyl sidearm for asymmetric Michael-type Friedel-Crafts reaction

English

  • Bisoxazoline ligands (box ligands), as one class of the privileged ligands [1, 2], have made a great contribution to the progress of asymmetric catalysis over the past three decades [3-6]. In the course of developing chiral box ligands, modifying the parent bisoxazoline ligands [7-9] has emerged as an efficient approach to seek novel chiral ligands. By using this strategy, a series of representative chiral bisoxazoline ligands were successfully designed and applied in the asymmetric radical addition reactions [10], Friedel-Crafts reactions [11], cyclopropanation reactions [12], aerobic oxidation reactions [13], and radical coupling reactions (Fig. 1a) [14]. The success of these representative chiral bisoxazoline ligands benefited from the C2 symmetry of the ligands [15] and the good coordination ability of the bisoxazolines, providing valuable experience for the development of novel ligands.

    Figure 1

    Figure 1.  Design and application of the chiral spiro-fused bisoxazoline ligands.

    Although plenty of representative chiral bisoxazoline ligands have been prepared to date, the development of novel chiral bisoxazoline ligands is still of central importance and highly desirable in the context of expanding the types of chiral ligands and solving the new synthetic challenges. Inspired by the structural features of the chiral spiro monooxazoline in the Spymox ligand (Fig. 1b, left) [16-18] and the deep chiral pocket in the WingPhos ligand (Fig. 1b, middle), [19] we designed a class of chiral spiro-fused bisoxazoline ligands based on the modification of the parent bisoxazoline ligands (Fig. 1b, right). In our proposal, the chiral spiro-fused bisoxazoline ligands can be constructed by connecting the chiral binaphthyl moieties to the oxazoline rings. The developed chiral spiro-fused bisoxazoline ligands possess three characteristics: (1) The rigidity of the spirocyclic units can decrease the flexibility of these bisoxazoline ligands and their related complexes, which facilitates the control of enantioselectivity in asymmetric catalysis. (2) The C2 symmetry of the chiral spiro-fused bisoxazoline ligands can reduce the number of possible transition states in asymmetric reactions [14]. (3) Two chiral binaphthyl sidearms can protrude directly forward to the coordination site of substrate, generating a deep chiral pocket. In this communication, we describe the preparation and application of the chiral spiro-fused bisoxazoline ligands.

    All-carbon quaternary stereocenters exist in numerous pharmaceuticals and bioactive natural products [20-22]. Although the chiral all-carbon quaternary stereocenters could be constructed effectively, the catalytically asymmetric synthesis of the trifluoromethylated all-carbon quaternary centers remains a synthetic challenge [23-27]. The asymmetric Michael-type Friedel-Crafts reaction of indole with the β-CF3-β-disubstituted nitroalkenes was a direct and efficient approach to prepare the active indole derivatives [28, 29] bearing a trifluoromethylated all-carbon quaternary stereocenter. As a pioneering work, Gao et al. utilized a Ni(ClO4)2-bisoxazoline complex as catalyst to prepare the active indole derivatives (up to 97% ee) [30, 31]. Similarly, Zhu et al. exploited a Ni(Ⅱ)/(imidazoline-oxazoline) complex as catalyst to produce the trifluoromethylated bis(3-indolyl)methanes (up to 94% ee) [32]. These two cases implied that the chiral catalysts based on the oxazoline ligands might have an advantage in this asymmetric Michael-type Friedel-Crafts reaction. Additionally, the studies from Akiyama group and Meggers group demonstrated that the chiral calcium BINOL bis(phosphate) complex [33] and the biscyclometalated iridium(Ⅲ) complex [34] could be employed as catalysts to construct the trifluoromethylated all-carbon quaternary stereocenters as well (up to 98% ee). In order to assess the catalytic performances of the chiral spiro-fused bisoxazoline ligands and improve the enantioselectivity in the construction of the trifluoromethylated all-carbon quaternary stereocenters, the chiral spiro-fused bisoxazoline ligands were employed in the Michael-type Friedel-Crafts reaction of indole with β-CF3-β-disubstituted nitroalkenes.

    The chiral spiro-fused bisoxazoline ligands were synthesized according to the procedures modified from the route to prepare the Spymox ligands in Shibatomi's reports. [16-18]. The synthesis started from an amidation of the unnatural chiral 1, 2-aminol [16, 35] with the disubstituted malonyl dichloride, which generated the desired bis-hydroxyamide compounds [9]. Subsequently, the hydroxy groups of the bis-hydroxyamide compounds were converted to good leaving groups via mesylation. Finally, the chiral spiro-fused bisoxazoline ligands (Fig. 2, L1L3), which were named as (R, R)-Bin-box ligands, were obtained under basic conditions (See the Section 2.1 in Supporting information for the synthetic details).

    Figure 2

    Figure 2.  Preparation of the chiral spiro-fused bisoxazoline ligands.

    With the chiral spiro-fused bisoxazoline ligands in hand, we studied their catalytic performances in the Michael-type Friedel-Crafts reaction using indole 1a and the (E)−1-phenyl-1-trifluoromethyl-2-nitroethene 2a as model substrates. Investigation of a series of parameters revealed that the Ni(ClO4)2·6H2O complex of the ligand L1 was the optimal catalyst. In presence of this catalyst, the reaction proceeded smoothly in toluene at 85 ℃, delivering the product 3a with excellent yield (99%) and excellent enantioselectivity (> 99% ee). The enantioselectivity was higher than that induced by the optimized ligand in Jia's report [30], due to the unique effects of ligand structure and the noncovalent interactions in the ligand L1 (Table 1, entry 1). Although the ligands L2 and L3 were similar to the ligand L1, they were not suitable for this asymmetric reaction because the unique chelate bite angles and the sizes of the bridging substituents of L2 and L3 might decrease the reactivities and enantioselectivities (Table 1, entries 2 and 3) [36, 37]. Notably, only poor yield and poor enantioselectivity were achieved in the presence of L4, which implied that C2 symmetric ligand L1 possessed great advantages over the chiral spiro monooxazoline ligand L4 in this asymmetric transformation (Table 1, entry 4). In addition, when Ni(OTf)2 was utilized in this reaction, the yield would decrease obviously (Table 1, entry 5). Moreover, the enantioselectivity would decline slightly when the amount of the 1a was increased to 2 equiv. or the loading of ligand L1 was raised to 12 mol% (Table 1, entries 6 and 7). At last, if the reaction temperature was raised to 100 ℃, the reaction time could be reduced to 12 h and the enantioselectivity only decreased slightly (98.6% ee), which indicated that higher temperature might be used as an alternative option in this asymmetric transformation.

    Table 1

    Table 1.  Optimization of the reaction conditions.a
    DownLoad: CSV

    To evaluate the generality of this asymmetric Michael-type Friedel-Crafts reaction, we investigated the scope of the substituted indoles and β-CF3-β-disubstituted nitroalkenes (Scheme 1). Initially, the indoles bearing a substituent at the C5 position were subjected to the reaction. The results demonstrated that indole with an electron-donating group (3b) possessed an advantage over that bearing an electron-withdrawing group (3c and 3d) both in the reactivity and enantioselectivity. Next, the effect of indoles bearing different functional groups at the C6 position was examined (3e3k). To our delight, the catalytic system was highly compatible with diverse functional groups at the C6 position of indole ring and the corresponding products were obtained in good to excellent yields (87%−99%) with excellent enantioselectivities (> 99%– > 99.9% ee). In particular, the highly enantiopure products (3f and 3k) were observed (> 99.9% ee) when 6-benzyloxyindole and methyl indole-6-carboxylate were subjected in this asymmetric reaction, respectively. Additionally, the 7-methoxy-1H-indole and 7-methylindole could deliver the products in excellent enantioselectivities as well, although the reactivities decreased because of the influence of steric hindrance (3i and 3m).

    Scheme 1

    Scheme 1.  Substrate scope of the asymmetric Michael-type Friedel-Crafts reaction. Standard reaction conditions: 1a (0.15 mmol), 2a (0.1 mmol), Ni(ClO4)2·6H2O (11 mol%), L1 (10 mol%) in anhydrous toluene (1 mL) under Ar atmosphere at 85 ℃ for 36 h. The yields were the isolated yields. The ee values were determined by HPLC. a At 100 ℃ for 96 h. b At 100 ℃ for 36 h.

    The scope of β-CF3-β-disubstituted nitroalkenes was also assessed under the optimal reaction conditions (3n3v). When the para-substituted (E)−1-phenyl-1-trifluoromethyl-2-nitroethenes were used in the reaction, different compatibilities were observed. The nitroalkene bearing an electron-donating group gave the product (3n) with excellent yield (99%) and excellent enantioselectivity (96% ee). On the other hand, the nitroalkenes bearing an electron-withdrawing group delivered the products (3o3q) with good to excellent enantioselectivities (78%−96% ee), indicating that a stronger electron-withdrawing property would lead to a lower enantioselectivity. After that, the effect of meta-substituted (E)−1-phenyl-1-trifluoromethyl-2-nitroethenes were studied. The excellent enantioselectivities were observed whether the nitroalkene possessed an electron-donating or electron-withdrawing substituent (3r and 3s), while the reactivity would decrease obviously in the presence of the electron-withdrawing substituent (3s). Subsequently, it was found that the β-CF3-β-naphthyl substituted nitroalkene could be converted to the desired product (3t) successfully at 100 ℃ in high yield (93%) with excellent enantioselectivity (99.7% ee). Furthermore, the β-CF3-β-alkyl-disubstituted nitroalkenes could also afford the corresponding products smoothly with moderate to excellent enantioselectivities (3u and 3v).

    In order to expand the types of substrates, we also investigated the reactive properties of pyrrole, β-CF2H-β-disubstituted nitroalkene and alkyl substituted nitroalkenes, respectively. To our delight, pyrrole could deliver the product 3w in a moderate yield (66%) with an excellent enantioselectivity (98% ee), which suggested that the chiral Ni(Ⅱ)-Bin-box complex has large potential in the enantioselective Friedel-Crafts alkylation of 4, 7-dihydroindoles [38]. Next, the β-CF2H-substituted nitroalkene was investigated in the asymmetric reaction, affording the product 3x with 90% yield and 98% ee. The high yield and excellent enantioselectivity implied that the noncovalent interactions of the fluorinated groups were beneficial for the asymmetric induction [33]. Additionally, other nitroacrylates and the methyl substituted nitroalkene could give the corresponding products (3y, 3z and 3aa) with obviously decreased enantioselectivities in the reaction, owing to the attenuation of noncovalent interactions between the substrates and chiral ligand. The outcomes of the 3y, 3z and 3aa suggested that the stereo control of the chiral spiro-fused bisoxazoline ligand was different with the chiral box ligands bearing two phenyl groups in the asymmetric alkylation of α-substituted β-nitroacrylates [39] and β-aryl β-alkyl nitroalkenes [40]. In the end, 1-methylindole was subjected to the Michael-type Friedel-Crafts reaction, whereas no target product 3ab was obtained, indicating that N-H bond of indole was necessary in this nickel(Ⅱ)-catalyzed asymmetric alkylation.

    The absolute configurations of the known indole derivatives 3a–3d, 3i, 3m, 3n, 3p, 3r, 3t, 3u and 3x were assigned as (R) configuration and the known products 3y, 3z and 3aa exhibited (S) configuration, which were identified by comparing with the data of optical rotation from previous studies [30, 33, 34].

    Apart from the nickel(Ⅱ) catalyst based on the chiral ligand L1, four types of chiral catalysts have been utilized in the Enanti–Oselective synthesis of 3a [18-21]. As shown in Fig. 3, chiral catalyst Ni(Ⅱ)/L1 had an advantage over other chiral catalysts in terms of the enantioselectivity (Fig. 3, entry 8). By analyzing the results related to the chiral oxazoline ligands at 100 ℃ (Fig. 3, entries 1–8), it was found that the ligand with a deeper chiral pocket could extremely enhance the confinement of substrate 2a, leading to a better enantioselectivity in this asymmetric reaction.

    Figure 3

    Figure 3.  Chiral catalysts utilized in the enantioselective synthesis of 3a.

    To understand the steric, electronic and stereochemical characteristics of the chiral spiro-fused bisoxazoline ligand L1, a model of chiral pocket was established according to the representation from Trost and Oslop [41-43], which might be effective in this asymmetric reaction. The naphthyl moieties with different stereochemical features played diverse roles, such as a flap or wall, forming a confined reaction space residing inside the deep chiral pocket (Figs. 4a and b) [44]. In the asymmetric addition of indole 1a with β-CF3-substituted nitroalkene 2a, the conformation of the activated 2a was stabilized by the noncovalent interactions [45], which were possibly enhanced by the fluorine-π interaction [46, 47] between the naphthyl moiety D and substrate 2a (Fig. 4c). Meanwhile, the deep chiral pocket reduced the number of possible trajectories of the incoming indole nucleophile, generating the desired product (R)−3a (Fig. 4d).

    Figure 4

    Figure 4.  Proposed model of the deep chiral pocket for asymmetric induction.

    In conclusion, we designed and synthesized the chiral spiro-fused bisoxazoline ligands, which have been successfully employed in the asymmetric Michael-type Friedel-Crafts reaction of substituted indole with β-CF3-β-disubstituted nitroalkenes. A range of indole derivatives bearing a trifluoromethylated all-carbon quaternary stereocenter were constructed with moderate to excellent yields and good to excellent enantioselectivities. Moreover, a proposed model of stereoinduction suggested that the deep chiral pocket could shield the Si-face of the β-CF3-β-disubstituted nitroalkenes effectively. The further applications of the chiral spiro-fused bisoxazoline ligands are ongoing in our lab.

    The authors report no declarations of interest and the authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    We are grateful for the financial support from the National Natural Science Foundation of China (No. 21871067), the Natural Science Foundation of Guangdong Province (Nos. 2018A030313038 and 2021A1515010190), the Shenzhen Fundamental Research Projects (No. JCYJ20180306171838187), the Harbin Institute of Technology (Shenzhen) (Talent Development Starting Fund from Shenzhen Government), the Open Project Program of State Key Laboratory of Elemento-Organic Chemistry (No. 202009) and the Guangdong Provincial Key Laboratory of Catalysis (No. 2020B121201002). Robert Li-Yuan Bao also thanks Dr. Rong Zhao for the helpful discussion about the manuscript, and Prof. Min Wang (College of Materials, Chemistry and Chemical Engineering, Hangzhou Normal University) for the assistance on the measurement of the HRMS.

    Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.cclet.2021.10.062.


    1. [1]

      Q.L. Zhou, in: Privileged Chiral Ligands and Catalysts, Wiley-VCH, Weinheim, 2011, pp. 171-212.

    2. [2]

      T.P. Yoon, E.N. Jacobsen, Science 299(2003) 1691-1693. doi: 10.1126/science.1083622

    3. [3]

      J.S. Johnson, D.A. Evans, Acc. Chem. Res. 33(2000) 325-335. doi: 10.1021/ar960062n

    4. [4]

      G. Desimoni, G. Faita, K.A. Jørgensen, Chem. Rev. 106(2006) 3561-3651. doi: 10.1021/cr0505324

    5. [5]

      G. Desimoni, G. Faita, K.A. Jørgensen, Chem. Rev. 111(2011) PR284-PR437. doi: 10.1021/cr100339a

    6. [6]

      S. Liao, X.L. Sun, Y. Tang, Acc. Chem. Res. 47(2014) 2260-2272. doi: 10.1021/ar800104y

    7. [7]

      R.E. Lowenthal, A. Abiko, S. Masamune, Tetrahedron Lett. 31(1990) 6005-6008. doi: 10.1016/S0040-4039(00)98014-6

    8. [8]

      D.A. Evans, K.A. Woerpel, M.M. Hinman, M.M. Faul, J. Am. Chem. Soc. 113(1991) 726-728. doi: 10.1021/ja00002a080

    9. [9]

      E.J. Corey, N. Imai, H.Y. Zhang, J. Am. Chem. Soc. 113(1991) 729-730. doi: 10.1021/ja00002a082

    10. [10]

      M.P. Sibi, J. Ji, J. Org. Chem. 62(1997) 3800-3801. doi: 10.1021/jo970558y

    11. [11]

      J. Zhou, Y. Tang, J. Am. Chem. Soc. 124(2002) 9030-9031. doi: 10.1021/ja026936k

    12. [12]

      S. Bellemin-Laponnaz, L.H. Gade, Angew. Chem. Int. Ed. 41(2002) 3473-3475. doi: 10.1002/1521-3773(20020916)41:18<3473::AID-ANIE3473>3.0.CO;2-N

    13. [13]

      W. Ding, L.Q. Lu, Q.Q. Zhou, et al., J. Am. Chem. Soc. 139(2017) 63-66. doi: 10.1021/jacs.6b11418

    14. [14]

      C. Jiang, L. Wang, H. Zhang, et al., Chem 6(2020) 2407-2419. doi: 10.1016/j.chempr.2020.07.003

    15. [15]

      J.K. Whitesell, Chem. Rev. 89(1989) 1581-1590. doi: 10.1021/cr00097a012

    16. [16]

      K. Shibatomi, T. Muto, Y. Sumikawa, A. Narayama, S. Iwasa, Synlett 2(2009) 241-244.

    17. [17]

      K. Shibatomi, A. Narayama, Y. Soga, T. Muto, S. Iwasa, Org. Lett. 13(2011) 2944-2947. doi: 10.1021/ol201007e

    18. [18]

      K. Shibatomi, Y. Soga, A. Narayama, I. Fujisawa, S. Iwasa, J. Am. Chem. Soc. 134(2012) 9836-9839. doi: 10.1021/ja304806j

    19. [19]

      G. Liu, X. Liu, Z. Cai, et al., Angew. Chem. Int. Ed. 52(2013) 4235-4238. doi: 10.1002/anie.201300646

    20. [20]

      Z.L. Song, C.A. Fan, Y.Q. Tu, Chem. Rev. 111(2011) 7523-7556. doi: 10.1021/cr200055g

    21. [21]

      R. Long, J. Huang, J. Gong, Z. Yang, Nat. Prod. Rep. 32(2015) 1584-1601. doi: 10.1039/C5NP00046G

    22. [22]

      Z. Wang, Org. Chem. Front. 7(2020) 3815-3841. doi: 10.1039/D0QO00763C

    23. [23]

      B.M. Trost, C. Jiang, Synthesis 3(2006) 369-396.

    24. [24]

      B. Wang, Y.Q. Tu, Acc. Chem. Res. 44(2011) 1207-1222. doi: 10.1021/ar200082p

    25. [25]

      K.W. Quasdorf, L.E. Overman, Nature 516(2014) 181-191. doi: 10.1038/nature14007

    26. [26]

      X.P. Zeng, Z.Y. Cao, Y.H. Wang, F. Zhou, J. Zhou, Chem. Rev. 116(2016) 7330-7396. doi: 10.1021/acs.chemrev.6b00094

    27. [27]

      J. Nie, H.C. Guo, D. Cahard, J.A. Ma, Chem. Rev. 111(2011) 455-529. doi: 10.1021/cr100166a

    28. [28]

      S. Lancianesi, A. Palmieri, M. Petrini, Chem. Rev. 114(2014) 7108-7149. doi: 10.1021/cr400676v

    29. [29]

      J.B. Chen, Y.X. Jia, Org. Biomol. Chem. 15(2017) 3550-3567. doi: 10.1039/C7OB00413C

    30. [30]

      J.R. Gao, H. Wu, B. Xiang, et al., J. Am. Chem. Soc. 135(2013) 2983-2986. doi: 10.1021/ja400650m

    31. [31]

      H. Wu, R.R. Liu, Y.X. Jia, Synlett 25(2014) 457-460.

    32. [32]

      W.J. Zhu, J.F. Gong, M.P. Song, J. Org. Chem. 85(2020) 9525-9537. doi: 10.1021/acs.joc.0c00336

    33. [33]

      I. Ibánz, M. Kaneko, Y. Kamei, et al., ACS Catal. 9(2019) 6903-6909. doi: 10.1021/acscatal.9b01811

    34. [34]

      W. Xu, X. Shen, Q. Ma, et al., ACS Catal. 6(2016) 7641-7646. doi: 10.1021/acscatal.6b02080

    35. [35]

      J.P. Mazaleyrat, A. Gaucher, M. Wakselmaa, L. Tchertanov, J. Guilhemb, Tetrahedron Lett. 37(1996) 2971-2974. doi: 10.1016/0040-4039(96)00502-3

    36. [36]

      I.W. Davies, L. Gerena, L. Castonguay, et al., Chem. Commun. 15(1996) 1753-1754.

    37. [37]

      S.E. Denmark, C.M. Stiff, J. Org. Chem. 65(2000) 5875-5878. doi: 10.1021/jo0007175

    38. [38]

      H. Wu, R.R. Liu, C. Shen, et al., Org. Chem. Front. 2(2015) 124-126. doi: 10.1039/C4QO00265B

    39. [39]

      J.Q. Weng, Q.M. Deng, L. Wu, et al., Org. Lett. 16(2014) 776-779. doi: 10.1021/ol403480v

    40. [40]

      H. Wu, W.J. Sheng, B. Chen, et al., Synlett 26(2015) 2817-2820. doi: 10.1055/s-0035-1560538

    41. [41]

      B.M. Trost, J.D. Oslop, J. Am. Chem. Soc. 121(1999) 3057-3064. doi: 10.1021/ja983617d

    42. [42]

      B.M. Trost, F.D. Toste, J. Am. Chem. Soc. 121(1999) 4545-4554. doi: 10.1021/ja9828713

    43. [43]

      B.M. Trost, M.R. Machacek, A. Aponick, Acc. Chem. Res. 39(2006) 747-760. doi: 10.1021/ar040063c

    44. [44]

      B. Mitschke, M. Turberg, B. List, Chem 6(2020) 2515-2532. doi: 10.1016/j.chempr.2020.09.007

    45. [45]

      A. Fanourakis, P.J. Docherty, P. Chuentragool, R.J. Phipps, ACS Catal. 10(2020) 10672-10714. doi: 10.1021/acscatal.0c02957

    46. [46]

      M.C. Sherman, M.R. Ams, K.D. Jordan, J. Phys. Chem. A 120(2016) 9292-9298. doi: 10.1021/acs.jpca.6b09193

    47. [47]

      P. Li, J.M. Maier, E.C. Vik, et al., Angew. Chem. Int. Ed. 56(2017) 7209-7212. doi: 10.1002/anie.201702950

  • Figure 1  Design and application of the chiral spiro-fused bisoxazoline ligands.

    Figure 2  Preparation of the chiral spiro-fused bisoxazoline ligands.

    Scheme 1  Substrate scope of the asymmetric Michael-type Friedel-Crafts reaction. Standard reaction conditions: 1a (0.15 mmol), 2a (0.1 mmol), Ni(ClO4)2·6H2O (11 mol%), L1 (10 mol%) in anhydrous toluene (1 mL) under Ar atmosphere at 85 ℃ for 36 h. The yields were the isolated yields. The ee values were determined by HPLC. a At 100 ℃ for 96 h. b At 100 ℃ for 36 h.

    Figure 3  Chiral catalysts utilized in the enantioselective synthesis of 3a.

    Figure 4  Proposed model of the deep chiral pocket for asymmetric induction.

    Table 1.  Optimization of the reaction conditions.a

    下载: 导出CSV
  • 加载中
计量
  • PDF下载量:  10
  • 文章访问数:  1267
  • HTML全文浏览量:  96
文章相关
  • 发布日期:  2022-05-15
  • 收稿日期:  2021-05-11
  • 接受日期:  2021-10-21
  • 修回日期:  2021-10-21
  • 网络出版日期:  2021-10-27
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

/

返回文章