Marion H.
Emmert‡
*a,
Cyndi Qixin
He‡
*b,
Akshay A.
Shah‡
c and
Stephanie
Felten‡
a
aProcess Research & Development, MRL, Merck & Co. Inc., 770 Sumneytown Pike, West Point, PA 19486, USA. E-mail: marion.emmert@merck.com
bComputational and Structural Chemistry, MRL, Merck & Co. Inc., 126 E Lincoln Ave, Rahway, NJ 07065, USA. E-mail: cyndi.he@merck.com
cDiscovery Chemistry, MRL, Merck & Co. Inc., 770 Sumneytown Pike, West Point, PA 19486, USA
First published on 5th February 2021
This manuscript reports the development of a mild, highly functional group tolerant and metal-free C–H aminoalkylation of azoles via a three-component coupling approach. This method enables the C–H functionalization of diverse azole substrates, such as oxazoles, benzoxazoles, thiazoles, benzothiazoles, imidazoles, and benzimidazoles. DFT calculations identify a key deprotonation equilibrium in the mechanism of the reaction. Using DFT as a predictive tool, the C–H aminoalkylation of initially unreactive substrates (imidazoles/benzimidazoles) can be enabled through an in situ protecting/activating group strategy. The DFT-supported mechanistic pathway proposes key interactions between the azole substrate and the Lewis acid/base pair TBSOTf/EtNiPr2 that lead to azole activation by deprotonation, followed by C–C bond formation between a carbene intermediate and an iminium electrophile. Two diverse approaches are demonstrated to explore the amine substrate scope: (i) a DFT-guided predictive analysis of amine components that relates reactivity to distortion of the iminium intermediates in the computed transition state structures; and (ii) a parallel medicinal chemistry workflow enabling synthesis and isolation of several diversified products at the same time. Overall, the presented work enables a metal-free approach to azole C–H functionalization via Lewis acid mediated azole C–H deprotonation, demonstrating the potential of a readily available, Si-based Lewis acid to mediate new C–C bond formations.
This strategy would also contribute to establishing more general and synthetically useful methods for azole C–H alkylation, which are generally challenging transformations. In fact, most azole C–H alkylation methods are promoted by transition metals and require forcing conditions, such as high temperatures, stoichiometric amounts of strong bases, and high catalyst loadings.3–5 Although various metal-free azole C–H alkylations are reported in the literature,2b,3,6 no examples of direct C–H aminoalkylations of azoles under mild reaction conditions have been described.
Addressing these limitations, this manuscript describes a Lewis acid promoted direct C–H functionalization of azole substrates, which proceeds under mild conditions (Scheme 2). Trapping the in situ generated iminium electrophile with an azole nucleophile produces a diversity of heterocyclic structures in a three-component coupling approach.
Scheme 2 Lewis-acid promoted, redox-neutral, one-pot, three-component azole C–H aminoalkylation & proposed mechanism. |
The mildness of the protocol is showcased by its ability to proceed in the presence of reactive functional groups such as amides, heterocycles, alkynes, and even enolizable ketones. Furthermore, its compatibility with a parallel medicinal chemistry workflow has been demonstrated, allowing access to compound arrays in parallel. Data from reaction optimization and DFT calculations support the need for Lewis acid and amine base to act in concert for successful C–H aminoalkylation, supporting the intermediacy of a silylated carbene as the key reactive species in these reactions.
A multidimensional optimization of reactant loadings (Scheme 5; 1.0 to 3.0 eq. TBSOTf, 1.0 to 2.5 eq. benzoxazole, 1.0 to 2.0 eq. TMS–NEt2) revealed that all varied reaction components showed a significant influence on the reaction outcome. The highest assay yields (65–66% AY; determined by calibrated UPLC analysis) were obtained under conditions with an excess of both basic and acidic reactants (2.5 or 3.0 eq. TBSOTf, 2.0 eq. TMS–NEt2, 1.5 eq. benzoxazole). These findings imply that high concentrations of acidic reagent (TBSOTf) and basic reagents promote high yields. As both basic components of the reaction mixture (TMS–NEt2, benzoxazole) contribute to the skeletal atoms of the desired product, reducing the required stoichiometry (ideally to 1:1) is highly desirable.
We reasoned that substituting the excess of basic reactants with an exogenous base may result in even better conversion and higher efficiency. To test this hypothesis, we surveyed the effects of five base additives (pyridine, EtNiPr2, DBU, NMO, 2,6-lutidine; Scheme 6) at two different loadings (0.5 eq. and 1.0 eq.) in the presence of two concentrations of TMS–NEt2 (1.25 eq. and 2.0 eq.). In this series, most of the tested conditions afforded yields close to 50% AY. The highest yield was obtained with the sterically bulky base EtNiPr2 (1.0 eq.) in combination with 1.25 eq. TMS–NEt2 (96% AY), while other common bases (pyridine, 2,6-lutidine) were detrimental to the reactivity under analogous conditions (22% and 17% AY, respectively; 1.0 eq. base, 1.25 eq. TMS–NEt2).
Scheme 6 High throughput screen: base optimization. AY = assay yield, as determined by calibrated UPLC analysis (210 nm), using biphenyl as internal standard. |
Based on this reasoning, our mechanistic hypothesis (Scheme 7) assigns important roles to EtNiPr2 and TBSOTf: we propose that the TBSOTf first reacts with the azole substrate 1 to form silylated intermediate 2, which then undergoes deprotonation. In a subsequent reaction step inspired by the mechanism of the Stetter reaction,13 the formed carbene 3 undergoes C–C bond formation with the pre-formed iminium electrophile 4. Hydrolysis of 5 during workup and/or isolation leads to the final product 6.
Preliminary experimental studies that react 1 with TBSOTf, iPr2NEt, and TfOH–D1 indeed detect mass spectroscopic evidence for 2 as well as deuterium incorporation in 1. These data (see ESI† for details) support the accessibility of a deprotonated species such as 3 under the tested conditions.
Scheme 8 (A) Calculated thermodynamics for deprotonation equilibria and C–C bond formation with benzoxazole and silylated benzoxazole. (B) HOMO of silylated benzoxazole. (C) Lowest-energy transition state of C–C bond formation. (D) Energy profile. See ESI† for details of computational methods. M06-2X/def2TZVPP/SMD(ether)//M06-2X/6-31G(d)/SMD(ether). |
Scheme 8B shows the electrostatic potential surface of the silylated azole 2 and the pictorial representation of the HOMO for the proposed carbene intermediate 3. While there are two electropositive regions in 2, the azole 2-position is spatially more accessible. Upon deprotonation, the large blue orbital lobe localized at the azole 2-position in 3 shows that this is the most nucleophilic site in the molecule.
Finally, the lowest energy transition states for the deprotonation and the C–C bond formation steps were computed (Scheme 8C and D). The proton transfer transition state was located with an activation barrier of 17.0 kcal mol−1 and is the rate-determining step of the reaction. The deprotonation barrier excluding the triflate counterion in the model is 2.1 kcal mol−1 lower. Upon formation of the carbene, the C–C bond forming TS between the silylated carbene intermediate (with TBS as silyl group and benzoxazole as azole) and the iminium electrophile [PhCHNEt2]+ is 13.6 kcal mol−1. The TS structure features a forming C–C bond distance of 2.5 Å and a favorable electrostatic interaction between the α-C–H of the TBS group and the phenyl group of the iminium reactant at 2.8 Å. Higher energy conformers are reported in the ESI.†
Scheme 9 Deprotonation (DP) equilibrium calculated for each azole substrate (shown in Scheme 10). |
Excitingly, a broad variety of azole substrates successfully underwent C–H aminoalkylation (Scheme 10), including benzoxazoles with electron-withdrawing and electron-donating substituents (7 to 10), oxazoles with various functionalized substituents (11 to 16), thiazoles (17, 18), and benzothiazoles (19, 20). Interestingly, imidazoles and benzimidazoles (22 to 28) afforded only traces of products, regardless of the substitution pattern at the N atom (H, alkyl, aryl, benzyl, allyl).
Scheme 10 Azole substrate scope. Successful substrates are shown in cyan (ΔGDP < 14 kcal mol−1); unsuccessful substrates are shown in red (ΔGDP > 14 kcal mol−1). ΔGDP = free energy of deprotonation (in kcal mol−1). Conditions: azole (1.0 eq.), TBSOTf (0.416 mL, 1.813 mmol, 3.0 eq.), TMS–NEt2 (1.25 eq.), PhCHO (1.0 eq.), EtNiPr2 (1.0 eq.), N2, 50 °C, 18 h. a Modification for substrates with acidic C–H or N–H bonds: TBSOTf (4.0 eq.), EtNiPr2 (3.0 eq.); workup for desilylation: 8 eq. KF, pyridine. b From N-tosyl imidazole; tosyl hydrolysis after reaction: pyridine (2.0 mL), water (0.5 mL), 50 °C, 3 h. See ESI† for details of computational methods. M06-2X/def2TZVPP/SMD(ether)//M06-2X/6-31G(d)/SMD(ether). |
When considering the calculated free energies of deprotonation ΔGDP for the respective silylated substrates in combination with the experimental findings, a clear pattern emerges: all azoles that show very low conversion (<2 LCAP; products 21 to 28) exhibit values for ΔGDP that are larger than 14 kcal mol−1. Furthermore, the deprotonation barrier of the least activated silylated 5-methylthiazole (product 21) is 4.1 kcal mol−1 higher than that of TBS–benzoxazole. This suggests that imidazoles and benzimidazoles as well as aryl- or alkyl-protected imidazoles are not reactive under the optimized conditions due to unfavorable deprotonation kinetics and thermodynamic equilibrium.
To circumvent this significant limitation in the azole scope, we considered activating protecting groups for imidazoles and benzimidazoles that may be able to shift the deprotonation equilibrium below the threshold energy of 14 kcal mol−1. Indeed, we discovered that urea derivative 29 (Scheme 10; ΔGDP 14.0 kcal mol−1) afforded reasonable crude yields of the desired coupling product (41 LCAP); however, isolation by column chromatography on silica led to decomposition of the product, resulting in a low isolated yield (3%). This observation suggests that electron-withdrawing protecting groups (such as acyl groups) can activate imidazole or benzimidazole substrates for C–H deprotonation and subsequent aminoalkylation.
In an extension of the activating protecting group strategy, subjecting N-tosyl imidazole (ΔGDP 10.5 kcal mol−1) to the standard conditions afforded an 88% isolated yield of the desired imidazole product 28. To enable analogous reactivity with more diverse imidazoles and benzimidazoles that might not be as easily available as tosylated precursors, we further developed a procedure (Scheme 11) combining in situ tosylation, C–H aminoalkylation, and detosylation in one pot. With this sequence, imidazole and benzimidazole coupling products 28 and 31 were synthesized, starting from unprotected imidazole and benzimidazole substrates in moderate to good yields (45% and 74%, respectively) without the need to isolate the N-tosyl intermediates. These results demonstrate that the use of electron-withdrawing protecting groups can be implemented successfully to enable the C–H aminoalkylation of less activated azoles.
We first calculated transition state energies for reactions between a series of differently substituted iminium electrophiles (32; Scheme 12A) and carbene intermediate 3 (Table S2†). The substitution patterns for 32 were chosen to reflect common amine cores in pharmaceutical chemistry and to vary between low and high steric demand. DFT calculations predict, accordingly, activation energies that span from 7.6 to 20.3 kcal mol−1 (with respect to the silylated carbene and iminium). Interestingly, the C–C bond forming step for the dicyclohexyl iminium electrophile is 12.6 kcal mol−1 higher than that for diethyl iminium, which is higher than the deprotonation barrier in the reaction pathway (Scheme 8D).
Scheme 12 Distortion–interaction analysis of reactions with iminium electrophiles. (A) Transition state energies ΔG‡ for series of iminium electrophiles. (B) Overlay of transition state structures. (C) Distortion energies ΔEdistversus ΔG‡ of C–C bond formation. See ESI† for details of computational methods. M06-2X/def2TZVPP/SMD(ether)//M06-2X/6-31G(d)/SMD(ether). |
A distortion–interaction analysis14 (Scheme 12B and C) reveals that the total distortion energies ΔEdist for the components in each reaction are unique and correlate with the activation free energies ΔG‡ of C–C bond formation. In contrast, the values of the interaction energy term between 3 and each of the iminium intermediates 32 lie within 0.5 kcal mol−1. Scheme 12B depicts an overlay of all computed transition state structures. It shows that carbene intermediate 3 remains in a very similar geometry in all computed transition state structures (ΔEdist,3 = 0.1–0.7 kcal mol−1). In contrast, the iminium electrophiles distort away from their preferred planar geometry to avoid potential steric clashes with 3 going from reactant to the transition state (see Fig. S11† for a breakdown of the total distortion energies). These data suggest that sterically congested iminium electrophiles such as the electrophile obtained by condensation of PhCHO and Cy2NH may not undergo C–C bond formation due to the high distortion energy necessary to achieve the transition state.
Scheme 13 Amine scope I. Conditions: see ESI.† |
This presented an additional challenge for reaction development: the majority of complex amines are only commerically available as their corresponding HCl salts and direct use of these salts did not afford any product under standard conditions. We hypothesized that chloride ions are detrimental to the reactivity, as chloride might react with TBSOTf and form less reactive TBSCl. Therefore, we developed a protocol that first removes chloride under non-aqueous conditions (Scheme 14, first step). After filtration, the resulting solution of the free amine can then be directly used for silylation and iminium ion formation. Using Et2NH2Cl as the amine with this protocol, 59% of the desired C–H aminoalkylation product 7 was isolated, demonstrating that employing HCl salts directly is feasible with this modified protocol.
Scheme 14 Non-aqueous free-basing conditions allowing use of amine HCl salts in azole C–H aminoalkylation. |
Ten complex amines that are commercially available as HCl salts were thus subjected to the PMC workflow under the modified conditions (Scheme 15). The amine structures depicted were chosen to reflect structures of interest for medicinal chemistry, as they are highly saturated15 and incorporate fluorinated moiteties.16 All employed monomers successfully formed product; nine out of ten of these products were isolated in ≥95% purity (suitable for primary assay screenings) by mass-guided, automated preparative HPLC. This is equivalent to an excellent success rate of 90%, as calculated by the number of amines yielding the intended product divided by the total number of amines employed. This result suggests that the established conditions will be applicable for synthesizing compound arrays of interest to the medicinal chemistry community.
Scheme 15 Amine scope II: azole C–H aminoalkylation in parallel. Conditions: see ESI.† All reactions were conducted on a 0.18 mmol scale. Products were isolated as TFA salts after automated preparative HPLC purification. |
Scheme 16 Aldehyde scope. Conditions: azole (1.0 eq.), TBSOTf (0.416 mL, 1.813 mmol, 3.0 eq.), TMS–NEt2 (1.25 eq.), PhCHO (1.0 equiv.), EtNiPr2 (1.0 equiv.), N2, 50 °C, 18 h. |
Excitingly, a thiophene-derived aldehyde was also employed successfully, forming product 62 in 82 LCAP (26% IY). Furthermore, pivaldehyde was successfully reacted to afford 63, showcasing that aliphatic aldehydes are generally also reactive under the established standard conditions. In contrast, aliphatic aldehydes with acidic C–H bonds in β-position to the in situ formed iminium intermediate do not form any product (64 to 66). This can be rationalized by the ability of the iminium intermediates to undergo deprotonation to the enamine products (shown in insert in Scheme 16). This would remove both the electrophile and the reactive base from the reaction solution and is thus expected to be detrimental for C–H aminoalkylation. Indeed, analysis of a crude reaction mixture (for details see ESI†) detects evidence for the formation of enamine side products.
Overall, this manuscript outlines a new strategy to diversify azoles via Lewis-acid mediated C–H functionalization under much milder conditions than those typically employed with transition metal catalysts.3–5 Moreover, the protocol employs both acidic and basic components, which is another distinguishing feature from more common, base-promoted azole C–H functionalization approaches.17 Importantly, this leads to broad functional group tolerance, which is typically inaccessible through other direct azole functionalizations approaches proceeding through catalytic or stoichiometric metalation.3–6
Footnotes |
† Electronic supplementary information (ESI) available: Detailed experimental procedures, characterization data for new compounds, details of DFT calculations. See DOI: 10.1039/d0sc06868c |
‡ The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. |
This journal is © The Royal Society of Chemistry 2021 |