Next Article in Journal
Novel Intensified Alternatives for Purification of Levulinic Acid Recovered from Lignocellulosic Biomass
Next Article in Special Issue
Digital Twin of mRNA-Based SARS-COVID-19 Vaccine Manufacturing towards Autonomous Operation for Improvements in Speed, Scale, Robustness, Flexibility and Real-Time Release Testing
Previous Article in Journal
Mechanical Response and Deformation Mechanisms of TB17 Titanium Alloy at High Strain Rates
Previous Article in Special Issue
PAT for Continuous Chromatography Integrated into Continuous Manufacturing of Biologics towards Autonomous Operation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Process Analytical Technology for Precipitation Process Integration into Biologics Manufacturing towards Autonomous Operation—mAb Case Study

Institute for Separation and Process Technology, Clausthal University of Technology, Leibnizstr. 15, 38678 Clausthal-Zellerfeld, Germany
*
Author to whom correspondence should be addressed.
Processes 2021, 9(3), 488; https://doi.org/10.3390/pr9030488
Submission received: 27 January 2021 / Revised: 15 February 2021 / Accepted: 25 February 2021 / Published: 9 March 2021
(This article belongs to the Special Issue Towards Autonomous Operation of Biologics and Botanicals)

Abstract

:
The integration of real time release testing into an advanced process control (APC) concept in combination with digital twins accelerates the process towards autonomous operation. In order to implement this, on the one hand, measurement technology is required that is capable of measuring relevant process data online, and on the other hand, a suitable model must be available to calculate new process parameters from this data, which are then used for process control. Therefore, the feasibility of online measurement techniques including Raman-spectroscopy, attenuated total reflection Fourier transformed infrared spectroscopy (ATR-FTIR), diode array detector (DAD) and fluorescence is demonstrated within the framework of the process analytical technology (PAT) initiative. The best result is achieved by Raman, which reliably detected mAb concentration (R2 of 0.93) and purity (R2 of 0.85) in real time, followed by DAD. Furthermore, the combination of DAD and Raman has been investigated, which provides a promising extension due to the orthogonal measurement methods and higher process robustness. The combination led to a prediction for concentration with a R2 of 0.90 ± 3.9% and for purity of 0.72 ± 4.9%. These data are used to run simulation studies to show the feasibility of process control with a suitable digital twin within the APC concept.

1. Introduction

Innovations in biologics manufacturing are accelerating and driving towards continuous operation to cope with increasing titers and the multitude of new entities such as antibody fragments, virus-like particles (VLPs), exosomes and messenger ribonucleic acid (mRNA) [1,2]. Due to the rising demand, old platform processes must be optimized, redesigned, or converted to novel products. Process control offers opportunities for optimization but requires a deep process comprehension. Achieving this, the way to process modeling is paved. The underlying process for this study is a monoclonal antibody manufacturing (mAb) in CHO cells (Chinese hamster ovary) [3,4,5,6], which is still the typical workhorse to present innovations.
Recent works have shown that selective precipitation [4,7,8,9,10,11,12,13,14,15] is a suitable alternative to protein A chromatography, which was identified to be the bottleneck in mAbs manufacturing. Precipitation can be achieved differently, e.g., with cold ethanol [14,16], polyethylene glycols (PEG) of different sizes (3350–8000 g/mol) [4,7,8], with a combination of PEG and salts [11,17] or even with salts alone [10].
Precipitation has been known for a long time and was formerly used for blood plasma fractionation, which has not changed to the current day [18,19]. The precipitation process is fast, does not require expensive equipment and is a concentration-independent process. Which means that it can cope with increasing or variating titers. Nonetheless, precipitation is not widespread in biologics manufacturing, which has begun to change in the past decade. Moreover, several working groups have shown that precipitation can be operated continuously by suitable operating modes such as a feed-and-bleed configuration [9,10,12].
Due to the rising demand of new entities, the FDA has established a general guideline (Quality by Design, QbD) to standardize process development and ensure uniform product quality. These guidelines include early-stage decisions on process layout and design such as the selection of stainless-steel equipment or disposables, but also the definition of critical quality attributes (CQAs). In the case of precipitation antibody concentration, yield and purity are defined as CQAs. Furthermore, the guideline includes a suitable control and regulation concept to ensure the compliance throughout the process. Occurring deviations of the CQAs must be monitored and compensated for a uniform product quality. To monitor this, mostly at-line or offline analytics are used that are time intensive and which only approximate the actual value in the process due to the delay between sampling and analysis. Moreover, it can only be partially automated. Process analytical technology (PAT), on the other hand, offers the possibility of recording real-time data during the process without taking samples. A huge advantage is also that data can be analyzed automatically [20,21,22,23,24,25].
Since precipitation is nevertheless a rarely used purification technique, not many works on online measurement technology in connection with precipitation have been published. PAT-related works have been accomplished for alcohol precipitation [26,27,28,29,30,31,32,33]; acid precipitation [34] and PEGylation [35]. Most of these works used NIR as spectroscopic method [26,27,28,29,30,31,32,33,34,35] and only a few have dealt with FTIR and Raman [36,37]. Furthermore, Focus beam reflectance measurement (FBRM) and particle vision measurement (PVM) have been used to study particles during precipitation regarding their number distribution and shape [38]. Furthermore, detectors such as Raman [39,40,41,42], FTIR (Fourier transform infrared spectroscopy) [43,44,45,46], DAD (diode array detector) [47] and fluorescence [48] have already been successfully used in other unit operations to monitor CQAs [49,50]. Which is why they are selected as promising candidates for precipitation and are tested in the present work.
To apply the online measured data to the process, a process model (digital twin) must be developed first to pave the way for process control. With the aid of the model and the online monitoring of CQAs, the process optimization can be conducted as well as the monitoring of aging processes. This is important especially for membranes.
Figure 1 shows a simplified flowsheet of one precipitation subunit, which consist of a stirred tank reactor where precipitation proceeds and a hollow fiber module where filtration and dissolution take place. After precipitation in the tank, precipitates are transferred to the filter module and are deposited onto the membrane through dead-end filtration. Hereby, valves are such positioned, that the permeate is discharged into the waste. Afterwards valves are closed to enable a closed loop for dissolution. Once the online measured redissolution of mAb is sufficient, valves are opened, and the intermediate product is filtered through the membrane and transferred to chromatography. The interconnection of several subunits such as the one shown in Figure 1 enable continuous processing of the incoming feed (light phase, LP) from aqueous-two phase extraction (ATPE).
The aim of this work is, firstly, a feasibility study on spectroscopic methods to monitor mAb concentration, yield, and purity throughout the precipitation process. For this purpose, the compatibility of the individual detectors is tested. Secondly, the combination of several detectors is tested, and a recommendation is given which detector is best suited for the precipitation unit. Thirdly, an advanced process control (APC) concept is proposed, which uses the online measured data for process optimization and monitoring the actual state of the process.

2. Materials and Methods

During this experimental series, three process cycles, starting with cultivation followed by extraction, precipitation, chromatography, and the UF/DF were performed at the Institute of Thermal Process Engineering. This paper only deals with the precipitation part. The overview of the whole study has been previously published [6]. Cell culture broth is obtained through cultivation of Chinese hamster ovary cells (CHO DG44) producing an industrially interesting monoclonal antibody (mAb) from the subclass IgG1. Detailed information about culture broth [1,3,4,5,41] and ATPE [4,51] have also been previously published.
The precipitation unit in this work consist of two parts: the precipitation step, where mAb precipitates from the feed solution and the dissolution, where resolution of these precipitates takes place. Operation procedure is conducted similar to a previous work described by Lohmann et al. [4]. In this work batch units are connected to a continuous process. Precipitation is accomplished after ATPE from the polymer rich phase (light phase, LP) using a Polyethylene glycol (PEG 4000) stock solution (40 wt %, pellets, from Sigma Aldrich, dissolved in 60 wt % HQ H2O). For complete precipitation, streams are adjusted to result in a PEG content of 12 wt % respective to light phase volume. Precipitation sampling required a workaround because precipitation is such a spontaneous process. In order to train the partial-least-square model (PLS) for spectral analysis several equilibrium stages are used to simulate the precipitation and display the changes in the spectra. Samples from 1 wt % to 17 wt % PEG (in 2 wt % PEG steps) are prepared.
Solid–liquid separation is performed using a hollow fiber module with a pore size of 0.2 µm (88 cm2, Midikros, Repligen Corporation, Rancho Dominguez, CA, USA). To minimize product loss, the module is operated in dead end mode to completely deposit the precipitates onto the membrane surface. Precipitates are filtered in dead-end mode at constant flux until the critical transmembrane pressure of 1800 mbar is reached. Thereafter, the feed stream is switched to the next filter module. Afterwards, precipitates are washed with PEG 4000 (12 wt %) solution to maintain precipitation conditions and prevent resolution of mAb. Subsequently, dissolution is performed by recycling a 50 mM sodium dihydrogen phosphate/disodium phosphate buffer at pH = 5.5 along the fibers for 14 h. Sampling rate in dissolution was every 10 mins during first two hours and then one sample every 60 mins. After dissolution, the pH is adjusted to 3.7 and held for 60 min for virus inactivation [52]. Investigation of suitable PAT detectors and sensors (pH, Turbidity, conductivity) is performed separately for precipitation and dissolution. The following PAT detectors have been used throughout the study: Raman (785 nm laser, 1.5 mW, Ocean Insight, Ostfildern, Germany, with an InPhotonics Raman probe), FTIR (Alpha II, Bruker, Billerica, MA, USA), DAD (Smartline DAD 2600, Knauer Wissenschaftliche Geräte GmbH, Berlin, Germany), and fluorescence (Jasco FP-2020 Fluorescence detector). With the exception that the fluorescence detector is not used in precipitation due to precipitates which are harmful for the flow cells of the fluorescence detector and they bias the measurement results. PLS models are built for three analytes: target component (TC), concentration of high molecular weight impurities (HMW) and light weight molecular impurities (LMW), investigating the change of concentration and purity over process time.

2.1. Analytics

2.1.1. Online

Data acquisition in Raman spectroscopy is set to automatically average a total of three spectra, with an integration time of 1s each. Data acquisition in FTIR spectroscopy is set to automatically average a total of 60 spectra, with an integration time of 1s each. Data acquisition in UV-Vis spectroscopy is set to continuously record at a sampling rate of 0.2 s, each of which have 32 msec of integration time. Data acquisition in fluorescence spectroscopy (Jasco FP-2020 Fluorescence detector) is conducted with a gain factor of 1. Spectra are measured with an excitation wavelength of 280 nm and the emission spectra is detected between 280 and 900 nm. For each sample at least three spectra are measured and averaged prior to PLS-regression.

2.1.2. Offline

The monoclonal antibody is quantified by Protein A chromatography (PA ID Sensor Cartridge, Applied Biosystems, Bedford, MA, USA). Dulbecco’s PBS buffer (Sigma-Aldrich, St. Louis, MO, USA) is used as loading buffer at pH of 7.4 and as elution buffer at pH of 2.6. The size exclusion chromatography (SEC) is performed with a YarraTM 3µm SEC 3000 column (Phenomenex Ltd., Aschaffenburg, Germany) utilizing 0.1 M Na2SO4, 0.1 M Na2HPO4, and 0.1 M NaH2PO4 (Merck KGaA, Darmstadt, Germany) as buffering system. The absorbance for both methods is monitored at 280 nm.

2.2. Preprocessing of Spectral Data

Spectral data are processed and analyzed using Unscrambler® X (Camo Analytics AS, Oslo, Norway). First, raw spectral data are trimmed to the fingerprint region for proteins which is for Raman 1800–400 cm−1 and for FTIR 1800–900 cm−1 respectively. DAD and Fluorescence spectra are not trimmed and used completely for analysis. The preprocessing strategy can vary for the obtained data and must be adjusted. To perform an adapted preprocessing, descriptive statistic is analyzed. Hereby, the scatter effects plot is used, which depicts the sample spectrum against the mean spectrum. The following scatter effects shown in Figure 2 can be present: (a) additive effects, (b) scatter effects, (c) additive and scatter effects and (d) complex effects.
Additive effects are eliminated by conducting a derivation or applying a baseline correction. Multiplicative scatter effects can be removed using the method of standard normal variate (SNV) or a multiplicative scatter correction (MSC). Combinations of additive and multiplicative effects are removed by a combination of baseline correction, and scatter correction and complex effects can be eliminated by using an extended multiplicative scatter correction (EMSC). Depending on the scatter effects, the pre-processing is adapted to extract a maximal information content. Pre-processed spectra are then used to correlate the changes in component concentration to the changes in spectral intensity at specific spectral regions by partial least squares regression. A maximum of four principal components was allowed in the regression to not overfit the data. The quality of the model was evaluated by means of spectral line loadings plot, score plot, and the plot of explained variance against number of principal components. If Pat inline would be feasible, another tool for APC would be a suitable digital twin based on process modeling as described in the following. Based on the precipitation unit proposed in a previous work [4], the process model for the APC concept is presented. The model is divided in to three sub models: Precipitation of mAb, filtration of precipitates and resolution of precipitated mAb which is schematically shown in Figure 3. Division into sub models is necessary for separation of effects within precipitation, filtration, and dissolution.

2.3. Precipitation Model

Precipitation takes place batchwise and is modeled as a continuous stirred tank reactor (CSTR). Precipitation with PEG is a spontaneous process. In literature fractal analysis is used to infer the rate limiting step in precipitation. Hereby, the fractal dimension is used to categorize in reaction limited and diffusion limited precipitation. Accordingly, a reaction-limited precipitation is present at low fractal dimensions of 1.7 and a diffusion-limited one at a fractal dimension of 2.4 [53,54]. Satzer et al. have demonstrated a methodology to 3D reconstruct precipitates using light microscopy in combination with advanced image evaluation [55]. The 3D reconstructs have been used to perform a fractal analysis to determine fractal dimension and fractality on micro- and nanoscale. Satzer et al. showed that PEG precipitates have an average fractal dimension of 2.4 [55].
To evaluate this, precipitates have been analyzed through light microscopy, which has shown that the particles consist of sphere-like structures that agglomerate and tend to monofractility (see Figure 4). Fractality is determined through the method published by Satzer et al. [55] and lead to comparable results of a fractality of 2.4. Therefore, it is expected to be a diffusion-limited precipitation, which is in good agreement with the experimental results described in Lohmann et al. [4] and the theory that IgG diffuses very slowly due to its size.
For this reason, concentration decrease during precipitation is described as a reactive crystallization due to the spontaneous and diffusion-limited process. The concentration change in the supernatant is calculated by applying a 1st order reaction kinetics [56] as can be seen in Equation (1):
r p r e c = k p r e c * ( c m A b c e q )
where r p r e c represents the reaction rate which is associated with the concentration change of the target component over process time, k p r e c the rate constant of precipitation, c m A b the instantaneous mAb concentration and c e q p r e c its equilibrium concentration. The precipitation constant was determined via turbidity in a previous work [4]. In literature precipitation is described as the change of chemical potential which leads to formation of dense precipitates as their status is more stable under new surrounding conditions [57]. This results in a shift of protein solubility (S) which is described by Equation (2) [58]. In the proposed model the maximal protein solubility equals the equilibrium concentration of mAb ( c e q ).
log ( S ) = β * ω + κ
Whereas, β stands for precipitation efficiency, ω for the weight percentage of PEG added to light phase and κ the intrinsic protein solubility [58]. Sim et al. derived a model to describe the precipitation efficiency in presence of PEG molecules based on attractive depletion and excluded volume [57,59]. β is determined as described in Sim et al. [57] Values for γ =   0.076 and δ = 0.045 are also adopted from [57]. The needed weight percentage of PEG for complete precipitation is determined through equilibrium experiments in a previous work [4].
β = γ * r h , P E G 0.211 * r h , P r o t e i n + δ * r h , P r o t e i n
Hydrodynamic radius for the target component is calculated with Equation (4). The hydrodynamic radius for PEG is obtained from literature [59].
r h , P r o t e i n   = ( 3 η M r , P r o t e i n 4 π   2.5   N A   ) 1 3

2.4. Filtration Model

The filtration of precipitates is performed with a hollow fiber module through dead-end filtration at constant flux. Hollow fiber filters are not designed for dead end operation, but their design has an advantage for subsequent redissolution. Modeling of dead-end filtration is commonly described as the occurring blocking mechanism, which can be calculated as a pressure change over time [60,61,62]. The commonly used blocking mechanisms are depicted in Figure 5 [62]. Complete blocking occurs when the particles to be filtered have approximately the same size as the pores in the membrane (Figure 5a). Standard blocking appears while particles are smaller than the pores and deposit inside the pores (Figure 5b). Cake filtration occurs when the particles form a dense layer on the membrane (Figure 5d) and intermediate blocking is a composed blocking mechanism of cake filtration and complete blocking (Figure 5c). A widely maintained opinion in literature is that description by only one blocking mechanism is not sufficient due to several superimposed effects [62,63,64,65].
In this case the filtration could be described best by a composite model of cake filtration (see Equation (5)) and standard blocking (see Equation (6)). First, cake filtration occurs which changes to standard blocking. p 0 , c a k e is the initial pressure at the beginning of the experiment and p 0 , s t a n d a r d means the current pressure at the switching point to standard blocking. J 0 is the transmembrane flux, which is kept constant throughout filtration and v stands for the membrane specific volume. Blocking constants K s and K c are determined for the specific filter with a calibration experiment in advance.
p = ( 1 + K c · J 0 · v ) · p 0 , c a k e
p = p 0 , s t a n d a r d ( 1 ( K s · v 2 ) ) 2

2.5. Dissolution Model

Model equations for dissolution are based on the same mechanism as precipitation, however, with a different mass transfer coefficient for dissolution. Product losses in dissolution may occur due to two cases: limitation through the feed concentration or the maximal protein solubility. The limitation by the feed concentration ( c m a x ) appears when the equilibrium concentration in the dissolution buffer is significantly higher than the mAb concentration in the feed. In this case the complete initially precipitated mass of mAb can be redissolved (see Equation (7)). If the depletion of the PEG or the volume ratio (DR) between feed and buffer is too low, limitation by equilibrium concentration (see Equation (8)) is the result due to a residual PEG content which should not exceed 3 wt % [13]. r d i s s describes the reaction kinetic dissolution of precipitates dependent on the circulation flowrate through the fibers. c m A b it the current mAb concentration dissolved in the dissolution buffer and c e q d i s s the equilibrium concentration of mAb in dissolution. As described before, equilibrium parameters do not change in dissolution, but the equilibrium concentration does due to reduction of the PEG content in solution. DR describes the ratio between feed and dissolution buffer volume.
r d i s s =   k d i s s · ( c e q d i s s c m A b )
The mass transfer from the membrane surface back into dissolution is dependent on the recirculation flow rate trough the fibers and is calculated as reversed concentration polarization. Material transfer in concentration polarization is described through the coefficient k e f f which is dependent on the flow regime. First, Reynolds number (Re) must be calculated for identification of the current flow regime. This relates inertial forces and viscous forces in the fluid [66].
R e = u e f f · ρ η
The effective velocity is referring to the empty tube velocity. V ˙ is the flow rate through one fiber during dissolution, A the cross-sectional area of all fibers and ϵ the porosity.
u e f f = V ˙ A · ϵ
Second, Schmid number Sc is determined for relation of diffusive impulse transport to diffusive mass transport. Whereas, D stands for diffusion coefficient of mAb, η the viscosity of the fluid and ρ for the density [67].
S c = η ρ · D
From Reynolds and Schmidt follows the Sherwood number which describes the ratio of the effective amount of substance transferred to the amount of substance transported by diffusion. The correlation by Leveque applies only for Reynolds number below 1800 [68]:
S h = 1.62 · R e 0.33 · S c 0.33 · ( d H L c h ) 0.33   Re < 1800
With dH as hydraulic diameter and Lch as characteristic length of the module. Using the Sherwood number, the effective mass transfer coefficient related to the surface of the membrane follows as:
k e f f = S h · D d H
The time dependent parameter k d i s s of the process model is obtained by multiplication of k e f f with the volume specific surface of the fibers.
k d i s s = a · k e f f

3. PAT Feasibility Study Results and Discussion

Due to the different experimental setup, the results for precipitation and dissolution are shown separately. Investigated components are high molecular weight impurities (HMWs), Target component (TC) and low molecular weight impurities (LMWs). Categorization of very good (1 < R2 < 0.9), moderate (0.9 < R2 < 0.8) and poor (R2 < 0.8) is based on the R2 of generated PLS models. These ranges are narrow due to the goal to replace standard offline analytics (HCPLs).

3.1. Raman Spectroscopy

In this section Raman spectral data are discussed, first for precipitation and following for dissolution. The results (mAb titer and purity) obtained in dissolution are the results of interest for the following unit operation. First, raw data are reduced to 1800–400 cm−1 Raman shift equally for precipitation and dissolution results. It becomes obvious, that spectral raw data for both steps are very distinct (Figure 6a,b). This can be explained by the different matrix of the samples. In precipitation the continuous phase is still the light phase obtained from ATPE, including the phase forming components phosphate and PEG 400 as well as the majority of impurities generated during USP.
The purity of light phase is usually between 15 ± 2 % and can be improved to 85   ± 5 % through precipitation and the subsequent filtration. Further, the surrounding medium is exchanged in dissolution with a sodium phosphate buffer. Figure 7 shows a SEC chromatogram to visualize the purity improvement gained through the precipitation unit. Hence, the described operation procedure does not only result in a significant gain in purity but also in the strong visible differences of the spectroscopic data between precipitation and dissolution in Figure 6.
Preprocessing is conducted due to evaluation of the scatter effects plots, which are depicted in Figure 8a for precipitation and in Figure 8c for dissolution. Parallel and superimposed scatter effects represent a linear offset in spectral raw data and are reduced with a first derivation using Savitzky Golay with three smoothing points. New scatter effects seem to run through one point with varying slopes and are depicted in Figure 8b,d. Therefore, SNV and MSC are tested. Neither SNV nor MSC gave satisfying result. Hence, data were not preprocessed any further.
In the case of Raman, preprocessing was equal for precipitation and dissolution. The preprocessed spectral data are presented in Figure 9a,c. PLS regression for precipitation was very good (R2 of 0.95) referring to LMWs for HMWs (R2 of 0.62) and the Target (R2 of 0.63) component only poor correlations are found. In dissolution convincing results are achieved for the target component (R2 of 0.93) only. Correlation for LMWs and HMWs could not be found.
Prediction of purity in precipitation and dissolution is very poor for LMWs (R2 of 0.46 and R2 of 0.22) and HMWs (R2 of 0.3 and R2 of 0.13). However, the purity prediction for the target component resulted in R2 of 0.66 for precipitation and R2 of 0.85 for dissolution which is acceptable. It can be concluded that Raman is a suitable detector in precipitation for LMWs and in dissolution for the target component.

3.2. Fourier-Transform Infrared Spectroscopy

The evaluation of the spectral data for the FTIR is carried out in accordance with the procedure in the previous Section 3.1 FTIR spectral raw data are also reduced to 1800–400 cm−1 wavenumber equally for precipitation and dissolution. Analyzing spectral raw data during precipitation, only slight differences are observed which is mainly the magnitude of absorbance. Further preprocessing is performed after evaluation of scatter effects shown Figure 10a,b. Evaluation of scatter effects for precipitation indicate that scatter effects run through one point with different slopes, wherefore, SNV using Savitzky Golay with three smoothing points is chosen as pretreatment. The new scatter effects plot changed slightly which indicates that the analyzed data does not contain a lot of overlaying interferences. In contrast, for dissolution larger deviations can be monitored between wavenumbers 1800–1400 cm−1. Scatter effects for dissolution do not show one of the effects mentioned above, wherefore, no further preprocessing is necessary in case of FTIR data in dissolution.
As an additional analyzation method, a PCA is performed for dissolution raw data, which indicates that most spectral information is located between 1800–850 cm−1. Therefore, spectral data are reduced further to this region (Figure 11b,e). Here again the observation is made, that in precipitation the PLS regression for LMWs is very good with a R2 of 0.94 but is very poor for HMWs (R2 of 0.43) and the target component (R2 of 0.42). In dissolution convincing results are achieved (R2 of 0.91) for the target component and only poor results for HMWs (R2 of 0.31) and LMWs (R2 of 0.54). PLS regression for purity did not lead to any satisfactory results in neither precipitation nor dissolution. Nevertheless, it can be concluded that FTIR is a suitable detector for LMWs during precipitation and for the target component in dissolution.

3.3. Diode-Array Detector

Spectral data of DAD are not reduced and the full spectrum is evaluated ranging from 190 nm to 520 nm. Analyzing spectral raw data for precipitation, arbitrary units decrease for a rising PEG content. This is convincing since less proteins are in solution with an increased PEG amount. A DAD detector can detect proteins reliably, but distinction between different proteins is difficult since most proteins have the highest absorbance at 280 nm. A similar trend can be observed in DAD raw data for dissolution (Figure 12b) but in the opposite direction. Arbitrary units increase steadily with increasing process time, as does the concentration of the re-dissolved antibody. This is to be expected as the absorbance of a DAD increases according to the law of Lambert–Beer.
Evaluation of scatter effects (see Figure 13) lead to a baseline correction in case of precipitation and to a first deviation in case of dissolution (Savitzky Golay with three smoothing points). The DAD results (see Figure 14) fit the previously described observations for FTIR and Raman quite well. The concentration of side components such as LMWs (R2 of 0.97) can be described appropriately during precipitation while spectral data in dissolution are suitable for the prediction of the target component (R2 of 0.93). Additionally, here the prediction of purity is not satisfactory with R2 of 0.57 for LMWs and R2 of 0.59 for TC during precipitation.
It can be concluded that DAD is a suitable detector to monitor LMWs in the precipitation process. Furthermore, DAD is appropriate to display the redissolution of the target component. Reliable determination of purity was not possible.

3.4. Fluorescence Detector

Fluorescence spectra are taken for dissolution only since the precipitates are harmful for the flow cell of the detector. In case of dissolution the entire spectrum ranging from 280 nm to 900 nm is used for analysis. Analyzing the raw data in Figure 15a a fluorescence maximum between 570 and 675 nm can be identified, which describes the concentration increase during dissolution adequately. At the beginning of the dissolution process, the antibody dissolves more rapidly and approaches an equilibrium value over process time. This can also be seen in the arbitrary units of the fluorescence; whose values change proportionally with the increase in concentration measured offline.
Evaluation of scatter effects, which are shown in Figure 16 does not show one of the effects described in Figure 2, hence, no preprocessing of spectral data is performed. The data are just trimmed to 500–900 nm, which is the region of interest. PLS regression resulted in a R2 of 0.90 for the concentration increase of the antibody. Purity could not be correlated to the spectral data. The regression results of the concentration prediction is depicted in Figure 15c.

3.5. Combination of Detectors

The combination of several detectors can lead to prediction accuracy and a higher robustness of the process. For this reason, the prediction quality of a combined PLS, which is trained with data from different measurement methods, is tested. Since the concentration and purity during dissolution are decisive for the next unit, only a combination for dissolution is shown. DAD and Raman have each given the best performance, which are summarized in Table 1 for both detectors.
Regression results of the generated PLS model can be seen in Figure 17a shows the regression results for the concentration of the target component and (b) results for the purity prediction. Both detectors in combination achieved a R2 of 0.90 for concentration and a R2 of 0.72 for purity. This is slightly worse than the results for the detectors alone, but still provides a satisfactory prediction for target component and the purity. Moreover, prediction of purity with a DAD alone is not sufficient. The quality of the data gains in any case because Raman and DAD are based on two different measurement methods and the combination of two orthogonal detectors improves the robustness of online analysis method.

3.6. Additional Process Data of Sensors

In this section further detectors are evaluated within the PAT concept. Sensors can be used to monitor changes during the manufacturing process. Investigation of pH, conductivity and turbidity sensors during precipitation are shown in Figure 18. pH value remains constant throughout precipitation which indicates that no process information can be correlated to this metrology. Light phase has a conductivity of approximately 12 mS/cm in contrast to the PEG solution, which is not conductive. Therefore, conductivity drops as soon as PEG is added and remains constant afterwards. Conductivity indicates the composition of the dispersion. Hence, conductivity can be used to monitor the component ratio between light phase and PEG 4000 online within a PAT concept to ensure optimal precipitation conditions. Turbidity has, as conductivity, a strong deflection during addition of PEG, because precipitate formation occurs instantaneous. For this reason, turbidity is a measure for the precipitation progress. However, turbidity is not specific to one component, but it gives an estimate on process duration.
Utilization of the different sensors during the dissolution step did not lead to any results. The pH value remained constant due to the buffering system, consisting of sodium dihydrogen phosphate and disodium phosphate. Overall conductivity in the system also corresponded to the conductivity of the dissolution buffer. Turbidity also gave no result, which is not surprising, since precipitates go back into solution. According to this study, Raman is recommended as a detector for precipitation, since it has shown good results in the precipitation and the subsequent dissolution with respect to the target component but also side components (LMWs). Additionally, the purity of the target component in the dissolution could be correlated. As an orthogonal measurement strategy, a DAD is recommended. Furthermore, a conductivity probe is recommended as PAT control strategy for detection of optimal precipitation conditions.

4. APC Proposal and Simulation Studies

Variations in product titer and purity of the target component are to be expected in bioprocesses. These may occur, e.g., due to a reduced growth rate of cells or through limitation/inhibition by the substrate. The reasons for titer and purity fluctuations are divers. The question is whether the precipitation unit presented can react to such fluctuations. Further, how online measurement techniques can be integrated to monitor the actual conditions present in the process. The digital twin is fed with the input data after ATPE and the online measured spectral data. The model estimates new process parameters and generates control values to ensure compliance of the critical quality attributes.
In this section an overview of process procedure of the precipitation unit is given as well as a simulation study on reaction possibilities to process deviation. Furthermore, the benefits of the proposed APC are presented.

4.1. Operation of Precipitation Unit

To begin with this section on APC and simulation studies, the function of the precipitation unit is explained for a better process comprehension. Basically, the precipitation unit consists of four subunits which has been defined in Figure 1. In the continuous set up, the individual precipitation tanks of the subunits are merged into one tank. In the following the operation of the continuous unit is explained. Throughout the process each subunit passes through different phases in a certain order, namely filtration of precipitates, followed by redissolution and, finally, a rinsing process to regenerate the membrane. Further, each subunit passes all phases multiple times, but with a time delay to ensure continuous processing of the incoming feed (light phase from ATPE) and to provide a constant output stream. Figure 19 shows the flow sheet of the complete unit including the proposed online measurement technique (online spectral analysis) as well as pressure sensors (P) and mass flow controllers (MFC).
All online measured values are fed back into the digital twin. Figure 20 shows a simplified illustration of the described APC concept. The blue dashed lines represent the signal connection between precipitation unit and digital twin. Hereby, the signal direction is indicated by the arrows. The adjacent boxes show the process parameters that are exchanged between the model and the precipitation unit. The online measurement train between ATPE and precipitation unit determines the initial parameters (IgG concentration, yield, and purity) for the precipitation unit. Further, the mass flow of the incoming feed is measured. Precipitation itself is a concentration independent process [52], therefore, the needed precipitant is calculated dependent on the incoming mass flow. Based on the MFC information the model calculates the mass flow rate of the precipitant (PEG 4000) and the speed of the PEG pump (E-200 in) is aligned. A conductivity probe in the precipitation tank is used to monitor the ratio between feed and precipitant. In the event of a deviation the speed of the PEG pump is adjusted to compensate the mass flow variation of the feed. The turbidity probe monitors the precipitation progress and ensures that the precipitation is completed, and no more antibody remains dissolved in the supernatant. Once the defined level in the precipitation tank is reached and the precipitation is finished (monitored trough turbidity), valves are opened, and filtration starts with the fist module (marked as Subunit 1 in). Filtration is performed in dead-end mode to avoid product loss and to deposit all precipitates on the membrane surface. Pressure sensors are used to measure the transmembrane pressure to prevent pressure peaks out of the operation range of the filters. At a critical pressure valves are switched, and the dispersion is directed to the next module.
The first loaded module moves on to the next phase, the dissolution of the precipitates. Here, the input variables from ATPE and the specification of the final concentration are used to determine model-based the required buffer volume in dissolution. Dissolution is performed in the membrane module by tangentially overflowing the surface area at a high flow rate. Dead-end filtration has compressed the precipitates, which leads to decelerated dissolution kinetics. The experiments have shown that a time of about 14 h is required to reach a yield of >90% antibody recovery after dissolution. The dissolution circuit also contains the measuring equipment described above, which controls the speed of the dissolution buffer pump as well as the valves to the product tank. As soon as the dissolution of the antibody has reached a stationary value, the dissolution is terminated. The module moves on to regeneration, to be ready for the next filtration phase.

4.2. Precipitation Subunit Modeling

The presented simulation study is accomplished with the digital twin presented for one subunit as shown in Figure 1. This study aims to demonstrate that the precipitation unit can react to process variations during manufacturing and is able to provide a constant output concentration. The simulation studies are executed in Aspen Custom Modeler™ (ACM).
In the initial simulation study feed titers ranging from 1 to 3 g/L have been simulated, since this is a realistic fluctuation for the described laboratory process. The final mAb concentration of 3 g/L is attempted to be achieved after dissolution. Figure 21a shows the concentration decrease in the precipitation step over process time for varying incoming feed concentrations. It becomes obvious that complete precipitation results within a few seconds after PEG addition. Thus, concentration variations do not affect the precipitation step if the precipitation conditions of 12 wt % PEG are maintained. Furthermore, a concentration adjustment in this step is not desirable since the aim is to fully precipitate the target and to prevent product loss. This occurs as soon as a quantity of the antibody remains dissolved in the supernatant. Since precipitation takes place in a stirred tank, temporary concentration fluctuations are mitigated and not directly passed on to filtration. An unstable feed concentration does influence the precipitated mass though. Varying mass of precipitates would lead to extended or shortened filtration time. Hence, concentration adjustments must be addressed in dissolution. To provide the constant output concentration the dissolution ratio is modified Figure 21b, a constant mAb concentration of around 2.7 g/L (± 0.05) in the dissolved output can be observed. This was achieved with the before mentioned adjustment of the dissolution ratio. This simulation study demonstrated that precipitation is capable to provide a constant concentration which is important for the following chromatography.

4.3. Reaction to Process Disturbances

To investigate the response capabilities of the unit operation to eventual process deviations, three theoretical process variations are simulated that are realistic variations for the described process. The simulation studies are labeled as scenario 1–3:
  • Concentration variations
  • Mass flow fluctuations
  • Purity changes
As mentioned before product titer and purity can variate in a bioprocess due to growth rates or other complex conditions. Mass flow deviation may occur through pump failures or can be caused during ATPE. The expected values for the process are a mAb titer of 2 g/L, an ATPE mass flow rate of 1 g/min as well as the purity after ATPE (15 ± 2%). All process variations are discussed for the precipitation step and the dissolution separately.
The first study within scenario one includes the proposed APC concept to simulate a process fluctuation during continuous operation. For this purpose, a periodic concentration change is used as input signal, which fluctuates from 2.5 to 1.5 g/L. A final concentration of 2 g/L is attempted to be achieved after dissolution. Figure 22a demonstrates the simulation results with and without the APC. Without APC the concentration fluctuations are still obvious after dissolution (grey line). But if the dissolution ratio is adjusted through APC the unit can provide a constant mAb concentration of 1.85 (±0.03 g/L) besides any feed fluctuations (red line). Figure 22b visualizes which parameters must be controlled to ensure a consistent concentration. By adjusting the dissolution ratio (brown columns), either slight concentration of the target component or dilution is achieved to result in a constant output concentration. Associated with this is the change in the residual PEG content (orange columns), which can become limiting due to the thermodynamic equilibrium. This is shown in Figure 22c. The grey line depicts the equilibrium concentration of mAb in presence of an increasing PEG content. The red line represents the solubility of IgG in the feed at different PEG concentrations. The blue shaded area indicates the thermodynamically given operating window, within which the laboratory process can be carried out without. Outside of this operation window limitation by PEG occur during redissolution, which leads to product loss. It becomes obvious that a PEG content above 5 wt % limits the resolution of the antibody (see Figure 22c). In an industrial environment, this window must be selected the more narrowly, due to higher titers that are achieved during Upstream processing. Then, the remaining PEG should not exceed a proportion of 3 wt % [13].
In the following section simulation studies regarding scenario two, which address mass flow fluctuations, are presented. Initial simulations have shown that volume fluctuations mainly affect the precipitation step. This occurs because the supernatant is removed by filtration and is, therefore, not passed on to the subsequent dissolution. The simulations have also shown that an elevated PEG content does not lead to yield loss, hence, here only the simulation results in case of an insufficient PEG content are shown. points out the potential failure mode with and without the integration of the proposed APC. Changes in the mass flow impact the present precipitation conditions in the precipitation tank, which are monitored with the conductivity probe. If this disturbance occurs over a longer period than a few seconds a wrong ratio between PEG and light phase would be present in the precipitation tank in relation to its residence time. This case is shown in Figure 23a. During operation, the required PEG content of 12 wt % can no longer be maintained (orange columns), which leads to an increase in the equilibrium concentration of IgG (light blue line). As a result, the target component no longer precipitates completely, and a residual portion remains dissolved which are potential yield losses. This dissolved residue would then be discarded through dead-end filtration which is not desirable. This loss is indicated by the decreased yield in the precipitation step (blue columns).
The second simulation study within scenario two addresses also the continuous operation with integrated APC (see Figure 23b). In this case the PEG mass flow is adjusted as soon as the PEG content is decreasing in the precipitation tank. This results in a brief increase of the equilibrium concentration in relation to the residence time (10 min). Afterwards the operating point with the correct equilibrium concentration is reached again, and thus no more product loss occurs. For this study the difference in yield loss is about 13%.
Scenario three was also investigated, but the compensation of a purity changes in the feed is not possible with the presented precipitation unit. The reason for this is that in the precipitation and dissolution step, an equilibrium is reached. It is not possible to influence the extent to which the side components precipitate or redissolve. Hence, no specific control parameter is available to adjust the purity. Nonetheless, the purity gain of this unit operation is enormous. Starting with an initial purity of 15 ± 2% after ATPE the purity is increased to an extent of 85 ± 5% after dissolution. This purity increase measured with SEC chromatography is previously shown in Figure 7.
Concluding it can be said that concentration and mass flow changes can both be compensated in the precipitation unit within a short time span of a few seconds. This is reached through the integration of the proposed APC. Nonetheless, both failure modes cannot be compensated at once, since a constant mAb output concentration is always accompanied by a fluctuation in the output mass flow. So, either concentration can be compensated or the mass flow. All simulation results are summarized in Table 2.

5. Conclusions

It has been shown that precipitation as a unit operation is very robust and achieves a high depletion of the side components due to the mode of operation. The purity is increased from 15 ± 2% after ATPE to 85 ± 5% after dissolution. The remaining impurities mainly consisting of LMWs and PEG, which can be separated easily in flow-through mode in the following chromatography (IEX). The dissolution buffer is at the same time the running buffer of the subsequent chromatography, which makes the transfer simple and no additional diafiltration is required. All in all, shows that precipitation is suitable for purification of mAb and can be ideally integrated into the alternative manufacturing process between extraction and chromatography. A further advantage is that no cost-intensive materials or equipment are needed. On the contrary, all used components and equipment devices are already standard in industry and can simply be scaled to other scales.
Moreover, the feasibility of using online measurement technology has been demonstrated within the framework of the PAT initiative. All detectors delivered very good results in the analysis of product concentration. However, the best result was achieved by Raman, which reliably detected product concentration and purity online. Furthermore, the combination of DAD and Raman has been tested, which is a recommended extension due to the orthogonal measurement methods, resulting in higher process robustness. Finally, in a simulation study, which has been performed with aid of a digital twin, the online measurement data have been successfully used to control the process. This led to minimized yield losses at potential feed flow or concentration fluctuations. Further studies will focus on continuous operation mode and transfer the process to other biologics such as fragments, peptides, pDNA and VLPs at the Institute.

Author Contributions

L.J.L. conceived, designed, and performed the experiments as well as wrote the paper; all authors interpreted the data; J.S. substantively revised the work and contributed the materials and analysis tools; J.S. is responsible for the conception and supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This research received funding from BMWi project “Traceless Plant Traceless Production”. We also acknowledge financial support by Open Access Publishing Fund of Clausthal University of Technology.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors would like to acknowledge their institute’s laboratory colleagues, especially Frank Steinhäuser and Volker Strohmeyer as well as Thomas Knebel. In addition, the authors would like to thank BMWI, especially M. Gahr, for project funding “Traceless Plant Traceless Production” and the whole TPTP consortium. The authors would also like to thank Heribert Helgers and Florian Vetter for fruitful discussion, Axel Schmidt for his support with ATPE and the discussion of experimentally obtained results, as well as Mourad Mouellef for his support with multivariate data analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kornecki, M.; Strube, J. Accelerating Biologics Manufacturing by Upstream Process Modelling. Processes 2019, 7, 166. [Google Scholar] [CrossRef] [Green Version]
  2. Subramanian, G. Continuous Biomanufacturing—Innovative Technologies and Methods; Wiley-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2017; ISBN 9783527699902. [Google Scholar]
  3. Zobel-Roos, S.; Schmidt, A.; Mestmäcker, F.; Mouellef, M.; Huter, M.; Uhlenbrock, L.; Kornecki, M.; Lohmann, L.; Ditz, R.; Strube, J. Accelerating Biologics Manufacturing by Modeling or: Is Approval under the QbD and PAT Approaches Demanded by Authorities Acceptable Without a Digital-Twin? Processes 2019, 7, 94. [Google Scholar] [CrossRef] [Green Version]
  4. Lohmann, L.J.; Strube, J. Accelerating Biologics Manufacturing by Modeling: Process Integration of Precipitation in mAb Downstream Processing. Processes 2020, 8, 58. [Google Scholar] [CrossRef] [Green Version]
  5. Kornecki, M.; Schmidt, A.; Lohmann, L.; Huter, M.; Mestmäcker, F.; Klepzig, L.; Strube, J. Accelerating Biomanufacturing by Modeling of Continuous Bioprocessing—Piloting Case Study of Monoclonal Antibody Manufacturing. Processes 2019, 7, 495. [Google Scholar] [CrossRef] [Green Version]
  6. Helgers, H.; Schmidt, A.; Lohmann, L.J.; Vetter, F.L.; Juckers, A.; Jensch, C.; Mouellef, M.; Zobel-Roos, S.; Strube, J. Towards Autonomous Operation by Advanced Process Control—Process Analytical Technology for Continuous Biologics Antibody Manufacturing. Processes 2021, 9, 172. [Google Scholar] [CrossRef]
  7. Oelmeier, S.A.; Ladd-Effio, C.; Hubbuch, J. Alternative separation steps for monoclonal antibody purification: Combination of centrifugal partitioning chromatography and precipitation. J. Chromatogr. A 2013, 1319, 118–126. [Google Scholar] [CrossRef]
  8. Hammerschmidt, N.; Hobiger, S.; Jungbauer, A. Continuous polyethylene glycol precipitation of recombinant antibodies: Sequential precipitation and resolubilization. Process. Biochem. 2016, 51, 325–332. [Google Scholar] [CrossRef] [Green Version]
  9. Burgstaller, D.; Jungbauer, A.; Satzer, P. Continuous integrated antibody precipitation with two-stage tangential flow microfiltration enables constant mass flow. Biotechnol. Bioeng. 2019, 116, 1053–1065. [Google Scholar] [CrossRef] [Green Version]
  10. Dutra, G.; Komuczki, D.; Jungbauer, A.; Satzer, P. Continuous capture of recombinant antibodies by ZnCl 2 precipitation without polyethylene glycol. Eng. Life Sci. 2020, 20, 265–274. [Google Scholar] [CrossRef] [Green Version]
  11. Gu, Q.; Li, Z.; Coffman, J.L.; Przybycien, T.M.; Zydney, A.L. High throughput solubility and redissolution screening for antibody purification via combined PEG and zinc chloride precipitation. Biotechnol. Prog. 2020, e3041. [Google Scholar] [CrossRef]
  12. Li, Z.; Chen, T.-H.; Andini, E.; Coffman, J.L.; Przybycien, T.; Zydney, A.L. Enhanced filtration performance using feed-and-bleed configuration for purification of antibody precipitates. Biotechnol. Prog. 2020, e3082. [Google Scholar] [CrossRef] [PubMed]
  13. Matheus, S.; Friess, W.; Schwartz, D.; Mahler, H.-C. Liquid high concentration IgG1 antibody formulations by precipitation. J. Pharm. Sci. 2009, 98, 3043–3057. [Google Scholar] [CrossRef]
  14. Tscheliessnig, A.; Satzer, P.; Hammerschmidt, N.; Schulz, H.; Helk, B.; Jungbauer, A. Ethanol precipitation for purification of recombinant antibodies. J. Biotechnol. 2014, 188, 17–28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Großhans, S.; Wang, G.; Fischer, C.; Hubbuch, J. An integrated precipitation and ion-exchange chromatography process for antibody manufacturing: Process development strategy and continuous chromatography exploration. J. Chromatogr. A 2018, 1533, 66–76. [Google Scholar] [CrossRef]
  16. Hammerschmidt, N.; Hintersteiner, B.; Lingg, N.; Jungbauer, A. Continuous precipitation of IgG from CHO cell culture supernatant in a tubular reactor. Biotechnol. J. 2015, 10, 1196–1205. [Google Scholar] [CrossRef] [PubMed]
  17. Li, Z.; Gu, Q.; Coffman, J.L.; Przybycien, T.; Zydney, A.L. Continuous precipitation for monoclonal antibody capture using countercurrent washing by microfiltration. Biotechnol. Prog. 2019, 35, e2886. [Google Scholar] [CrossRef] [PubMed]
  18. Cohn, E.J.; Strong, L.E. Preparation and properties of serum and plasma proteins; a system for the separation into fractions of the protein and lipoprotein components of biological tissues and fluids. J. Am. Chem. Soc. 1946, 68, 459–475. [Google Scholar] [CrossRef]
  19. Bertolini, J. The purification of plasma proteins for therapeutic use. In Rossi’s Principles of Transfusion Medicine; John Wiley & Sons: Hoboken, NJ, USA, 2016; pp. 302–320. [Google Scholar] [CrossRef]
  20. Rathore, A.S.; Winkle, H. Quality by Design for Biopharmaceuticals. Nat. Biotechnol. 2009, 27, 26–34. [Google Scholar] [CrossRef]
  21. ICH. ICH Q11 Guideline: Development and Manufacture of Drug Substances (Chemical Entities and Biotechnological/Biological Entities): Questions and Answers. 2017. Available online: https://database.ich.org/sites/default/files/Q11_Q%26As_Q%26As.pdf (accessed on 2 March 2021).
  22. ICH. ICH Guideline Q11 on Development and Manufacture of Drug Substances (Chemical Entities and Biotechnological/Biological Entities). 2011. Available online: https://www.ema.europa.eu/en/documents/scientific-guideline/draft-ich-guideline-q11-development-manufacture-drug-substances-chemical-entities-biotechnological/biological-entities_en.pdf (accessed on 2 March 2021).
  23. ICH. ICH-Endorsed Guide for ICH Q8/Q9/Q10 Implementation. 2006. Available online: https://database.ich.org/sites/default/files/Q8_Q9_Q10_Q%26As_R4_Points_to_Consider_0.pdf (accessed on 2 March 2021).
  24. ICH. Q12 Technical and Regulatory Considerations for Pharmaceutical Product Lifecycle Mamanagement Annexes. Available online: https://www.fda.gov/regulatory-information/search-fda-guidance-documents/q12-technical-and-regulatory-considerations-pharmaceutical-product-lifecycle-management-annex (accessed on 2 March 2021).
  25. EMA. EMA-FDA Pilot Program for Parallel Assessment of Quality-by-Design Applications: Lessons Learnt and Q&A Resulting from the First Parallel Assessment. 2013. Available online: https://www.ema.europa.eu/en/documents/other/european-medicines-agency-food-drug-administration-pilot-programme-parallel-assessment-quality_en.pdf (accessed on 2 March 2021).
  26. Jin, Y.; Wu, Z.; Liu, X.; Wu, Y. Near infrared spectroscopy in combination with chemometrics as a process analytical technology (PAT) tool for on-line quantitative monitoring of alcohol precipitation. J. Pharm. Biomed. Anal. 2013, 77, 32–39. [Google Scholar] [CrossRef]
  27. Wu, Z.; Xu, B.; Du, M.; Sui, C.; Shi, X.; Qiao, Y. Validation of a NIR quantification method for the determination of chlorogenic acid in Lonicera japonica solution in ethanol precipitation process. J. Pharm. Biomed. Anal. 2012, 62, 1–6. [Google Scholar] [CrossRef]
  28. Huang, H.; Qu, H. In-line monitoring of alcohol precipitation by near-infrared spectroscopy in conjunction with multivariate batch modeling. Anal. Chim. Acta 2011, 707, 47–56. [Google Scholar] [CrossRef] [PubMed]
  29. Sun, Z.; Wang, J.; Nie, L.; Li, L.; Cao, D.; Fan, J.; Wang, H.; Liu, R.; Zhang, Y.; Zang, H. Calibration transfer of near infrared spectrometers for the assessment of plasma ethanol precipitation process. Chemom. Intell. Lab. Syst. 2018, 181, 64–71. [Google Scholar] [CrossRef]
  30. Xu, B.; Wu, Z.; Lin, Z.; Sui, C.; Shi, X.; Qiao, Y. NIR analysis for batch process of ethanol precipitation coupled with a new calibration model updating strategy. Anal. Chim. Acta 2012, 720, 22–28. [Google Scholar] [CrossRef] [PubMed]
  31. Li, L.; Ding, B.; Yang, Q.; Chen, S.; Ren, H.; Wang, J.; Zang, H.; Wang, F.; Zang, L. The relevance study of effective information between near infrared spectroscopy and chondroitin sulfate in ethanol precipitation process. J. Innov. Opt. Health Sci. 2014, 7, 1450022. [Google Scholar] [CrossRef] [Green Version]
  32. Huang, H.; Qu, H. In-situ monitoring of saccharides removal of alcohol precipitation using near-infrared spectroscopy. J. Innov. Opt. Health Sci. 2018, 11, 1850027. [Google Scholar] [CrossRef] [Green Version]
  33. Xu, B.; Lin, Z.; Shi, X.; Qiao, Y.; Du, M.; Wu, Z. NIR Determination of Three Critical Quality Attributes in Alcohol Precipitation Process of Lonicerae Japonicae with Uncertainty Analysis. In Proceedings of the IEEE International Conference on Biomedical Engineering and Biotehnology, Macau, China, 28–30 May 2012; pp. 1566–1571. [Google Scholar] [CrossRef]
  34. Sun, Z.; Fan, J.; Wang, J.; Wang, F.; Nie, L.; Li, L.; Dong, Q.; Li, C.; Du, R.; Quan, S.; et al. Assessment of the human albumin in acid precipitation process using NIRS and multi-variable selection methods combined with SPA. J. Mol. Struct. 2020, 1199, 126942. [Google Scholar] [CrossRef]
  35. Hebbi, V.; Thakur, G.; Rathore, A.S. Process analytical technology application for protein PEGylation using near infrared spectroscopy: G-CSF as a case study. J. Biotechnol. 2020. [Google Scholar] [CrossRef] [PubMed]
  36. Meyer-Kirschner, J.; Kather, M.; Pich, A.; Engel, D.; Marquardt, W.; Viell, J.; Mitsos, A. In-line Monitoring of Monomer and Polymer Content During Microgel Synthesis Using Precipitation Polymerization via Raman Spectroscopy and Indirect Hard Modeling. Appl. Spectrosc. 2016, 70, 416–426. [Google Scholar] [CrossRef]
  37. Schöll, J.; Lindenberg, C.; Vicum, L.; Brozio, J.; Mazzotti, M. Precipitation of alpha L-glutamic acid: Determination of growth kinetics. Faraday Discuss. 2007, 136, 247–264; discussion 309–328. [Google Scholar] [CrossRef]
  38. Wu, H.; White, M.; Khan, M.A. Quality-by-Design (QbD): An integrated process analytical technology (PAT) approach for a dynamic pharmaceutical co-precipitation process characterization and process design space development. Int. J. Pharm. 2011, 405, 63–78. [Google Scholar] [CrossRef]
  39. Feidl, F.; Garbellini, S.; Vogg, S.; Sokolov, M.; Souquet, J.; Broly, H.; Butté, A.; Morbidelli, M. A new flow cell and chemometric protocol for implementing in-line Raman spectroscopy in chromatography. Biotechnol. Prog. 2019, 35, e2847. [Google Scholar] [CrossRef] [PubMed]
  40. Kornecki, M.; Schmidt, A.; Strube, J. Pat as Key-Enabling Technology for Qbd in Pharmaceutical Manufacturing—A Conceptual Review on Upstream and Downstream Processing. Chem. Today 2018, 44–48. [Google Scholar]
  41. Kornecki, M.; Strube, J. Process Analytical Technology for Advanced Process Control in Biologics Manufacturing with the Aid of Macroscopic Kinetic Modeling. Bioengineering 2018, 5, 25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Buckley, K.; Ryder, A.G. Applications of Raman Spectroscopy in Biopharmaceutical Manufacturing: A Short Review. Appl. Spectrosc. 2017, 71, 1085–1116. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, Y.; Boysen, R.I.; Wood, B.R.; Kansiz, M.; McNaughton, D.; Hearn, M.T.W. Determination of the secondary structure of proteins in different environments by FTIR-ATR spectroscopy and PLS regression. Biopolymers 2008, 89, 895–905. [Google Scholar] [CrossRef] [PubMed]
  44. Boulet-Audet, M.; Kazarian, S.G.; Byrne, B. In-column ATR-FTIR spectroscopy to monitor affinity chromatography purification of monoclonal antibodies. Sci. Rep. 2016, 6, 30526. [Google Scholar] [CrossRef] [Green Version]
  45. Walther, C.; Mayer, S.; Jungbauer, A.; Dürauer, A. Getting ready for PAT: Scale up and inline monitoring of protein refolding of Npro fusion proteins. Process. Biochem. 2014, 49, 1113–1121. [Google Scholar] [CrossRef]
  46. Großhans, S.; Rüdt, M.; Sanden, A.; Brestrich, N.; Morgenstern, J.; Heissler, S.; Hubbuch, J. In-line Fourier-transform infrared spectroscopy as a versatile process analytical technology for preparative protein chromatography. J. Chromatogr. A 2018, 1547, 37–44. [Google Scholar] [CrossRef] [PubMed]
  47. Brestrich, N.; Rüdt, M.; Büchler, D.; Hubbuch, J. Selective protein quantification for preparative chromatography using variable pathlength UV/Vis spectroscopy and partial least squares regression. Chem. Eng. Sci. 2018, 176, 157–164. [Google Scholar] [CrossRef]
  48. Sauer, D.G.; Melcher, M.; Mosor, M.; Walch, N.; Berkemeyer, M.; Scharl-Hirsch, T.; Leisch, F.; Jungbauer, A.; Dürauer, A. Real-time monitoring and model-based prediction of purity and quantity during a chromatographic capture of fibroblast growth factor 2. Biotechnol. Bioeng. 2019, 116, 1999–2009. [Google Scholar] [CrossRef]
  49. Zobel-Roos, S.; Mouellef, M.; Siemers, C.; Strube, J. Process Analytical Approach towards Quality Controlled Process Automation for the Downstream of Protein Mixtures by Inline Concentration Measurements Based on Ultraviolet/Visible Light (UV/VIS) Spectral Analysis. Antibodies 2017, 6, 24. [Google Scholar] [CrossRef] [Green Version]
  50. Brestrich, N.; Briskot, T.; Osberghaus, A.; Hubbuch, J. A tool for selective inline quantification of co-eluting proteins in chromatography using spectral analysis and partial least squares regression. Biotechnol. Bioeng. 2014, 111, 1365–1373. [Google Scholar] [CrossRef] [PubMed]
  51. Schmidt, A.; Richter, M.; Rudolph, F.; Strube, J. Integration of Aqueous Two-Phase Extraction as Cell Harvest and Capture Operation in the Manufacturing Process of Monoclonal Antibodies. Antibodies 2017, 6, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Martins, D.L.; Sencar, J.; Hammerschmidt, N.; Flicker, A.; Kindermann, J.; Kreil, T.R.; Jungbauer, A. Truly continuous low pH viral inactivation for biopharmaceutical process integration. Biotechnol. Bioeng. 2020, 117, 1406–1417. [Google Scholar] [CrossRef]
  53. Robinson; Earnshaw. Experimental study of colloidal aggregation in two dimensions. I. Structural aspects. Phys. Rev. A 1992, 46, 2045–2054. [Google Scholar] [CrossRef] [PubMed]
  54. Neil, C.W.; Ray, J.R.; Lee, B.; Jun, Y.-S. Fractal aggregation and disaggregation of newly formed iron(iii ) (hydr)oxide nanoparticles in the presence of natural organic matter and arsenic. Environ. Sci. Nano 2016, 3, 647–656. [Google Scholar] [CrossRef] [Green Version]
  55. Satzer, P.; Burgstaller, D.; Krepper, W.; Jungbauer, A. Fractal dimension of antibody-PEG precipitate: Light microscopy for the reconstruction of 3D precipitate structures. Eng. Life Sci. 2020, 20, 67–78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Müller-Erlwein, E. (Ed.) Chemische Reaktionstechnik; Teubner: Wiesbaden, Germany, 2007; ISBN 978-3-8351-0187-6. [Google Scholar]
  57. Sim, S.-L.; He, T.; Tscheliessnig, A.; Mueller, M.; Tan, R.B.H.; Jungbauer, A. Protein precipitation by polyethylene glycol: A generalized model based on hydrodynamic radius. J. Biotechnol. 2012, 157, 315–319. [Google Scholar] [CrossRef]
  58. Juckles, I. Fractionation of proteins and viruses with polyethylene glycol. Biochim. Biophys. Acta 1971, 229, 535–546. [Google Scholar] [CrossRef]
  59. Atha, D.H.; Inghamm, K.C. Mechanism of precipitation of Proteins by Polyethylen Glycols. J. Biol. Chem. 1981, 256, 12108–12117. [Google Scholar] [CrossRef]
  60. Thakur, G.; Hebbi, V.; Parida, S.; Rathore, A.S. Automation of Dead End Filtration: An Enabler for Continuous Processing of Biotherapeutics. Front. Bioeng. Biotechnol. 2020, 8, 758. [Google Scholar] [CrossRef] [PubMed]
  61. Van Reis, R.; Zydney, A. Bioprocess membrane technology. J. Membr. Sci. 2007, 297, 16–50. [Google Scholar] [CrossRef]
  62. Iritani, E.; Katagiri, N. Developments of Blocking Filtration Model in Membrane Filtration. KONA 2016, 33, 179–202. [Google Scholar] [CrossRef] [Green Version]
  63. Sampath, M.; Shukla, A.; Rathore, A.S. Modeling of Filtration Processes-Microfiltration and Depth Filtration for Harvest of a Therapeutic Protein Expressed in Pichia pastoris at Constant Pressure. Bioengineering 2014, 1, 260–277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Ho, C.C.; Zydney, A.L. A Combined Pore Blockage and Cake Filtration Model for Protein Fouling during Microfiltration. J. Colloid Interface Sci. 2000, 232, 389–399. [Google Scholar] [CrossRef]
  65. Goldrick, S.; Joseph, A.; Mollet, M.; Turner, R.; Gruber, D.; Farid, S.S.; Titchener-Hooker, N.J. Predicting performance of constant flow depth filtration using constant pressure filtration data. J. Membr. Sci. 2017, 531, 138–147. [Google Scholar] [CrossRef]
  66. Cheryan, M.; Strauss, S. Ultrafiltration and Microfiltration Handbook, 2nd ed.; Chapman and Hall/CRC: Boca Raton, FL, USA, 1998; ISBN 1566765986. [Google Scholar]
  67. Mersmann, A.; Kind, M.; Stichlmair, J. Thermische Verfahrenstechnik: Grundlagen und Methoden; Springer: Berlin, Germany, 2005; ISBN 3-540-23648-1. [Google Scholar]
  68. Leveque, M.A. Les Lois de la Transmission de Chaleur par Convection; Dunod: Berlin, Germany, 1928. [Google Scholar]
Figure 1. Simplified illustration of a precipitation subunit with the tested online measurement.
Figure 1. Simplified illustration of a precipitation subunit with the tested online measurement.
Processes 09 00488 g001
Figure 2. Exemplary scatter effects plots: (a) shows additive effects, (b) scatter effects, (c) a combination of additive and scatter effects and (d) shows complex effects [6]. The different colors represent the scatter effect of different spectral data.
Figure 2. Exemplary scatter effects plots: (a) shows additive effects, (b) scatter effects, (c) a combination of additive and scatter effects and (d) shows complex effects [6]. The different colors represent the scatter effect of different spectral data.
Processes 09 00488 g002
Figure 3. Representation of the areas that were used for the material balance in the developed process model. Red represents the balance room for precipitation, yellow for filtration and green for dissolution.
Figure 3. Representation of the areas that were used for the material balance in the developed process model. Red represents the balance room for precipitation, yellow for filtration and green for dissolution.
Processes 09 00488 g003
Figure 4. Light microscopy images of precipitates. (a) shows a magnification of 10-fold and (b) of 50-fold.
Figure 4. Light microscopy images of precipitates. (a) shows a magnification of 10-fold and (b) of 50-fold.
Processes 09 00488 g004
Figure 5. Common blocking mechanisms for dead-end filtration: (a) complete blocking, (b) standard blocking, (c) intermediate blocking, (d) cake filtration.
Figure 5. Common blocking mechanisms for dead-end filtration: (a) complete blocking, (b) standard blocking, (c) intermediate blocking, (d) cake filtration.
Processes 09 00488 g005
Figure 6. Raman spectral raw data for precipitation (a) and dissolution (b).
Figure 6. Raman spectral raw data for precipitation (a) and dissolution (b).
Processes 09 00488 g006
Figure 7. Size exclusion chromatography (SEC) chromatogram of the feed obtained after aqueous-two phase extraction (ATPE) (gray). Supernatant after precipitation is eliminated through dead end filtration. Redissolved mAb in sodium phosphate buffer I depicted in green.
Figure 7. Size exclusion chromatography (SEC) chromatogram of the feed obtained after aqueous-two phase extraction (ATPE) (gray). Supernatant after precipitation is eliminated through dead end filtration. Redissolved mAb in sodium phosphate buffer I depicted in green.
Processes 09 00488 g007
Figure 8. Raman scatter effects plot for evaluation of needed preprocessing of spectral data. (a,c) Show scatter effects of raw spectral data and (b,d) show scatter effects after preprocessing. The sample order is sorted from red to blue.
Figure 8. Raman scatter effects plot for evaluation of needed preprocessing of spectral data. (a,c) Show scatter effects of raw spectral data and (b,d) show scatter effects after preprocessing. The sample order is sorted from red to blue.
Processes 09 00488 g008
Figure 9. Raman spectral data for precipitation (a,b) and dissolution (c,d). Preprocessed spectral data are presented in (a,c). (b,d) show the PLS regression results of the prediction of the investigated CQAs. Different colors represent the spectral data for different samples and are sorted from red to blue. Sample names are equal to Figure 6. Spectral data (c,d) have been previously published [6].
Figure 9. Raman spectral data for precipitation (a,b) and dissolution (c,d). Preprocessed spectral data are presented in (a,c). (b,d) show the PLS regression results of the prediction of the investigated CQAs. Different colors represent the spectral data for different samples and are sorted from red to blue. Sample names are equal to Figure 6. Spectral data (c,d) have been previously published [6].
Processes 09 00488 g009
Figure 10. FTIR scatter effects plot for evaluation of needed preprocessing of spectral data. Scatter effects of raw spectral data (a) and preprocessed data (b) of precipitation are shown as well as the scatter effects of raw data in dissolution (c). No further preprocessing was needed for FTIR raw data for dissolution (d).
Figure 10. FTIR scatter effects plot for evaluation of needed preprocessing of spectral data. Scatter effects of raw spectral data (a) and preprocessed data (b) of precipitation are shown as well as the scatter effects of raw data in dissolution (c). No further preprocessing was needed for FTIR raw data for dissolution (d).
Processes 09 00488 g010
Figure 11. FTIR spectral data for precipitation (ac) and dissolution (df). (a,c) raw data, (b,e) preprocessed data and (c,f) show the PLS regression results. The sample order is sorted from red to blue. Sample names are equal in each row and are always shown in (a,d). Spectral data (df) have been previously published [6].
Figure 11. FTIR spectral data for precipitation (ac) and dissolution (df). (a,c) raw data, (b,e) preprocessed data and (c,f) show the PLS regression results. The sample order is sorted from red to blue. Sample names are equal in each row and are always shown in (a,d). Spectral data (df) have been previously published [6].
Processes 09 00488 g011
Figure 12. Spectral raw data of DAD for precipitation (a) and dissolution (b).
Figure 12. Spectral raw data of DAD for precipitation (a) and dissolution (b).
Processes 09 00488 g012
Figure 13. Diode array detector (DAD) scatter effects plot for evaluation of needed preprocessing of spectral data. (a,c) show scatter effects of raw spectral data and (b,d) show scatter effects after preprocessing.
Figure 13. Diode array detector (DAD) scatter effects plot for evaluation of needed preprocessing of spectral data. (a,c) show scatter effects of raw spectral data and (b,d) show scatter effects after preprocessing.
Processes 09 00488 g013
Figure 14. DAD spectral data for precipitation (a,b) and dissolution (c,d). Preprocessed spectral data are presented in (a,c). (b,d) show the PLS regression results of the prediction of the investigated CQAs. Different colors represent the spectral data for different samples and are sorted from red to blue. Sample names are equal to Figure 12. Spectral data (c,d) have been previously published [6].
Figure 14. DAD spectral data for precipitation (a,b) and dissolution (c,d). Preprocessed spectral data are presented in (a,c). (b,d) show the PLS regression results of the prediction of the investigated CQAs. Different colors represent the spectral data for different samples and are sorted from red to blue. Sample names are equal to Figure 12. Spectral data (c,d) have been previously published [6].
Processes 09 00488 g014
Figure 15. Fluorescence spectral data for dissolution. (a) raw data, (b) data are trimmed to the region of interest and (c) show the PLS regression results. The sample order is sorted from red to blue. Sample names are equal in each row and are shown in (a,b). Spectral data (a) has been previously published [6].
Figure 15. Fluorescence spectral data for dissolution. (a) raw data, (b) data are trimmed to the region of interest and (c) show the PLS regression results. The sample order is sorted from red to blue. Sample names are equal in each row and are shown in (a,b). Spectral data (a) has been previously published [6].
Processes 09 00488 g015
Figure 16. Fluorescence spectral data for dissolution. (a) raw data, (b) preprocessed data. In this case no further preprocessing was necessary due to evaluation of scatter effects. Spectral data (a) and (b) have been previously published [6].
Figure 16. Fluorescence spectral data for dissolution. (a) raw data, (b) preprocessed data. In this case no further preprocessing was necessary due to evaluation of scatter effects. Spectral data (a) and (b) have been previously published [6].
Processes 09 00488 g016
Figure 17. Prediction of process data using combination of DAD and Raman. IgG concentration (a), purity of the target component (b).
Figure 17. Prediction of process data using combination of DAD and Raman. IgG concentration (a), purity of the target component (b).
Processes 09 00488 g017
Figure 18. Course of conductivity, and turbidity probes during precipitation.
Figure 18. Course of conductivity, and turbidity probes during precipitation.
Processes 09 00488 g018
Figure 19. Process control strategy for the precipitation unit including online spectral analysis measuring points. Abbreviation: MFC = mass flow controller; Cond = conductivity probe; Turb = turbidity probe; R = Raman; FTIR = Fourier transform infrared spectroscopy; DAD = diode array detector; F = fluorescence; P = pressure sensor; L = level sensor. The presented APC concept has shortly been published previously [6].
Figure 19. Process control strategy for the precipitation unit including online spectral analysis measuring points. Abbreviation: MFC = mass flow controller; Cond = conductivity probe; Turb = turbidity probe; R = Raman; FTIR = Fourier transform infrared spectroscopy; DAD = diode array detector; F = fluorescence; P = pressure sensor; L = level sensor. The presented APC concept has shortly been published previously [6].
Processes 09 00488 g019
Figure 20. Simplified illustration of proposed advanced process control (APC) for the precipitation unit. Precipitation tanks are merged into one tank. Therefore, subunits mean the filter modules.
Figure 20. Simplified illustration of proposed advanced process control (APC) for the precipitation unit. Precipitation tanks are merged into one tank. Therefore, subunits mean the filter modules.
Processes 09 00488 g020
Figure 21. Simulation results of scenario one (concentration fluctuations) are presented for precipitation (a) and for dissolution (b). These results have shortly been published previously [6].
Figure 21. Simulation results of scenario one (concentration fluctuations) are presented for precipitation (a) and for dissolution (b). These results have shortly been published previously [6].
Processes 09 00488 g021
Figure 22. Second simulation study results of scenario one (concentration fluctuations). (a) shows the comparison of incoming concentration and the adjusted IgG concentration after dissolution because of the integrated APC. (b) shows the parameter that must be controlled to achieve the constant output concentration of 2 g/L and (c) depicts the equilibrium concentration of IgG in presence of PEG. The blue zone in (c) marks the operating window of residual PEG content.
Figure 22. Second simulation study results of scenario one (concentration fluctuations). (a) shows the comparison of incoming concentration and the adjusted IgG concentration after dissolution because of the integrated APC. (b) shows the parameter that must be controlled to achieve the constant output concentration of 2 g/L and (c) depicts the equilibrium concentration of IgG in presence of PEG. The blue zone in (c) marks the operating window of residual PEG content.
Processes 09 00488 g022
Figure 23. Simulation studies regarding volume deviation during the process. (a) shows the effect of volume deviations and its impact on process yield. In (b) the same scenario is shown, but with the integration of the proposed APC.
Figure 23. Simulation studies regarding volume deviation during the process. (a) shows the effect of volume deviations and its impact on process yield. In (b) the same scenario is shown, but with the integration of the proposed APC.
Processes 09 00488 g023
Table 1. PLS results of process analytical technology (PAT) feasibility study for Raman and DAD referring to mAb concentration and purity.
Table 1. PLS results of process analytical technology (PAT) feasibility study for Raman and DAD referring to mAb concentration and purity.
ConcentrationPurity
Raman0.93 ± 3.9%0.85 ± 4.9%
DAD0.93 ± 2.8%0.59 ± 2.6%
Combination0.90 ± 3.9%0.72 ± 4.9%
Table 2. Summary of simulations results.
Table 2. Summary of simulations results.
Scenario
(Deviation ± 50%)
PrecipitationDissolution
Concentration
(a) high
(b) low
(a + b) Concentration deviation does not have impact on precipitation.
→ Precipitant mass flow is calculated based on Feed mass flow
(a) Volume dissolution buffer ↑
(b) Volume dissolution buffer ↓
→ Solubility limits volume reduction of dissolution buffer downwards
→ Control Parameter: mass flow
→ Control Strategy: pump speed
Feed mass flow
(a) Increased
(b) decreased
(a) PEG Stream ↑
(b) PEG Stream ↓
→ Control Parameter: Conductivity
→ Control Strategy: pump speed
(a) Volume dissolution buffer ↓
(b) Volume dissolution buffer ↑
→ Control Parameter: Conductivity
→ Control Strategy: pump speed
Purity of IgG
(a) Increased
(b) decreased
Precipitation does not have an impact on purity of the targetDissolution does not have an impact on purity of the target
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lohmann, L.J.; Strube, J. Process Analytical Technology for Precipitation Process Integration into Biologics Manufacturing towards Autonomous Operation—mAb Case Study. Processes 2021, 9, 488. https://doi.org/10.3390/pr9030488

AMA Style

Lohmann LJ, Strube J. Process Analytical Technology for Precipitation Process Integration into Biologics Manufacturing towards Autonomous Operation—mAb Case Study. Processes. 2021; 9(3):488. https://doi.org/10.3390/pr9030488

Chicago/Turabian Style

Lohmann, Lara Julia, and Jochen Strube. 2021. "Process Analytical Technology for Precipitation Process Integration into Biologics Manufacturing towards Autonomous Operation—mAb Case Study" Processes 9, no. 3: 488. https://doi.org/10.3390/pr9030488

APA Style

Lohmann, L. J., & Strube, J. (2021). Process Analytical Technology for Precipitation Process Integration into Biologics Manufacturing towards Autonomous Operation—mAb Case Study. Processes, 9(3), 488. https://doi.org/10.3390/pr9030488

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop