Nnamdi
Akporji
a,
Ruchita R.
Thakore
a,
Margery
Cortes-Clerget
b,
Joel
Andersen
c,
Evan
Landstrom
a,
Donald H.
Aue
*a,
Fabrice
Gallou
b and
Bruce H.
Lipshutz
*a
aDepartment of Chemistry & Biochemistry, University of California Santa Barbara, Santa Barbara, CA 93106, USA. E-mail: lipshutz@chem.ucsb.edu
bNovartis Pharma AG, Basel, Switzerland
cDepartment of Chemistry, University of Cincinnati, Cincinnati, OH 45221, USA
First published on 7th May 2020
A new biaryl phosphine-containing ligand from an active palladium catalyst for ppm level Suzuki–Miyaura couplings, enabled by an aqueous micellar reaction medium. A wide array of functionalized substrates including aryl/heteroaryl bromides are amenable, as are, notably, chlorides. The catalytic system is both general and highly effective at low palladium loadings (1000–2500 ppm or 0.10–0.25 mol%). Density functional theory calculations suggest that greater steric congestion in N2Phos induces increased steric crowding around the Pd center, helping to destabilize the 2:1 ligand–Pd(0) complex more for N2Phos than for EvanPhos (and less bulky ligands), and thereby favoring formation of the 1:1 ligand–Pdo complex that is more reactive in oxidative addition to aryl chlorides.
Scheme 1 HandaPhos, EvanPhos, and next generation ligand, N2Phos (L1), for SMCs (top). Ligands prepared and screened based on the EvanPhos biaryl skeleton (bottom). |
Typically, ligand design for Pd-catalyzed cross coupling reactions takes, in large measure, steric and electronic factors into consideration, as they can weigh heavily on the activity of the active palladium catalyst.22–25 While these factors are surely at play as well under micellar catalysis conditions, the switch from a traditional organic solvent medium to water presents new rules that must also be taken into account (e.g., ligand lipophilicity).26 Thus, further derivatization of the biaryl skeletal of EvanPhos (Scheme 1), with the goal of increasing the activity of the resulting ligated Pd complex, focused on modifications ranging from increasing steric bulk around the phosphine moiety by altering the alkoxy residues (L2–L4) to enhancing steric interactions based on additional substitution on the naphthyl ring. Results from the former changes showed no significant rate increase, while replacement of an alkoxy group at C-2 in the naphthyl ring showed significant improvements in reaction rates relative to those observed upon chelation of Pd with EvanPhos. Additional steric crowding by placement of two dialkylamino groups at the C-2 and C-7 locations afforded a catalyst, N2Phos, with the highest activity. The final ratio of N2Phos (L1) to Pd(OAc)2 of 1.8:1 led to an active pre-catalyst capable of mediating SMCs at 1000 ppm of Pd under mild, aqueous conditions.
Initial screening of the catalyst derived from N2Phos complexed with Pd(OAc)2 (1.8:1) versus that formed using EvanPhos under otherwise identical micellar conditions showed a significant increase in reactivity for the former in all substrates tested (Scheme 3). At 45 °C, the new ligand system led to good isolated yields of biaryl products, while the extent of conversion, and therefore, isolated yields of the same reactions, using EvanPhos was ca. 25% lower over the same time period. EvanPhos, however, showed comparable yields to N2Phos when given longer reaction times (16–24 h) for the substrates tested. One notable advantage of the N2Phos/Pd(OAc)2-derived catalyst is that pre-activation using commercially available DIBAL in toluene is no longer needed, as is the case with the corresponding EvanPhos-derived catalyst.15
The same pre-catalyst combination of N2Phos/Pd(OAc)2 could also be applied to SMCs at the 1000 ppm level, run under traditional conditions involving an organic solvent such as dioxane28 (Scheme 4). Relative to the corresponding reaction run under micellar catalysis conditions, these tended to require longer reaction times to reach completion (e.g., formation of biaryl 2; 8 h vs. 12 h; see Scheme 3vs.Scheme 4). When compared to commercially available SPhos,29a tri-t-butylphosphine (P(t-Bu)3), and tri(o-tolyl)phosphine (P(o-tol)3), N2Phos shows similar reactivity in dioxane for the substrates tested. Hence, under such conditions in organic media, there appears to be no real benefit to using this ligand over those currently readily available.
In water, however, the choice of a 2 wt% aqueous solution of designer surfactant TPGS-750-M21 was based on prior efforts that showed it to be an especially enabling and recyclable medium for several types of reactions, such as peptide couplings,32,33 SNAr,34 bio-catalytic processes,35 as well as a wide variety of metal-catalyzed cross coupling reactions. In the case of SMCs, high reactivity can be achieved using this new catalytic system (L1/Pd(OAc)2), with inclusion of K3PO4·H2O (1.5 equiv). This salt exhibited superior performance compared to other commonly employed bases, such as Et3N or K2CO3. Toluene was selected as co-solvent (10% by volume),36 given the poor solubility of some of the chosen substrates. Reactions run in the presence of this co-solvent appeared as nicely stirring emulsions. Broad functional group tolerance associated with either reaction partner is apparent, leading to good yields of coupled products. Hence, an assortment of challenging bromide partners can be coupled using this system at a loading of only 1000 ppm Pd (Scheme 5) including those bearing an ester (12), aldehyde (4), nitro (3, 8), lactam (13), and sulfonamide (7) residues. Heteroaromatic bromides and boronic acids are also amenable, such as those containing pyridine, pyrimidine, benzofuran, indoles, furans, and benzothiophenes. Not surprisingly, aryl iodides are also converted quickly to the corresponding biaryls in excellent yields. Alternative boron-containing partners beyond boronic acids, including a Bpin derivative (leading to product 10), a potassium trifluoroborate30 (affording product 15), and a MIDA boronate31 (giving product 16), could be used as well.
Although EvanPhos was designed as a quickly-synthesized, alternative ligand to HandaPhos for Pd-catalyzed SMCs, it has shown limited, if any, applicability to aryl chlorides, especially at the ppm level of precious metal.14,15 It is known to exist as a four-coordinate complex: (EvanPhos)2PdX2, based on both an X-ray crystal structure and computational work.15 However, N2Phos-complexed palladium can be used for couplings with aryl chlorides regardless of electronic influences on the ring. This may be a result of the differential modes of complexation, where EvanPhos can spatially form 2:1 complexes with Pd. As shown, however, by modeling and calculations (vide infra), the far greater size of N2Phos makes 2:1 complexation significantly less favorable. Thus, aryl chloride precursors reflecting both electron-rich, electron-poor, and neutral educts afford product biaryls 20, 22, and 23, respectively, although an increase in Pd loading to 2500 ppm (0.25 mol%) was needed (Scheme 6, top). Comparisons in each of these three cases were also made with EvanPhos, SPhos, P(t-Bu)3, and P(o-tol)3, the results from which indicate that N2Phos affords a catalyst that is equal to or better than these others in terms of reactivity at this loading of Pd under micellar catalysis conditions. The difference in conversion of aryl chlorides is quite noticeable, leading, e.g., to biaryl 20. Thus, while N2Phos gave 88% of product 20 after 8 hours, EvanPhos led to only 7% conversion after 8 h, and only 9% after 16 h (96% with N2Phos), after which time there was both unreacted aryl chloride present, as well as the product of protodeborylation.
The low yields with aryl chlorides appears, then, to be the result of a slow oxidative addition with EvanPhos compared to N2Phos. How these two ligands differ in their behavior and the possible role of steric effects is explored computationally (vide infra).
A similar outcome was noted in the case of the ester of indomethacin (Scheme 6, bottom), where a better result was observed in comparisons with catalysts derived from either SPhos or BI-DIME.24,37 Moreover, that the results between N2Phos and HandaPhos were within experimental error further suggests that the derived catalysts in each case are roughly comparable in activity imparted to their 1:1 complexes with Pd. The biaryls prepared and isolated from a more general study on the scope of these couplings with various aryl/heteroaryl chlorides are illustrated in Scheme 7. Several heterocyclic arrays, such as those in products 24–27 were formed under the same coupling conditions used for analogous bromides, although somewhat increased reaction times and catalyst loadings were needed to achieve high levels of conversion, and hence, good isolated yields.
Modified known targets that coupled smoothly included fenofibrate (leading to 28), used to treat high cholesterol, and previously mentioned indomethacin (giving products 29a-b), an anti-inflammatory drug. Very highly functionalized and especially challenging chloroarenes also participated, including rivaroxaban (a medication used to prevent blood clots) affording product 30, and glibenclamide (used to treat type 2 diabetes) leading to biaryl 31 (Scheme 8). Palladium loadings of 5000 ppm (0.50 mol%) together with a reaction temperature of 60 °C were required, along with water-miscible DMSO (10%) as co-solvent. Direct comparisons with the most commonly used, state-of-the-art Pd/ligand combination (i.e., SPhos),29b clearly indicated that under otherwise identical micellar conditions and ppm levels of Pd, (N2Phos)Pd gave far better results.
Scheme 8 Comparison of N2Phos and SPhos with highly functionalized aryl/heteroaryl chlorides. Isolated yields. |
Recycling of the 2 wt% aqueous surfactant medium could be readily achieved as a means of reducing the amount of organic waste produced, as measured by E factors (Scheme 9).38,39 Each recycle was carried out starting with an in-flask extraction of the aqueous reaction with ethyl acetate, followed by introduction of additional catalyst (i.e., an additional 1000 ppm Pd) into the reaction medium. Based on this representative biaryl coupling, an E factor of only 4.7 is indicative of the limited amounts of solvent waste generated. This value is far below those typically seen based on usage of organic solvents as both reaction medium and for extraction purposes.
Levels of residual Pd found in the products using ppm loadings of Pd catalysts enabled by aqueous nanomicellar technology are typically below the 10 ppm FDA allowed limit. This is yet another major advantage associated with this chemistry in water which avoids the time and expense of scavenging metal impurities. In this case, analyses of three products via ICP-MS led to observed levels all within the targeted range (Fig. 1).7
A representative 3-step, 1-pot sequence illustrative of the infinite opportunities now available for combining chemo- and bio-catalysis in water35 is illustrated in Scheme 10. Thus, following an enzymatic reduction of p-bromoacetophenone in an aqueous buffered medium (step A), the product (without isolation) is then used for a 1000 ppm N2Phos Pd-catalyzed Suzuki–Miyaura coupling (step B). Again, without processing the newly formed biaryl, nitro group reduction utilizing carbonyl iron powder (CIP)40 (step C) ultimately affords nonracemic aminoalcohol 32 in 81% overall yield.
Scheme 10 3-Step, 1-pot sequence involving an initial enzymatic reduction followed by a SMC, and then a nitro group reduction. |
Density functional theory (DFT) calculations and an X-ray crystal structure were obtained to gain insight into the structural factors that contribute to the greater reactivity of the N2Phos/Pd(OAc)2-derived species formed in solution relative to the corresponding catalyst derived from EvanPhos.
The X-ray crystal structure of N2Phos is shown in Fig. 2a. Two notable features include the dihedral angle between the biaryl groups of 79.47° and the orientation of the two N,N-dibenzyl moieties, the steric requirements for each forcing the aromatic rings of the biaryl unit to be almost perpendicular. Quantum-calculated geometry optimizations of the free ligand N2Phos at the B3LYP/6-31G(d) level were completed for seventeen of the most reasonable conformations. The next most stable, by 0.60 kcal mol−1, of these calculated structures had the same basic conformation as in the X-ray crystal structure. At the B3LYPD3/6-31+G(d,p) level with D3 empirical dispersion corrections and a larger basis set and at the M06/6-31+G(d,p) level, this structure was not always the lowest in energy but was within about 2 kcal mol−1 of being the lowest energy. Small computational inconsistencies and/or small crystal-packing effects could account for these differences between experiment and theory. Structural parameters for the X-ray structure and for the calculated structures at the three levels for the same basic conformation gave a close correlation between experiment and theory at all levels for bond distances and angles. Average errors were about 0.010 Å in some selected distances and 2.0–2.1° for selected angles. For selected dihedral angles that determine the exact conformation, the differences were larger, averaging 10–14°, as expected since these bond rotations have shallow energy wells. The B3LYPD3/6-31+G(d,p) level of theory gave slightly smaller geometry errors than the other two levels of theory. These geometry comparisons between theory and experiment suggest that the levels of theory chosen might also be expected to give reliable reaction energies for interconversion of the conformers (see ESI-2 for details†).
The structures and relative energies of the pre-reductive elimination intermediates, Ph2PdL, for SMCs between two simple phenyl rings to form biphenyl were calculated at the B3LYP/6-31G(d)-SDD(Pd) level using a conformation for the ligand close to the experimental conformation of the free ligand from the X-ray structure for N2Phos and for a low-energy conformer for EvanPhos.15 The calculations reveal that intermediates from both ligands form a square planar complex, as expected for Pd(II). Key distinguishing features, however, show a C–Pd–C angle of 82.17°, along with a bond length of 2.453 Å between palladium and the methoxy oxygen on the resorcinol ring of N2Phos, well in the range of the van der Wall's radii of the two elements indicative of a weak interaction (Fig. 1b). The former observation highlights a notable difference as compared with the corresponding EvanPhos-containing intermediate possessing a greater C–Pd–C angle of 87.68°. The enhanced proximity of the two phenyl rings could translate into an increased rate of reductive elimination for the N2Phos-containing intermediate compared with that with EvanPhos in this complex.41 The explanation can be attributed to the N,N-dibenzyl moieties of the naphthyl ring, such that the –NBn2 residues significantly crowd the available space around palladium, thereby forcing the two phenyl rings into closer proximity and increasing their rate of reductive elimination. No such phenomenon is in play with EvanPhos. The steric requirements of the N2Phos ligand make the formation of a 2:1 complex with diphenylpalladium impossible, with no energy minimum found for the Ph2Pd(N2Phos)2 complex, though less hindered versions lacking dibenzylamino substituents were able to form such 2:1 complexes.
More important, however, than the reductive elimination step in our analysis as to why the N2Phos ligand gives improved yields over EvanPhos with aryl chlorides is the possibly rate-limiting oxidative addition step.42 Oxidative additions of ligated Pdo (LPdo and L2Pdo; Fig. 2b) to chlorobenzene forming PhPd(Cl)L and PhPd(Cl)L2, respectively, were studied with both EvanPhos and N2Phos ligands. Oxidative addition was found to be much more downhill in free energy when carried out with monoligated species LPdo in accord with literature expectations.43,44 The oxidative addition reaction summarized in Scheme 11 the computational results with M06D3/6-31+G(d,p)-SDD(Pd) level free energies. The free energies for the oxidative addition are 8 kcal mol−1 more downhill for the monoligated species LPdo for EvanPhos than for N2Phos. This seems to contradict our presumption that catalysis by the Pd-complexing N2Phos ligand would be favored, although it does follow based on steric arguments. The equilibria for further ligation of LPdo to form L2Pdo have free energies of reaction of −36 and −16 kcal mol−1 for the EvanPhos and N2Phos ligands, respectively. Thus, the di-ligated L2Pdo is the predominant form of the Pdo catalyst, and the overall reaction to form PhPd(Cl)L from L2Pdo is predicted to be downhill by 3 kcal mol−1 for N2Phos and uphill by 9 kcal mol−1 for EvanPhos, leading to a strong 12 kcal mol−1 overall preference for the N2Phos reaction. Inclusion of solvation by toluene computationally using a continuum model makes little difference, showing the same 12 kcal mol−1 preference (ESI for details†).
Scheme 11 Reaction scheme for oxidative addition with free energies at 298 K in kcal mol−1 with EvanPhos (and N2Phos energies in parentheses) ligands at the M06D3/6-31+G(d,p)-SDD(Pd) level of theory. |
Electron-donating effects of the dibenzylamino nitrogens were also considered as a possible influence on the electron density in the vicinity of the phosphine group and perhaps, catalyst reactivity. Natural population analysis calculations at several levels of theory on the free ligands and diphenylpalladium complexes for EvanPhos, N2Phos, and an N2Phos with the nitrogens replaced with CH groups showed little or no regular variation of charge densities at phosphorus or palladium. The fact that the dihedral angles between the planes of the biaryl groups in the ligands is near 80° for the free ligands and palladium complexes would certainly be expected to seriously diminish any putative π-donating electronic effect. This suggests that the efficacy of N2Phos is likely the result of steric, rather than any significant electronic effects.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0sc00968g |
This journal is © The Royal Society of Chemistry 2020 |