Richard S.
Massey
a,
Jacob
Murray
a,
Christopher J.
Collett
b,
Jiayun
Zhu
a,
Andrew D.
Smith
b and
AnnMarie C.
O'Donoghue
*a
aDepartment of Chemistry, Durham University, South Road, Durham DH1 3LE, UK. E-mail: annmarie.odonoghue@durham.ac.uk
bEaStCHEM, School of Chemistry, University of St. Andrews, North Haugh, St Andrews KY16 9SY, UK
First published on 14th December 2020
Steady-state kinetic and structure–activity studies of a series of six triazolium-ion pre-catalysts 2a–2f were investigated for the benzoin condensation. These data provide quantitative insight into the role of triazolium N-aryl substitution under synthetically relevant catalytic conditions in a polar solvent environment. Kinetic behaviour was significantly different to that previously reported for a related thiazolium-ion pre-catalyst 1, with the observed levelling of initial rate constants to νmax at high aldehyde concentrations for all triazolium catalysts. Values for νmax for 2a–2f increase with electron withdrawing N-aryl substituents, in agreement with reported optimal synthetic outcomes under catalytic conditions, and vary by 75-fold across the series. The levelling of rate constants supports a change in rate-limiting step and evidence supports the assignment of the Breslow-intermediate forming step to the plateau region. Correlation of νmax reaction data yielded a positive Hammett ρ-value (ρ = +1.66) supporting the build up of electron density adjacent to the triazolium N-Ar in the rate-limiting step favoured by electron withdrawing N-aryl substituents. At lower concentrations of aldehyde, both Breslow-intermediate and benzoin formation are partially rate-limiting.
Fig. 1 Previous mechanistic studies of the NHC-catalysed benzoin condensation and the catalysts studied in this work. |
In terms of triazolium catalysis, our groups have reported the isolation and structural characterisation of a series of triazolium-derived adducts (III; X = Y = N) and demonstrated that (III) form reversibly from free aldehyde and triazolium pre-catalyst. Kinetic analysis of the H/D-exchange reactions at the α-position enabled structural effects on Breslow intermediate formation from (III) to be directly assessed (Fig. 1). Stoichiometric studies provided access to rate constants for initial adduct (III) formation and decay from triazolium pre-catalyst.6 An earlier study by Huskey and Jordan reported a related kinetic analysis of thiazolium-aldehyde adducts and the determination of a C(α)–H pKa of 15.7 for one example.7 To date, there has been no report of experimental adduct (III) C(α)–H pKa values for the triazolyl series.
Regarding detailed kinetic analysis of the overall NHC-catalysed benzoin reaction, integrated-rate analyses from various studies gave different and often ambiguous results, with all studies mainly focusing on thiazolium catalysis (I; X = S, Y = C).2b,8 Key issues in these investigations include the observation that the reaction order changes with the concentration of benzaldehyde ([PhCHO]); catalyst concentration is not necessarily constant due to degradation; concentration of base can be affected by formation of benzoic acid side-product over time; and finally, the assumption of irreversibility is not always appropriate. For enzymatic processes, which suffer from similar problems, kineticists employ an initial rate method (IRM) at low % product conversion to circumvent these issues. Importantly, Leeper and co-workers reported an IRM screening of thiazolium catalyst 1, under true catalytic conditions (Fig. 1), reporting a first-order dependence on [PhCHO], with a pseudo-first order rate constant, kcat = 6 × 10−6 s−1. Three rate-limiting scenarios were considered: (i) adduct (III) formation; (ii) Breslow intermediate (IV) formation; (iii) attack of IV to a second PhCHO (Scheme 1). Further stoichiometric studies and deuterium labelling determined that all three stages of the process must be in-part rate-limiting to account for all the observations. At higher concentrations of benzaldehyde, rate is more limited by deprotonation to form IV (kBI); conversely, at reduced concentrations of benzaldehyde, the initial adduct and product-forming steps are considered to be more rate-limiting.9
Leeper's IRM and stoichiometric studies of the thiazolium-catalysed mechanism highlight the importance of rigorous understanding of the kinetics of catalytic processes. Knowledge of the rate-limiting processes within a reaction allow for targeted development of catalysts for improved methods. Moreover, as the push towards aqueous organocatalysis increases, detailed mechanistic understanding of key rate-limiting factors will aid catalytic development for a more synthetically challenging environment.1b,c 1,2,4-Triazolium ions are often the preferred NHC catalysts due to increased propensities for functionalization and therefore stereoselectivity, but also increased activities.1a pKa values of triazoliums are 1–3 units lower than their thiazolium relatives, correlating to enhanced rates within the benzoin reaction.6a,10
Despite the widespread application of 1,2,4-triazolium ions in contemporary organocatalysis, quantitative kinetic insights into the triazolium-catalysed benzoin condensation are limited. Herein, we report an IRM-study of the triazolium catalysed reaction, under true catalytic conditions, and a structure–activity study of reaction data. A polar, protic MeOH environment was chosen for kinetic analysis, to enable comparison with the previous thiazolium ion study by Leeper. Our results build upon our previous stoichiometric triazolium catalyst studies providing quantitative insight of the interplay of structure and kinetic activity for these commonly-used organocatalyst systems.9
ν = kcat[PhCHO]n | (1) |
Using identical reaction conditions as employed for thiazolium salt 1 (triethylamine-buffered MeOH at 50 °C), triazolium pre-catalysts 2a–f showed distinctly different kinetic behaviour. Triazolium pre-catalysts 2a–f clearly demonstrate saturation kinetic profiles at higher benzaldehyde concentrations (e.g.Fig. 2 for 2c; ESI, section 2, Fig. S14–18† for 2a–b, 2d–f). By contrast, saturation was not observed for plots of ν versus [PhCHO]0 for 1 at 12, 24 and 30 mM thiazolium catalyst (ESI, Fig. S13†).
Fig. 2 Triazolium catalysis of the benzoin condensation under catalytic conditions at 50 °C: initial rate (ν, M s−1) plot for 2c at 24 (red), 12 (green) and 6 (blue) mM catalyst loading. The solid lines show the fit of the reaction data to eqn (2). |
Initial rate data for 2a–f were fit to a steady-state rate equation (eqn (2)), derived for the Breslow intermediate IV in Scheme 2, with the assumption that the concentration of adduct ([III]) was equal to that of catalyst ([I]). This assumption is validated by our previous kinetic NMR studies with stoichiometric aryl aldehyde and triazolium salt concentrations in the same triethylamine-buffered methanol medium, which showed that adduct (III) forms rapidly and reversibly from these reactants at 25 °C.6b These previous studies for a range of aryl aldehydes and triazolium catalysts included the stoichiometric kinetic analysis of 2a–d and benzaldehyde, permitting access to rate and equilibrium constants for formation of (III).6b Half-lives (t1/2) for adduct III formation at 25 °C could be calculated as ∼670 s and ∼16 s for the lowest and highest [PhCHO]0, which would be predicted to be a minimum of ∼10-fold lower again at 50 °C, and significantly smaller than the kinetic experiment measurement times employed in the initial rate analysis (ESI, section S4†). Using the experimental equilibrium constants for adduct formation determined under stoichiometric conditions, values for the percentage of free pre-catalyst (%I) at equilibrium can also be estimated (Table S14†). Under the catalytic conditions of the initial rate study with a 10–100 fold excess of benzaldehyde relative to pre-catalyst, the equilibrium significantly favours adduct. Additionally, the use of the IRM approach permits the assumption that the reverse reaction of benzoin is not significant under the reaction conditions.‡
(2) |
The experimental data under catalytic conditions fit well to eqn (2) and, at high [PhCHO]0, the values for the initial rates plateau and approach νmax. There is some deviation of datapoints at the lowest [PhCHO]0 perhaps owing to a more significant % of pre-catalyst I still present initially at the lowest excess of aldehyde. Uncertainties in νmax were obtained from the kinetic fits and, with the exception of 2e, are relatively small. The uncertainty for the νmax value for 2e is larger due to this system not reaching a final plateau at the highest benzaldehyde concentration employed. Table 1 summarises the kinetic constants obtained from the fit to eqn (2) for catalysts 2a–f (ESI, section S2† for kinetic fits for 2a,b,d–f). The absolute values of these kinetic parameters are unique to the triethylamine buffer system employed in methanol as solvent, however, their comparison allows the effects of structural changes of catalyst on individual reaction steps to be determined. Values for νmax vary greatly depending on the catalyst, with a large 75-fold difference between the slowest 2a and fastest 2f triazolium catalysts, and the largest value observed for a N-pentafluorophenyl substituent. This highlights the importance of targeted tuning of the electronic nature of the N-aryl substituent for efficient catalysis.
Catalyst | [I]/mM | ν max/×10−7 M s−1 | k BI/s−1c | k BI (rel)d | k p/k−BI |
---|---|---|---|---|---|
a Values determined using an initial rate method (IRM) to a maximum of 10% conversion to benzoin product with initial benzaldehyde concentrations ranging from 0.16–1.6 M, in 2:1 0.16 M Et3N/Et3NH+Cl− buffered MeOH solutions at 50 ± 0.01 °C. b Reaction in the presence of catalyst 2f was particularly fast, thus requiring lower loadings to enable study by the IRM. c Uncertainties in kBI can be approximated as similar to those in νmax, as kBI = νmax/[I], with uncertainties in [I] unlikely to have significant impact. d k BI(rel) determined as kBI(2a–f)/kBI (2c) at the same catalyst concentration. | |||||
2a | 12 | 4.08 ± 0.53 | 3.40 × 10−5 | 0.37 | 5.7 |
6 | 2.27 ± 0.30 | 3.78 × 10−5 | 0.41 | 5.0 | |
2b | 12 | 6.52 ± 0.49 | 5.44 × 10−5 | 0.59 | 4.8 |
2c | 24 | 19.2 ± 3.9 | 8.02 × 10−5 | 1.6 | |
12 | 11.1 ± 2.5 | 9.24 × 10−5 | 1.0 | 2.0 | |
6 | 5.52 ± 1.18 | 9.21 × 10−5 | 1.0 | 2.4 | |
2d | 12 | 16.2 ± 2.1 | 1.35 × 10−4 | 1.5 | 4.2 |
2e | 12 | 143 ± 43 | 1.19 × 10−3 | 13 | 0.94 |
2f | 6b | 171 ± 3.8 | 2.85 × 10−3 | 31 | 4.0 |
The plateauing of initial rate plots with respect to [PhCHO]0 confirms a change in order of reaction and therefore change in rate-limiting step. Using eqn (2), at high [PhCHO]0 when kp[PhCHO] ≫ k−BI, values for νmax approach kBI[I] i.e. deprotonation to form the Breslow intermediate (IV) becomes rate-limiting for triazolium salts 2a–f at high concentrations of aldehyde. By contrast, as a saturation point is not witnessed within the IRM for thiazolium salt 1, the formation of Breslow intermediate IV does not become cleanly rate-limiting under catalytic conditions in this case. Catalysts 2a and 2c were explored at different catalyst concentrations ([I]). Values for the pseudo-first order kinetic constant kBI, determined as νmax/[I], were consistent within <10% for 2a and <1% for 2c, supporting the present kinetic and mechanistic analysis.
Electron withdrawing N-aryl substituents are expected to favour formation of Breslow intermediate IV from adduct III by stabilising the negative charge formed upon deprotonation and thereby increasing the kinetic acidity at the α-position. The observed νmax values are highest for more electron-withdrawing substituents, thus supporting the proposal of rate-limiting formation of IV at high [PhCHO]0 and assignment of kBI as the rate constant for this step.
Magnitudes of kBI determined herein are in good agreement with our reported values of kex for H/D-exchange of O-methylated adducts 3 of an ortho-substituted aldehyde with the same series of 1,2,4-triazolium catalysts (Scheme 3).11O-Methylation prevents both reverse equilibration to free aldehyde and catalyst, in addition to onward reaction to benzoin product, thus permitting focus just on the Breslow-intermediate forming step. Although this H/D-exchange study was performed in a predominantly aqueous D2O medium at a lower temperature of 25 °C, the relative values of kex and kBI for a given triazolium pre-catalyst are closely similar. Again, this supports our conclusion that the plateau region of the initial rate plots corresponds to the Breslow intermediate (IV) forming step and highlights the similarity of triazolium substituent effects in methanol and water media.
Scheme 3 Previous studies of H/D-exchange of O-methylated hydroxyaryl adducts 3 derived from 2b, 2c and 2f.11 Relative values for kex (s−1) are similar to kBI (rel) determined herein. |
Synthetic observations for the benzoin reaction with catalysts 2a–2f report increased yields of products with more electron withdrawing N-aryl substituents.1a,10b,12 Values for νmax (and kBI) increase with more electron withdrawing N-aryl substituents, thus substituent effects on νmax under catalytic conditions correlate with synthetic outcome. An excellent linear Hammett correlation is observed for para-substituted catalysts 2a and 2c–e with ρ = +1.66 ± 0.05 (Fig. 3).13 The same Hammett ρ-value is obtained from either log(νmax) or log(kBI) correlations given the common catalyst concentration. Our experimental Hammett ρ-value is intermediary between reference values of ρ = +1.44 and ρ = +3.32 for the acid dissociations of benzoic acids and protonated anilines in methanol, respectively (Fig. 4).14 This significant sensitivity to substituent change supports the formation of an enaminol-like transition state for the conversion of III to IV with significant charge neutralisation on the N-aryl nitrogen of catalyst in the rate-limiting step (cf.IV-b). A significantly smaller Hammett ρ value would be expected for a more localised carbanion-like transition state with charge change more remote from the N-Ar of catalyst (cf.IV-a). The solid state characterisation of Breslow-like intermediates achieved to date (e.g.Fig. 1, IV′, IV′′, IV′′′) also support enaminol-like structures, thus the present Hammett analysis permits unification with behaviour in solution.
Fig. 4 Comparison of the Hammett ρ values for the reference benzoic acid dissociation (+1.44), the dissociation of anilines (+3.32) and the experimental value obtained herein (+1.66) in methanol {values in braces indicate reference values in water}.14 Assuming νmax reflects rate-limiting deprotonation of III to IV, the experimental value for ρ suggests significant increase in electron density near the N-aryl substituent in the transition state suggestive of a greater contribution from an enaminol-like resonance canonical structure IV-b rather than a more localised carbanion IV-a. |
Data for ortho-substituted triazolium salts 2b and 2f were not included in the correlation. Three literature Hammett σ-values have been reported for the pentafluoro-substituent of 2f, however, there is no reported value for the N-mesityl group of 2b. As is frequently observed owing to the presence of additional steric and electronic factors in ortho-substituted cases, the datapoint for 2f significantly deviates from the correlation line for the para-substituted triazolium salts (section S5, ESI†). Irrespective of the Hammett σ-value employed, the log (νmax) value for 2f falls significantly below the correlation line.
Values of log(kBI/[I]) may also be correlated with pKas of the triazolium ion pre-catalysts 2a–2f in both water and DMSO (Fig. S24, ESI†).6a,10c,d Although strictly not Brønsted plots, which would require correlation with α-H pKa values of III, these plots serve to demonstrate the increase in kBI with more acidic triazolium pre-catalysts. Harper et al. have reported Hammett analyses of pKa (C3–H) values for three different series of triazolium salts in both water and DMSO.10d A Hammett ρ-value of −0.93 ± 0.05 was obtained for the aqueous pKas of the same family of triazolium catalysts as utilised in the present study with a decrease in pKa observed for more electron-withdrawing N-aryl substituents. Thus, electron-withdrawing N-aryl substituents favour acid dissociation of triazolium pre-catalyst in addition to the deprotonation step from adduct III to Breslow intermediate IV.§
The effect of N-aryl substituent on product partitioning of the Breslow intermediate IV may also be evaluated by comparison of values for kp/k−BI obtained from the fit to eqn (2) (Table 1). With the exception of catalyst 2e, all ratios are greater than 1, suggesting greater contribution from the product benzoin forming step. Values for kp/k−BI vary in the range 0.94–5.7, with the reliance on two separate kinetic constants being the likely reason for the absence of a clear trend with the electronic nature of the substituent. Interestingly, the highest product ratio is obtained for the most electron-donating N-aryl substituent in 2a, suggesting that the nucleophilicity of IV possibly dominates in the benzoin forming step.
Values for νmax vary greatly depending on the catalyst, with a large 75-fold difference between the slowest 2a and fastest 2f triazolium catalysts. These data provide quantitative insight into the role of catalyst N-aryl substitution on the rate-limiting step under synthetically-relevant catalytic conditions in a polar solvent environment. In particular, Hammett structure–activity analysis provides evidence for the formation of an enaminol-like transition state IV-b for the Breslow intermediate-forming step, permitting the alignment of previous solid state structural characterisation of IV with solution state behaviour. With a drive to further the organocatalysis field to more aqueous media, thorough mechanistic analysis in such environments is essential. This work expands our knowledge of the well-explored benzoin condensation, providing mechanistic insight into synthetic observations.
Footnotes |
† Electronic supplementary information (ESI) available: Synthetic and kinetic methods; kinetic data fitting and analysis. See DOI: 10.1039/d0ob02207a |
‡ The IRM was chosen, with the reaction followed to <10% product formation, to minimize any contribution of the reverse benzoin condensation. In our previous stoichiometric studies involving 1:1 triazolium:aldehyde, which were followed to high % product, we investigated the contribution of the reverse reaction for the 2,4,6-trichlorophenyltriazolium tetrafluoroborate salt.6b Using two representative substituted benzoins as starting material, <12% retro benzoin reaction was detectable under these conditions over long timescales. |
§ Rate constants for proton transfer from 2a–f to HO− are in the range 106–108 M−1 s−1 and reprotonation of the triazolyl carbene by water has been demonstrated to be close to the diffusional limit for dielectric relaxation of water (1011 s−1).6a As protonation of diphenyl carbene by methanol has also been shown to be limited by dielectric relaxation of solvent,15 the proton transfer step between triazolium ions 2a–f and the corresponding triazolyl carbenes is very unlikely to be rate-limiting in methanol solvent. |
This journal is © The Royal Society of Chemistry 2021 |