Open Access Article
Peter
Sagmeister
ab,
Jason D.
Williams
ab,
Christopher A.
Hone
ab and
C. Oliver
Kappe
*ab
aCenter for Continuous Synthesis and Processing (CCFLOW), Research Center Pharmaceutical Engineering (RCPE), Inffeldgasse 13, 8010 Graz, Austria
bInstitute of Chemistry, University of Graz, NAWI Graz, Heinrichstrasse 28, A-8010 Graz, Austria. E-mail: oliver.kappe@uni-graz.at
First published on 8th April 2019
Currently, the monitoring of multistep continuous flow processes by multiple analytical sources is still seen as a resource intensive and specialized activity. In this article, the coupling of a modular microreactor platform with real-time monitoring by inline IR and NMR, in addition to online UPLC, is described. Using this platform, we rapidly generated experimental data (17 iterations in under 2 hours) to access information on the different chemical species at multiple points within the reactor and to generate process understanding. We highlight the application of the platform in the optimization of a multistep organolithium transformation. The optimized continuous flow conditions were demonstrated in a scale-out experiment with in-process monitoring to afford the desired product in 70% isolated yield and provided a throughput of 4.2 g h−1.
There is a current paradigm shift in the pharmaceutical industry from traditional batch manufacturing to continuous processing for the preparation of active pharmaceutical ingredients (APIs),5 which has been supported by the increased availability and implementation of PAT.6 Recent reviews have highlighted the successful implementation of PAT within continuous flow environments.7 PAT has been utilized within continuous flow systems for enabling feedback loops for process control,8 black box optimization such as Design of Experiments (DoE) and automated self-optimization,9 kinetic model discrimination and parameter estimation,10 and to guide drug discovery programmes.11 The most commonly used real-time inline PAT tools are UV/Vis,12 Raman,13 IR,14 and NMR.15 Recently, their use has become more popular, partly due to increased commercial availability of benchtop devices and flow-through cells. In addition, the utilization of integrated chromatographic analysis techniques, such as high/ultra performance liquid chromatography (HPLC/UPLC),16 and GC,17 within continuous flow platforms has also been described.18 However, in these instances additional sample preparation is often required before sample injection. It is still commonplace to perform all of these analytical techniques offline, within both academic and industrial laboratories.19 In these cases, time is wasted with sample preparation and the analyzed sample is not always representative of the reaction performance. Processes which utilize real-time analysis will lead to a faster and often more reliable process optimization as compared to cases using only offline analysis.
The aforementioned examples demonstrate that the integration of a single analytical instrument within a continuous flow system for analysis of a single-step reaction is now well established within the continuous processing community. However, a key benefit of flow chemistry is the ability to perform complex multistep transformations within an integrated continuous flow system without manual intervention.20 Currently, there is limited precedent for combining different reactor modules and analytics without interruption into a fully integrated continuous process. Continuous flow technology is often inflexible and compatibility issues arise between different steps. Reaction analysis at multiple points within a continuous flow system by different techniques is vital for expeditious parameter optimization and effective control over process performance to fully realize multistep transformations. Existing studies of multistep continuous flow transformations generally use different PAT instruments in separated continuous flow steps.21
Carefully designed continuous flow reactors have been shown to bestow significant benefits for improving reaction performance in a variety of systems; particularly those which are mixing sensitive, highly exothermic, or involve short-lived intermediates.22 Accordingly, modular reactor systems consisting of multiple discrete units are most suited to these applications, since each “module” can be optimized to suit the required chemistry and reconfigured to fit a myriad of different reaction systems.23
Significant improvements can be observed for reactions using organometallic reagents performed within microreactor systems due to their enhanced heat and mass transfer characteristics,24 but organometallic reagents can pose significant processing challenges caused by moisture sensitivity and solid formation within the reactor channels.25 In a batch reactor, this type of reaction is generally conducted at cryogenic temperatures, but flow processing can be tolerant of higher temperatures.26 Our current research interests led us to investigate the integration of multiple PAT instruments (IR, NMR and UPLC) within a modular flow reactor system for the monitoring of a challenging multistep organometallic transformation.
The continuous flow setup utilized a Modular MicroReaction System (MMRS) equipped with a Lonza FlowPlate® Lab reactor, manufactured by Ehrfeld Mikrotechnik.28 The system is designed to be reconfigured for use in multiple distinct applications, but can also be scaled up in a facile manner.28e,f Initially, a FlowPlate reactor with a minimum diameter of 0.2 mm and total volume 0.44 mL was used in this study. In the flow setup (Fig. 1b) four syringe pumps (HiTec Zang SyrDos2) were used to introduce the feed solutions. Feed solutions were comprised of: (1) 0.4 M solution of substrate 1 and internal standard (biphenyl) in THF, (2) 0.44 M solution of LDA in THF,29 (3) 0.48 M solution of 4-fluorobenzaldehyde 3 in THF, and (4) water. Substrate 1 and LDA solutions were mixed first to form the enolate 2 within the FlowPlate, with a residence time of approximately 5 seconds. Subsequently, the aldehyde 3 solution was introduced into the FlowPlate, reacting for an additional 80 seconds prior to the quench. The residence time for the quench corresponded to approximately 18 seconds. The effluent from the FlowPlate was passed through an IR flow-through cell (Ehrfeld Mikrotechnik) for inline IR analysis (ReactIR 15 DiComp probe, Mettler-Toledo). The reaction mixture was then passed through an inline flow cell (0.8 mL internal volume) for monitoring by a benchtop NMR spectrometer (Spinsolve Ultra 43 MHz, Magritek). Water was then introduced to quench the reaction mixture and the outlet mixture was collected. UPLC analysis was conducted online after the introduction of water and offline by sampling the fractionated reactor output. In the case of online analysis, a subsample of process stream was diverted and then diluted by H2O/MeCN. Subsequently an internal 6-port valve (2.6 μL injection volume) injected aliquots of the diluted stream onto the LC column. The pumps and sensors within the continuous flow setup were connected to a HiTec Zang LabManager unit and controlled by LabVision software.
Initial tests confirmed the expected suitability of the analytical device for this purpose, where the tert-butyl propionate starting material 1 was observed (C
O stretch at 1730 cm−1). Introduction of the LDA stream induced complete deprotonation, even at the shortest attempted residence time of 3.9 seconds (ESI† Table S1). A concentration vs. response curve of propionate 1 revealed its limit of quantification to be around 10% of the initial concentration used, implying a degree of confidence in the successful deprotonation of >90% starting material, throughout the course of reaction monitoring. Although the enolate 2 could also be observed (1644 cm−1, Fig. 2a), its quantification in the acquisition rate of 15 seconds allows for timely parameter correction in the event of any process disturbance.
We observed that the electrophile, 4-fluorobenzaldehyde 3, had a similar C
O stretch at 1700 cm−1 (Fig. 2b), showing partial overlap with the propionate 1 peak. Nevertheless, it was possible to treat these in isolation and obtain a concentration of both components, by comparison of the peak height versus a two-point baseline. Following reaction with electrophile 3, no additional carbonyl stretch was observed, indicating that the deprotonated intermediate was lithium enolate 4, rather than the analogous lithium alkoxide species. This had the implication that quantification of 4 could not be achieved at this stage due to multiple overlapping signals, but instead required an alternative analytical technique.
Instead, 1H NMR was utilized, with 4 scans per data point, which resulted in a relatively fast acquisition time of 43 seconds. Herein, the aldehyde peak of the electrophile 3 (9.7 ppm) was sufficiently separated from all other reaction components. Furthermore, one set of aryl peaks from the aldehyde 3 could be distinguished from those of intermediate 4, meaning that product quantification was also possible in this manner (Fig. 3), and was found to provide similar results compared to those based on the aldehyde proton (ESI† Fig. S18).
To enable sample injections in quick succession, an isocratic elution method was developed, obviating any column equilibration between injections.31 With a trifluoroacetic acid-modified solvent system, clear separation of the aldehyde starting material 3 and desired product 5 could be achieved within 1 minute. Indeed, this method was even found to be sufficient to give some separation between the two product diastereomers, bestowing an additional level of analytical detail (Fig. 4). Within the present study, the diastereomeric ratio was observed to be 1
:
1.2 for all experiments, but this analytical method could also be used in studies aiming to achieve a different ratio. The analysis time was increased to 2 minutes to allow the inclusion of an internal standard (biphenyl), detection of a less polar side product 6, and to ensure that there was no carryover between injections.
The higher sensitivity of UPLC means that it also tolerates a far lower analyte concentration compared to traditional HPLC methods.31 Despite the use of a small volume sample loop (2.6 μL, ESI† Fig. S26), significant further dilution was required to reach a suitable concentration. Method development using the autosampler showed that injections of 3 nmol were acceptable, implying that a dilution factor of ∼100 from the process stream was required. Multiple solutions to this challenge were explored (ESI† section 2.9), but the most reliable was deemed to be the use of a low volume piston (HPLC) pump to continuously sub-sample the process stream, followed by an additional dilution pump (Fig. 5). When paired with an internal standard, this method enabled accurate analysis, providing results every 2 minutes which were found to be equivalent to those obtained by offline analysis (ESI† Fig. S35). Due to the low flow rate of this subsample stream (0.02 mL min−1), the relatively small volume between the process stream and UPLC injection caused a major delay in obtaining representative data. Thus, a “flush” script was incorporated between injections (ESI† Fig. S23), to increase the flow rate for a short time and ensure that sample injections were representative of the process stream.
![]() | ||
| Fig. 5 A schematic representation of the continuous subsampling and dilution setup used for online UPLC analysis, via automated sample injections. | ||
Furthermore, analysis of pressure readings taken during operation revealed an unstable profile (ESI† Fig. S7), whereby the pressure was observed in some cases to increase by over 20 bar, due to partial channel blockages. We propose that the blockages were caused by lithium salt precipitation. In most cases, the resulting increase in pressure relieved the blockage after some time, restoring the system pressure to a normal level. Nevertheless, partial blockages are detrimental to long-term reaction stability, accelerating pump wear, but also restricting reagent flow and inhibiting mixing. Eventually, complete blockages will occur and require temporary shutdown of the process. In a manufacturing environment this type of failure could require diversion of vast amounts of material to waste (particularly during the process restart phase), resulting in supply chain complications and loss of capital.
For all subsequent experiments a reactor plate with wider channels (0.5 mm minimum diameter) and larger volume (1.63 mL) was installed (ESI† Fig. S6). Within this reactor configuration, the residence time for the deprotonation changed to approximately 11 seconds and there was an additional 99 seconds after the addition of the electrophile prior to the quench. The residence time for the quench was approximately 18 seconds. With this new reactor plate, no further blockage issues were observed, resulting in a stable pressure profile.
The process was found to operate most favorably at 20 °C. Warming to 40 °C was anticipated to improve solubility of lithium salts and reduce the likelihood of precipitation, but had a detrimental effect on product 5 formation. Operation of organometallic reactions, which usually require cryogenic conditions in batch, at closer to ambient temperature has been demonstrated to have the potential to save significant quantities of energy when considered on an industrial scale.32
Subsequently, the responses from the optimization experiments were fitted to polynomials models using a statistical experimental design software package (Modde v11). Models were fitted from the NMR and UPLC data for both aldehyde 3 and desired product 5 by using multiple linear regression (ESI† section 3). Models were generated by including all main, square and interaction terms and then non-significant terms were removed. A good fit was achieved for all the models with R2 values of 0.74 (3, NMR), 0.89 (3, UPLC), 0.93 (5, NMR) and 0.84 (5, UPLC). These models also showed a moderate to good level of predictability with Q2 values of 0.50 (3, NMR), 0.79 (3, UPLC), 0.66 (5, NMR) and 0.66 (5, UPLC). Aldehyde 3 equivalents were demonstrated to have no influence on the desired product 5 yield over the range of equivalents explored. Increasing LDA equivalents was shown to have positive influence on the product 5 yield. A negative squared term for the effect of temperature was observed, describing the increase in yield of 5 from 0 to 20 °C and the decrease from 20 to 40 °C. A model was also successfully fitted for side product formation based on using the UPLC data (ESI† Fig. S39). The models were subsequently used to explore the experimental design space (Fig. 6).
![]() | ||
| Fig. 6 Contour plot of predicted product yields from model fitted using the experimental data from offline UPLC analysis. | ||
Complete deprotonation of the ester 1 was observed across almost the entire reaction period, implying that reagent delivery was consistent with negligible influence of any flow rate variation (Fig. 7b, red line). When two brief disturbances were observed (at 29 and 32 min), a 5% decrease in the extent of ester 1 deprotonation caused a 10% reduction in yield for a period of almost 10 minutes. This significant effect caused by such a minor deviation emphasizes the sensitivity of the process to minor changes, and the power of real-time reaction monitoring.
UPLC analysis showed a 70–85% yield of desired product 5 throughout the duration of the experiment (Fig. 7b, green line), corresponding well with the predicted values from the previously generated polynomial models. NMR monitoring, however, provided a consistently higher value of 85–90% for enolate 4. This discrepancy is thought to arise from the poor differentiation between desired intermediate 4versus side product 6 by NMR. This side product arises from addition of the product enolate 4 to a second molecule of electrophile (ESI† section 4.2), and its increase throughout the reaction corresponds with a decrease in aldehyde 3. This highlights the power of UPLC as a quantitative analytical technique, for its resolution between numerous chemical species. Despite the discrepancy in absolute values, the observed trends by both techniques were essentially identical, providing further confidence in the acquired data.
The scaled-out continuous flow process was successfully operated over a 70 min time period. 4.9 g of desired product was isolated after purification by column chromatography, corresponding to a 70% yield and productivity of 4.2 g h−1.
O stretch could be clearly quantified. Subsequently, progress of lithium enolate addition to the electrophile 3 was monitored by inline NMR, using the aldehyde proton and aryl protons as distinct markers. The final reaction performance was quantified by online UPLC, which was enabled by a custom-built continuous subsampling-dilution system.
The large volume of data was exploited for rapid exploration of reaction parameters, mapping the experimental space. A good fit was obtained for all polynomial models generated, allowing an estimation of the impact that small variations in each parameter will have upon the overall reaction performance. The process was scaled out, using the integrated sensors and inline/online analytics to monitor each reaction step. After 70 min processing time, 4.9 g of material was isolated, with a productivity of 4.2 g h−1.
The developed reactor platform with integrated PAT supports a data rich laboratory environment for real-time multistep reaction monitoring; of critical importance in fulfilling the ethos of Industry 4.0 for continuous processing and the laboratories of tomorrow within the pharmaceutical and fine chemical industries. Ongoing work within our research group will further improve and expand the platform.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/c9re00087a |
| This journal is © The Royal Society of Chemistry 2019 |