Yining
Ji
*,
Hongming
Li
,
Alan M.
Hyde
,
Qinghao
Chen
,
Kevin M.
Belyk
,
Katrina W.
Lexa
,
Jingjun
Yin
,
Edward C.
Sherer
,
R. Thomas
Williamson
,
Andrew
Brunskill
,
Sumei
Ren
,
Louis-Charles
Campeau
,
Ian W.
Davies
and
Rebecca T.
Ruck
Department of Process Research & Development, Merck & Co., Inc., Rahway, New Jersey 07065, USA. E-mail: yining.jichen@merck.com
First published on 19th January 2017
Significant catalyst loading reduction and increased reaction robustness have been achieved for a Pd-catalyzed asymmetric intramolecular C–N coupling through comprehensive mechanistic studies. Detailed kinetic, spectroscopic, and crystallographic analyses revealed that the mono-oxidation of the bis-phosphine ligand is critical for a successful transformation. 31P NMR studies provided an understanding of the inefficient activation of the Pd(OAc)2/(R,R)-QuinoxP* pre-catalyst to form the active bis-phosphine mono-oxide–Pd(0) catalyst with competitive formation of a less active (R,R)-QuinoxP*·PdBr2 complex. Based on these detailed mechanistic studies, a new series of bis-phosphine mono-oxides (BPMO)-ligated Pd(II) pre-catalysts have been rationally developed that allow for reliable and complete catalyst activation which should have general utility in academic and industrial settings.
We recently reported9 a novel, enantioselective hemiaminal synthesis via a palladium-catalyzed C–N coupling using chiral bis-phosphine ligands (Scheme 1) that has been used in the preparation of ruzasvir, a clinically relevant hepatitis C virus NS5A inhibitor.10
In this Pd-catalyzed process, we also observed that the in situ mono-oxidation of the bis-phosphine ligand (R,R)-QuinoxP* (abbreviated as L*) to access the active catalyst L*(O)Pd(0) B from the pre-catalyst L*Pd(OAc)2A was critical for effective catalysis. However, different ligands displayed diverse behaviors as the bis-phosphine versus their mono-oxide (Fig. 1). While (R,R)-QuinoxP* and the parent mono-oxide gave similar reaction conversions under screening conditions, (R,R)-Et-DuPhos proved to be inferior to its respective mono-oxide derivative and finally, fully reversed reactivity was observed for (S)-BINAP/(S)-BINAP(O) system. The inefficient in situ oxidation of (R,R)-Et-DuPhos to form its mono-oxide counterpart and a possible over-oxidation of BINAP(O) leading to the inactive bis-oxide5,8 during the reduction of palladium(II) likely account for the differential reaction performance.11
Fig. 1 Comparison of reaction conversion and enantioselectivity of different ligands with respect to their parent mono-oxidized form for the reaction in Scheme 1. Note: in screening mode, the pre-catalyst is formed in situ from a mixture of Pd(OAc)2 and the corresponding bis-phosphine ligand. |
Experimental nuances, such as “age” time3e during the pre-catalyst preparation, order of addition of the reagents, ligand/metal ratio7a,8b and the reaction conditions, in general, could all affect the efficiency of the in situ oxidation and significantly impact both reaction performance and outcome. Therefore, interpreting the actual intrinsic reactivity of these ligands based solely on the screening results is potentially misleading. In theory, a detailed investigation of the efficiency of pre-catalyst activation and mechanism would be required prior to evaluation of any results, which would undermine the inherent value of high throughput experimentation (HTE).
In this article, we report a systematic kinetic and spectroscopic investigation into the asymmetric C–N coupling reaction shown in Scheme 1. We gained insight into the efficiency of pre-catalyst activation and potential deactivated catalyst states. This deeper understanding enabled the successful development of more effective palladium–BPMO pre-catalysts and an optimized reaction with significantly reduced catalyst loading and increased reaction robustness. Based on these studies, we anticipate that use of these pre-catalysts will also enable a reliable assessment of the innate reactivity of BPMOs and facilitate widespread employment of this under-utilized class of ligands.
Fig. 2 Temporal concentration profiles monitored by 1H NMR spectroscopy for the reaction of Scheme 1. [1]0 as noted; [A]0 = 3.25 mM; K3PO4 (5.5 equiv. based on [1]0 = 0.0650 M and 11.0 equiv. based on [1]0 = 0.0325 M). |
The current reaction system (Scheme 1) is heterogeneous due to the low solubility of K3PO4 in toluene. While we see rate acceleration with increasing amount of phosphate (see ESI†), the fact that the reaction exhibited positive first order kinetics in concentration of L*·Pd(OAc)2 precursor A (see ESI†) indicates that we are not operating under mass transfer limited conditions. The base appears to play two main roles in the system: (1) in combination with water, it facilitates reduction of Pd(II) to Pd(0) with concomitant phosphine oxidation13 and (2) it serves to deprotonate substrate 1 at the benzylic position allowing turnover of the catalytic cycle. The observed rate increase with more base present could be rationalized by either higher surface area and/or increased water content.14
Fig. 3 Temporal concentration profiles monitored by 1H NMR spectroscopy for the reaction of Scheme 1. [1]0 as noted; [A]0 = 3.25 mM; K3PO4 (5.5 equiv. based on [1]0 = 0.0650 M and 11.0 equiv. based on [1]0 = 0.0325 M). Product addition conditions: [1]0 = 0.0325 M, [2]0 = 0.030 M. |
In our previous work, we reported that the analogous substrate 3 (Scheme 3) gave the same enantiomer of the corresponding product when subjected to the identical reaction conditions.9 Treatment of ArBr 3 with L*(O) and Pd[P((o-tolyl)3)2] also afforded a single pair of characteristic 31P NMR resonances at δ 48.6 and 34.9 ppm that correspond to the phosphine oxide and phosphine moieties of the BPMO, respectively (see ESI†). Gratifyingly, the isolated complex D allowed the growth of X-ray quality single crystals by slow evaporation from toluene. The crystal structure obtained (Fig. 5a) corroborated the structure obtained from NMR spectroscopic and computational modeling analysis (Fig. 5b). In light of these results, we were confident to apply the combination of NMR and DFT modeling techniques to characterize complex C, which proved to possess an almost identical structure as complex D (Fig. 5c).
Fig. 5 (a) X-ray diffraction structure of complex D; (b) overlay between the NMR-DFT modeling (blue) and the X-ray structure (red) of complex D; (c) NMR-DFT modeling structure of complex C. |
These complexes feature the Pd atom in a square planar arrangement, with the arene and bromide in a trans-relationship to each other. The bis-phosphine mono-oxide is directly ligated to the metal center through the phosphorus atom on one side, leaving the phosphine-oxide oxygen 3.125 Å from the Pd center. The imine-type nitrogen is ligated to the metal center, forming a concave boat conformation and, consequently, forcing the nearby rings into proximity. This binding mode facilitates the aromatic group at the hemiaminal position moving to the axial position to avoid the steric repulsion from the bulky ligand. In both complexes C and D, the observed absolute hemiaminal stereochemistry was consistent with the desired enantiomer of the product 2 (Scheme 4).
Fig. 6 Temporal concentration profiles monitored by 1H NMR spectroscopy for the reaction of Scheme 1. [1]0 = 0.0650 M; [C]0 = 3.25 mM; [A]0 = 3.25 mM; 5.5 equiv. of K3PO4; the amounts of AcOH and PivOH as noted. |
The superior catalytic performance of complex C (in the presence of an appropriate amount of carboxylate) could be an indication of an inefficient activation of pre-catalyst A. To investigate this hypothesis, 31P NMR spectra for both catalyst systems were compared over the course of the two reactions. As shown in Fig. 7, for complex C, a clean 31P NMR spectra shows the two characteristic 31P resonances were maintained throughout the reaction course. Good superposition from the time-adjusted point for two experiments that differ in initial concentration of substrate but share the same catalyst concentration suggests that neither catalyst deactivation nor product inhibition processes are significant under these conditions (Fig. 8). The kinetic behavior of the catalyst appears to be robust, consistent with the catalyst remaining viable throughout the reaction. It is worth noting that, under these conditions, the reaction appeared to exhibit saturation kinetics in [1] with the initial rate of the two reactions only matching up to 20% conversion (see ESI†).
Fig. 7 202 MHz 31P NMR spectra of aliquots taken at noted times during the reaction of Scheme 1 (complex C used as the catalyst, 10 mol% AcOH added). All spectra were recorded by diluting an aliquot of the reaction mixture in toluene-d8. |
Fig. 8 Temporal concentration profiles monitored by 1H NMR spectroscopy for the reaction of Scheme 1 using 5 mol% of complex C as the catalyst. [1]0 as noted; [C]0 = 3.25 mM; 10 mol% PivOH; K3PO4 (5.5 equiv. based on [1]0 = 0.0650 M and 11.0 equiv. based on [1]0 = 0.0325 M). |
In contrast, we found that pre-catalyst A, L*·Pd(OAc)2, was only partially activated under standard reactions conditions. In the 31P NMR spectrum over the course of the reaction, we observed a minor species showing phosphorus signals at δ 48.6 and 34.9 ppm, corresponding to the oxidative addition complex C (Fig. 9). A major phosphorus resonance was also observed at 52.6 ppm, very close to the 52.3 ppm resonance in pre-catalyst A. This new species showed very high stability under the reaction conditions and was even present at the end of the reaction (t = 180 min).17
Fig. 9 202 MHz 31P NMR spectra of aliquots taken at noted times during the reaction of Scheme 1. All spectra were recorded by diluting an aliquot of the reaction mixture in toluene-d8. |
Because of the ambiguous 31P NMR spectra, we turned to the 1H NMR (Fig. 10c) to deconvolute the species. We identified the doublet at 1.16 ppm belonging to the tert-butyl group of the minor catalytic species, complex C. However, the tert-butyl resonances of the major palladium species showed different chemical shift (0.9 ppm) vs. the pre-catalyst A (1.11 ppm). Furthermore, the doublet corresponding to the methyl group of the ligand moiety at 1.94 ppm in the pre-catalyst A was not present at the same chemical shift in the spectrum of the reaction mixture.
Fig. 10 500 MHz 1H NMR spectra of (a) pre-catalyst A; (b) complex C, (c) aliquot taken at t = 15 min of reaction at Scheme 1. All spectra were recorded in toluene-d8. |
The high stability allowed us to isolate this persistent species, and the structure was unambiguously established by X-ray diffraction analysis as L*PdBr2E,18 originating from an anion exchange between acetate and bromide. The newly discovered complex E was found to match the major species present in the reaction mixture in both 31P (52.5 vs. 52.6 ppm) and 1H NMR spectra (Fig. 11).
Fig. 11 500 MHz 1H NMR spectra of (a) L*PdBr2, E; (b) aliquot taken at t = 15 min of reaction at Scheme 1. All spectra were recorded in toluene-d8. |
Based on quantitative 1H and 31P NMR analyses of a standard reaction at ca. 20% conversion (see ESI†), the observed species account for >90% of the total phosphorus-containing Pd species and, of that, only 27% is complex C. The prevalence of complex E likely accounts for the lack of overlay observed between the experiments that were designed to test the catalyst robustness (vide supra).
The lower reaction rate obtained with pre-catalyst A compared to oxidative addition complex C suggests that the formation of complex E is detrimental to the reaction. This hypothesis was confirmed by comparing the kinetic profiles obtained using pre-catalyst E and A (Fig. 12, green diamonds vs. blue circles, respectively). The lower reaction rate provided by complex E is attributable to its inferior reactivity toward activation to form the active BPMO–Pd(0) complex than pre-catalyst A. Notably, when using complex E as the catalyst, the addition of an excess of CsOAc (1 equiv. with respect to 1) to drive the equilibrium from E to A increased the reaction rate four-fold (green triangles). 31P NMR analysis at 15 min indicated a higher concentration of the oxidative addition complex C, likely associated with a shift of the equilibrium toward pre-catalyst A by CsOAc (Fig. 13).
Fig. 12 Comparison of temporal concentration profiles monitored by 1H NMR spectroscopy for the reaction of Scheme 1 using different catalyst systems, complex Aversus complex E (with/without CsOAc). [1]0 = 0.065 M; [A]0 = 3.25 mM; [E]0 = 3.25 mM; 5.5 equiv. of K3PO4; CsOAc as noted. |
Fig. 13 202 MHz 31P NMR spectra of aliquots taken at t = 15 min for the reaction of Scheme 1: (a) 5 mol% complex E; 1 equivalent of CsOAc (b) 5 mol% complex E; (c) 5 mol% of pre-catalyst A. All spectra were recorded by diluting an aliquot of the reaction mixture in toluene-d8. |
31P NMR analysis of the reaction mixture under these conditions showed that species C is not present, while ligand L*(O) and a new species identified by mass spectroscopy as (L*)2Pd(0) (F) were formed (Fig. 14a). Free L* readily displaces L*(O) to produce inactive species F (Scheme 5) that could also be formed upon simple treatment of L* and Pd[P(o-tolyl)3]2 (Fig. 14b).
Fig. 14 202 MHz 31P NMR spectra of (a) aliquot taken at t = 180 min of Scheme 1 in the presence of excess 5 mol% of L*; (b) Pd[P(o-tolyl)3]2:L* (1:1). All spectra were recorded in toluene-d8. |
Fig. 16 Temporal enantioselectivity of the remaining substrate 1versus reaction conversion for the reaction in Scheme 1 using pre-catalyst A or complex C as catalyst. |
Computational studies21 of the formation of complex C were carried out using M06/6-31+G(d) with SDD for Pd and Br in Gaussian 09. Solvation in toluene was accounted for with the SMD solvent model (Fig. 17). The ΔΔG‡ (ca. 1.5 kcal mol−1) favors the formation of the desired complex C (S configuration at the hemiaminal position) over the undesired diastereomer (R-C). Complex C is also favored thermodynamically over the undesired diastereomer R-C by 3.67 kcal mol−1. The R-C substrate must undergo re-puckering at the benzylic methylene in order to alleviate steric clashes during formation of the oxidative addition complex. As a result, the TS-R is a later transition state than TS-S, and can be described by longer distances and wider angles between the reacting atoms. The Pd–C bond being formed is 2.14 Å in the TS-R and 2.05 Å in TS-S, and the Pd–Br bond being formed is 3.58 Å in TS-R and only 3.07 Å in TS-S. Also, the Br–C–Pd angle is sharper in the TS-S than the TS-R, with angles of 90.9° and 109.2° respectively.
Based on the high diastereoselectivity observed experimentally, correspondingly high product ee (>98%) would be expected. However, 91% ee is obtained under catalytic conditions. Further, reacting a single diastereomer of complex C with 5.5 equivalents of K3PO4 at the standard reaction temperature (50 °C, Scheme 8), gave a slight product ee erosion (from 95% ee22 to 91% ee at ca. 10% and 98% conversions, respectively) (Fig. 18). Notably, a very similar ee reduction was also observed under standard catalytic conditions (see ESI†). Product 2 was shown to be configurationally stable to base at the reaction temperature (see ESI†). Therefore, the ee erosion should derive from some downstream catalytic step following the formation of complex C.
Fig. 18 Temporal reaction ee monitored over time for the stoichiometric reaction of complex C with 5.5 equivalents of K3PO4. |
Two plausible intermediates could be derived from the deprotonation of complex C (Scheme 10), namely: (a) complex H with retention of the stereochemistry at the hemiaminal position, from which stereochemical integrity is maintained and/or; (b) complex G in which the stereochemical information at the hemiaminal center is lost through ring opening. Along the deprotonation pathway, DFT calculations indicated a stable intermediate, complex I (Scheme 9), could be formed where C–H bond cleavage has occurred but the bromide remains bound. Calculations of interconversion through the open form G indicated both C and R-C were accessible, with I being thermodynamically preferred relative to R-I when closing from the open form G (Table 1). Furthermore, our simulations showed that while complex G is stable and energetically accessible, interconversion between the C and R-C form through the zwitterion J is less likely due to the high thermodynamic penalty (Table 1) and dramatic structural rearrangement required.
Complex interconversion | ΔΔGr (kcal mol−1) gas phase | ΔΔGr (kcal mol−1) toluene |
---|---|---|
H → J | +20.85 | +17.27 |
I → G | +10.27 | +11.60 |
I → R-I | +3.71 | +4.95 |
I(K+) → G(K+) | −1.43 | −6.49 |
In order to examine if complex G could be a competent intermediate, we treated the open form 4 with Pd[P(o-tolyl)3]2 and L*(O) (Scheme 10). The reaction provided the desired product in 90% ee. Presumably, the transformation occurred through the intermediacy of complex G. To account for the high observed enantioselectivity, the ring closure of complex G leading to H to allow the reductive elimination to take place should occur with high diastereoselectivity, presumably aided by the stereochemical information relay from the chiral ligand. Notably, since no proton abstraction was required, the reaction occurred instantaneously at room temperature.
Based on the above results and in consonance with the slight reaction ee change as the reaction progresses, two different enantio-determining steps, namely, the diastereoselective formation of complex C and the ring closure of complex G to form complex H, appear operational.
While the formation of C is expected to be highly diastereoselective, the possible competing side reactions representing the formation of G, either directly from C or through trapping of B by the open form 4, likely explain the observed lower ee. Our computational studies further supported the feasibility of both pathways.
Scheme 11 Proposed mechanism for the reaction of Scheme 1 highlighting critical species and catalyst poisons. |
Active catalyst B could be formed from a mixture of free ligand L*, Pd(OAc)2, water and base, but the complex proved unstable and quickly decomposed. Alternatively, trapping with ArBr 1 was complicated by cyclization of complex C in the presence of base preventing its isolation. Based on this knowledge, we envisioned that we could use simple aryl halides to trap active catalyst B as stable oxidative addition complexes that would only release active L*(O)–Pd(0) catalyst B after a promoted reductive elimination.
We identified a series of aryl halides with relatively high reactivity toward oxidative addition,23 including iodobenzene, bromobenzene and p-bromobenzonitrile, that cleanly afforded the desired oxidative complexes K24 (Scheme 12). These complexes were shown to be air-stable and could be purified by silica gel chromatography. Chlorobenzene failed to undergo oxidative addition. The introduction of an electron-withdrawing substituent (i.e. nitrile) increased the reactivity with respect to bromobenzene and afforded the corresponding oxidative addition complex (K-3) in an increased 86% yield on 10 g scale. The installation of a para-nitrile group also allowed for direct isolation by crystallization from toluene, a highly desired feature from a practical perspective.
With these newly synthesized oxidative addition complexes in hand, we evaluated their performance as pre-catalysts in the standard C–N coupling reaction (Scheme 13). As expected, utilizing iodobenzene complex K-1 as catalyst, in the absence of additive, failed to produce any product. Using triethylsilane as a reductant, 15% conversion was obtained at 1.5 mol% catalyst loading, providing proof of concept for our activation hypothesis. Higher conversion (32%) was observed on introduction of bis-pinacolborane, and employing of phenylboronic acid proved optimal, giving 98% conversion. The observation of biphenyl formation (1.5 mol%, in accordance with the catalyst loading) in the latter case confirmed activation proceeded via Suzuki–Miyaura coupling. In accord with our hypothesis, the use of these pre-catalysts did not alter the enantioselectivity of the reaction, and we observed the expected oxidative addition complex C during the course of the reaction.
Unfortunately, the use of phenylboronic acid as a reductant led to low reaction conversion at lower catalyst loading (11% conversion at 0.3 mol%). Therefore, we screened a range of arylboronic acids with the p-bromobenzonitrile complex K-3 to identify a superior additive (see ESI†). Both steric and electronic properties of the boronic acid were shown to play critical roles in the efficiency of the activation. o-Tolylboronic acids bearing electron-donating groups at the para-position to increase nucleophilicity of the boronic acids25 were found to be optimal. Utilizing 0.3 mol% of the pre-catalyst K-3, 84% conversion to product was achieved after 28 h in the presence of 4-methoxy-2-tolylboronic acid.
The use of these new series of Pd(II) pre-catalysts activated with addition of arylboronic acids provides the highest catalytic performance, delivering maximal on-cycle palladium in the active catalyst form L*(O)–Pd(0) B. Notably, as was expected from our previous kinetic data, PivOH addition was also required to achieve the desired reactivity. Using pre-catalyst K-3 and 4-methoxy-2-tolylboronic acid as additive, 100% conversion to the desired product was achieved with as low as 0.5 mol% catalyst loading, a 10-fold reduction in catalyst loading while maintaining the required high product enantioselectivity (Scheme 14).
Footnote |
† Electronic supplementary information (ESI) available: Detailed experimental procedures and characterization data for all new compounds. CCDC 1512973–1512975. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c6sc05472b |
This journal is © The Royal Society of Chemistry 2017 |