Caio M. Pachecoa,
Wesley Limaa,
Fernanda A. Limaa,
Mauro R. B. P. Gomeza,
Isabela G. da Silvaa,
Leandro S. M. Mirandaa,
Pierre M. Estevesb,
Ivaldo Itabaiana Jrc,
Robert Wojcieszakde,
Raquel A. C. Leãoa and
Rodrigo O. M. A. de Souza*a
aBiocatalysis and Organic Synthesis Group, Chemistry Institute, Federal University of Rio de Janeiro, Rio de Janeiro, 21941910, Brazil
bInstituto de Química, Universidade Federal do Rio de Janeiro, Av. Athos da Silveira Ramos, 149, CT, A-622, Cid. Univ., Rio de Janeiro, RJ 21941-909, Brazil
cDepartment of Biochemical Engineering, School of Chemistry, Federal University of Rio de Janeiro, Rio de Janeiro, 21941909, Brazil
dUniv. Lille, CNRS, Centrale Lille, Univ. Artois, UMR 8181 – UCCS – Unité de Catalyse et Chimie du Solide, Lille, France
eUniversité de Lorraine, CNRS, L2CM UMR 7053, Nancy, F-54000, France. E-mail: souzarod21@gmail.com
First published on 29th October 2024
The global imperative to shift towards renewable and sustainable resources has spurred significant interest in exploring and utilizing platform chemicals derived from renewable feedstocks. Among these, levoglucosenone (LGO) and Cyrene™ have emerged as promising candidates. LGO, derived from the pyrolysis of cellulose and hemicellulose, exhibits structural versatility, making it an attractive starting material for various valuable products. Its chemical transformations can yield a diverse array of derivatives, including levulinic acid, furan derivatives, and intermediates for pharmaceutical and agrochemical synthesis, as well as bio-based materials such as bioplastics and resins. Cyrene™, produced through the hydrogenation of LGO, serves as a renewable, biodegradable, and non-toxic dipolar aprotic solvent, offering sustainability advantages for green chemistry applications. Herein we report our results on the continuous-flow cascade transformation of LGO into Cyrene™ and then (R)-γ-carboxy-γ-butyrolactone in good yields with an additional mechanistic investigation for the Baeyer–Villiger oxidation of Cyrene™.
Among the possibilities, the synthesis of (R)-γ-carboxy-γ-butyrolactone (4), which is an intermediate for the synthesis of substituted lactones with potential biological activity is of growing interest.19–22 Herein we report our efforts on (R)-γ-carboxy-γ-butyrolactone (4) synthesis starting from LGO (1). Our group has envisioned a strategic approach that begins with continuous-flow hydrogenation of levoglucosenone (LGO, 1), resulting in the production of Cyrene™ (2). Subsequently, a cascade continuous-flow Baeyer–Villiger/oxidation protocol is employed, culminating in the production of (R)-γ-carboxy-γ-butyrolactone (4) (Scheme 1). The Baeyer–Villiger reaction of Cyrene™ (2) has been previously documented by Bonneau and colleagues.18 However, the unusual product from this transformation prompted us to undertake a theoretical investigation of the reaction mechanism.18 This study was conducted to gain a deeper understanding of this reaction for further developments.
Entry | T (°C) | Rt (sec) | Conv. (%)b | Selectivity (%)b | |
---|---|---|---|---|---|
Cyrene™ (2) | 5 | ||||
a Reaction conditions: A solution of levoglucosenone (LGO, 1) in water (0.1 M), 250 mg of Pd/C 10% catalyst cartridge (Thales Nano), 200 μL internal volume, 15–30 bar H2.b Conversion (%) and selectivity (%) values based on GC/FID absolute area analysis (chromatographic profile: Fig. S4 to S9, see ESI). | |||||
1 | 40 | 100 | >99 | >99 | — |
2 | 50 | >99 | >99 | — | |
3 | 60 | 100 | >99 | 95 | 5 |
4 | 50 | >99 | 87 | 13 | |
5 | 80 | 100 | >99 | 40 | 60 |
6 | 50 | >99 | 50 | 50 |
The condition with maximized Cyrene™ (2) production were observed at 40 °C and a flow rate of 1.0 mL min−1 (entry 2, Table 1). Under these conditions, Cyrene™ (2) was produced with a selectivity of 99% after only 50 seconds (entry 2, Table 1). When the temperature was raised to 60 and 80 °C, the formation of a single side-product was observed, and its appearance was greatly affected by increasing temperatures (entries 3–6, Table 1). The single side product observed was characterized as the 1,4-addition of water product 5, previously reported by other groups,23 and NMR spectroscopy (Fig. S37–S40 in the ESI†) has confirmed the chemical structure.
Further reduction of the residence time at 40 °C did not permit to reach full conversion limiting the reaction to a 50 seconds reaction. The limiting concentration that could be used for this reaction has also been evaluated and reproducible results (>99% conversion and selectivity) were achieved until 1.5 M at 40 °C and 100 seconds residence time.
With Cyrene™ (2) in hands we proceeded to the Baeyer–Villiger protocol developed by Bonneau and co-workers.18 While the Baeyer–Villiger protocol has been efficiently optimized under batch conditions, its translation to a continuous-flow system necessitates the adaptation of specific reaction parameters, especially reaction temperature, to fully leverage the potential benefits of process intensification. We initially focused on increasing reaction temperature over different residence times to observe the effect on reaction conversion. Then we tested the limiting concentration of the starting material (Table 2). In batch conditions, the Baeyer–Villiger protocol employs a three-temperature setup, beginning with the addition of hydrogen peroxide at 0 °C, followed by a temperature increase to 50 °C for lactone formation followed by another temperature increase to 90 °C for the decomposition of any remaining H2O2 and deformylation of the final intermediate towards 2H-HBO (3) (see batch conditions Section 1.1 in the ESI†). Due to the rapid heat transfer characteristic of continuous-flow reactors, we opted to forego the three-temperature approach, as depicted on Table 2.
Entry | Rt (min) | T (°C) | Conv. (%)b |
---|---|---|---|
a Reaction conditions: solution A: Cyrene™ (4.61 mL, 45 mmol) in 2.84 mL of H2O (6.04 M), solution B: H2O2 30% (v/v) (5.57 mL, 54 mmol, 1.2 eq) in 0.715 mL of H2O (8 M). The solutions were pumped with a Syrris® syringe pump, mixed in a T connector, and passed through a 12 mL PTFE coil reactor with 1/16′ inner diameter, the line was pressurized with a back pressure regulator set at 4.0 bar.b Analysis by CG/FID absolute area (chromatographic profile: Fig. S11–S18, at ESI). | |||
1 | 30 | 100 | 99 |
2 | 15 | 99 | |
3 | 7.5 | 94 | |
4 | 3.75 | 72 | |
5 | 15 | 150 | 99 |
6 | 7.5 | 99 | |
7 | 3.75 | 99 | |
8 | 5 | 99 |
Table 2 demonstrates that at 100 °C, 7.5 and 3.75 minutes of residence times can result in moderate to good conversions (entries 3 and 4). A further increase in reaction temperature can enhance the reaction conversion at 3.75 minutes of residence time and increase reaction concentration to 6.03 M without sacrificing efficiency (entries 7 and 8). Product 3 can be easily isolated without compromising reaction yield. At 100 °C, where we found 72% conversion, we could identify the final formylated intermediate (FBO) as a single by-product. While the FBO intermediate is known in the literature,18 no mechanistic investigation has been conducted to shed light on this unusual Baeyer–Villiger reaction.
The next step involves the oxidation of 2H-HBO (3) to form a carboxylic acid derivative 4 (see batch conditions Section 1.2 in the ESI†). To optimize this step, we chose to use the crude reaction product from the previous step to mimic the reaction contents as closely as possible, facilitating the implementation of a cascade reaction. First, we decided to explore the residence time based on a standard protocol that uses TEMPO (0.075 eq), NaBr (0.8 eq), and NaOCl (1.5 eq) in order to optimize the production of 4. Results are presented in Table 3.
Entry | Rt (min) | Conv. 4 (%)b |
---|---|---|
a Reaction conditions: solution A: 2H-HBO (3) (610 μL, 6.46 mmol) in 2 mL of H2O (3.2 M); solution B: TEMPO (0.075 g, 0.075 eq), NaBr (0.531 g, 0.8 eq), NaClO (650 μL, 1.5 eq), 40 μL of NaOH solution 1 M in 2 mL of H2O (pH: 10). The solutions were pumped with a Syrris® syringe pump, mixed in a T connector and passed through a 12 mL coil reactor with 1/16′ inner diameter a temperature of 25 °C, the line was pressurized with back pressure regulator set at 4.0 bar.b Analysis by CG/FID absolute area (chromatographic profile: Fig. S20–S24, at ESI). | ||
1 | 60 | >99 |
2 | 30 | >99 |
3 | 20 | 84 |
4 | 15 | 68 |
5 | 5 | 31 |
Based on this initial evaluation, we could observe that short residence times could be achieved by standard TEMPO oxidation protocol, where after only 20 minutes, a good conversion could still be obtained (entries 1–3, Table 3). Further decrease in the residence time has a dramatic effect on reaction conversion (entries 4 and 5, Table 3). Although 30 min of residence time provided excellent conversion, we decided to optimize the reaction conditions, keeping the residence time at 20 minutes and improving the reaction efficiency through this operational window. The oxidation step was optimized based on the number of equivalents of TEMPO, NaBr, and NaOCl. Results are presented in Table 4.
Entry | TEMPO (eq) | NaBr (eq) | NaOCl (eq) | Conv. 4 (%)b |
---|---|---|---|---|
a Reaction conditions: solution A: 2H-HBO (3) (610 μL, 6.46 mmol) in 2 mL of H2O (3.2 M); solution B: TEMPO (0.075 g, 0.075 eq), NaBr (0.531 g, 0.8 eq), NaClO (866 μL, 2.0 eq), 40 μL of NaOH solution 1 M in 2 mL of H2O (pH: 10). The solutions were pumped with a Syrris® syringe pump, mixed in a T connector and passed through a 12 mL PTFE coil reactor with 1/16′ inner diameter a temperature of 25 °C with a residence time of 20 minutes, the line was pressurized with a back pressure regulator set at 4.0 bar.b Analysis by CG/FID absolute area (chromatographic profile: Fig. S25–S31, at ESI). | ||||
1 | 0.075 | 0.8 | 1.5 | 84 |
2 | 0.075 | 0.8 | 2.0 | >99 |
3 | 0.11 | 0.8 | 1.5 | 93 |
4 | 0.04 | 0.8 | 1.5 | 84 |
5 | 0.075 | 1.2 | 1.5 | 66 |
6 | 0.075 | 0.8 | 0.7 | 68 |
7 | 0.05 | 1 | 2 | 78 |
As observed in Table 4, decreasing the amount of TEMPO or/and increasing the number of equivalents of NaBr has a deleterious effect on reaction conversion (entries 4, 5 and 7, Table 4), as well as the reduction in NaOCl equivalents (entry 6, Table 4). In our experiments, we were able to note that a slightly increase in the amount of NaOCl, from 1.5 to 2.0 equivalents, led us from 84% of conversion to >99% (entries 1 and 2, Table 4). It is important to note that the increased amount of TEMPO did not have the same effect where since conversion has raised to 93% (entry 3, Table 4). Isolation of the product to establish the isolated yield was performed as mentioned in the experimental section (see Experimental section). The optimized conditions obtained for the oxidation step were then used to couple both steps on a cascade protocol (Scheme 2).
We initiated the process by producing Cyrene™ (2) via hydrogenation of LGO (1) under mild conditions (using water as a solvent at 40 °C) and a very short residence time. Under these conditions, we achieved complete selectivity for the desired product. However, the first step is limited by the use of 0.1 M of LGO (1), whereas steps 2 (Baeyer–Villiger) and 3 (oxidation) can be conducted at significantly higher concentrations, up to 6.03 M. This disparity directly impacts the space-time-yield of the process. Subsequently, a 6.03 M solution of the synthesized Cyrene™ (2) was pumped through the cascade reactor. Initially, it reacted with H2O2 at 150 °C via the Baeyer–Villiger protocol to produce 2H-HBO (3), which then was cooled down in a heat exchange coil to react with the oxidation reagents and deliver carboxylic acid 4 in quantitative isolated yields, the whole cascade process led to a productivity of 5.11 g h−1.
Aiming to better understand this unusual Baeyer–Villiger reaction, we investigated the reaction mechanism both by DFT calculations (IEFPCM(H2O)/M06-2X/6-311++G(d,p)) and some control experiments (see Computational details for further information). In addition to the unusual Baeyer–Villiger rearrangement, the Cyrene™ (2) was also reported to participate in an abnormal Beckmann rearrangement, yielding a nitrile as the major product, instead of the usual amide.24 In the case of this abnormal Beckmann rearrangement, such behavior is reported when a carbocation stabilizing group is present anti to the departing leaving group. We hypothesized that such a stabilizing effect may also be the reason for the observed lactone formation in the BV reaction of Cyrene™ (2). This mechanistic rationale is depicted in Scheme 3.
Scheme 3 Pictorial representation of a mechanistic proposal for the Baeyer–Villiger rearrangement of Cyrene™ (2). |
In this mechanistic hypothesis, after the classical Baeyer–Villiger step, the presence of an acetal adjacent to the ester would stabilize the carbocation, which leads to the carboxylate intermediate 10, resulting then in the observed product after a 5-exo-tet cyclization of 10.
To theoretically evaluate this hypothesis, the geometries of reactant, intermediates, transition states, and products were optimized. Neutral and protonated reactants were considered in the calculations, as well as the explicit solvation of water molecules.
According to the theoretical model used to evaluate the reaction coordinate, under neutral conditions, the product of the addition of hydrogen peroxide to Cyrene™ is a slightly exergonic process when compared to the pre-association complex 6 between hydrogen peroxide and Cyrene™. Interestingly, intermediate 7 leads to intermediate carboxylate 10 in a highly exergonic process that does not pass through intermediate 8 but through 9. Such a sequence of events represents a Grob-like fragmentation24,25 of the addition intermediate 7 instead of the classical Bayer–Villiger mechanism which would lead to 8. From the kinetic point of view, the most significant barrier is the one related to the ring-opening step with a reaction barrier (ΔGTS1≠) of 34.2 kcal mol−1. In this transition state (TS2) the interatomic distance between the carbocation carbon and the carbonyl carbon is 2.009 Å (Fig. 3). The observed lactone product 11 is then derived from intermediate 10 through a 5-exo-tet cyclization (evidence by NMR 1H spectra in Fig. S44†), kinetically favored over the concurrent 6-exo-tet cyclized product 12 (Fig. 1).
Fig. 1 Reaction mechanism profile, with optimized geometries at IEFPCM(H2O)/M06-2X/6-311++G(d,p) level, for the species shown in Scheme 3. |
The influence of acid catalysis and explicit solvent assistance effects on the fragmentation step was then investigated. The presence of an acid leads to the elimination of water instead of the hydroxide anion, dramatically decreasing the reaction barrier, with a reaction barrier (ΔG≠TS1) of only 1.1 kcal mol−1, with a much earlier transition state, Fig. 2b. The presence of explicit solvation water molecules also impacts this step, however, in different directions in the acid-catalyzed and non-catalyzed reactions. The presence of explicit solvent molecules in the non-catalyzed process decreases the reaction barrier of TS1 to ΔGTS1≠ = 30.9 kcal mol−1), probably through solvation (and stabilization) of the hydroxide leaving group, which increases its nucleofugacity, Fig. 2a. The interaction of the explicit solvent molecules also impacts the acid-catalyzed reaction, where the reaction barrier is now increased from ΔGTS1≠ = 25.9 kcal mol−1 to ΔGTS1≠ of 12.7 kcal mol−1, a result of the increased stabilization of the intermediate 3 due to solvation, Fig. 2c. These reaction barriers in the presence of explicit solvent molecules represent a more realistic scenario once a reaction barrier of only 1.1 kcal mol−1 is highly unlikely.
The IBO analysis26,27 of the rearrangement step shows that the C–C bond density flows in order to displace the leaving group, confirming the BV nature of the transformation. However, the reactive intermediate, an oxygen-stabilized carbocation (or oxonium ion), is formed and can undergo internal rotations that afford the different products, formed by an intramolecular substitution, Fig. 3.
Fig. 3 IBO analysis of the BV step, showing the progressive flow of electron density from the CC bond (blue) to the oxygen group from the leaving OH. |
In conclusion, we have optimized two crucial steps in the synthesis of (R)-γ-carboxy-γ-butyrolactone (4), beginning with the hydrogenation of levoglucosenone (LGO, 1) to produce Cyrene™ (2), followed by the Baeyer–Villiger (BV) protocol and subsequent oxidation steps. The hydrogenation step, employing continuous-flow conditions with the H-Cube Mini reactor achieved a remarkable 99% selectivity towards the desired product after only 50 seconds. Upscale results were attained to 1.5 M concentration at 40 °C and 50 seconds residence time, providing valuable insights for scalability. Regarding the Baeyer–Villiger protocol, increasing reaction temperature to expedite ester hydrolysis, led us to achieve moderate to good conversions at 100 °C with varying residence times and concentrations. In the final oxidation step, we optimized conditions based on the TEMPO oxidation protocol, adjusting the equivalents of reagents to maximize conversion. By balancing reagent ratios, we achieved over 99% conversion, highlighting the importance of fine-tuning reaction parameters for optimal outcomes.
Furthermore, our mechanistic investigations shed light on the unusual Baeyer–Villiger rearrangement observed during 2H-HBO (3) synthesis, providing valuable insights into the mechanism for product formation. Through computational modeling, we proposed a mechanism for the formation of lactone product 3, underscoring the significance of understanding reaction mechanisms for rational design and optimization.
Flow synthesis was performed using a 3-module setup: a syringe pump utilizing a set of 1 and 0.5 mL syringes; a PTFE coil reactor with 12 mL volume and 1/8′′ OD; and a back pressure regulator.
NMR spectra were recorded on a Bruker 400 MHz spectrometer in a deuterated solvent. The values of chemical shifts (δ) are expressed in ppm with reference to tetramethylsilane (TMS) for 1H NMR and chloroform signal for 13C NMR. Coupling constants (J) are expressed in hertz (Hz). Gas chromatography analyses were carried out in GC-MS using Shimadzu GC 2010 (gas chromatography/flame-ionization detection/helium was used as carrier gas) and GC/FID using Shimadzu GC 2010 – DB5 HT column, 30 m, with automatic injector.
Scheme 4 Geometries of the species involved in the reaction, computed at IEFPCM(H2O)/M06-2X/6-311++G(d,p) level. Selected geometric parameters are shown, when appropriate. |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ra06403h |
This journal is © The Royal Society of Chemistry 2024 |