Using Compact 1 H NMR , NIR , and Raman Spectroscopy Combined with Multivariate Data Analysis to Monitor a Biocatalyzed Reaction in a Microreaction System

Process analytical technology aims at process knowledge and process improvement, efficiency, and sustainability. A prerequisite is process monitoring. .e combination of microreaction systems and spectroscopy proved suitable due to dimension and compound reduction and real-time monitoring capabilities. Compact 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 real-time. To facilitate the identification of signals suitable for the extraction of concentrationtime (c-t) 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. .e data obtained were interpreted in terms of a firstorder kinetic model, and corresponding reaction rate constants were extracted. An alternative, elegant, and fit-for-automation approach for the kinetic analysis of the spectra was demonstrated in using multivariate curve resolution (MCR). .e results of the univariate and multivariate approaches were comparable with regard to reaction rates and concentrations. While the manual integration of the 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.


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 fast-growing research area [1][2][3].It has established itself as an important element of industrial production processes [4].us, process analytical instruments are located in or nearby large-scale processes and reactors.Process understanding results from their application, which in fact is process monitoring.e process knowledge obtained aims at ensuring a constant final product quality.It further enables process control and allows process optimization.e 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 near-infrared (NIR) play major roles [6][7][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.
erefore, 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 in-line, on-line, at-line or-less desirable-off-line.For in-line monitoring, a probe is inserted into the reaction vessel.When the measurement is conducted on-line, a bypass is used.For at-line analysis, a sample is taken from the reaction and treated before analysis.Hence, in-, on-, and at-line monitoring types usually provide results in real-time [1,2,7,11].
In order to compensate for the low resolution of, e.g., NIR or compact bench-top 1 H NMR spectrometers, suitable multivariate data analysis procedures have been devised and applied where univariate data processing did not prove sufficient [11][12][13].Another strategy for enhancing spectral resolution has been developed by using covariance data processing [14][15][16].
is approach has been derived from statistical mathematics.Hence, different spectra that possess a common so-called perturbation dimension are correlated, where the spectral intensity represents a function of two independent variables such as wavenumber or frequency.
is type of two-dimensional correlation spectroscopy has been introduced by Noda et al. for IR and also later for other vibrational and optical absorption methods [16][17][18].For NMR spectroscopy, Brueschweiler et al. and Eads and Noda have reported on the use of covariance transformations to correlate frequency domains [19][20][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.
Concentration-time (c-t) profiles from reaction monitoring can be obtained either by univariate methods or by multivariate methods such as the multivariate curve resolution alternating least squares (MCR-ALS) [24,25].Using this self-modelling curve resolution method, reaction kinetics can be determined directly from the recorded spectra series [26][27][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 real-time 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. e recorded spectra were transformed into 2Dheterocovariance correlation maps for the identification of specific resonances for univariate kinetic analysis.e results were compared to the kinetic parameters obtained from multivariate MCR-ALS analysis.

Multivariate Curve Resolution-Alternating Least Squares (MCR-ALS). Principal component nalysis (PCA)
and partial least square regression (PLS-R) 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. is is evenly true for the MCR-ALS algorithm.Physical-chemical 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 X 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, X represents spectral data.
e data matrix can be decomposed into the small matrices S and C which consist of pure spectra and the associated concentration profiles and the residual matrix E as follows [25,[34][35][36]: e matrix X 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.
e transposed matrix S T of dimensions a × s contains the spectra of the pure components as rows a.In matrix C of dimensions z × a, the contributions of the pure spectra in columns a correspond to the concentrations in the corresponding spectrum.e residual matrix E 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 C and S can be computed by MCR-ALS.ey are iteratively calculated from a first estimate for C or S in an alternating least-squares method as shown below [24,35,37]: where ∧ denotes the estimated value and + a pseudoinverse.e decomposition of X is achieved by iterative least-squares minimization of ‖E‖ 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 physical-chemical knowledge of the system as possible into the constraints [34,38].

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 one-dimensional spectra as responses to a quantitative, measurable perturbation applied to the sample in a systematic way [16,18].is perturbation may be temperature, pressure, concentration, stress, electrical field, and most importantly a chemical reaction [17,41].A correlation spectrum or map, X(] 1 , ] 2 ), is defined according to the following equation: where  y(] 1 , t) is the spectrum affected by the perturbation, called the dynamic spectrum, and ] 1 and ] 2 are two different spectral variables at a given time or a fixed interval of the external variable t. e symbol 〈 〉 denotes a suitable correlation function.e dynamic spectrum is computed by subtracting the average spectrum from each spectrum of a series of spectra.e 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 X(] 1 , ] 2 ) can be considered as the covariance matrix, when the series of dynamic spectra  y(] 1 , t) is represented by a matrix.Assuming that X(] 1 , ] 2 ) 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].e proof is based essentially on Parseval's theorem.
e 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.is is referred to as heterocovariance spectroscopy [21,41].A heterocorrelation map will result according to the following equation: where μ 1 and ] 2 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 F is the 2D spectrum, C is the covariance matrix, and S and G are the two data matrices, such as a series of spectra with a common perturbation dimension.Cross-peaks 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.

Materials and Methods
3.1.Spectral Recording and Processing.NMR spectra were recorded on-line using a compact NMR picoSpin80 spectrometer with a proton Larmor frequency of 82 MHz ( ermo Fisher Scientific GmbH, Dreieich, Germany).It had a built-in 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.e sample cell possessed a total volume of 40 µL and an active volume of 40 nL.e operation temperature was 36 °C.e spectrometer was controlled by ermo Fisher PicoSpin software 0.9.3 running on the spectrometer control board and accessed via a web interface from a laptop computer.e free induction decay (FID) was recorded as 4092 data points that were zero filled to 8 k points prior to the Fourier transformation.e 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.e 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. e spectrum size after zero-filling and Fourier transformation (FT) resulted in 64 k.Pretreatment methods were used for spectra optimization.e 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 ™ Smar- tRaman spectrometer ( ermo Fisher Scientific GmbH, Dreieich, Germany) and processed using the Omnic software, version 9.2 ( ermo 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. e probe was positioned at its focal point of 5 mm above the solution in the reaction vessel.Sixty-four spectra comprising a spectral range of 250-3500 cm −1 with a resolution of 5 cm −1 were acquired at each reaction time point.e laser was operated at 785 nm and 150 mW.For processing, the Omnic software Scheme 1: Reaction scheme for the biocatalyzed hydrolysis of acetic anhydride followed by esterification with isoamyl alcohol to isoamyl acetate.Immobilized CALB catalyses both reaction steps as one-pot reaction in the microreaction system [29].
9.2 ( ermo 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.e 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 FT-NIR Analyzer using the software Omnic, version 8.0 ( ermo Fisher Scientific GmbH, Dreieich, Germany).e 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. e 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.

Microreaction Assembly.
Stainless-steel 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.e 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.e immobilized enzyme was then transferred into the cartridge of the Ehrfeld reactor F 200.
e filled volume amounted to 1.8 of the 2 mL cartridge volume.e reaction was initiated through heating the reactor to 35 °C.e residence time of a solution fraction, approximately 2 mL, in the reactor amounted to 2 min.e solution, then containing product, was recirculated into the storage vessel.e microreaction assembly was controlled through the software LabVision 2.10 (HiTec Zang GmbH, Herzogenrath, Germany).
e 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.e 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).e reaction was conducted over a period of 1380 min.

Covariance and Multivariate Data Analysis.
e 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. e data files were transformed into homo-or heterocorrelation maps and displayed as contour or color plots using appropriate threshold levels.
For MCR-ALS analysis, NMR, Raman, and NIR spectra were automatically imported into MatLab, version 2017b, wherein the MCR-ALS toolbox 2 was applied [35,47].ree latent variables or components were selected on the basis of visual inspection of the computational results.
e 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 .e agreement between calculated and experimental concentration-time profile was taken as goodness of the obtained kinetic parameters.

Qualitative Reaction Monitoring.
e enzyme catalysed acetic anhydride hydrolysis, and subsequent esterification with isoamyl alcohol was conducted in a microreaction platform.
e enzyme was immobilized and located in a cartridge reactor as a heterogeneous catalyst.e process was monitored using on-line 1D 1 H NMR and in-line 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 c-t 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).e 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.e 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.e 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.e 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 well-separated signals could be attributed.is is a prerequisite for the application of univariate data analysis.Hence, c-t curves were readily obtained for the reaction, (see Section 4.3).
While high resolution is intrinsic to high-field NMR spectrometers, the much smaller bench-top spectrometer operating at 82 MHz proton Larmor frequency provided sufficient resolution to obtain relevant educt and product bands.
e following paragraph illustrates how the information gained from NMR analysis can be quickly transferred to Raman and NIR analyses.

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.Heterocorrelation spectra do not contain diagonal peaks [16,48].us, 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.e product's methylene resonance at 3.8 ppm displayed only a very weak in-phase or positive correlation to the stretching C�O vibration from acetic acid at 1700 cm −1 .e positive sign is indicative for a product signal.e analogue applies for the correlations at 2.0 ppm, 2980 cm −1 and 1.8 ppm, 2980 cm −1 : e negative (blue) signal indicates an educt-product relationship, the positive (red) cross-peak a product-product one.In the range of 810-910 cm −1 , weak C-C stretching vibrations of acetates occurred.is 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 low-concentration byproducts [49].
e enzyme-catalysed 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.e 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.is 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 C-H 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 NMR-NIR 2D heterocorrelation (cf. Figure 3(b)), signals reflecting the hydrolysis were observed, but a detailed assignment could not be achieved.Yet, from the 2D NMR-NIR correlation spectrum, the overtone and combination vibrations of C-H and O-H 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.e 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.

Quantitative Reaction Monitoring.
Process analytics generally aims for the reception of concentration-time (c-t) data.eir extraction from spectra proceeded straightforward using the in-line and on-line monitoring at the microreaction assembly.rough 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 c-t plots was applied.Nevertheless, series of spectra with the fully recorded spectral range were subjected to multivariate data analysis in the form of MCR-ALS.
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.e spectra, in which the product signal was superimposed by the water shift, were discarded from kinetic evaluations.Since the area-underthe-curve 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).
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.e 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.8ppm and the corresponding Raman band at 800 cm −1 .An acceleration of the reaction might be initiated by choosing a somewhat higher temperature.e esterification might occur eventually more prominent.
e univariate evaluation of the signals from   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, MCR-ALS analysis is seemed to be an attractive method as identification of suitable bands would not  be required.e MCR-ALS analyses of NMR and Raman spectra are shown in Figure 6.Before computing c-t curves, the first-order reaction mechanism, A ⟶ B, an initial guess of the rate constants, and the educt starting concentration were defined as boundary conditions [53].e algorithm transforms spectral data into c-t data according to the predefined model, i.e., the first-order reaction.
e correspondence of the extracted data to the model is visualized by the c-t 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)).e correspondence between c-t data and the first-order 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 c-t data.Neither the second-order reaction A + B ⟶ C + D nor the follow-up reaction A ⟶ B ⟶ C displayed data fits with evenly good agreement as the first-order model.However, the low yield of the isoamyl acetate and therefore the resulting low signal-to-noise 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.
e results of the reaction monitoring of the enzymaticcatalysed hydrolysis reaction are summarized for univariate and multivariate analyses via MCR-ALS in Table 1.
All spectroscopic methods used revealed rate constants of a first-order chemical kinetic model.Because of enzyme catalysis, the reaction is likely to follow a pseudo-first-order mechanism, which would not be discernible from a first-order 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 MCR-ALS algorithm was not able to identify isoamyl signals and-in case of Raman-bands that possessed a sufficient signal-tonoise ratio.e signal intensity of the product proved too low for the determination of a reasonable c-t profile through the MCR-ALS algorithm, whereas a trained analyst was able to extract consistent data.us, the univariate data analysis offered in this case the advantage to quantify extremely small signals, which also requires a greater effort for the evaluation.
e enzyme is also known to suffer from acid pH which was a consequence of the formation of acetic acid.e 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   In summary, two approaches may be followed to obtain c-t data, which proved fit for process analytical technology: univariate analysis and MCR-ALS 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 c-t 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.is knowledge may then be used to automate process control.

Conclusions
Analytical bench-top instruments, i.e., a compact 1 H NMR operating at 82 MHz, NIR and a Raman spectrometer, were successfully combined with a stainless-steel microreaction assembly.e 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.e c-t data obtained from NMR, Raman, and NIR proved of sufficient quality to test a first-order or pseudofirst-order kinetic model yielding consistent results.
e MCR-ALS algorithm was applied successfully for a biocatalyzed reaction.Yet, an improved signal-to-noise ratio, e.g., through increasing the number of acquisitions but thus prolongation of the measurement time, will enhance the recognition of low-abundant species.Furthermore, the MCR-ALS algorithm may be simultaneously applied to the data of the three analytical techniques.
e 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, MCR-ALS was judged better suitable than univariate data treatment.Educt feed, reaction temperature, and flow rate should be optimized for best yields.
ese process parameters may then be adjusted automatically by detecting deviations from the corresponding reaction kinetics and counteracting through a feedback loop.

Figure 1 :Figure 2 :
Figure 1: (a) Process flow-chart of the microreaction assembly: 1, Raman and NIR probe; 2, storage vessel containing the majority of the reaction mixture; 3, piston pump; 4, heat exchanger; 5, cartridge reactor; 6, temperature sensor; and then a compact 1 H NMR spectrometer is located in a bypass.(b) Picture of the microreaction platform (from left to right): inlet module, heat exchanger, cartridge reactor, T-sensor, optical flow cell, and outlet module.Raman probe and transflection NIR immersion probe (front).

Figure 4 :
Figure 4: Raman (a), NIR (b), and1 H NMR (c) spectra from on-line reaction monitoring of the enzyme catalysed reaction in the microreaction assembly.Spectra were recorded from the neat reaction mixture during 23 hours.NMR spectra, 16 accumulations, were acquired at 82 MHz and T � 36 °C.NIR spectra were recorded using an immersion probe.Raman spectra, 64 accumulations, were acquired with laser excitation at 785 nm.Color coding refers to corresponding reaction times among the spectroscopic methods.e reaction was followed over a period of 1380 min.

Figure 5 :
Figure5: c-t plots as area-under-the-curve vs. time of the enzyme catalysed esterification from univariate data analysis and corresponding fits computed according to first-order reaction based on NMR (a) and Raman (b) spectra; acetic acid (blue), acetic anhydride (orange), and isoamyl acetate (yellow).

Table 1 :
Reaction order and rate constants according to 1 H NMR, Raman, and NIR spectroscopic methods for the enzyme-catalysed hydrolysis of acetic anhydride using univariate and multivariate spectral analyses.two different bands.us, the actual concentrations or concentration ratios of products and educts could not be derived from the Raman bands contrast to the NMR resonances.