Integration Approaches to Model Bioreactor Hydrodynamics and Cellular Kinetics for Advancing Bioprocess Optimisation

Large-scale bioprocesses are increasing globally to cater to the larger market demands for biological products. As fermenter volumes increase, the efficiency of mixing decreases, and environmental gradients become more pronounced compared to smaller scales. Consequently, the cells experience gradients in process parameters, which in turn affects the efficiency and profitability of the process. Computational fluid dynamics (CFD) simulations are being widely embraced for their ability to simulate bioprocess performance, facilitate bioprocess upscaling, downsizing, and process optimisation. Recently, CFD approaches have been integrated with dynamic Cell reaction kinetic (CRK) modelling to generate valuable information about the cellular response to fluctuating hydrodynamic parameters inside large production processes. Such coupled approaches have the potential to facilitate informed decision-making in intelligent biomanufacturing, aligning with the principles of "Industry 4.0" concerning digitalisation and automation. In this review, we discuss the benefits of utilising integrated CFD-CRK models and the different approaches to integrating CFD-based bioreactor hydrodynamic models with cellular kinetic models. We also highlight the suitability of different coupling approaches for bioprocess modelling in the purview of associated computational loads.

An Automatic Method for Generation of CFD-Based 3D Compartment Models: Towards Real-Time Mixing Simulations

This article aims to develop a method to automatically generate CFD-based compartment models. This effort to simplify mixing models aims at capturing the interactions between material transport and chemical/biochemical conversions in large-scale reactors. The proposed method converts the CFD results into a system of mass balance equations for each defined component. The compartmentalization method is applied to two bioreactor geometries and was able to replicate tracer mixing profiles observed in CFD simulations. The generated compartment models were successfully coupled with, a simple Monod-type biokinetic model describing microbial growth, substrate consumption and product formation. The coupled model was used to simulate a four-hour fermentation in a 190 L reactor and a 10 m3 reactor. Resolving the substrate gradients had a clear impact on the biokinetics, increasing with the scale of the reactor. Moreover, the coupled model could simulate the fermentation faster than real-time. Having a real-time-solvable model is essential for implementations in digital twins and other real-time applications using the models as predictive tools.

Milligrams to kilograms: making microbes work at scale

If biomanufacturing can become a sustainable route for producing chemicals, it will provide a critical step in reducing greenhouse gas emissions to fight climate change. However, efforts to industrialize microbial synthesis of chemicals have met with varied success, due, in part, to challenges in translating laboratory successes to industrial scale. With a particular focus on Escherichia coli, this review examines the lessons learned when studying microbial physiology and metabolism under conditions that simulate large-scale bioreactors and methods to minimize cellular waste through reduction of maintenance energy, optimizing the stress response and minimizing culture heterogeneity. With general strategies to overcome these challenges, biomanufacturing process scale-up could be de-risked and the time and cost of bringing promising syntheses to market could be reduced.

Kinetic studies and CFD-based reaction modeling for insights into the scalability of ADC conjugation reactions

The manufacturing of antibody-drug conjugates (ADCs) involves the addition of a cytotoxic small-molecule linker-drug (= payload) to a solution of functionalized antibodies. For the development of robust conjugation processes, initially small-scale reaction tubes are used which requires a lot of manual handling. Scale-up to larger reaction vessels is often knowledge-driven and scale-comparability is solely assessed based on final product quality which does not account for the dynamics of the reaction. In addition, information about the influence of process parameters, such as stirrer speed, temperature, or payload addition rates, is limited due to high material costs. Given these limitations, there is a need for a modeling-based approach to investigate conjugation scale-up. In this work, both experimental kinetic studies and computational fluid dynamics (CFD) conjugation simulations were performed to understand the influence of scale and mixing parameters. In the experimental part, conjugation kinetics in small-scale reaction tubes with different mixing types were investigated for two ADC systems and compared to larger bench-scale reactions. It was demonstrated that more robust kinetics can be achieved through internal stirrer mixing instead of external mixing devices, such as orbital shakers. In the simulation part, 3D-reactor models were created by coupling CFD-models for three large-scale reaction vessels with a kinetic model for a site-specific conjugation reaction. This enabled to study the kinetics in different vessels, as well as the effect of process parameter variations in silico. Overall, it was found that for this conjugation type sufficient mixing can be achieved at all scales and the studied parameters cause only deviations during the payload addition period. An additional time-scale analysis demonstrated to aid the assessment of mixing effects during ADC process scale-up when mixing times and kinetic rates are known. In summary, this work highlights the benefit of kinetic models for enhanced conjugation process understanding without the need for large-scale experiments.

Microfluidic single-cell scale-down bioreactors: A proof-of-concept for the growth of Corynebacterium glutamicum at oscillating pH values

In large-scale bioreactors, gradients in cultivation parameter such as oxygen, substrate and pH result in fluctuating cell environments. pH fluctuations were identified as a critical parameter for bioprocess performance. Traditionally, scale-down systems at the laboratory scale are used to analyze the effects of fluctuating pH values on strain and thus process performance. Here, we demonstrate the application of dynamic microfluidic single-cell cultivation (dMSCC) as a novel scale-down system for the characterization of Corynebacterium glutamicum growth using oscillating pH conditions as a model stress-factor. A detailed comparison between two-compartment reactor (two-CR) scale-down experiments and dMSCC was performed for one specific pH oscillation between reference pH 7 (~ 8 min) and disturbed pH 6 (~2 min). Similar reductions in growth rates were observed in both systems (dMSCC 21% and two-CR 27%) compared to undisturbed cultivation at pH 7. Afterwards, systematic experiments at symmetric and asymmetric pH oscillations between pH ranges of 4-6 and 8-11 and different intervals from 1 minute to 20 minutes, were performed to demonstrate the unique application range and throughput of the dMSCC system. Finally, the strength of the dMSCC application was demonstrated by mimicking fluctuating environmental conditions of a putative large-scale bioprocesse, which is difficult to conduct using two-CRs. This article is protected by copyright. All rights reserved.

Digital models in biotechnology: Towards multi-scale integration and implementation

Industrial biotechnology encompasses a large area of multi-scale and multi-disciplinary research activities. With the recent megatrend of digitalization sweeping across all industries, there is an increased focus in the biotechnology industry on developing, integrating and applying digital models to improve all aspects of industrial biotechnology. Given the rapid development of this field, we systematically classify the state-of-art modelling concepts applied at different scales in industrial biotechnology and critically discuss their current usage, advantages and limitations. Further, we critically analyzed current strategies to couple cell models with computational fluid dynamics to study the performance of industrial microorganisms in large-scale bioprocesses, which is of crucial importance for the bio-based production industries. One of the most challenging aspects in this context is gathering intracellular data under industrially relevant conditions. Towards comprehensive models, we discuss how different scale-down concepts combined with appropriate analytical tools can capture intracellular states of single cells. We finally illustrated how the efforts could be used to develop digitals models suitable for both cell factory design and process optimization at industrial scales in the future.

Microfluidic Reproduction of Dynamic Bioreactor Environment Based on Computational Lifelines

The biotechnological production of fine chemicals, proteins and pharmaceuticals is usually hampered by loss of microbial performance during scale-up. This challenge is mainly caused by discrepancies between homogeneous environmental conditions at laboratory scale, where bioprocesses are optimized, and inhomogeneous conditions in large-scale bioreactors, where production takes place. Therefore, to improve strain selection and process development, it is of great interest to characterize these fluctuating conditions at large-scale and to study their effects on microbial cells. In this paper, we demonstrate the potential of computational fluid dynamics (CFD) simulation of large-scale bioreactors combined with dynamic microfluidic single-cell cultivation (dMSCC). Environmental conditions in a 200 L bioreactor were characterized with CFD simulations. Computational lifelines were determined by combining simulated turbulent multiphase flow, mass transport and particle tracing. Glucose availability for Corynebacterium glutamicum cells was determined. The reactor was simulated with average glucose concentrations of 6 g m−3, 10 g m−3 and 16 g m−3. The resulting computational lifelines, discretized into starvation and abundance regimes, were used as feed profiles for the dMSCC to investigate how varying glucose concentration affects cell physiology and growth rate. In this study, each colony in the dMSCC device represents a single cell as it travels through the reactor. Under oscillating conditions reproduced in the dMSCC device, a decrease in growth rate of about 40% was observed compared to continuous supply with the same average glucose availability. The presented approach provides insights into environmental conditions observed by microorganisms in large-scale bioreactors. It also paves the way for an improved understanding of how inhomogeneous environmental conditions influence cellular physiology, growth and production.

Stochastic parcel tracking in an Euler-Lagrange compartment model for fast simulation of fermentation processes

The compartment model (CM) is a well‐known approach for computationally affordable, spatially resolved hydrodynamic modeling of unit operations. Recent implementations use flow profiles based on Computational Fluid Dynamics (CFD) simulations, and several authors included microbial kinetics to simulate gradients in bioreactors. However, these studies relied on black‐box kinetics that do not account for intracellular changes and cell population dynamics in response to heterogeneous environments. In this paper, we report the implementation of a Lagrangian reaction model, where the microbial phase is tracked as a set of biomass‐parcels, each linked with an intracellular composition vector and a structured reaction model describing their intracellular response to extracellular variations. A stochastic parcel tracking approach is adopted, in contrast to the resolved trajectories used in CFD implementations. A penicillin production process is used as a case study. We show good performance of the model compared with full CFD simulations, both regarding the extracellular gradients and intracellular pool response, using the mixing time as a matching criterion and taking into account that the mixing time is sensitive to the number of compartments. The sensitivity of the model output towards some of the inputs is explored. The coarsest representative CM requires a few minutes to solve 80 h of flow time, compared with approximately 2 weeks for a full Euler–Lagrange CFD simulation of the same case. This alleviates one of the major bottlenecks for the application of such CFD simulations towards the analysis and optimization of industrial fermentation processes.

Microbial technologies for biotherapeutics production: Key tools for advanced biopharmaceutical process development and control

Current trends in the biopharmaceutical market such as the diversification of therapies as well as the increasing time-to-market pressure will trigger the rethinking of bioprocess development and production approaches. Thereby, the importance of development time and manufacturing costs will increase, especially for microbial production. In the present review, we investigate three technological approaches which, to our opinion, will play a key role in the future of biopharmaceutical production. The first cornerstone of process development is the generation and effective utilization of platform knowledge. Building processes on well understood microbial and technological platforms allows to accelerate early-stage bioprocess development and to better condense this knowledge into multi-purpose technologies and applicable mathematical models. Second, the application of verified scale down systems and in silico models for process design and characterization will reduce the required number of large scale batches before dossier submission. Third, the broader availability of mathematical process models and the improvement of process analytical technologies will increase the applicability and acceptance of advanced control and process automation in the manufacturing scale. This will reduce process failure rates and subsequently cost of goods. Along these three aspects we give an overview of recently developed key tools and their potential integration into bioprocess development strategies.

Potential of Integrating Model-Based Design of Experiments Approaches and Process Analytical Technologies for Bioprocess Scale-Down

Typically, bioprocesses on an industrial scale are dynamic systems with a certain degree of variability, system inhomogeneities, and even population heterogeneities. Therefore, the scaling of such processes from laboratory to industrial scale and vice versa is not a trivial task. Traditional scale-down methodologies consider several technical parameters, so that systems on the laboratory scale tend to qualitatively reflect large-scale effects, but not the dynamic situation in an industrial bioreactor over the entire process, from the perspective of a cell. Supported by the enormous increase in computing power, the latest scientific focus is on the application of dynamic models, in combination with computational fluid dynamics to quantitatively describe cell behavior. These models allow the description of possible cellular lifelines which in turn can be used to derive a regime analysis for scale-down experiments. However, the approaches described so far, which were for a very few process examples, are very labor- and time-intensive and cannot be validated easily. In parallel, alternatives have been developed based on the description of the industrial process with hybrid process models, which describe a process mechanistically as far as possible in order to determine the essential process parameters with their respective variances. On-line analytical methods allow the characterization of population heterogeneity directly in the process. This detailed information from the industrial process can be used in laboratory screening systems to select relevant conditions in which the cell and process related parameters reflect the situation in the industrial scale. In our opinion, these technologies, which are available in research for modeling biological systems, in combination with process analytical techniques are so far developed that they can be implemented in industrial routines for faster development of new processes and optimization of existing ones.

Understanding gradients in industrial bioreactors

Gradients in industrial bioreactors have attracted substantial research attention since exposure to fluctuating environmental conditions has been shown to lead to changes in the metabolome, transcriptome as well as population heterogeneity in industrially relevant microorganisms. Such changes have also been found to impact key process parameters like the yield on substrate and the productivity. Hence, understanding gradients is important from both the academic and industrial perspectives. In this review the causes of gradients are outlined, along with their impact on microbial physiology. Quantifying the impact of gradients requires a detailed understanding of both fluid flow inside industrial equipment and microbial physiology. This review critically examines approaches used to investigate gradients including large-scale experimental work, computational methods and scale-down approaches. Avenues for future work have been highlighted, particularly the need for further coordinated development of both in silico and experimental tools which can be used to further the current understanding of gradients in industrial equipment.

A CFD model for predicting protein aggregation in low-pH virial inactivation for mAb production

Significant amounts of soluble product aggregates were observed in the low-pH viral inactivation (VI) operation during an initial scale-up run for an immunoglobulin-G 4 (IgG4) monoclonal antibody (mAb IgG4-N1). Being earlier in development, a scale-down model did not exist, nor was it practical to use costly Protein A eluate (PAE) for testing the VI process at scale, thus, a computational fluid dynamics (CFD)-based high-molecular weight (HMW) prediction model was developed for troubleshooting and risk mitigation. It was previously reported that the IgG4-N1 molecules upon exposure to low pH tend to change into transient and partially unfolded monomers during VI acidification (i.e., VIA) and form aggregates after neutralization (i.e., VIN). Therefore, the CFD model reported here focuses on the VIA step. The model mimics the continuous addition of acid to PAE and tracks acid distribution during VIA. Based on the simulated low-pH zone (≤pH 3.3) profiles and PAE properties, the integrated low-pH zone (ILPZ) value was obtained to predict HMW level at the VI step. The simulations were performed to examine the operating parameters, such as agitation speed, acid addition rate, and protein concentration of PAE, of the pilot scale (50-200 L) runs. The conditions with predictions of no product aggregation risk were recommended to the real scale-up runs, resulted in 100% success rate of the consecutive 12 pilot-scale runs. This study demonstrated that the CFD-based HMW prediction model could be used as a tool to facilitate the scale up of the low-pH VI process directly from bench to pilot/production scale.

Scale-down simulators for mammalian cell culture as tools to access the impact of inhomogeneities occurring in large-scale bioreactors

During the scale-up of a bioprocess, not all characteristics of the process can be kept constant throughout the different scales. This typically results in increased mixing times with increasing reactor volumes. The poor mixing leads in turn to the formation of concentration gradients throughout the reactor and exposes cells to varying external conditions based on their location in the bioreactor. This can affect process performance and complicate process scale-up. Scale-down simulators, which aim at replicating the large-scale environment, expose the cells to changing environmental conditions. This has the potential to reveal adaptation mechanisms, which cells are using to adjust to rapidly fluctuating environmental conditions and can identify possible root causes for difficulties maintaining similar process performance at different scales. This understanding is of utmost importance in process validation. Additionally, these simulators also have the potential to be used for selecting cells, which are most robust when encountering changing extracellular conditions. The aim of this review is to summarize recent work in this interesting and promising area with the focus on mammalian bioprocesses, since microbial processes have been extensively reviewed.

Bacterial Flow Cytometry and Imaging as Potential Process Monitoring Tools for Industrial Biotechnology

Minimizing process variations by early identification of deviations is one approach to make industrial production processes robust. Cell morphology is a direct representation of the physiological state and an important factor for the cell’s survival in harsh environments as encountered during industrial processing. The adverse effects of fluctuating process parameters on cells were studied using flow cytometry and imaging. Results showed that altered pH caused a shift in cell size distribution from a heterogeneous mix of elongated and short cells to a homogenous population of short cells. Staining based on membrane integrity revealed a dynamics in the pattern of cluster formation during fermentation. Contradictory findings from forward scatter and imaging highlight the need for use of complementary techniques that provide visual confirmation to interpret changes. An atline flow cytometry or imaging capable of identifying subtle population deviations serves as a powerful monitoring tool for industrial biotechnology.

Coupled metabolic-hydrodynamic modeling enabling rational scale-up of industrial bioprocesses

Metabolomics aims to address what and how regulatory mechanisms are coordinated to achieve flux optimality, different metabolic objectives as well as appropriate adaptations to dynamic nutrient availability. Recent decades have witnessed that the integration of metabolomics and fluxomics within the goal of synthetic biology has arrived at generating the desired bioproducts with improved bioconversion efficiency. Absolute metabolite quantification by isotope dilution mass spectrometry represents a functional readout of cellular biochemistry and contributes to the establishment of metabolic (structured) models required in systems metabolic engineering. In industrial practices, population heterogeneity arising from fluctuating nutrient availability frequently leads to performance losses, that is reduced commercial metrics (titer, rate, and yield). Hence, the development of more stable producers and more predictable bioprocesses can benefit from a quantitative understanding of spatial and temporal cell-to-cell heterogeneity within industrial bioprocesses. Quantitative metabolomics analysis and metabolic modeling applied in computational fluid dynamics (CFD)-assisted scale-down simulators that mimic industrial heterogeneity such as fluctuations in nutrients, dissolved gases, and other stresses can procure informative clues for coping with issues during bioprocessing scale-up. In previous studies, only limited insights into the hydrodynamic conditions inside the industrial-scale bioreactor have been obtained, which makes case-by-case scale-up far from straightforward. Tracking the flow paths of cells circulating in large-scale bioreactors is a highly valuable tool for evaluating cellular performance in production tanks. The "lifelines" or "trajectories" of cells in industrial-scale bioreactors can be captured using Euler-Lagrange CFD simulation. This novel methodology can be further coupled with metabolic (structured) models to provide not only a statistical analysis of cell lifelines triggered by the environmental fluctuations but also a global assessment of the metabolic response to heterogeneity inside an industrial bioreactor. For the future, the industrial design should be dependent on the computational framework, and this integration work will allow bioprocess scale-up to the industrial scale with an end in mind.