Sebastian
Arndt‡
a,
Birgit
Grill‡
b,
Helmut
Schwab
bc,
Georg
Steinkellner
bd,
Urška
Pogorevčnik
b,
Dominik
Weis
a,
Alexander M.
Nauth
a,
Karl
Gruber
b,
Till
Opatz
a,
Kai
Donsbach
e,
Siegfried R.
Waldvogel
*a and
Margit
Winkler
*bc
aDepartment of Chemistry, Johannes Gutenberg University Mainz, Duesbergweg 10-14, 55128 Mainz, Germany. E-mail: waldvogel@uni-mainz.de
bacib – Austrian Centre of Industrial Biotechnology, Krenngasse 37, A-8010 Graz, Austria
cInstitute of Molecular Biotechnology, Graz University of Technology, Petersgasse 14, 8010 Graz, Austria. E-mail: margitwinkler@acib.at
dInnophore GmbH, Am Eisernen Tor 3, 8010 Graz, Austria
ePharmaZell GmbH, Hochstrass-Süd 7, 83064 Raubling, Germany
First published on 8th December 2020
Levetiracetam is an active pharmaceutical ingredient widely used to treat epilepsy. We describe a new synthesis of levetiracetam by a dynamic kinetic resolution and a ruthenium-catalysed ex-cell anodic oxidation. For the enzymatic resolution, we tailored a high throughput screening method to identify Comamonas testosteroni nitrile hydratase variants with high (S)-selectivity and activity. Racemic nitrile was applied in a fed-batch reaction and was hydrated to (S)-(pyrrolidine-1-yl)butaneamide. For the subsequent oxidation to levetiracetam, we developed a ligand-free ruthenium-catalysed method at a low catalyst loading. The oxidant was electrochemically generated in 86% yield. This route provides a significantly more sustainable access to levetiracetam than existing routes.
Scheme 1 Technical and auxiliary-based approaches to levetiracetam (selected key-intermediates). X = Cl, Br. |
Synthetic route via | AE [%] | EF [kg kg−1] | Overall yield [%] |
---|---|---|---|
Comparability might be limited by scale- and optimisation effects, and by the choice of the starting material. TM = transition metal, AE = atom efficiency, EF = E-factor. | |||
Strecker synthesis | 12 | 52 | 21 |
α-Halogenation | 5 | 95 | 22 |
Evans’ auxiliary | 5 | 465 | 37 |
Sulfinimine auxiliary | 0.8 | 572 | 18 |
Ugi reaction | 38 | 1544 | 5 |
Ru catalysis | 3 | 5110 | 25 |
TM hydrogenation | 39 | 151 | 18 |
Our approach | 44 | 21 | 44 |
Academic strategies relied upon the use of chiral auxiliaries, such as oxazolidinones (Scheme 1, III) or sulfinimines.12 Those routes gave indeed higher yields, but the removal of the auxiliary and the large number of synthetic steps again resulted in a low overall efficiency. For this reason, those methods are economically prohibitive.
A particularly short method employed the Ugi three-component reaction that, given by the low number of synthetic steps, seemed very efficient at first sight (Scheme 1, IV).13 Revising the resolution,14 however, a very poor overall yield of 5% was obtained accompanied by an uneconomic use of chemicals. Moreover, toxic 2,2,2-trifluoroethanol was used as the solvent. The syntheses of the auxiliaries themselves were not taken into account, which would have impaired the metrics further.15 Catalytic consecutive double-bond migration and oxidative cleavage of (S)-2 under ruthenium catalysis was reported (Scheme 2a), wherein levetiracetam was obtained in 25% over six steps.16 The extensive use of chemicals resulted in by far the lowest sustainability of this route within our comparison. Furthermore, an unacceptable co-solvent according to pharmaceutical guidelines was used (CCl4).17 Bandichhor separately reported the oxidation of (S)-3 and obtained levetiracetam in 42% after amidation of (S)-4.18 Rhodium or cobalt-catalysed approaches were reported by Shevlin and Chirik for the enantioselective hydrogenation of dehydro-levetiracetam (5).6,19–21 Reviewing this pathway,20,21 levetiracetam was obtained in 18% yield, with 39% AE, and with an EF of 151, which still is poor for a technical scale (Scheme 2b). It should be pointed out that N-heterocyclic carbene or chiral phosphine ligands are expensive. Their molar price may be higher than that of the noble transition metal itself. An alternative electrochemical oxidation of (S)-3 was developed by Stahl et al. using a TEMPO-related mediator (Scheme 2a, step II).22 A patent elaborates on an enzymatic approach to levetiracetam and involved an (S)-selective nitrile hydratase (NHase) in a classical enzymatic resolution of (rac)-6, leading to a maximum theoretical yield of 50% after a low yielding two-step chemical synthesis of the oxo-nitrile precursor (Scheme 2c).23
Scheme 2 (a–d) Transition metal-catalysed or electrochemical approaches to levetiracetam (selected intermediates). Ox. = oxidation, sel. = selective, ee = enantiomeric excess. |
NHases exhibit a broad substrate scope for catalytic enzymatic transformations of nitriles to the corresponding amides.24 The perfectly atom economic reaction proceeds in water – which is both the reagent and solvent, under ambient conditions. Nitrile hydratases are applied on industrial scale e.g. for the production of the bulk chemical acrylamide.25 Likewise, electrolyses are highly atom-efficient, inherently safe, and environmentally benign as they do not require primary oxidants. Electricity is inexpensive, readily available, and sustainable if produced from renewable energy sources.26 Combining these two technologies offers new avenues for sustainable production processes.
In conclusion, technical routes to levetiracetam are more sustainable than ‘academic’ routes, but still suffer from low overall yields. Yet half of the material is lost in a resolution if a recycling of the wrong enantiomer is impossible. Chiral auxiliary-based routes indeed circumvent the resolution, but these methods are accompanied by a large step-count. Unfavourably, the synthesis and the removal of the auxiliary are extremely uneconomic. Catalytic routes ought to feature higher atom efficiencies and a better waste management in theory, but the elaborate syntheses of starting materials, the utilisation of chiral auxiliaries, and the employment of expensive ligands eliminate those advantages. Furthermore, they suffer from high overall costs that impede the application on a technical scale.
The aim of this work was to find both a sustainable and economical route for the synthesis of levetiracetam. From a retro-synthetic point of view, we envisioned a Strecker reaction with the benefit of low-priced and readily available starting materials. As opposed to previous routes, the pyrrolidine moiety shall be directly installed saving one synthetic step. The Strecker reaction is telescoped by an enzymatic dynamic kinetic resolution of 2-(pyrrolidine-1-yl)butanenitrile (rac)-7 to the corresponding amide (S)-8. Finally, we envisioned an electrochemical stereo-conservative and regiospecific oxidation of (S)-8 in position α to the amine (Scheme 2d).
Screening a panel of known and new NHases,30 we identified the thermotolerant cobalt-dependent Comamonas testosteroni NHase (CtNHase) as the most promising enzyme candidate (ESI Table S1†).31 It showed activity from pH 6.0 to 9.5 (ESI Fig. S2†), hydrolysed 7 in the temperature range between 25 and 50 °C (ESI Fig. S3†), and showed the highest (S)-selectivity among the candidates with 89% ee. At 10 mM (rac)-7 substrate concentration the product was detected in 1.1 mM concentration (ESI Table S1†). This small amount of detected product indicated a strong need for a deeper understanding of the limitations of the reaction system.
We investigated potential product inhibition and found that up to 25 mM, no inhibition by (rac)-8 was obvious. Aminonitrile 7 is in equilibrium with its three components propanal, pyrrolidine and hydrocyanic acid in aqueous solution, which is the key requirement for a dynamic kinetic resolution. The equilibrium is pH dependent, with high pH values favouring the formation of 7. On the one hand, the amounts of pyrrolidine, propanal and cyanide in the reaction solution depend on the pH, on the other hand, they also influence the pH. Each of these components may inhibit CtNHase, and especially cyanide has been reported as an NHase inhibitor before.32 Cyanide indeed inhibited CtNHase (ESI Fig. S4†) whereas pyrrolidine and propanal did not (ESI Fig. S5†). In successive optimisation of reaction conditions including increased substrate concentrations, we aimed for minimal cyanide concentrations by pH control and concomitantly adapted buffer capacities in analytical scale reactions. In doing so, product levels could be increased to more than 40% at 50 mM substrate load, and 15% at 200 mM substrate concentration while product ee ranged from 70–90% (ESI Fig. S6†). Time resolved monitoring of 7 hydration showed that the reaction was very fast, but ceased within a short period of time (ESI, Fig. S7†). The reason was not fast deactivation of CtNHase,30,31 but the formation of a side product: propanal and cyanide form 2-hydroxybutanenitrile – a competitive substrate for CtNHase which gives 2-hydroxybutanamide upon hydration, as confirmed by HPLC-MS (ESI Fig. S8†). Formation of (S)-8, enantioselectivity and the substrate selectivity needed to be increased to create a potent biocatalyst suited for industrial application. We therefore focused on CtNHase on the molecular level, docked (S)-7 into the active site of a model (ESI Fig. S11†) and identified amino acid residues within 4 Å around the docked compound as potential targets to shape the active site for improved (S)-7 binding. Putative Cobalt-binding residues, strictly conserved residues in the protein family and Arg52 in the β-subunit (a residue likely involved in proton transfer33) were excluded.
We created site saturation libraries for ten positions: Q93, W120, P126, K131 and R169 of the α-subunit and M34, F37, L48, F51 and Y68 of the β-subunit. To screen for improved variants, a colorimetric assay was established. Specifically, amides can be transacylated by amidases in the presence of hydroxylamine to give hydroxamic acids, which can then be visualised as iron-complex (Scheme 3).34 Using an (S)-selective amidase in the screening assay provides additional selection for (S)-8 and strictly (S)-selective Rhodococcus erythropolis amidase35 (ReAmd) was chosen for this step. The 10 libraries were screened in 96 deep-well plates. As an example, the results for position β-L48 are shown in the ESI (Fig. S12†).
Cell suspensions of the most promising variants from this screening were used for biotransformation reactions (ESI chapter 5.1 and 5.2†) to simultaneously determine concentration of 8 and the corresponding ee by chiral HPLC. In this re-screening step, false positives generated by potentially formed 2-hydroxybutaneamide were unambiguously eliminated. Whereas variants in the α-subunit showed little effect, several mutants from the β-subunit gave increased formation of 8 and/or ee. Most promising hits were analysed by sequencing (for a selected example see Table S3†). Hydrophobic amino acids instead of Phe in position 51 increased product levels up to twofold with ee of up to 92% at 50 mM of (rac)-7.
Random engineering by error prone PCR was used to complement the rational protein engineering approach.36 To increase the throughput, the above described assay was established on colony level (ESI chapter 8.2†). Colonies are attached to a membrane and separated from their growth medium. First, the membranes are treated with substrate solution. In the second phase they are transferred to a solution of amidase and hydroxylammonium chloride. Finally, the red complex is developed by transfer of the membrane attached colonies to acidic FeSO4 solution. Instead of the full length sequence, four short stretches in CtNHase lining the active site were defined to decrease the library size and therefore reduce the number of variants that must be screened. Random libraries were constructed for each stretch (coloured elements in Fig. 1). Specifically, the regions were amino acid 70–110 (α1) and 120–175 (α2) in the α-, and 30–71 (β1) and 124–170 (β2) in the β-subunit, respectively. The β1 library was based on wild-type CtNHase, whereas the other three libraries were generated on β1 variant βF51L. At least 11000 clones per library were screened and 900–1700 variants per library were used for re-screening in the liquid format as described for the site-saturation libraries. Hit positions were evaluated on the 500 μL scale to determine concentration of 8 and ee by chiral HPLC. Library β1 revealed a high number of improved clones, whereby position F51 was found in 16 of 26 sequenced clones, with substitutions to either Ile, Leu or Val (Table 2, entries 3–5) alone or in combination with other amino acid exchanges or silent mutations. Also, substitution of βG54 was found 9 times (entries 6 and 7). The highest ee was found for variant βL48R (96.1%). In contrast, the α1 library consisted of many parents and only pos 110 carried an Ile in 7 of 16 clones (entry 8). This position significantly increased product levels as compared to parent βF51L and slightly the product ee. Similarly, library α2 revealed a substitution of P121 by Ser, Val or Thr beneficial for the formation of 8, whereas ees were typically lower than that of the parent βF51L (entry 9). The hits of library β2 contained mutations in numerous positions, however, only one combination with improved ee was found (entry 10). Associating most impact on CtNHase stretch β1, we screened a focused library that aimed for combinations of the three critical positions 48, 51 and 54, respectively (ESI chapter 9.3†). Substitutions in pos. 48 lead to ees up to 99%, however, product formation were strongly compromised (Table 2, entry 11). Around 60% of amide 8 were obtained (based on rac-7 addition) with βF51V variants combined with Ile, Arg or Val in pos. 54 (entry 12). Finally, combining the most promising mutations from all libraries revealed 10 variants with ees >98.5%. Most of them combined mutations of L48 and G54 (entry 14). Mutations in position 51 alone lead to both, an increase in selectivity and activity, but in combination with other mutations, activity decreased below wild-type level or ees did not reach beyond 94% (data not shown). The best variants and relevant control strains were compared in reactions which were supplemented with 1 equivalent of propanal (Fig. 2) to push the equilibrium towards (rac)-7 formation and compensate for the volatility of propanal. All variants displayed ee >98%. The double mutant αP121T/ βL48R showed two-fold increased product formation in comparison to the wild-type. Three double mutants based on the βL48P exchange reached ee 99.8% or higher. Clearly, position 48 in the β-subunit was responsible for the most significant increase in enantioselectivity, potentially due to improved contact of the amino acid side chain to (S)-7.
Entry | CtNHase | Library | Anal. yielda [%] | ee(S) [%] |
---|---|---|---|---|
a (rac)-7 (100 mM), resting cells (8.5 mg mL−1), Tris-HCl buffer (500 mM pH 7.0), 25 °C, 700 rpm, 2 h. b Cells cultivated in shake flasks. c Libraries based on mutant βF51L. | ||||
1 | Wild-typeb | — | 24 | 85 |
2 | βL48R | Random β1 | 26 | 96 |
3 | βF51I | Random β1 | 62 | 90 |
4 | βF51L | Random β1 | 30 | 91 |
5 | βF51V/ βE70L | Random β1 | 50 | 91 |
6 | βH53L/ βG54V | Random β1 | 24 | 92 |
7 | βN43I/ βG54C | Random β1 | 45 | 90 |
8 | αV110I/ βF51L | Random α1c | 57 | 93 |
9 | αP121T/ βF51L | Random α2c | 65 | 90 |
10 | βF51L/ βH146L/ βF167Y | Random β2c | 50 | 95 |
11 | βL48R/ βF51I/ βG54I | β1-Focused | 9 | 98 |
12 | βF51V/ βG54V | β1-Focused | 66 | 94 |
13 | βP121T/ βL48Rb | Combined | 42 | 98 |
14 | βL48R/ βG54Vb | Combined | 31 | 99 |
On the preparative scale, (rac)-7 was hydrated to (S)-8 in pH controlled stirred 50 mL volume fed-batch reactions, using whole cell biocatalysts (ESI chapter 10†). Variant P121T/L48R showed the highest concentration of the desired product on analytical scale and hydrated (rac)-7 (1.3 g substrate in total) to (S)-8 (1.1 g, 95% ee) on a 50 mL scale. The more selective variant L48R/G54V gave (S)-8 in 98% ee (0.78 g at 1.21 g substrate load). The yields of (S)-8 were defined by the degree of cyanohydrin hydration (Scheme 4) and reached 73.3% analytical yield (ESI Fig. S14 and 15†) in case of variant P121T/L48R and 61.7% in case of variant L48R/G54V, respectively.
For the final step towards levetiracetam, a direct electrochemical oxidation was envisaged first. Pyrrolidines have been electrolysed in Shono-type oxidations,37 whereby pyrrolidones may be obtained in water as reaction medium.38 Stereo-conservative Shono oxidations have been reported for sulphonamides.39 Cyclic voltammetry indicated an irreversible oxidation of (S)-8 in the desired α-position of the amine. Up to 29% of 1 was determined by GC in optimised batch-electrolyses (detailed in experimental section in ESI, chapter 17†).40 However, chiral HPLC revealed a racemisation of the product during electrolysis. We reasoned that an N-radical intermediate must be formed during the electrolysis causing the racemisation. This prompted us to pursue an ex-cell electrochemical oxidation.
The ex-cell approach was investigated utilising ruthenium catalysts. Oxo ruthenium-catalysed α-oxidation of N-protected pyrrolidines generally feature mild reaction conditions, short reaction times and high selectivity.41 More importantly, RuO4 is known to perform two consecutive concerted oxidation steps that should result in the retention of the absolute configuration.42 Reported favourable conditions are solvent mixtures of CCl4/MeCN/H2O, a low to neutral pH, and low temperatures. Periodate is often used as terminal oxidant.43 The catalysis was performed by dissolving the hydrated ruthenium(IV) oxide (RuO2) and the desired oxidant. The mixture was suspended until the ruthenium tetroxide (RuO4) was formed as indicated by the appearance of a pale-yellow colour. Subsequently, (S)-8 was added as a solid and the reaction progress was monitored by GC. After complete reaction, the yield was determined either by GC using caffeine as an internal standard or the product was isolated by flash column chromatography. Best results were obtained for hydrated RuCl3 or RuO2 pre-catalysts and sodium metaperiodate. The catalytic reactions were performed in mixtures of aqueous and organic solvent such as dimethylformamide, acetonitrile, or acetone (Table 3, entries 2–4). Ethers, benzene, alcohols and tetrachloromethane were not investigated for safety- and regulatory reasons.44 A neutral pH was critical, and an excess of periodate (≥2.6 eq.) as well as a low temperature (0 °C, Table 3, entries 5 and 6) were beneficial. The product (S)-1 was formed in 66% yield and it was isolated in 49% yield and 99.6% ee. A substantial loss of material was observed in the gas chromatograms and in the mass balance after flash column chromatography. To rule out misleading GC data, the intermediate 9 and the potential over-oxidation product (S)-2-(2,5-dioxopyrrolidin-1-yl)butanamide were isolated or synthesised, respectively, but the clear appearance of their GC signals indicated a sufficient GC-stability and response factor. Accordingly, ring opening- and/or polymerisation pathways of the unstable hemi-aminal 9 intermediate were hypothesised.45 To prevent those side reactions, a biphasic solvent system with ethyl acetate/water was investigated similar to the protocols of Hamed et al. and Sashida et al.46 Up to 76% of the product was obtained even at room temperature and at a low amount of organic solvent (Table 3, entry 7), which is favourable from a technical point of view. We reasoned that the RuO4 is transferred into the organic layer where, due to the low water content, the intermediate might be protected from side reactions. The more polar product in turn is transferred to the aqueous phase where it might be protected from over-oxidation. The reaction was scaled up 10-fold obtaining a reproducible yield of 74% [Table 3, entry 8, 1.0 g (S)-8]. The RuO2 was efficiently recovered by filtration over aluminium oxide and the iodate was recovered quantitatively by crystallisation with methanol (up to 96% isolated yield).
Entry | Solvent | [v%] | T [°C] | GC yield [%] | |||
---|---|---|---|---|---|---|---|
8 | 9a | 1 | ∑ | ||||
a Relative intensity vs. standard, eq. = equivalents, T = temperature, DMF = dimethylformamide, isol. = isolated yield. Screenings were performed on a 100 mg scale of (S)-8. | |||||||
1 | H2O | — | 20 | 5 | 9 | 31 | 45 |
2 | DMF | 33 | 20 | 16 | 18 | 37 | 72 |
3 | MeCN | 33 | 20 | 1 | 4 | 57 | 62 |
4 | Acetone | 33 | 20 | 1 | 4 | 60 | 65 |
5 | MeCN | 50 | 0 | 1 | 2 | 65 | 68 |
6 | MeCN | 70 | 0 | 0 | 0 | 66 | 66 |
7 | EtOAc | 33 | 20 | — | — | 76 | Isol. |
8 | EtOAc | 33 | 20 | — | — | 74 | Isol. |
Next, an electrochemical recycling method was developed for the re-oxidation of the iodate generated in the ruthenium catalysis. High-grade periodate is expensive and generates substantial amounts of waste. Therefore, the recycling decreases the E-factor and the costs of the process. The direct oxidation of common iodide was recently published by our group.47 Based on our previous research, the iodate was oxidised in caustic soda. Water is generally considered as an abundant, environmentally benign and non-hazardous solvent.48 A divided electrolysis cell was used equipped with a Nafion membrane, a stainless steel cathode, and a boron-doped diamond (BDD) anode.40 Notably, BDD anodes are innovative and high performance metal-free electrodes, which are considered to be sustainable since they are made from methane.49 In anodic electro-organic synthesis they proved to be superior and technically viable as compared to other electrode materials.50 The conditions were optimised by a statistical screening approach,51 whereby a robust process was found. Optimum conditions were found for an applied charge of Q = 3 F, a current density of j = 10 mA cm−2, and a hydroxide/iodate concentration of 1 M/0.21 M (Table 4, entries 1–3). Paraperiodate (H2IO63−) was obtained in 86% yield and was converted to metaperiodate (IO4−) by acidic recrystallisation (up to 71% isolated yield).52 The electrochemically obtained periodate was tested in the Ru-catalysed oxidation to give reproducible results. Finally, the iodate electrolysis was scaled up in flow electrolysis.53 The electrolysis conditions were adjusted to an increased current density of j = 100 mA cm−2 (Table 4, entries 4–7) and an applied charge of Q = 4 F (Table 4, entries 8–10). The product was obtained in 78% yield, which corresponds to 48 g of paraperiodate.
Entry | V [mL] | Q [F] | j [mA cm−2] | C(NaOH) [M] | LC-PDA yield [%] | ||
---|---|---|---|---|---|---|---|
IO3− | IO4− | ∑ | |||||
BBD = boron-doped diamond, VA = stainless steel, rt = room temperature. Reactions were performed with an initial concentration of iodate of 0.21 M. | |||||||
1 | 6 | 3 | 2 | 3 | 20 | 83 | 103 |
2 | 6 | 3 | 10 | 1 | 16 | 83 | 99 |
3 | 6 | 3 | 10 | 1 | 11 | 86 | 97 |
4 | 50 | 3 | 10 | 1 | 20 | 85 | 105 |
5 | 50 | 3 | 50 | 1 | 26 | 73 | 99 |
6 | 50 | 3 | 100 | 1 | 34 | 63 | 97 |
7 | 50 | 3 | 300 | 1 | 63 | 37 | 100 |
8 | 500 | 3 | 100 | 1 | 29 | 70 | 99 |
9 | 1000 | 3 | 100 | 1 | 37 | 72 | 109 |
10 | 1000 | 4 | 100 | 1 | 18 | 78 | 96 |
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0gc03358h |
‡ Contributed equally. |
This journal is © The Royal Society of Chemistry 2021 |