

Control of meta-selectivity in the Ir-catalyzed aromatic C-H borylation directed by hydrogen bond interaction: A combined computational and experimental study
English
Control of meta-selectivity in the Ir-catalyzed aromatic C-H borylation directed by hydrogen bond interaction: A combined computational and experimental study
-
Key words:
- Meta-C-H borylation
- / Regioselectivity
- / DFT calculation
- / Ligand design
- / Hydrogen bond
-
Regioselective C—H bond activation of arenes provides an efficient approach for obtaining functionalized aromatics in organic synthesis [1-3]. The directed ortho-selective C—H functionalization of arenes is well established through transition-metal catalysis [4-7]. Compared with proximal ortho-C-H bonds, the meta- and para-C-H bonds are away from the directing group but are adjacent to each other, leading to more difficulties in achieving the high regioselectivity [8,9]. In recent years, quite a few elegant methods are developed to control remote selective C—H activation, such as σ-bond activation-assisted functionalization [10-13], template strategy [14-22], noncovalent interaction direction [23-33], Pd(Ⅱ)/norbornene cooperative catalysis [34-36], traceless directing group strategy [37], and catalyst or reagent controlled strategy [38-43]. Recently, noncovalent interaction directed iridium-catalyzed borylation becomes a robust approach [44-47] because the borylated products can be served as versatile synthetic building blocks in chemical and industrial synthesis [48-51].
In 2015, Kanai and coworkers reported a meta-selective C—H borylation of arenes using a bipyridine-derived ligand with a urea moiety [23]. Aromatic amides, phosphonates and phosphine oxides are suitable for this reaction, but the ratios of meta-selectivity highly depend on the substrates, ranging from 30:1 to 0.5:1. The authors claimed that hydrogen bond interaction between ligand and substrate placed the iridium center in close proximity to the meta-C-H bond of substrates and controlled the regioselectivity (Fig. 1A). The mechanism of Ir-catalyzed C—H borylation of arenes has been established by Sakaki [52], Hartwig [53], Houk [54] and Ke [55] groups, in which the rate-determining step is C—H bond oxidative addition (Fig. 1B). Sunoj and coworkers reported a computational study on Kanai's work, and they also proved that the aromatic C—H activation is the regioselectivity-determining step. In addition, they demonstrated that the observed high meta-selectivity was predominantly attributed to a good number of noncovalent interactions between catalyst/ligand and substrate [56]. However, through a combined computational and experimental study, we found that the origin of meta-selectivity is much more complex than expected.
Figure 1
Figure 1. (A) Hydrogen bond interaction directed meta-selective C-H borylation proposed by Kanai and coworkers. (B) Catalytic cycle for the Ir-catalyzed borylation of arenes using B2pin2.We first explored the electronic effect of substituents in the urea moiety of ligand to validate the important role of hydrogen bond interaction in controlling regioselectivity. We used 2,2′-bipyridine (bpy) as ligand to evaluate the intrinsic borylation selectivity of benzamide 1a, which gave a meta to para ratio of 0.7:1 (Fig. 2). Meanwhile, ligands L1-L5 were synthesized and tested in the iridium-catalyzed borylation of benzamide 1a. As shown in Fig. 2, for L1 and L2, with electron-donating groups (OMe and Me), the ratios of meta to para were 2.9:1 and 3.2:1, respectively. For ligand L3, a slightly increased meta-selectivity (3.9:1) was observed. However, ligands L4 and L5 with electron-withdrawing groups (Cl and CF3) exhibited even lower meta-selectivities (1.4:1 and 1.1:1). Especially for L5, the selectivity of borylated products was very close to that using bpy ligand without hydrogen bond interaction. It is well known that ureas bearing electron-withdrawing groups are better hydrogen bond donors than those with electron-donating groups, and it should enhance the hydrogen bond interaction between ligand and substrate [57]. With stronger hydrogen bonds, the meta-selectivity would have been further improved, rather than becoming much worse. What is the deeper reason behind this unusual experimental phenomenon?
Figure 2
To answer this question, DFT calculations were conducted using L5 as ligand (Fig. 3). According to Kanai's hydrogen bond recognition model shown in Fig. 1A, we located the meta- and para-C-H activation transition states TS1-meta and TS1-para directed by hydrogen bond interaction. Computational results showed that the free energy of TS1-meta is 3.6 kcal/mol lower than that of TS1-para, implying that L5 favors meta-selectivity greatly. This result was obviously inconsistent with our experimental data, in which only 1.1:1 of meta to para selectivity was obtained. We then conducted extensive conformational search and discovered another two non-directed transition states TS1-meta-ND and TS1-para-ND with even lower energies, which were not reported in previous study by Sunoj and coworkers [56]. In the non-directed pathway, the urea moiety in ligand recognizes the O atom in Bpin bound to iridium center via N—H⋯O hydrogen bonding, rather than the C=O group in substrate. Taking these four transition states into account (Fig. 3), it is now clear that the non-directed TS1-meta-ND and TS1-para-ND control the ratio of the meta- and para-borylated products. As their free energy difference is only 0.6 kcal/mol, predicting an approximate 1:1 mixture of meta- and para-borylated products, this is in good agreement with the experimental result (L5, m: p = 1.1, Fig. 2). These calculations indicated that the previouly overlooked non-directed pathway played a critical role in determining regioselectivity, even overriding the directed pathway.
Figure 3
Figure 3. DFT-computed geometries and Gibbs free energies for the C-H activation transition states of 1a-Me using B2pin2 and L5 as ligand, computed with SMD(p-xylene)-ωB97X-D/6–311+G(d, p)[SDD for Ir]//M06/6–31G(d)[SDD for Ir]. For improved clarity, most hydrogen atoms are omitted. All distances are in Å.As shown in Scheme 1, Kanai and coworkers reported that, in the presence of B2pin2, the borylation of 1a gave a 13:1 ratio of meta to para using ligand L6 [23]. We speculated that the non-directed pathway would also exist in this case, and that it would become predominant when the B2pin2 was replaced by a smaller diboron reagent, which favors the hydrogen bond interaction between the boron group and the urea moiety of ligand, leading to deteriorated meta-selectivity. To test the steric effect of the boron group, DFT calculations were carried out using the less bulky diboron reagent DB-1 (Scheme 1), and four transition states were located (Fig. 4). As expected, both non-directed transition states TS2-meta-ND and TS2-para-ND had lower free energies than two directed transition states TS2-meta and TS2-para, indicating that the regioselectivity was again determined by the non-directed pathway. The negligible energy difference of 0.2 kcal/mol predicted a poor regioselectivity. To confirm this prediction, DB-1 was synthesized and tested using the same ligand L6. Indeed, it delivered meta- and para-borylated products as a 1:1 isomeric mixture (Scheme 1). On the other hand, the diboron reagent with bulkier substituent was also investigated. B2Epin2 (DB-2), a more stable boronic ester than B2pin2, was synthesized via two steps from 3-pentanone according to a known procedure [58]. Unfortunately, B2Epin2 exhibited poor reactivity in the borylation of 1a using L6 (conv. < 5%). Therefore, ligand evolution may be a more promising approach to improve the meta-selectivity.
Scheme 1
Scheme 1. Influence of diboron reagents on regioselectivity. a data from Ref. [23]. b Conversion of 1a. c At 50 ℃, m:p = 0.5:1, 90% conversion of 1a.Figure 4
Figure 4. DFT-computed geometries and Gibbs free energies for the C—H activation transition states of 1a-Me using DB-1 and L6 as ligand, computed with SMD(p-xylene)-ωB97X-D/6–311+G(d, p)[SDD for Ir]//M06/6–31G(d)[SDD for Ir]. For improved clarity, most hydrogen atoms are omitted. All distances are in Å.With the computational and experimental results in hand, we then focus on how to utilize our findings to improve the previously reported suboptimal results. Kanai and coworkers reported the borylation of heteroaromatic amide 2a using L6 as ligand, which gave a moderate regioselectivity (m: p = 6.4:1, Scheme 2). We located the meta- and para-C-H activation transition states of this reaction (Fig. 5A). Among them, TS3-meta and TS3-para-ND were responsible for generating the meta- and para-borylated product. The preference of meta- over para-C-H borylation is 1.2 kcal/mol, in accordance with the reported selectivity (m: p = 6.4:1). Due to its more electron-withdrawing aromatic ring, amide 2a is a poorer hydrogen bond acceptor than benzamide. As a result, its carbonyl group has less advantage against the O atom of Bpin to interact with the urea moiety of ligand, further favoring the non-directed pathway. We envisioned that, through increasing the steric repulsion between the urea moiety of ligand and the Bpin group bound to iridium center, the non-directed pathway would be inhibited. Therefore, we changed the cyclohexyl to a much bulkier 2-adamantyl and designed ligand L7. Further DFT calculations showed that the steric repulsion between the urea moiety and Bpin increased the energies of non-directed transiton states TS4-meta-ND and TS4-para-ND by 3–5 kcal/mol, impeding the non-directed pathway effectively (Fig. 5B). Meanwhile, as the free energy of directed transition state TS4-meta is much lower than those of three other transition states, a higher meta-selectivity is anticipated. When the 2-adamantyl-bearing ligand L7 was synthesized and used in the C—H borylation reaction, it indeed gave a higher selectivity (m: p = 12:1, Scheme 2). This realized the control of meta-selectivity through inhibiting the non-directed para-C-H borylation by inceasing the size of substituent on the urea moiety.
Scheme 2
Scheme 2. Steric effect of ligands on the meta-selective borylation of 2a. * data from Ref. [23].Figure 5
Figure 5. DFT-computed geometries and Gibbs free energies for the C—H activation transition states of 2a using B2pin2 and L6 or L7 as ligand, computed with SMD(p-xylene)-ωB97X-D/6–311+G(d, p)[SDD for Ir]//M06/6–31G(d)[SDD for Ir]. For improved clarity, most hydrogen atoms are omitted. All distances are in Å.Encouraged by this positive result, ligand L7 was used for the substrates with poor meta-selectivities in Kanai's work (Scheme 3) [23]. In the borylation of benzamide 3a, a 6.9:1 ratio of meta to para was reported by using L6, while an enhanced selectivity (m: p = 12:1) was obtained by using our ligand L7. Previously, N-methylisoindolin-1-one 4a showed a meta-regioselectivity of only 3.3:1. This selectivity was as high as 15:1 when L7 was used as ligand. Finally, for aryl phosphonate 5a, the para-borylated product was generated as major product (m: p = 0.5:1) using L6, which was similar to the result using dtbpy, a typical ligand exhibiting intrinsic selectivity of substrates in the Ir-catalyzed borylation [59,60]. Significantly, our ligand L7 enabled the meta-selective borylation of 5a, delivering the desired meta-product with a high selectivity (m: p = 12:1).
Scheme 3
Scheme 3. The comparison of meta-selectivity in the C—H borylation of arenes by using ligands L6 and L7. * data from Ref. [23].In summary, this work combines calculations and experiments to elucidate the origin of regioselectivity in hydrogen bond interaction directed meta-selective C—H borylation of benzamides. The non-directed pathway, in which the ligand recognizes Bpin instead of substrate, is disclosed. We find that the non-directed para-borylation pathway plays an important role in controlling regioselectivity through competing with the meta-borylation pathways. This study demonstrates that the electronic and steric effects of ligand and the size of diboron reagent affect the regioselectivity by tuning the barriers of directed and non-directed pathways. Based on these findings, the ligand is improved and applied into some previously reported unsuccessful arenes, giving excellent performance. This work is a cornerstone of our ligand design in the remote C—H activation, which leads to the tunable meta- and para-selective C—H borylation of arenes bearing a variety of multi-transformable directing groups [61,62].
Declaration of competing interest
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.
Acknowledgments
We are grateful for financial support from the Fundamental Research Funds for the Central Universities (Nos. 020514380253, 020514380277), the Natural Science Foundation of Jiangsu Province (No. BK20211555), and the Jiangsu Innovation & Entrepreneurship Talents Plan. We thank the High Performance Computing Center (HPCC) of Nanjing University for doing the numerical calculations in this paper on its blade cluster system.
-
-
[1]
G. Meng, N.Y.S. Lam, E.L. Lucas, et al., J. Am. Chem. Soc. 142 (2020) 10571–10591. doi: 10.1021/jacs.0c04074
-
[2]
T.W. Lyons, M.S. Sanford, Chem. Rev. 110 (2010) 1147–1169. doi: 10.1021/cr900184e
-
[3]
R. Bisht, C. Haldar, M.M.M. Hassan, et al., Chem. Soc. Rev. 51 (2022) 5042–5100. doi: 10.1039/d1cs01012c
-
[4]
V. Snieckus, Chem. Rev. 90 (1990) 879–933. doi: 10.1021/cr00104a001
-
[5]
S.R. Neufeldt, M.S. Sanford, Acc. Chem. Res. 45 (2012) 936–946. doi: 10.1021/ar300014f
-
[6]
Z. Chen, B. Wang, J. Zhang, et al., Org. Chem. Front. 2 (2015) 1107–1295. doi: 10.1039/C5QO00004A
-
[7]
Z.R. Wang, P.P. Xie, Y.Z. Xia, Chin. Chem. Lett. 29 (2018) 47–53. doi: 10.11648/j.ajmsp.20180303.12
-
[8]
M.T. Mihai, G.R. Genov, R.J. Phipps, Chem. Soc. Rev. 47 (2018) 149–171. doi: 10.1039/C7CS00637C
-
[9]
U. Dutta, S. Maiti, T. Bhattacharya, D. Maiti, Science 372 (2021) eabd5992. doi: 10.1126/science.abd5992
-
[10]
J.A. Leitch, C.G. Frost, Chem. Soc. Rev. 46 (2017) 7145–7153. doi: 10.1039/C7CS00496F
-
[11]
N. Hofmann, L. Ackermann, J. Am. Chem. Soc. 135 (2013) 5877–5884. doi: 10.1021/ja401466y
-
[12]
C.J. Teskey, A.Y.W. Lui, M.F. Greaney, Angew. Chem. Int. Ed. 54 (2015) 11677–11680. doi: 10.1002/anie.201504390
-
[13]
C.C. Yuan, L. Zhu, C.P. Chen, et al., Nat. Commun. 9 (2018) 1189–1198. doi: 10.1038/s41467-018-03341-6
-
[14]
A. Dey, S.K. Sinha, T.K. Achar, D. Maiti, Angew. Chem. Int. Ed. 58 (2019) 10820–10843. doi: 10.1002/anie.201812116
-
[15]
D. Leow, G. Li, T.S. Mei, J.Q. Yu, Nature 486 (2012) 518–522. doi: 10.1038/nature11158
-
[16]
R.Y. Tang, G. Li, J.Q. Yu, Nature 507 (2014) 215–220. doi: 10.1038/nature12963
-
[17]
S. Lee, H. Lee, K.L. Tan, J. Am. Chem. Soc. 135 (2013) 18778–18781. doi: 10.1021/ja4107034
-
[18]
Y.F. Yang, G.J. Cheng, P. Liu, et al., J. Am. Chem. Soc. 136 (2014) 344–355. doi: 10.1021/ja410485g
-
[19]
S. Bag, T. Patra, A. Modak, et al., J. Am. Chem. Soc. 137 (2015) 11888–11891. doi: 10.1021/jacs.5b06793
-
[20]
M. Bera, A. Maji, S.K. Sahoo, D. Maiti, Angew. Chem. Int. Ed. 54 (2015) 8515–8519. doi: 10.1002/anie.201503112
-
[21]
Z.P. Zhang, K. Tanaka, J.Q. Yu, Nature 543 (2017) 538–542. doi: 10.1038/nature21418
-
[22]
M. Liu, L.J. Li, J. Zhang, et al., Chin. Chem. Lett. 31 (2020) 1301–1304. doi: 10.1016/j.cclet.2019.09.057
-
[23]
Y. Kuninobu, H. Ida, M. Nishi, M. Kanai, Nat. Chem. 7 (2015) 712–717. doi: 10.1038/nchem.2322
-
[24]
R. Bisht, B. Chattopadhyay, J. Am. Chem. Soc. 138 (2016) 84–87. doi: 10.1021/jacs.5b11683
-
[25]
H.J. Davis, M.T. Mihai, R.J. Phipps, J. Am. Chem. Soc. 138 (2016) 12759–12762. doi: 10.1021/jacs.6b08164
-
[26]
H.J. Davis, G.R. Genov, R.J. Phipps, Angew. Chem. Int. Ed. 56 (2017) 13351–13355. doi: 10.1002/anie.201708967
-
[27]
M.E. Hoque, R. Bisht, C. Haldar, B. Chattopadhyay, J. Am. Chem. Soc. 139 (2017) 7745–7748. doi: 10.1021/jacs.7b04490
-
[28]
L. Yang, K. Semba, Y. Nakao, Angew. Chem. Int. Ed. 56 (2017) 4853–4857. doi: 10.1002/anie.201701238
-
[29]
R. Bisht, M.E. Hoque, B. Chattopadhyay, Angew. Chem. Int. Ed. 57 (2018) 15762–15766. doi: 10.1002/anie.201809929
-
[30]
L.C. Yang, N. Uemura, Y. Nakao, J. Am. Chem. Soc. 141 (2019) 7972–7979. doi: 10.1021/jacs.9b03138
-
[31]
G. Liao, Y.J. Wu, B.F. Shi, Acta Chim. Sin. 78 (2020) 289–298. doi: 10.6023/a20020027
-
[32]
J. Trouvé, P. Zardi, S. Al-Shehimy, T. Roisnel, R. Gramage-Doria, Angew. Chem. Int. Ed. 60 (2021) 18006–18013. doi: 10.1002/anie.202101997
-
[33]
S. Lu, T. Zheng, J. Ma, et al., Angew. Chem. Int. Ed. 61 (2022) e202201285.
-
[34]
X.C. Wang, W. Gong, L.Z. Fang, et al., Nature 519 (2015) 334–338. doi: 10.1038/nature14214
-
[35]
H. Shi, A.N. Herron, Y. Shao, et al., Nature 558 (2018) 581–585. doi: 10.1038/s41586-018-0220-1
-
[36]
J.C. Wang, G.B. Dong, Chem. Rev. 119 (2019) 7478–7528. doi: 10.1021/acs.chemrev.9b00079
-
[37]
J. Cornella, M. Righi, I. Larrosa, Angew. Chem. Int. Ed. 50 (2011) 9429–9432. doi: 10.1002/anie.201103720
-
[38]
R.J. Phipps, M.J. Gaunt, Science 323 (2009) 1593–1597. doi: 10.1126/science.1169975
-
[39]
Y. Nakao, Y. Yamada, N. Kashihara, T. Hiyama, J. Am. Chem. Soc. 132 (2010) 13666–13668. doi: 10.1021/ja106514b
-
[40]
C.L. Ciana, R.J. Phipps, J.R. Brandt, et al., Angew. Chem. Int. Ed. 50 (2011) 458–462. doi: 10.1002/anie.201004703
-
[41]
N.A. Romero, K.A. Margrey, N.E. Tay, D.A. Nicewicz, Science 349 (2015) 1326–1330. doi: 10.1126/science.aac9895
-
[42]
Y. Saito, Y. Segawa, K. Itami, J. Am. Chem. Soc. 137 (2015) 5193–5198. doi: 10.1021/jacs.5b02052
-
[43]
F. Berger, M.B. Plutschack, J. Riegger, et al., Nature 567 (2019) 223–228. doi: 10.1038/s41586-019-0982-0
-
[44]
C. Haldar, M.E. Hogue, R. Bisht, B. Chattopadhyay, Tetrahedron Lett. 59 (2018) 1269–1277. doi: 10.1016/j.tetlet.2018.01.098
-
[45]
Y. Kuroda, Y. Nakao, Chem. Lett. 48 (2019) 1092–1100. doi: 10.1246/cl.190372
-
[46]
X.L. Zou, S.M. Xu, Chin. J. Org. Chem. 41 (2021) 2610–2620. doi: 10.6023/cjoc202103020
-
[47]
S. Pandit, S. Maiti, D. Maiti, Org. Chem. Front. 8 (2021) 4349–4358. doi: 10.1039/d1qo00452b
-
[48]
I.A.I. Mkhalid, J.H. Barnard, T.B. Marder, et al., Chem. Rev. 110 (2010) 890–931. doi: 10.1021/cr900206p
-
[49]
J.W.B. Fyfe, A.J.B. Watson, Chem 3 (2017) 31–55. doi: 10.1016/j.chempr.2017.05.008
-
[50]
J. Hu, M. Ferger, Z. Shi, T.B. Marder, Chem. Soc. Rev. 50 (2021) 13129–13188. doi: 10.1039/d0cs00843e
-
[51]
J. Lv, B. Zhao, Y. Han, et al., Chin. Chem. Lett. 32 (2021) 691–694. doi: 10.1016/j.cclet.2020.06.028
-
[52]
H. Tamura, H. Yamazaki, H. Sato, S. Sakaki, J. Am. Chem. Soc. 125 (2003) 16114–16126. doi: 10.1021/ja0302937
-
[53]
T.M. Boller, J.M. Murphy, M. Hapke, et al., J. Am. Chem. Soc. 127 (2005) 14263–14278. doi: 10.1021/ja053433g
-
[54]
A.G. Green, P. Liu, C.A. Merlic, K.N. Houk, J. Am. Chem. Soc. 136 (2014) 4575–4583. doi: 10.1021/ja411699u
-
[55]
X. Zou, H. Zhao, Y. Li, et al., J. Am. Chem. Soc. 141 (2019) 5334–5342. doi: 10.1021/jacs.8b13756
-
[56]
A. Unnikrishnan, R.B. Sunoj, Chem. Sci. 10 (2019) 3826–3835. doi: 10.1039/C8SC05335A
-
[57]
G. Pupo, V. Gouverneur, J. Am. Chem. Soc. 144 (2022) 5200–5213. doi: 10.1021/jacs.2c00190
-
[58]
N. Oka, T. Yamada, H. Sajiki, et al., Org. Lett. 24 (2022) 3510–3514. doi: 10.1021/acs.orglett.2c01174
-
[59]
T. Ishiyama, J. Takagi, J.F. Hartwig, N. Miyaura, Angew. Chem. Int. Ed. 41 (2002) 3056–3058. doi: 10.1002/1521-3773(20020816)41:16<3056::AID-ANIE3056>3.0.CO;2-#
-
[60]
T. Ishiyama, J. Takagi, K. Ishida, et al., J. Am. Chem. Soc. 124 (2002) 390–391. doi: 10.1021/ja0173019
-
[61]
W. Chang, Y. Chen, S. Lu, et al., Chem 8 (2022) 1775–1788. doi: 10.1016/j.chempr.2022.04.025
-
[62]
Y. Wang, W. Chang, S. Qin, et al., Angew. Chem. Int. Ed. 61 (2022) e202206797.
-
[1]
-
Figure 3 DFT-computed geometries and Gibbs free energies for the C-H activation transition states of 1a-Me using B2pin2 and L5 as ligand, computed with SMD(p-xylene)-ωB97X-D/6–311+G(d, p)[SDD for Ir]//M06/6–31G(d)[SDD for Ir]. For improved clarity, most hydrogen atoms are omitted. All distances are in Å.
Scheme 1 Influence of diboron reagents on regioselectivity. a data from Ref. [23]. b Conversion of 1a. c At 50 ℃, m:p = 0.5:1, 90% conversion of 1a.
Figure 4 DFT-computed geometries and Gibbs free energies for the C—H activation transition states of 1a-Me using DB-1 and L6 as ligand, computed with SMD(p-xylene)-ωB97X-D/6–311+G(d, p)[SDD for Ir]//M06/6–31G(d)[SDD for Ir]. For improved clarity, most hydrogen atoms are omitted. All distances are in Å.
Scheme 2 Steric effect of ligands on the meta-selective borylation of 2a. * data from Ref. [23].
Figure 5 DFT-computed geometries and Gibbs free energies for the C—H activation transition states of 2a using B2pin2 and L6 or L7 as ligand, computed with SMD(p-xylene)-ωB97X-D/6–311+G(d, p)[SDD for Ir]//M06/6–31G(d)[SDD for Ir]. For improved clarity, most hydrogen atoms are omitted. All distances are in Å.
Scheme 3 The comparison of meta-selectivity in the C—H borylation of arenes by using ligands L6 and L7. * data from Ref. [23].
-

计量
- PDF下载量: 9
- 文章访问数: 1117
- HTML全文浏览量: 48