Research Article  Open Access
Robin Legner, Alexander Wirtz, Martin Jaeger, "Using Compact ^{1}H NMR, NIR, and Raman Spectroscopy Combined with Multivariate Data Analysis to Monitor a Biocatalyzed Reaction in a Microreaction System", Journal of Spectroscopy, vol. 2018, Article ID 5120789, 11 pages, 2018. https://doi.org/10.1155/2018/5120789
Using Compact ^{1}H NMR, NIR, and Raman Spectroscopy Combined with Multivariate Data Analysis to Monitor a Biocatalyzed Reaction in a Microreaction System
Abstract
Process analytical technology aims at process knowledge and process improvement, efficiency, and sustainability. A prerequisite is process monitoring. The combination of microreaction systems and spectroscopy proved suitable due to dimension and compound reduction and realtime monitoring capabilities. Compact ^{1}H NMR, NIR, and Raman spectroscopy were used to monitor the biocatalyzed hydrolysis and esterification of acetic anhydride to isoamyl acetate using immobilized Candida antarctica lipase B (CALB) in a microreaction system in realtime. To facilitate the identification of signals suitable for the extraction of concentrationtime (ct) graphs, 2D heterocorrelation spectra were generated through covariance transformations applied to 1D Raman, NIR, and NMR data. By means of this purely mathematical statistical procedure, the relevant signals of the process media were assigned to educts and products and thus made applicable for univariate data evaluation. The data obtained were interpreted in terms of a firstorder kinetic model, and corresponding reaction rate constants were extracted. An alternative, elegant, and fitforautomation approach for the kinetic analysis of the spectra was demonstrated in using multivariate curve resolution (MCR). The results of the univariate and multivariate approaches were comparable with regard to reaction rates and concentrations. While the manual integration of the ^{1}H NMR spectra followed by univariate analysis allowed to establish a concentration profile of the final product isoamyl acetate hence revealing more details, multivariate analysis was found more suitable for process automation.
1. Introduction
Since the joint initiative of regulatory authorities and the pharmaceutical industry, process analytical technology (PAT) has evolved from a mere production support discipline into a fastgrowing research area [1–3]. It has established itself as an important element of industrial production processes [4]. Thus, process analytical instruments are located in or nearby largescale processes and reactors. Process understanding results from their application, which in fact is process monitoring. The process knowledge obtained aims at ensuring a constant final product quality. It further enables process control and allows process optimization. The process may be improved towards efficiency, hence cost reduction, sustainability, and safety [1, 2]. Meeting the customer demand for reproducible product quality is an important prerequisite for the competitiveness of a company, which requires using facilities, resources, and energy in the most economical way [5].
Today’s process analytical tools comprise a variety of methodologies and take advantage of their inherent features [2]. Spectroscopy, spectrometry, and chromatography represent highly specific methods [4]. Among them, Raman and nearinfrared (NIR) play major roles [6–8]. More recently, nuclear magnetic resonance (NMR) spectroscopy has grown into focus [9, 10]. During process monitoring, the information sought is most often the quantity of a compound at a given time. Therefore, spectroscopic techniques providing even relatively lower resolution and sensitivity than those commonly required for structure elucidation purposes satisfy the requirements for process monitoring. Process monitoring can proceed inline, online, atline or—less desirable—offline. For inline monitoring, a probe is inserted into the reaction vessel. When the measurement is conducted online, a bypass is used. For atline analysis, a sample is taken from the reaction and treated before analysis. Hence, in, on, and atline monitoring types usually provide results in realtime [1, 2, 7, 11].
In order to compensate for the low resolution of, e.g., NIR or compact benchtop ^{1}H NMR spectrometers, suitable multivariate data analysis procedures have been devised and applied where univariate data processing did not prove sufficient [11–13]. Another strategy for enhancing spectral resolution has been developed by using covariance data processing [14–16]. This approach has been derived from statistical mathematics. Hence, different spectra that possess a common socalled perturbation dimension are correlated, where the spectral intensity represents a function of two independent variables such as wavenumber or frequency. This type of twodimensional correlation spectroscopy has been introduced by Noda et al. for IR and also later for other vibrational and optical absorption methods [16–18]. For NMR spectroscopy, Brueschweiler et al. and Eads and Noda have reported on the use of covariance transformations to correlate frequency domains [19–21]. A 2D homocorrelation map is generated by transformation of a series of 1D spectra or a 2D spectrum from one spectroscopic technique. When two series of 1D spectra or two spectra of different spectroscopic techniques are correlated, a 2D heterocorrelation map results [22, 23]. Covariance methods are especially useful for the identification of spectral resonances suitable for process analysis.
Concentrationtime (ct) profiles from reaction monitoring can be obtained either by univariate methods or by multivariate methods such as the multivariate curve resolution alternating least squares (MCRALS) [24, 25]. Using this selfmodelling curve resolution method, reaction kinetics can be determined directly from the recorded spectra series [26–28].
Isoamyl acetate belongs to flavoring substances, its odor being similar to both banana and pear. It is the main constituent of banana oil. Due to its odor, isoamyl acetate is used as artificial flavor [29, 30]. Its application as solvent and carrier, once known as aircraft dope, has vanished except for historical reproductions and scale models since airplanes are made of metal or carbon composites [31].
In this study, compact ^{1}H NMR, NIR, and Raman spectrometers were applied to monitor in realtime the biocatalyzed esterification to isoamyl acetate via the hydrolysis of acetic anhydride using immobilized Candida antarctica lipase B (CALB) in a microreaction system, cf. Scheme 1. The recorded spectra were transformed into 2Dheterocovariance correlation maps for the identification of specific resonances for univariate kinetic analysis. The results were compared to the kinetic parameters obtained from multivariate MCRALS analysis.
2. Theory
2.1. Multivariate Curve ResolutionAlternating Least Squares (MCRALS)
Principal component nalysis (PCA) and partial least square regression (PLSR) are among the most often used chemometric tools for process data analysis [32, 33]. Principal components and latent variables usually do not represent chemical and physical quantities. This is evenly true for the MCRALS algorithm. Physicalchemical knowledge on the process can be integrated into the chemometric analysis, as the number of components may be defined through the number of reactants and products [24, 25, 34]. From a data matrix containing all spectral information as superpositions due to the reaction mixture, the information about the pure components, e.g., the spectra, can be extracted. In case of reaction monitoring, represents spectral data. The data matrix can be decomposed into the small matrices and which consist of pure spectra and the associated concentration profiles and the residual matrix as follows [25, 34–36]:
The matrix of dimensions z × s contains the recorded NMR spectra as rows with the respective intensity values. From reaction monitoring, spectral series represent a function of time. The transposed matrix of dimensions a × s contains the spectra of the pure components as rows a. In matrix of dimensions z × a, the contributions of the pure spectra in columns a correspond to the concentrations in the corresponding spectrum. The residual matrix of dimensions z × s contains the spectral residues, which cannot be explained by the pure components. Residuals can represent random noise or systematic disturbances in the spectrum.
Matrices and can be computed by MCRALS. They are iteratively calculated from a first estimate for or in an alternating leastsquares method as shown below [24, 35, 37]:where denotes the estimated value and + a pseudoinverse. The decomposition of is achieved by iterative leastsquares minimization of under suitable constraints, i.e., a chemical kinetic model, nonnegativity in spectral, and concentration profiles as well as an initial guess of reagent and product concentrations. Since the uniqueness of the solution is determined by the secondary conditions, such as the chemical model, reaction order, and pure educt or product spectra, it is important to introduce as much physicalchemical knowledge of the system as possible into the constraints [34, 38].
2.2. Homo and Heterocorrelation Spectroscopy
In contrast to NMR spectroscopy, vibrational spectroscopic techniques are not based on frequency pulse sequences that allow us to encode two frequency domains. While NMR spectroscopy is intimately associated with the Fourier transformation to generate 2D correlation spectra, a different approach was suggested by Eads and Noda and Frasinski et al. [21, 22, 39, 40] Covariance transformations were used to create a correlation map by recording onedimensional spectra as responses to a quantitative, measurable perturbation applied to the sample in a systematic way [16, 18]. This perturbation may be temperature, pressure, concentration, stress, electrical field, and most importantly a chemical reaction [17, 41]. A correlation spectrum or map, , is defined according to the following equation:where is the spectrum affected by the perturbation, called the dynamic spectrum, and and are two different spectral variables at a given time or a fixed interval of the external variable . The symbol denotes a suitable correlation function. The dynamic spectrum is computed by subtracting the average spectrum from each spectrum of a series of spectra. The quantitative comparison of spectral intensity variations can then be performed by means of statistical theory. Variations from a mean value are represented by covariances. Hence, the correlation map can be considered as the covariance matrix, when the series of dynamic spectra is represented by a matrix. Assuming that is a complex number, it can be expressed according to the following equation consisting of the real part and the imaginary part :
For spectroscopic purposes, the relation between the covariance matrix and a spectrum, in particular 2D NMR spectra, was thoroughly proven [19, 20, 42]. The proof is based essentially on Parseval’s theorem. The numerical calculation of the correlation maps has been described elsewhere [16, 20].
Due to the statistical nature of the approach, one series of spectra may be compared to another series of spectra originating from a different spectroscopic technique. This is referred to as heterocovariance spectroscopy [21, 41]. A heterocorrelation map will result according to the following equation:where and are the spectral variables from different spectroscopic techniques. Equation (6) can be rewritten in matrix notation, cf. Equation (7), where the square root is a consequence of Parseval’s theorem [19, 20, 43]:where is the 2D spectrum, is the covariance matrix, and and are the two data matrices, such as a series of spectra with a common perturbation dimension. Crosspeaks in such correlation spectra have been reported as a help to disentangle crowded spectral regions [16, 43]. In this study, 2D heterocovariance correlation maps are computed from series of 1D ^{1}H NMR, Raman, and NIR spectra recorded during reaction monitoring.
3. Materials and Methods
3.1. Spectral Recording and Processing
NMR spectra were recorded online using a compact NMR picoSpin80 spectrometer with a proton Larmor frequency of 82 MHz (Thermo Fisher Scientific GmbH, Dreieich, Germany). It had a builtin flow cell and an electrical lock such that neat liquid process media were measured without further sample preparation or addition of deuterated solvents or reference standards. The sample cell possessed a total volume of 40 µL and an active volume of 40 nL. The operation temperature was 36°C. The spectrometer was controlled by Thermo Fisher PicoSpin software 0.9.3 running on the spectrometer control board and accessed via a web interface from a laptop computer. The free induction decay (FID) was recorded as 4092 data points that were zero filled to 8 k points prior to the Fourier transformation. The number of scans amounted to 16 and the pulse length was set to 60 µs corresponding to a 90° pulse. Spectra were recorded from 0.3 to 7 ppm. The bandwidth amounted to 4 kHz. All spectral processing was achieved through MestReNova 12.0.0 (Mestrelab Research S.L., Santiago de Compostela, Spain) on a laptop computer under Windows 7. The spectrum size after zerofilling and Fourier transformation (FT) resulted in 64 k. Pretreatment methods were used for spectra optimization. The spectra were binned with 0.004 ppm and normalized to the largest peak. Baseline correction was achieved using a Bernstein polynomial fit of 3rd order.
Raman spectra were recorded using the DXR™ SmartRaman spectrometer (Thermo Fisher Scientific GmbH, Dreieich, Germany) and processed using the Omnic software, version 9.2 (Thermo Fisher Scientific GmbH, Dreieich, Germany) running on a desktop PC under Windows 7. A Raman immersion probe from InPhotonics (InPhotonics, Inc., Norwood, MA, USA) was connected via a special module for fiber optics, allowing a distance between spectrometer and reaction vessel of up to 5 m. The probe was positioned at its focal point of 5 mm above the solution in the reaction vessel. Sixtyfour spectra comprising a spectral range of 250–3500 cm^{−1} with a resolution of 5 cm^{−1} were acquired at each reaction time point. The laser was operated at 785 nm and 150 mW. For processing, the Omnic software 9.2 (Thermo Fisher Scientific GmbH, Dreieich, Germany) was used. During 1380 min, 50 spectra were recorded. All spectra were baseline corrected using an asymmetric 2nd order Huber function. The spectra were normalized to the highest band. For baseline correction and normalization, MatLab R 2017b (MathWorks, Inc., Natick, MA, USA) was used.
NIR measurements were carried out with the Antaris II FTNIR Analyzer using the software Omnic, version 8.0 (Thermo Fisher Scientific GmbH, Dreieich, Germany). The spectral range was chosen from 4000–10000 cm^{−1} at a resolution of 4 cm^{−1} and the number of scans was set to 64. The spectrometer was equipped with a NIR transflection immersion probe Falcata (Hellma GmbH & Co. KG, Muellheim, Germany). All spectra were normalized after the baseline correction using an asymmetric Huber function 3rd order, again under MatLab R 2017b.
3.2. Microreaction Assembly
Stainlesssteel microreaction components (Ehrfeld Mikroreaktionstechnik, Wendelsheim, Germany) were assembled according to Figure 1 [44, 45]. In an open storage vessel, the reactants, 20 mL of acetic anhydride (>98%, Merck KGaA, Darmstadt, Germany), 40 mL of isoamyl alcohol (>98%, Merck KGaA, Darmstadt, Germany), and traces of distilled water for the benefit of the enzyme, were mixed and transferred via a peristaltic pump Smartline 100 (Knauer GmbH, Berlin, Germany) into the microreaction assembly at a rate of 5 mL/min. The excess of isoamyl alcohol ensured complete conversion of the anhydride. Immobilization of the enzyme Candida antarctica lipase B (Novozymes A/S, Bagsvaerd, Denmark) was achieved by mixing the enzyme with the epoxy acrylate resin Purolite® ECR8205F (Purolite GmbH, Ratingen, Germany) in 50 mM phosphate buffer and stirring for 20 h at ambient temperature. The immobilized enzyme was then transferred into the cartridge of the Ehrfeld reactor F 200. The filled volume amounted to 1.8 of the 2 mL cartridge volume. The reaction was initiated through heating the reactor to 35°C. The residence time of a solution fraction, approximately 2 mL, in the reactor amounted to 2 min. The solution, then containing product, was recirculated into the storage vessel. The microreaction assembly was controlled through the software LabVision 2.10 (HiTec Zang GmbH, Herzogenrath, Germany).
(a)
(b)
The reaction mixture was diverted to the NMR spectrometer sample cell via a valve (Figure 1). Setting the valve, the flow was stopped after five seconds such that 16 spectra could be accumulated per sample. The following sample was introduced into the spectrometer via the bypass, transferring the previous one back to the vessel at the same time (cf. Figure 1). For Raman and NIR spectral recording, corresponding probes were attached a few millimeters above the surface of the reaction liquid in the storage vessel or was immersed into it (cf. Figure 2). The reaction was conducted over a period of 1380 min.
3.3. Covariance and Multivariate Data Analysis
The 2D correlation maps were generated using the program 2Dshige running on a personal computer under Microsoft Windows 7 [46]. Spectral data were exported as ASCII files from the spectral processing software listed above. To reduce the size of the data sets, 1D spectra were reduced by binning to 1 k data points for NMR, 732 data points for Raman, and 1 k data points for NIR. The data files were transformed into homo or heterocorrelation maps and displayed as contour or color plots using appropriate threshold levels.
For MCRALS analysis, NMR, Raman, and NIR spectra were automatically imported into MatLab, version 2017b, wherein the MCRALS toolbox 2 was applied [35, 47]. Three latent variables or components were selected on the basis of visual inspection of the computational results. The nonnegativity conditions for spectra and concentration profiles were set as constraints. Furthermore, a kinetic model A ⟶ B was assumed. Other models tested were A ⟶ B ⟶ C and A + B ⟶ C + D. It is noteworthy that the number of variables or components is not interchangeable with the number of real chemical components of the process system. An initial guess of the educt concentration ranged between 3.0 and 3.5 mol·L^{−1}. As the starting value, the reaction rate constant was assumed to be 2.5·10^{−3} min^{−1}. The agreement between calculated and experimental concentrationtime profile was taken as goodness of the obtained kinetic parameters.
4. Results and Discussion
4.1. Qualitative Reaction Monitoring
The enzyme catalysed acetic anhydride hydrolysis, and subsequent esterification with isoamyl alcohol was conducted in a microreaction platform. The enzyme was immobilized and located in a cartridge reactor as a heterogeneous catalyst. The process was monitored using online 1D ^{1}H NMR and inline Raman and NIR spectroscopy. Spectra were recorded during the course of the reaction over 1380 min. One aim of process analysis, process understanding, is represented here by generating ct plots and their mathematical description using chemical kinetic theory. It is straightforward to seek educt or product specific signals in spectra, which can be interpreted in terms of concentrations. Starting with the interpretation of the NMR spectra, the usefulness of covariance spectroscopy, especially heterocorrelation spectroscopy shall be demonstrated. Yet, analysts more familiar with vibrational spectroscopy might start from Raman or NIR spectra as well. Relevant educt and product signals could be identified in ^{1}H NMR spectra (Figure 2).
The hydrolysis of acetic anhydride to free acetic acid appeared in the NMR spectrum through the methyl resonances at 2.0 and 1.8 ppm where the latter signal stems from acetic anhydride and the former from acetic acid. The following esterification reaction of free acetic acid with isoamyl alcohol to isoamyl acetate could be recognized by the very weak methyl resonance at 3.8 ppm. The water resonance, which was observed at 4.1 ppm at the beginning of the reaction, shifted in the course of the reaction and was found at 6.5 ppm by the end of the reaction. The signal migrated due to the change in the pH induced by the liberation of acetic acid. Due to the high selectivity of NMR spectroscopy, individual and wellseparated signals could be attributed. This is a prerequisite for the application of univariate data analysis. Hence, ct curves were readily obtained for the reaction, (see Section 4.3).
While high resolution is intrinsic to highfield NMR spectrometers, the much smaller benchtop spectrometer operating at 82 MHz proton Larmor frequency provided sufficient resolution to obtain relevant educt and product bands. The following paragraph illustrates how the information gained from NMR analysis can be quickly transferred to Raman and NIR analyses.
4.2. 2D Correlation Spectroscopy to Support Reaction Monitoring
To easily discover signals or bands suitable for reaction monitoring, heterocorrelation maps were computed from NMR, Raman, and NIR spectra, as shown in Figure 3. Since spectral differences along with the reaction progress were most easily recognized and explained in ^{1}H NMR spectra, the gained information was transferred to the Raman and NIR spectra. While in other cases Raman or NIR spectrum prove more readily interpretable, the inverse approach works evenly well.
(a)
(b)
Heterocorrelation spectra do not contain diagonal peaks [16, 48]. Thus, all correlations observed related signal changes at a given NMR chemical shift to intensity changes at wavenumbers in Raman spectra. It is therefore by the sign of a correlation peak that educt and product signals can be attributed. Using selected successfully attributed NMR signals, i.e., 2.0 ppm for acetic anhydride, 3.8 ppm for isoamyl acetate, and 1.8 ppm for acetic acid, correlations with the bands at 800, 1700 and 1740, and 3000 cm^{−1} were observed (cf. Figure 3). In the region between 1740 and 1790 cm^{−1}, asymmetric and symmetric stretching vibrations of the two C=O groups of carboxylic acid anhydrides were detected. The product’s methylene resonance at 3.8 ppm displayed only a very weak inphase or positive correlation to the stretching C=O vibration from acetic acid at 1700 cm^{−1}. The positive sign is indicative for a product signal. The analogue applies for the correlations at 2.0 ppm, 2980 cm^{−1} and 1.8 ppm, 2980 cm^{−1}: The negative (blue) signal indicates an eductproduct relationship, the positive (red) crosspeak a productproduct one. In the range of 810–910 cm^{−1}, weak CC stretching vibrations of acetates occurred. This signal was not directly apparent in the 1D Raman spectrum. Its counterpart however was visible in the NMR, which suggests that 2D correlation spectroscopy is also suitable for the identification of lowconcentration byproducts [49].
The enzymecatalysed hydrolysis of acetic anhydride could be followed through monitoring the C=O stretching vibration at 1740 cm^{−1} in Raman spectra and the methylene resonance in the range of 2.0 ppm in NMR spectra. The C=O stretching vibration displayed a negative correlation with the methylene resonance of acetic acid at 1.8 ppm, while the two educt and product signals run in opposite directions simultaneously. This is reflected by the positive correlations between 1740 cm^{−1} and 2.0 ppm, and between 1700 cm^{−1} and 1.8 ppm for the free acetic acid. Apart from a strong correlation in the range of the CH vibrations at 3000 cm^{−1}, a moderate correlation in the range of 800 cm^{−1} from the substances involved in the hydrolysis could be recognized. Firstly, a negative correlation of the acetic anhydride at 2.0 ppm resulted. Secondly, a positive correlation with the free acetic acid at 1.8 ppm was found [50].
In the ^{1}H NMRNIR 2D heterocorrelation (cf. Figure 3(b)), signals reflecting the hydrolysis were observed, but a detailed assignment could not be achieved. Yet, from the 2D NMRNIR correlation spectrum, the overtone and combination vibrations of CH and OH at 6250–7150 cm^{−1} and 4540–5000 cm^{−1} could be attributed by means of the respective positive and negative correlations to the methylene resonances at 2.0 and 1.8 ppm. It could also be recognized that a univariate evaluation would falsify the kinetic evaluation since observed bands represent superpositions of components [50].
In summary, the 2D heterocorrelation spectra helped to confirm the assignment of signals to an educt or a product by directly transferring the information from one spectroscopic technique to another. The 2D maps indicated the origin of vibrational bands from a single compound or as superposition. Heterocorrelation spectra were thus found to facilitate the identification and analysis of signals suitable for process monitoring. A selection of the 58 spectra recorded is shown in Figure 4 as stack plots for each spectroscopic method.
(a)
(b)
(c)
4.3. Quantitative Reaction Monitoring
Process analytics generally aims for the reception of concentrationtime (ct) data. Their extraction from spectra proceeded straightforward using the inline and online monitoring at the microreaction assembly. Through spectral analysis using heterocovariance correlation plots, the signals corresponding to the educts and the product were identified as described above. Based on distinct signals, univariate analysis leading to ct plots was applied. Nevertheless, series of spectra with the fully recorded spectral range were subjected to multivariate data analysis in the form of MCRALS.
Inspection of the NMR spectra (cf. Figure 4(a)) revealed that the water signal originating from the buffer system of the enzyme was well discernible. The spectra, in which the product signal was superimposed by the water shift, were discarded from kinetic evaluations. Since the areaunderthecurve of an NMR signal is proportional to the number of nuclear spins causing the signal, the actual relative concentration was determined by integrating the signals at 2.0, 1.8, and 3.8 ppm and by correcting them for the number of nuclei. To determine absolute concentrations, the internal standard method would need to be applied [51, 52]. Provided a sufficient linearity of the NMR amplifier, the absolute concentrations could be determined during the reaction since the initial concentrations of the reactants were known [51]. Due to nonlinearity of the amplifier, the relative composition of the mixture as obtained from direct integration of specific educt and product signals, hence univariate analysis, was plotted against reaction time (Figure 5).
(a)
(b)
While NMR spectra displayed a weak but interpretable signal of isoamyl acetate at 3.8 ppm, Raman and NIR spectra did not provide a corresponding band of sufficient intensity. Due to the choice of the process parameters, only a small amount of product was formed during the reaction period monitoring such that isoamyl ester product signals appeared weakly as seen by the NMR signal and its heterocorrelation at 3.8 ppm in the corresponding spectra. The first catalytic step, i.e., the hydrolysis of acetic anhydride to acetic acid, was monitored by the signals in the range of 2.0–1.8 ppm and the corresponding Raman band at 800 cm^{−1}. An acceleration of the reaction might be initiated by choosing a somewhat higher temperature. The esterification might occur eventually more prominent. The univariate evaluation of the signals from acetic acid and acetic anhydride yielded comparable reaction rate constants from Raman and NMR spectroscopic monitoring, showing the values of 2.1 and 2.2·10^{−3} min^{−1} for the transformation from acetic anhydride to acetic acid (cf. Table 1). Within Raman spectra, the bands at 1740 and 1700 cm^{−1} were used for educt and product monitoring, respectively. A trial to use univariate analysis on NIR data did not give satisfactory results. While envisaging automation for process control, MCRALS analysis is seemed to be an attractive method as identification of suitable bands would not be required. The MCRALS analyses of NMR and Raman spectra are shown in Figure 6. Before computing the ct curves, the firstorder reaction mechanism, A ⟶ B, an initial guess of the rate constants, and the educt starting concentration were defined as boundary conditions [53].

(a)
(b)
(c)
The algorithm transforms spectral data into ct data according to the predefined model, i.e., the firstorder reaction. The correspondence of the extracted data to the model is visualized by the ct curve. From Figure 6, it can be recognized that for Raman and NMR data, the agreement was very good. Only the NIR data did not prove of sufficient quality to yield good agreement to the model and thus minor scattering (cf. Figure 6(c)). The correspondence between ct data and the firstorder model was observed for all three spectroscopic techniques. Outliers, such as observed for the NIR analysis, may indicate disturbances, such as laboratory temperature change during sampling or measurement. Other reaction mechanisms were tested against the ct data. Neither the secondorder reaction A + B ⟶ C + D nor the followup reaction A ⟶ B ⟶ C displayed data fits with evenly good agreement as the firstorder model. However, the low yield of the isoamyl acetate and therefore the resulting low signaltonoise ratio influenced the performance of the algorithm so that no superior agreement with the models could be found.
Furthermore, the comparison of univariate to multivariate analyses showed a remarkable consistency with respect to the reaction order and the obtained rate constant. The results of the reaction monitoring of the enzymaticcatalysed hydrolysis reaction are summarized for univariate and multivariate analyses via MCRALS in Table 1.
All spectroscopic methods used revealed rate constants of a firstorder chemical kinetic model. Because of enzyme catalysis, the reaction is likely to follow a pseudofirstorder mechanism, which would not be discernible from a firstorder mechanism under the reaction conditions applied. Apparently, the catalysed hydrolysis proceeded faster or was preferred to the isoamyl ester formation. From univariate analysis of the NMR data (cf. Figure 5), a corresponding reaction constant of 1.3 10^{−3} min^{−1} was obtained. No multivariate description of this subsequent esterification could be achieved since isoamyl acetate was observed only in the NMR spectrum after about 2 hours of reaction time, provided careful visual inspection. It is assumed that the MCRALS algorithm was not able to identify isoamyl signals and—in case of Raman—bands that possessed a sufficient signaltonoise ratio. The signal intensity of the product proved too low for the determination of a reasonable ct profile through the MCRALS algorithm, whereas a trained analyst was able to extract consistent data. Thus, the univariate data analysis offered in this case the advantage to quantify extremely small signals, which also requires a greater effort for the evaluation.
The enzyme is also known to suffer from acid pH which was a consequence of the formation of acetic acid. The pH of the process solution dropped below the working range of the biocatalyst of pH 5–9. It was hence assumed that the CALB lost activity and ceased to produce isoamyl acetate. As a consequence, the subsequent reaction would not significantly affect the determined rate constants.
It should be noted that the intensity of vibrational bands depends on oscillator strength and polarizability, which is in general different for two different bands. Thus, the actual concentrations or concentration ratios of products and educts could not be derived from the Raman bands in contrast to the NMR resonances.
In summary, two approaches may be followed to obtain ct data, which proved fit for process analytical technology: univariate analysis and MCRALS as a multivariate approach. Covariance spectroscopy, here in the form of heterocorrelation spectroscopy, was suitable to help identify specific signals of educts and products. Fundamental process understanding could be achieved through describing the ct curves using kinetic models. Starting from the mathematical description of the data, it becomes possible to predict the course of the reaction and to quickly recognize deviations from the expected process. This knowledge may then be used to automate process control.
5. Conclusions
Analytical benchtop instruments, i.e., a compact ^{1}H NMR operating at 82 MHz, NIR and a Raman spectrometer, were successfully combined with a stainlesssteel microreaction assembly. The CALB catalysed hydrolysis, and subsequent esterification of acetic anhydride to isoamyl acetate served as model for a heterogeneous biocatalysis reaction. NMR spectra provided distinct signals of all educts and the product, and NIR and Raman bands could not easily be assigned on visual inspection of the spectra. To facilitate the assignments, 2D heterocorrelation maps were computed from series of 1D NMR and 1D vibrational spectroscopic data using covariance transformations.
Quantitative analysis of the spectroscopic data was achieved using univariate and multivariate methods. The ct data obtained from NMR, Raman, and NIR proved of sufficient quality to test a firstorder or pseudofirstorder kinetic model yielding consistent results. The MCRALS algorithm was applied successfully for a biocatalyzed reaction. Yet, an improved signaltonoise ratio, e.g., through increasing the number of acquisitions but thus prolongation of the measurement time, will enhance the recognition of lowabundant species. Furthermore, the MCRALS algorithm may be simultaneously applied to the data of the three analytical techniques. The knowledge on the process obtained by spectroscopic monitoring, heterocorrelationbased spectral interpretation, and univariate or multivariate kinetic analysis may allow future process control. For automation purposes, MCRALS was judged better suitable than univariate data treatment. Educt feed, reaction temperature, and flow rate should be optimized for best yields. These process parameters may then be adjusted automatically by detecting deviations from the corresponding reaction kinetics and counteracting through a feedback loop.
Data Availability
Data were recorded on NMR, Raman, and NIR spectrometers from Thermo Fisher Scientific. The raw data can be provided in the spectrometer file format. Data have been analyzed using MatLab. MatLab files can be obtained upon request from the corresponding author. Upon request, we can also include a link to download from the cloud Sciebo.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Acknowledgments
The authors are grateful to Thermo Fisher Scientific for the vibrational and nuclear magnetic resonance spectrometers within a SEED project and A. NickischHartfiel and P. Naderwitz for technical and biotechnological expertise. R.L is very grateful for a grant from the German Academic Scholarship Foundation. The authors thank the Niederrhein University of Applied Sciences (HSNR) for financial support.
References
 United States Food and Drug Administration, Guidance for Industry: PAT—A Framework for Innovative Pharmaceutical Development, Manufacturing, and Quality Assurance, U.S. Department of Health and Human Services, Washington, DC, USA, 2004.
 European Medicines Agency, Reflection Paper: Chemical, Pharmaceutical and Biological Information to be Included in Dossiers when Process Analytical Technology (PAT) is Employed, European Medicines Agency, Canary Wharf, London, UK, 2006.
 K. A. Bakeev, Process Analytical Technology, John Weiley & Sons Ltd, Chichester, England, 2nd edition, 2010.
 W. Chew and P. Sharratt, “Trends in process analytical technology,” Analytical Methods, vol. 2, no. 10, pp. 1412–1438, 2010. View at: Publisher Site  Google Scholar
 A. Chanda, A. M. Daly, D. A. Foley et al., “Industry perspectives on process analytical technology: tools and applications in API development,” Organic Process Research & Development, vol. 19, no. 1, pp. 63–83, 2014. View at: Publisher Site  Google Scholar
 R. J. A. Do Nascimento, G. R. De Macedo, E. S. Dos Santos, and J. A. De Oliveira, “Real time and in situ nearinfrared spectroscopy (NIRS) for quantitative monitoring of biomass, glucose, ethanol and glycerine concentrations in an alcoholic fermentation,” Brazilian Journal of Chemical Engineering, vol. 34, no. 2, pp. 459–468, 2017. View at: Publisher Site  Google Scholar
 R. Schalk, F. Braun, R. Frank et al., “Noncontact Raman spectroscopy for inline monitoring of glucose and ethanol during yeast fermentations,” Bioprocess and Biosystems Engineering, vol. 40, no. 10, pp. 1519–1527, 2017. View at: Publisher Site  Google Scholar
 K. Dégardin, A. Guillemain, and Y. Roggo, “Comprehensive study of a handheld Raman spectrometer for the analysis of counterfeit medicines,” Journal of Spectroscopy, vol. 2017, Article ID 3154035, 13 pages, 2017. View at: Publisher Site  Google Scholar
 B. Blümich and K. Singh, “Desktop NMR and its applications from materials science to organic chemistry,” Angewandte Chemie International Edition, vol. 57, no. 24, pp. 2–17, 2017. View at: Publisher Site  Google Scholar
 M. F. IsaacLam, “Determination of alcohol content in alcoholic beverages using 45 MHz benchtop NMR spectrometer,” International Journal of Spectroscopy, vol. 2016, Article ID 2526946, 8 pages, 2016. View at: Publisher Site  Google Scholar
 S. Kern, K. Meyer, S. Guhl et al., “Online lowfield NMR spectroscopy for process control of an industrial lithiation reaction—automated data analysis,” Analytical and Bioanalytical Chemistry, vol. 410, no. 14, pp. 3349–3360, 2018. View at: Publisher Site  Google Scholar
 F. Dalitz, M. Cudaj, M. Maiwald, and G. Guthausen, “Process and reaction monitoring by lowfield NMR spectroscopy,” Progress in Nuclear Magnetic Resonance Spectroscopy, vol. 60, pp. 52–70, 2012. View at: Publisher Site  Google Scholar
 N. Zientek, C. Laurain, K. Meyer, M. Kraume, G. Guthausen, and M. Maiwald, “Simultaneous 19F–1H medium resolution NMR spectroscopy for online reaction monitoring,” Journal of Magnetic Resonance, vol. 249, pp. 53–62, 2014. View at: Publisher Site  Google Scholar
 D. A. Snyder, R. Bruschweiler, G.E. Morris, and J.W. Emsley, “Multidimensional correlation spectroscopy by covariance NMR,” in Multidimensional NMR Methods for the Solution State, pp. 97–105, John Wiley & Sons Ltd, Chichester, UK, 2009. View at: Publisher Site  Google Scholar
 M. Jaeger and R. L. E. G. Aspers, “Covariance NMR and Small Molecule Applications,” Annual Reports on NMR Spectroscopy, vol. 83, 2014. View at: Google Scholar
 I. Noda, A. E. Dowrey, C. Marcott, G. M. Story, and Y. Ozaki, “Generalized twodimensional correlation spectroscopy,” Applied Spectroscopy, vol. 54, no. 7, pp. 236A–248A, 2000. View at: Publisher Site  Google Scholar
 I. Noda, “Generalized twodimensional correlation method applicable to infrared, Raman, and other types of spectroscopy,” Applied Spectroscopy, vol. 47, no. 9, pp. 1329–1336, 1993. View at: Publisher Site  Google Scholar
 I. Noda, “Twodimensional infrared spectroscopy,” Journal of the American Chemical Society, vol. 111, no. 21, pp. 8116–8118, 1989. View at: Publisher Site  Google Scholar
 N. Trbovic, S. Smirnov, F. Zhang, and R. Brüschweiler, “Covariance NMR spectroscopy by singular value decomposition,” Journal of Magnetic Resonance, vol. 171, no. 2, pp. 277–283, 2004. View at: Publisher Site  Google Scholar
 R. Brüschweiler, “Theory of covariance nuclear magnetic resonance spectroscopy,” Journal of Chemical Physics, vol. 121, no. 1, pp. 409–414, 2004. View at: Publisher Site  Google Scholar
 C. D. Eads and I. Noda, “Generalized correlation NMR spectroscopy,” Journal of the American Chemical Society, vol. 124, no. 6, pp. 1111–1118, 2002. View at: Publisher Site  Google Scholar
 L. J. Frasinski, K. Codling, and P. A. Hatherly, “Covariance mapping: a correlation method applied to multiphoton multiple ionization,” Science, vol. 246, no. 4933, pp. 1029–1031, 1989. View at: Publisher Site  Google Scholar
 I. Noda, “Modified twodimensional correlation spectra for streamlined determination of sequential order of intensity variations,” Journal of Molecular Structure, vol. 1124, pp. 197–206, 2016. View at: Publisher Site  Google Scholar
 L. R. Terra, M. N. Catrinck, and R. F. Teófilo, “MCRALS applied to the quantification of the 5hydroxymethylfurfural using UV spectra: study of catalytic process employing experimental design,” Chemometrics and Intelligent Laboratory Systems, vol. 167, pp. 132–138, 2017. View at: Publisher Site  Google Scholar
 W. Kessler and R. W. Kessler, “Multivariate curve resolution: a method of evaluating the kinetics of biotechnological reactions,” Analytical and Bioanalytical Chemistry, vol. 384, no. 5, pp. 1087–1095, 2006. View at: Publisher Site  Google Scholar
 Y. Le Dréau, N. Dupuy, J. Artaud, D. Ollivier, and J. Kister, “Infrared study of aging of edible oils by oxidative spectroscopic index and MCRALS chemometric method,” Talanta, vol. 77, no. 5, pp. 1748–1756, 2009. View at: Publisher Site  Google Scholar
 R. Tauler, “Multivariate curve resolution applied to second order data,” Chemometrics and Intelligent Laboratory Systems, vol. 30, no. 1, pp. 133–146, 1995. View at: Publisher Site  Google Scholar
 R. J. Pell and X. Chen, “Multivariate curve resolution for understanding complex reactions,” Journal of Chemometrics, vol. 28, no. 5, pp. 411–419, 2014. View at: Publisher Site  Google Scholar
 A. Güvenç, N. Kapucu, and Ü. Mehmetoǧlu, “The production of isoamyl acetate using immobilized lipases in a solventfree system,” Process Biochemistry, vol. 38, no. 3, pp. 379–386, 2002. View at: Publisher Site  Google Scholar
 F. W. Welsh, W. D. Murray, R. E. Williams, and I. Katz, “Microbiological and enzymatic production of flavor and fragrance chemicals,” Critical Reviews in Biotechnology, vol. 9, no. 2, pp. 105–169, 1989. View at: Publisher Site  Google Scholar
 C. H. Dickinson, J. Watson, and B. Wallace, “An impression method for examining epiphytic microorganisms and its application to phylloplane studies,” Transactions of the British Mycological Society, vol. 63, no. 3, pp. 616–619, 1974. View at: Publisher Site  Google Scholar
 Z. P. Chen, J. Morris, A. Borissova et al., “Online monitoring of batch cooling crystallization of organic compounds using ATRFTIR spectroscopy coupled with an advanced calibration method,” Chemometrics and Intelligent Laboratory Systems, vol. 96, no. 1, pp. 49–58, 2009. View at: Publisher Site  Google Scholar
 G. Mazarevica, J. Diewok, J. R. Baena, E. Rosenberg, and B. Lendl, “Online fermentation monitoring by midinfrared spectroscopy,” Applied Spectroscopy, vol. 58, no. 7, pp. 804–810, 2004. View at: Publisher Site  Google Scholar
 W. Kessler and R. W. Kessler, “Multivariate curve resolution—integration von Wissen in chemometrische modelle,” Chemie Ingenieur Technik, vol. 82, no. 4, pp. 441–451, 2010. View at: Publisher Site  Google Scholar
 J. Jaumot, N. Escaja, R. Gargallo, C. González, E. Pedroso, and R. Tauler, “Multivariate curve resolution: a powerful tool for the analysis of conformational transitions in nucleic acids,” Nucleic Acids Research, vol. 30, no. 17, p. 92e, 2002. View at: Publisher Site  Google Scholar
 R. Tauler, “Multivariate curve Resolution : theory and applications,” Trends in Analytical Chemistry, vol. 23, pp. 70–79, 2004. View at: Google Scholar
 W. Kessler, Multivariate Datenanalyse für die Pharma, Bio und Prozessanalytik, WileyVCH, Weinheim, Germany, 2006.
 A. L. Pomerantsev and O. Y. Rodionova, “Process analytical technology: a critical view of the chemometricians,” Journal of Chemometrics, vol. 26, no. 6, pp. 299–310, 2012. View at: Publisher Site  Google Scholar
 L. J. Frasinski, A. J. Giles, P. A. Hatherly, J. H. Posthumus, M. R. Thompson, and K. Codling, “Covariance mapping and triple coincidence techniques applied to multielectron dissociative ionization,” Journal of Electron Spectroscopy and Related Phenomena, vol. 79, pp. 367–371, 1996. View at: Publisher Site  Google Scholar
 R. L. E. G. Aspers, P. E. T. J. Geutjes, M. Honing, and M. Jaeger, “Using indirect covariance processing for structure elucidation of small molecules in cases of spectral crowding,” Magnetic Resonance in Chemistry, vol. 49, no. 7, pp. 425–436, 2011. View at: Publisher Site  Google Scholar
 M. Jaeger and R. Legner, “Homo and heterocovariance NMR spectroscopy and applications to process analytical technology,” in Spectroscopic Analyses—Developments and Applications, E. Sharmin and F. Zafar, Eds., pp. 55–80, 2017. View at: Google Scholar
 B. W. Hu, P. Zhou, I. Noda, and G. Z. Zhao, “An NMR approach applicable to biomolecular structure characterization,” Analytical Chemistry, vol. 77, no. 23, pp. 7534–7538, 2005. View at: Publisher Site  Google Scholar
 B. Hu, J. P. Amoureux, and J. Trebosc, “Indirect covariance NMR spectroscopy of throughbond homonuclear correlations for quadrupolar nuclei in solids under highresolution,” Solid State Nuclear Magnetic Resonance, vol. 31, no. 4, pp. 163–168, 2007. View at: Publisher Site  Google Scholar
 W. Ehrfeld, V. Hessel, and H. Loewe, “Microreactors—new technology for modern chemistry,” Organic Process Research & Development, vol. 5, no. 1, p. 89, 2001. View at: Publisher Site  Google Scholar
 R. Legner, S. Häfner, M. Voigt et al., “Micro process analytical technology—monitoring chemical processes in microreactors using miniature spectrometers,” GIT, vol. 34, pp. 38–40, 2016. View at: Google Scholar
 S. Morita, 2Dshige (c), KwanseiGakuin University, Nishinomiya, Japan, 20042005.
 J. Jaumot, A. de Juan, and R. Tauler, “MCRALS GUI 2.0: new features and applications,” Chemometrics and Intelligent Laboratory Systems, vol. 140, pp. 1–12, 2015. View at: Publisher Site  Google Scholar
 Y. Mee Jung and I. Noda, “New approaches to generalized two―dimensional correlation spectroscopy and its applications,” Applied Spectroscopy Reviews, vol. 41, no. 5, pp. 515–547, 2006. View at: Publisher Site  Google Scholar
 I. Noda, “Progress in twodimensional (2D) correlation spectroscopy,” Journal of Molecular Structure, vol. 799, no. 1–3, pp. 2–15, 2006. View at: Publisher Site  Google Scholar
 G. Socrates, Infrared and Raman Characteristic Group Frequencies, John Wiley & Sons Ltd., Chichester, UK, 3rd edition, 2004.
 G. Wider and L. Dreier, “Measuring protein concentrations by NMR spectroscopy,” Journal of the American Chemical Society, vol. 128, no. 8, pp. 2571–2576, 2006. View at: Publisher Site  Google Scholar
 F. Malz and H. Jancke, “Validation of quantitative NMR,” Journal of Pharmaceutical and Biomedical Analysis, vol. 38, no. 5, pp. 813–823, 2005. View at: Publisher Site  Google Scholar
 N. Viereck, L. Nørgaard, R. Bro, and S. B. Engelsen, “Chemometric analysis of NMR data,” Modern Magnetic Resonance, pp. 1833–1843, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Robin Legner et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.