The shortcomings of accurate rate estimations in cultivation processes and a solution for precise and robust process modeling

Analyzing and understanding the robustness of bioprocesses

The robustness of bioprocesses is becoming increasingly important. The main driving forces of this development are, in particular, increasing demands on product purities as well as economic aspects. In general, bioprocesses exhibit extremely high complexity and variability. Biological systems often have a much higher intrinsic variability compared with chemical processes, which makes the development and characterization of robust processes tedious task. To predict and control robustness, a clear understanding of interactions between input and output variables is necessary. Robust bioprocesses can be realized, for example, by using advanced control strategies for the different unit operations. In this review, we discuss the different biological, technical, and mathematical tools for the analysis and control of bioprocess robustness.

Comparison of mechanistic and hybrid modeling approaches for characterization of a CHO cultivation process: Requirements, pitfalls and solution paths

Despite the advantages of mathematical bioprocess modeling, successful model implementation already starts with experimental planning and accordingly can fail at this early stage. For this study, two different modeling approaches (mechanistic and hybrid) based on a four-dimensional antibody-producing CHO fed-batch process are compared. Overall, 33 experiments are performed in the fractional factorial four-dimensional design space and separated into four different complex data partitions subsequently used for model comparison and evaluation. The mechanistic model demonstrates the advantage of prior knowledge (i.e., known equations) to get informative value relatively independently of the utilized data partition. The hybrid approach displayes a higher data dependency but simultaneously yielded a higher accuracy on all data partitions. Furthermore, our results demonstrate that independent of the chosen modeling framework, a smart selection of only four initial experiments can already yield a very good representation of a full design space independent of the chosen modeling structure. Academic and industry researchers are recommended to pay more attention to experimental planning to maximize the process understanding obtained from mathematical modeling.

Online 2D Fluorescence Monitoring in Microtiter Plates Allows Prediction of Cultivation Parameters and Considerable Reduction in Sampling Efforts for Parallel Cultivations of Hansenula polymorpha

Multi-wavelength (2D) fluorescence spectroscopy represents an important step towards exploiting the monitoring potential of microtiter plates (MTPs) during early-stage bioprocess development. In combination with multivariate data analysis (MVDA), important process information can be obtained, while repetitive, cost-intensive sample analytics can be reduced. This study provides a comprehensive experimental dataset of online and offline measurements for batch cultures of Hansenula polymorpha. In the first step, principal component analysis (PCA) was used to assess spectral data quality. Secondly, partial least-squares (PLS) regression models were generated, based on spectral data of two cultivation conditions and offline samples for glycerol, cell dry weight, and pH value. Thereby, the time-wise resolution increased 12-fold compared to the offline sampling interval of 6 h. The PLS models were validated using offline samples of a shorter sampling interval. Very good model transferability was shown during the PLS model application to the spectral data of cultures with six varying initial cultivation conditions. For all the predicted variables, a relative root-mean-square error (RMSE) below 6% was obtained. Based on the findings, the initial experimental strategy was re-evaluated and a more practical approach with minimised sampling effort and elevated experimental throughput was proposed. In conclusion, the study underlines the high potential of multi-wavelength (2D) fluorescence spectroscopy and provides an evaluation workflow for PLS modelling in microtiter plates.

Model Transferability and Reduced Experimental Burden in Cell Culture Process Development Facilitated by Hybrid Modeling and Intensified Design of Experiments

Reliable process development is accompanied by intense experimental effort. The utilization of an intensified design of experiments (iDoE) (intra-experimental critical process parameter (CPP) shifts combined) with hybrid modeling potentially reduces process development burden. The iDoE can provide more process response information in less overall process time, whereas hybrid modeling serves as a commodity to describe this behavior the best way. Therefore, a combination of both approaches appears beneficial for faster design screening and is especially of interest at larger scales where the costs per experiment rise significantly. Ideally, profound process knowledge is gathered at a small scale and only complemented with few validation experiments on a larger scale, saving valuable resources. In this work, the transferability of hybrid modeling for Chinese hamster ovary cell bioprocess development along process scales was investigated. A two-dimensional DoE was fully characterized in shake flask duplicates (300 ml), containing three different levels for the cultivation temperature and the glucose concentration in the feed. Based on these data, a hybrid model was developed, and its performance was assessed by estimating the viable cell concentration and product titer in 15 L bioprocesses with the same DoE settings. To challenge the modeling approach, 15 L bioprocesses also comprised iDoE runs with intra-experimental CPP shifts, impacting specific cell rates such as growth, consumption, and formation. Subsequently, the applicability of the iDoE cultivations to estimate static cultivations was also investigated. The shaker-scale hybrid model proved suitable for application to a 15 L scale (1:50), estimating the viable cell concentration and the product titer with an NRMSE of 10.92% and 17.79%, respectively. Additionally, the iDoE hybrid model performed comparably, displaying NRMSE values of 13.75% and 21.13%. The low errors when transferring the models from shaker to reactor and between the DoE and the iDoE approach highlight the suitability of hybrid modeling for mammalian cell culture bioprocess development and the potential of iDoE to accelerate process characterization and to improve process understanding.

Monte Carlo-Based Error Propagation for a More Reliable Regression Analysis across Specific Rates in Bioprocesses

During process development, bioprocess data need to be converted into applicable knowledge. Therefore, it is crucial to evaluate the obtained data under the usage of transparent and reliable data reduction and correlation techniques. Within this contribution, we show a generic Monte Carlo error propagation and regression approach applied to two different, industrially relevant cultivation processes. Based on measurement uncertainties, errors for cell-specific growth, uptake, and production rates were determined across an evaluation chain, with interlinked inputs and outputs. These uncertainties were subsequently included in regression analysis to derive the covariance of the regression coefficients and the confidence bounds for prediction. The usefulness of the approach is shown within two case studies, based on the relations across biomass-specific rate control limits to guarantee high productivities in E. coli, and low lactate formation in a CHO cell fed-batch could be established. Besides the possibility to determine realistic errors on the evaluated process data, the presented approach helps to differentiate between reliable and unreliable correlations and prevents the wrong interpretations of relations based on uncertain data.

Generic and specific recurrent neural network models: Applications for large and small scale biopharmaceutical upstream processes

The calculation of temporally varying upstream process outcomes is a challenging task. Over the last years, several parametric, semi-parametric as well as non-parametric approaches were developed to provide reliable estimates for key process parameters. We present generic and product-specific recurrent neural network (RNN) models for the computation and study of growth and metabolite-related upstream process parameters as well as their temporal evolution. Our approach can be used for the control and study of single product-specific large-scale manufacturing runs as well as generic small-scale evaluations for combined processes and products at development stage. The computational results for the product titer as well as various major upstream outcomes in addition to relevant process parameters show a high degree of accuracy when compared to experimental data and, accordingly, a reasonable predictive capability of the RNN models. The calculated values for the root-mean squared errors of prediction are significantly smaller than the experimental standard deviation for the considered process run ensembles, which highlights the broad applicability of our approach. As a specific benefit for platform processes, the generic RNN model is also used to simulate process outcomes for different temperatures in good agreement with experimental results. The high level of accuracy and the straightforward usage of the approach without sophisticated parameterization and recalibration procedures highlight the benefits of the RNN models, which can be regarded as promising alternatives to existing parametric and semi-parametric methods.

Novel Strategy for the Calorimetry-Based Control of Fed-Batch Cultivations of Saccharomyces cerevisiae

Typical controllers for fed-batch cultivations are based on the estimation and control of the specific growth rate in real time. Biocalorimetry allows one to measure a heat signal proportional to the substrate consumed by cells. The derivative of this heat signal is usually used to evaluate the specific growth rate, introducing noise to the resulting estimate. To avoid this, this study investigated a novel controller based directly on the heat signal. Time trajectories of the heat signal setpoint were modelled for different specific growth rates, and the controller was set to follow this dynamic setpoint. The developed controller successfully followed the setpoint during aerobic cultivations of Saccharomyces cerevisiae, preventing the Crabtree effect by maintaining low glucose concentrations. With this new method, fed-batch cultivations of S. cerevisiae could be reliably controlled at specific growth rates between 0.075 h−1 and 0.20 h−1, with average root mean square errors of 15 ± 3%.

Error propagation in constraint-based modeling of Chinese hamster ovary cells

Chinese hamster ovary (CHO) cells are the most popular mammalian cell factories for the production of glycosylated biopharmaceuticals. To further increase titer and productivity and ensure product quality, rational system-level engineering strategies based on constraint-based metabolic modeling, such as flux balance analysis (FBA), have gained strong interest. However, the quality of FBA predictions depends on the accuracy of the experimental input data, especially on the exchange rates of extracellular metabolites. Yet, it is not standard practice to devote sufficient attention to the accurate determination of these rates. In this work, we investigated to what degree the sampling frequency during a batch culture and the measurement errors of metabolite concentrations influence the accuracy of the calculated exchange rates and further, how this error then propagates into FBA predictions of growth rates. We determined that accurate measurements of essential amino acids with low uptake rates are crucial for the accuracy of FBA predictions, followed by a sufficient number of analyzed time points. We observed that the measured difference in growth rates of two cell lines can only be reliably predicted when both high measurement accuracy and sampling frequency are ensured.

Time Integrated Flux Analysis: Exploiting the Concentration Measurements Directly for Cost-Effective Metabolic Network Flux Analysis

Background: Flux analyses, such as Metabolic Flux Analysis (MFA), Flux Balance Analysis (FBA), Flux Variability Analysis (FVA) or similar methods, can provide insights into the cellular metabolism, especially in combination with experimental data. The most common integration of extracellular concentration data requires the estimation of the specific fluxes (/rates) from the measured concentrations. This is a time-consuming, mathematically ill-conditioned inverse problem, raising high requirements for the quality and quantity of data. Method: In this contribution, a time integrated flux analysis approach is proposed which avoids the error-prone estimation of specific flux values. The approach is adopted for a Metabolic time integrated Flux Analysis and (sparse) time integrated Flux Balance/Variability Analysis. The proposed approach is applied to three case studies: (1) a simulated bioprocess case studying the impact of the number of samples (experimental points) and measurements’ noise on the performance; (2) a simulation case to understand the impact of network redundancies and reaction irreversibility; and (3) an experimental bioprocess case study, showing its relevance for practical applications. Results: It is observed that this method can successfully estimate the time integrated flux values, even with relatively low numbers of samples and significant noise levels. In addition, the method allows the integration of additional constraints (e.g., bounds on the estimated concentrations) and since it eliminates the need for estimating fluxes from measured concentrations, it significantly reduces the workload while providing about the same level of insight into the metabolism as classic flux analysis methods.