Cristiana
Margarita
,
Piret
Villo
,
Hernando
Tuñon
,
Oscar
Dalla-Santa
,
David
Camaj‡
,
Robin
Carlsson‡
,
Malin
Lill‡
,
Anja
Ramström‡
and
Helena
Lundberg
*
Department of Chemistry, KTH Royal Institute of Technology, Teknikringen 30, S-10044 Stockholm, Sweden. E-mail: hellundb@kth.se
First published on 24th September 2021
Kinetic analysis was used as a tool for rational optimization of a catalytic, direct substitution of alcohols to enable the selective formation of unsymmetrical ethers, thioethers, and Friedel–Crafts alkylation products using a moisture-tolerant and commercially available zirconium complex (2 to 8 mol%). Operating in air and in the absence of dehydration techniques, the protocol furnished a variety of products in high yields, including glycosylated alcohols and sterically hindered ethers. In addition, the kinetic studies provided mechanistic insight into the network of parallel transformations that take place in the reaction, and helped to elucidate the nature of the operating catalyst.
In direct etherification, two alcohols are joined with a loss of water. To enable the formation of unsymmetrical ethers from two different alcohols, the transformation typically relies on one alcohol being more electrophilic than the other, e.g. by π-activation. For this reason, secondary or highly electron rich benzylic alcohols are commonly employed as electrophiles in this type of transformation,22,47,61–65 with the more challenging direct etherification of unsubstituted benzyl alcohol being less prevalent.66–68 In addition, symmetric side-products are recurrently observed for such unsymmetrical intermolecular etherification protocols and can significantly decrease the yields of these processes.21,69–72
During the last decades, user-friendly methodologies have been developed for intuitive visual kinetic analysis of organic reactions, known as reaction progress kinetic analysis (RPKA) and variable time normalization analysis (VTNA).73–78 These methodologies enable facile extraction of kinetic data by utilization of full reaction profiles for a minimal number of experiments under synthetically relevant conditions. Typically, these kinetic tools have been used for elucidation of operating mechanisms for organic transformations once a protocol has been established. Less commonly, these methods can be integrated with synthetic method development to enable streamlined reaction optimisation.79–83 Herein, we describe how the use of kinetic analysis during the early method development phase of a zirconium-catalysed protocol for direct substitution of alcohols enabled rational optimization of reaction conditions and enhanced yields by suppression of the symmetric side-product formation.
Firstly, different solvents were assessed (Fig. 1A). Similarly to the observed behaviour of Zr-catalysed esterification,84 non-polar solvents were found to furnish the unsymmetrical product 2b in higher yields compared to more polar solvents. In addition to the results in Fig. 1A, acetonitrile was assessed as solvent and resulted in <5% yield 2b. The formation of 2b in sulfolane proceeded with a high rate, but the yield was limited by a mixture of unidentified by-products. The highest yields for 2b were achieved in heptane and benzotrifluoride (BTF) with similar reaction rates. While the former solvent may present benefits in terms of environmental impact,85,86 the latter was selected for further studies due to its ability to dissolve a wide range of organic substrates. An assessment of reaction temperature indicated that heating above 80 °C is required to achieve reasonable reaction rates (Fig. 1B). For all solvents, the symmetric side-product 2a was formed along with the desired unsymmetrical product 2b, thereby limiting the yield of the latter.
To increase the reaction selectivity, a more efficient incorporation of benzyl alcohol (1a) into product 2b was required, and a kinetic study was undertaken using the tools of RPKA and VTNA.73–78 To assess the driving forces for the formation of the two ethers, a set of different excess experiments were conducted to estimate the order in a given component. This was done by comparison of reactions where the initial concentration of the given component was varied while the others were kept constant.73,74 As shown in Fig. 2A and B the formation of 2a and 2b demonstrated a near first order rate dependence on [Zr]. The average order77 of [1a] for the formation of 2b was experimentally found to be 0.8 at low concentrations of 1a, whereas it was reduced to 0.35 at higher concentrations (Fig. 2C and D). In contrast, the rate dependence on [1a] for the formation of side-product 2a was found to be of first order for the whole concentration range examined (Fig. 2E). A slight positive rate dependence on [1b] was observed in the low concentration range for 2b formation (Fig. 2F). At higher [1b] a negative average order of −0.35 was observed for d[2b]/dt (Fig. 2G), whereas a more pronounced negative average order in [1b] was observed for the formation of side-product 2a for the whole concentration range.
As a result of the different responses to concentration changes for 2a and 2b formation, the combined data in Fig. 2 suggested that higher selectivity for the unsymmetrical product 2b could be obtained if the concentration of 1a was kept low and 1b high throughout the reaction. Indeed, sequential addition of 1a in five portions over four hours to a BTF-solution of the catalyst and 1 M 1b resulted in maintained yield even with an increase in [1a] to a 1:
1 molar ratio to [1b], all other concentrations kept constant (Fig. 3). Slow addition of 1a with a syringe pump (rate 10 μl h−1) in the equimolar reaction gave slightly lower yield for 2b compared with sequential addition, and was not explored further. Sequential addition of 1a using a 2
:
1 molar ratio of 1b to 1a resulted in 84% isolated yield of 2b after 24 hours, representing an increase in yield with 14 percent units compared to a setup where all material was added at the outset of the reaction (“batch addition”) using the same reactant stoichiometries (Fig. 3). As expected, the same yield of 84% of 2b was also obtained by dilution of [1a] while the concentrations of catalyst and 1b were kept constant (Fig. 3). Under these optimized “batch” conditions, the ratio of desired mixed ether 2b to symmetrical side-product 2a was measured to be 12
:
1 by 1H NMR spectroscopy of the crude reaction mixture.
While the sequential addition enables the use of less reagents, and effectively a lower catalyst loading with respect to the limiting substrate, batch addition using low [1a] offers operational simplicity. From a synthetic standpoint, the two strategies for improving the selectivity of the reaction are complementary. Interestingly, a good isolated yield (69%) could still be obtained when lowering the catalyst loading (0.01 M, 2 mol%) in neat conditions, using a 2:
1 molar ratio of 1b to 1a. However, the ratio for the resulting ethers 2b
:
2a was 4
:
1 in this case. Higher concentrations of catalyst and reactants did not improve yields further after 24 hours. No reaction took place in the absence of catalyst.
A set of same excess experiments, mimicking 25% conversion in the absence and presence of the corresponding amount of water indicated that insignificant catalyst deactivation or product inhibition occurs during the reaction (Fig. 4).
Apparent Lewis acid catalysis has occasionally been reported to be due to a Brønsted acid, formed by hydrolysis of the former.87,88 In the context of metallocene triflates, Hollis et al. reported that the titanium analogue of the herein employed catalyst may undergo hydrolytic decomposition to release trifluoromethanesulfonic (triflic) acid.89 To probe whether this hydrolysis product was the operating catalyst in our protocol, a set of experiments were carried out. While addition of di-tert-butyl-4-methylpyridine to the reaction mixture containing either zirconocene triflate or triflic acid as catalyst proved inconclusive (see ESI†),87,88 an assessment of the two catalysts in the presence of different amounts of water made it clear that the two catalysts were not displaying identical behaviours (Fig. 5A). While both catalysts appear to require a sufficient amount of water, as is evident from the negligible yield in the presence of molecular sieves, the yield of 2b in the triflic acid-catalysed reaction drops gradually as a function of water addition. This behaviour is not observed for the zirconocene complex, where the 24 h yields stayed around 70% for the reactions with ratios up to 25:
1 H2O
:
catalyst (1
:
1 ratio H2O
:
1a). In the 125
:
1 H2O
:
catalyst experiment (5
:
1 ratio H2O
:
1a) the yield of 2b did not exceed 2% in the presence of the Zr catalyst, whereas triflic acid catalysis resulted in inconsistent yields between 0 and 63%, 10% being the median value (Fig. 5A and ESI†). These observations suggest that while the zirconium complex is irreversibly deactivated at high concentrations of water, the inhibition of triflic acid is reversible. Furthermore, a set of triflic acid catalysed experiments with different [1b] were carried out, indicating an order of zero in [1b] for both d[2b]/dt and d[2a]/dt in the concentration range 1–2 M 1b (Fig. 5B and ESI†). These results stand in stark contrast to the observed negative orders in [1b] for both product and side-product formation in the zirconium-catalysed reactions (Fig. 2G and H). In addition, the observed selectivities 2b
:
2a were lower for reactions carried out in the presence of triflic acid. Thus, the different reaction behaviours observed when using either triflic acid or the zirconocene complex suggest that different active catalysts are operating in the two sets of reactions. By NMR-spectroscopy, a downfield shift was observed for 13C-labelled 1a in the presence of the Zr complex in d8-toluene, suggesting interaction between the metal and substrate (see ESI†). These combined observations support involvement of the zirconium centre in the mechanism.
The relative stability of the zirconium catalyst towards water, and the inhibition of catalytic activity in the presence of molecular sieves have previously been observed for zirconocene triflate.84 While the titanium analogue was reported to readily form aqua complexes with labile ligand structure,28 zirconocene with fluorinated sulfonate ligands can form dimers in the presence of water.90 Hence, the observed decelerating effect of molecular sieves may suggest that aqua complexes and/or multinuclear complexes are the catalytically active species in this system, where the former may or may not act as Brønsted acids.91 The need for a sufficient amount of water for efficient catalysis to occur has previously been reported for different types of Lewis acidic catalysts.92–95
A faster reaction rate, compared to benzyl alcohol 1a, was observed when the para-position was substituted with an electron-donating methoxy group, whereas an electron-withdrawing para-trifluoromethyl group displayed a considerably slower rate (Fig. 6). This behaviour is consistent with a mechanism involving a carbocationic intermediate, a notion further supported by observed racemization in ether 2c, formed from enantiopure (S)-1-phenylethanol and 1b (Scheme 1).
![]() | ||
Fig. 6 Reaction rate as a function of electrophile substitution. Standard conditions: electrophile (0.5 M), 1b (1 M), Zr(Cp)2(CF3SO3)2·THF (0.02 M), BTF, 100 °C. |
The observed first order rate dependence on [Zr] and the positive orders on [1a] for the formation of both product 2b and side-product 2a (Fig. 2) are compatible with a mechanism in which the two ethers are formed via a single carbocationic intermediate, with C–OH activation constituting a turnover limiting step of the catalytic cycle. The average order in [1a] for the formation of [2b] was reduced from 0.8 to 0.35 at higher concentrations of 1a, whereas the average order in [1a] for the formation of 2a remained at 1. This finding explains the observed selectivity enhancement with respect to the unsymmetrical product 2b as [1a] is decreased, since the rate of 2b formation will be less disfavoured compared to that of the symmetric side-product. It is interesting to note that the strategy of maintaining reagent concentrations low for enhanced selectivity has been successful for various transformations, despite having quite different mechanistic origins.96–99 The average order in [1b] was increased from −0.35 to −0.2 at higher [1a] (see ESI†). These findings may be interpreted as if the contribution of the nucleophilic addition to the carbocationic intermediate to the rate is more pronounced under these conditions. The negative orders in [1b] for both product and side-product formation (Fig. 2) suggest that this compound is involved in the formation of off-cycle species that reduces the concentration of active catalyst in the reaction and, in turn, the reaction rate. Inhibition of catalytic activity in the presence of alcohols has previously been observed for zirconium chloride complexes, and the corresponding zirconocene dichloride complex did not yield the desired product 2b under our conditions (see ESI†).100,101 A tentative mechanism for the direct etherification and the rate laws for product and side-product formation under synthetically relevant conditions are shown in Fig. 7. The viability of the proposed mechanism was supported by kinetic modelling (see ESI†). Depending on the nature of the substrate it can be assumed that a carbocationic intermediate may form in solution as a more or less tight ion pair.102 Consequently, the degree of SN1-character of the mechanism may vary between substrates.
In addition to direct etherification, 2b can form by transetherification of 2a with 1b under the catalytic conditions (see ESI†).47 The rate of this process was found to be positively correlated to [2a] and [Zr] with a negative rate dependence on [1b]. A comparison between direct etherification and transetherification demonstrated that the latter was considerably slower than the former even at high concentrations of 2a. Thus, transetherification of 2a was judged to constitute a minor pathway for the formation of 2b, mainly operating towards the end of the reaction when [1a] has been depleted and [2a] has reached its highest levels. It is to be noted that the integrated use of RPKA and VTNA in the synthetic optimization provided this mechanistic information already at the early stages of the project, thus obviating the need for a separate work-intensive and time-consuming study using initial rate kinetics.68
Using the optimized conditions, a substrate evaluation was carried out to form a range of differently substituted ethers (Scheme 2). Higher yields for the unsymmetrical ethers were observed for the majority of the substrates using low electrophile concentration, either by sequential addition of the electrophile or by dilution under batch conditions. For operational simplicity, dilute batch conditions were chosen as standard procedure for the rest of substrate evaluation, unless otherwise indicated. Remarkable selectivity towards the formation of unsymmetrical products was observed for most benzylic substrates (ratios for formed ethers were determined by 1H NMR spectroscopy of the crude reaction mixtures, and are reported in ESI† where applicable). Whereas benzyl alcohol 1a required a reaction temperature of 100 °C to efficiently convert to the product, benzylic alcohols with electron donating groups in para position were amenable to lower reaction temperatures and/or reduced reaction times (2bvs.2j and 2k). This observation is likely due to stabilization of the proposed carbocationic intermediate. Similarly, increased degree of substitution at the benzylic position (2d–f, 2i) or extension of the aromatic system (2l, 2p–s) enabled good to excellent yields of the unsymmetrical ethers at lower reaction temperatures and/or shorter reaction times compared to the formation of 2b. Gratifyingly, sterically hindered ethers 2d–e that can be synthetically challenging, formed in good yields, indicating that the current methodology can serve as a complementary strategy to decarboxylative routes and C–H functionalization strategies.103,104 In this context, in can be noted that compound 2i was isolated in good yield, while it was previously described to undergo elimination in the presence of TfOH.105
![]() | ||
Scheme 2 Substrate evaluation. Standard conditions: electrophile (0.25 M, 0.25 mmol scale unless otherwise specified in ESI†), nucleophile (1 M), Zr(Cp)2(CF3SO3)2·THF (0.02 M), BTF. Standard conditions (sequential addition): electrophile (0.5 mmol, added in five portions of 0.1 mmol over 4 hours), nucleophile (1 mmol, 1 M), Zr(Cp)2(CF3SO3)2·THF (0.02 M), BTF. aNMR yield (1,3,5-trimethoxybenzene as internal standard) belectrophile![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Various functional groups were found to be compatible with the catalytic conditions, including alkynes (2f–g), alkenes (2h, 2l), aryl bromides (2m), and carbocycles of various size (2d, 2r–s). While thiophene was tolerated under the reaction conditions (2o), substrates with basic nitrogen-containing heterocycles were not (see ESI†). Such groups have previously been shown to inhibit catalysis using the same zirconium complex in direct esterification.84 The selectivity of 2g:
2a and 2h
:
2a was slightly lower compared to the benchmark substrate 2b (see ESI†), which may be attributed to increased steric hindrance by the longer aliphatic chain that affects the selectivity determining step. In addition, the reduced yields of 2g and 2o can be explained by observed unidentified by-product formation from the nucleophilic alcohol under the catalytic conditions.
Secondary aliphatic alcohols performed well as nucleophiles (2n, 2q), albeit with a drop in ee between 2.5 and 14.5 percent units depending on reaction temperature and time. Glycosylation of 1b at 60 °C with 2,3,4,6-tetra-O-benzyl-D-glucopyranose was efficiently carried out to form 2t in good yield, demonstrating that benzyl ethers are stable as O-protecting groups under the applied conditions. Similarly, disaccharide 2u was formed from the same glucopyranose in good yield, highlighting the compatibility of cyclic acetals with the protocol. The reaction between β-cholestanol and p- methoxybenzyl (PMB) alcohol proceeded smoothly to form the PMB-ether 2w. Similarly, PMB-protected Boc-Ser-OMe 2v formed under the catalytic conditions, albeit in moderate yield. The use of secondary and tertiary allylic alcohols as electrophiles afforded the corresponding 2-phenethyl ethers 2l and 2x in good yields with isomerization of the double bond (Scheme 2B), while using primary allylic alcohol prenol as electrophile furnished 2x only in a moderate yield of 34% due to competing elimination. Phenols have been reported to serve as both O- and C-nucleophiles in Lewis acid catalysed direct substitution of alcohols.52,106 With our conditions, phenolic OH-groups were unreactive as O-nucleophiles, which allowed for mono-etherification of diols with full selectivity for the aliphatic alcohol (Scheme 2C). Due to low solubility of the diols in BTF, these reactions were carried out in THF.
With the phenolic OH-group in para position, ether 2y was isolated in 78% whereas the ortho substituted analogue 2y′ was isolated in 34%, likely due to the increased steric hindrance and/or increased degree of catalyst chelation. The meta-substituted analogue 2y′′ was not detected. These observations are in accordance with expected substituent effects of Hammett-type and are as such consistent with the proposed carbocationic intermediate.107
While the inductively electron-withdrawing character of the hydroxyl group governs the reactivity of the meta-substituted substrate, the benzylic carbon can be expected to benefit from resonance stabilization by the functionality when situated in ortho and para positions. Similarly, other electron withdrawing substituents on the aromatic ring were found to reduce reaction rates and yields of ether products (see ESI†). The catalytic methodology was smoothly extended to include thiols as nucleophiles, affording linear as well as sterically hindered thioethers in good to excellent yields (Scheme 3).
Similar to the observations made for ether formation, the yields of thioethers 3a and 3b increased when the electrophile concentration was reduced. Thiophenol acted as S-nucleophile to form thioethers (3b, Scheme 3), whereas phenols were unreactive as O-nucleophiles under the tested conditions (Scheme 2C).
However, in the absence of aliphatic O-nucleophiles, phenol was active as C-nucleophile in Friedel–Crafts alkylation (Scheme 4). When reacted with di- and triphenylmethanol, phenol selectively reacts in para position to form tertiary and quaternary carbon centres in the sterically hindered phenols 4a and 4b. The bulky 4-tritylphenol 4b and its substituted analogues constitute an interesting compound class of relevance for rotaxanes.108 The preference of phenol to react as a C- rather than O-nucleophile with electrophilic alcohols under the applied conditions stands in contrast to its reactivity under Pd catalysis.109 Friedel–Crafts alkylation was also feasible using N-methyl indole as nucleophile (4c, Scheme 4).
![]() | ||
Scheme 4 Friedel–Crafts alkylation of arylmethanols. Conditions: electrophile (0.25 M), nucleophile (1 M), Zr(Cp)2(CF3SO3)2·THF (0.02 M), BTF, 60 °C. |
BTF | Benzotrifluoride |
THF | Tetrahydrofuran |
Cp | Cyclopentadienyl |
ee | Enantiomeric excess |
NMR | Nuclear magnetic resonance |
HPLC | High performance liquid chromatography |
Footnotes |
† Electronic supplementary information (ESI) available: Experimental details, kinetic data, NMR (1H and 13C), HPLC chromatograms can be found in the ESI. See DOI: 10.1039/d1cy01219c |
‡ These authors contributed equally. |
This journal is © The Royal Society of Chemistry 2021 |