Investigating the interactions between physical and biological heterogeneities in bioreactors using compartment, population balance and metabolic models

A perspective-driven and technical evaluation of machine learning in bioreactor scale-up: A case-study for potential model developments

Bioreactor scale-up and scale-down have always been a topical issue for the biopharmaceutical industry and despite considerable effort, the identification of a fail-safe strategy for bioprocess development across scales remains a challenge. With the ubiquitous growth of digital transformation technologies, new scaling methods based on computer models may enable more effective scaling. This study aimed to evaluate the potential application of machine learning (ML) algorithms for bioreactor scale-up, with a specific focus on the prediction of scaling parameters. Factors critical to the development of such models were identified and data for bioreactor scale-up studies involving CHO cell-generated mAb products collated from the literature and public sources for the development of unsupervised and supervised ML models. Comparison of bioreactor performance across scales identified similarities between the different processes and primary differences between small- and large-scale bioreactors. A series of three case studies were developed to assess the relationship between cell growth and scale-sensitive bioreactor features. An embedding layer improved the capability of artificial neural network models to predict cell growth at a large-scale, as this approach captured similarities between the processes. Further models constructed to predict scaling parameters demonstrated how ML models may be applied to assist the scaling process. The development of data sets that include more characterization data with greater variability under different gassing and agitation regimes will also assist the future development of ML tools for bioreactor scaling.

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.

A comprehensive study of glucose and oxygen gradients in a scaled-down model of recombinant HuGM-CSF production in thermoinduced Escherichia coli fed-batch cultures

The effect of gradients of elevated glucose and low dissolved oxygen in the addition zone of fed-batch E. coli thermoinduced recombinant high cell density cultures can be evaluated through two-compartment scale-down models. Here, glucose was fed in the inlet of a plug flow bioreactor (PFB) connected to a stirred tank bioreactor (STB). E. coli cells diminished growth from 48.2 ± 2.2 g/L in the stage of RP production if compared to control (STB) with STB-PFB experiments, when residence time inside the PFB was 25 s (34.1 ± 3.5 g/L) and 40 s (25.6 ± 5.1 g/L), respectively. The recombinant granulocyte-macrophage colony-stimulating factor (rHuGM-CSF) production decreased from 34 ± 7% of RP in inclusion bodies (IB) in control cultures to 21 ± 8%, and 7 ± 4% during the thermoinduction production phase when increasing residence time inside the PFB to 25 s and 40 s, respectively. This, along with the accumulation of acetic and formic acid (up to 4 g/L), indicates metabolic redirection of central carbon routes through metabolic flow and mixed acid fermentation. Special care must be taken when producing a recombinant protein in heat-induced E. coli, because the yield and productivity of the protein decreases as the size of the bioreactors increases, especially if they are carried at high cell density.

Modeling large-scale bioreactors with diffusion equations. Part II: Characterizing substrate, oxygen, temperature, carbon dioxide, and pH profiles

Large-scale fermentation processes involve complex dynamic interactions between mixing, reaction, mass transfer, and the suspended biomass. Empirical correlations or case-specific computational simulations are usually used to predict and estimate the performance of large-scale bioreactors based on data acquired at bench scale. In this two-part-study, one-dimensional axial diffusion equations were studied as a general and predictive model of large-scale bioreactors. This second part focused on typical fed-batch operations where substrate gradients are known to occur, and characterized the profiles of substrate, pH, oxygen, carbon dioxide, and temperature. The physically grounded steady-state axial diffusion equations with first- and zeroth-order kinetics yielded analytical solutions to the relevant variables. The results were compared with large-scale Escherichia coli and Saccharomyces cerevisiae experiments and simulations from the literature, and good agreement was found in substrate profiles. The analytical profiles obtained for dissolved oxygen, temperature, pH, andCO 2 ${\text{CO}}_{2}$ were also consistent with the available data. Distribution functions for the substrate were defined, and efficiency factors for biomass growth and oxygen uptake rate were derived. In conclusion, this study demonstrated that axial diffusion equations can be used to model the effects of mixing and reaction on the relevant variables of typical large-scale fed-batch fermentations.

Cell-culture process optimization via model-based predictions of metabolism and protein glycosylation

In order to meet the rising demand for biologics and become competitive on the developing biosimilar market, there is a need for process intensification of biomanufacturing processes. Process development of biologics has historically relied on extensive experimentation to develop and optimize biopharmaceutical manufacturing. Experimentation to optimize media formulations, feeding schedules, bioreactor operations and bioreactor scale up is expensive, labor intensive and time consuming. Mathematical modeling frameworks have the potential to enable process intensification while reducing the experimental burden. This review focuses on mathematical modeling of cellular metabolism and N-linked glycosylation as applied to upstream manufacturing of biologics. We review developments in the field of modeling cellular metabolism of mammalian cells using kinetic and stoichiometric modeling frameworks along with their applications to simulate, optimize and improve mechanistic understanding of the process. Interest in modeling N-linked glycosylation has led to the creation of various types of parametric and non-parametric models. Most published studies on mammalian cell metabolism have performed experiments in shake flasks where the pH and dissolved oxygen cannot be controlled. Efforts to understand and model the effect of bioreactor-specific parameters such as pH, dissolved oxygen, temperature, and bioreactor heterogeneity are critically reviewed. Most modeling efforts have focused on the Chinese Hamster Ovary (CHO) cells, which are most commonly used to produce monoclonal antibodies (mAbs). However, these modeling approaches can be generalized and applied to any mammalian cell-based manufacturing platform. Current and potential future applications of these models for Vero cell-based vaccine manufacturing, CAR-T cell therapies, and viral vector manufacturing are also discussed. We offer specific recommendations for improving the applicability of these models to industrially relevant processes.

From Spatial-Temporal Multiscale Modeling to Application: Bridging the Valley of Death in Industrial Biotechnology

The Valley of Death confronts industrial biotechnology with a significant challenge to the commercialization of products. Fortunately, with the integration of computation, automation and artificial intelligence (AI) technology, the industrial biotechnology accelerates to cross the Valley of Death. The Fourth Industrial Revolution (Industry 4.0) has spurred advanced development of intelligent biomanufacturing, which has evolved the industrial structures in line with the worldwide trend. To achieve this, intelligent biomanufacturing can be structured into three main parts that comprise digitalization, modeling and intellectualization, with modeling forming a crucial link between the other two components. This paper provides an overview of mechanistic models, data-driven models and their applications in bioprocess development. We provide a detailed elaboration of the hybrid model and its applications in bioprocess engineering, including strain design, process control and optimization, as well as bioreactor scale-up. Finally, the challenges and opportunities of biomanufacturing towards Industry 4.0 are also discussed.

High-throughput genetic engineering of nonmodel and undomesticated bacteria via iterative site-specific genome integration

Efficient genome engineering is critical to understand and use microbial functions. Despite recent development of tools such as CRISPR-Cas gene editing, efficient integration of exogenous DNA with well-characterized functions remains limited to model bacteria. Here, we describe serine recombinase-assisted genome engineering, or SAGE, an easy-to-use, highly efficient, and extensible technology that enables selection marker-free, site-specific genome integration of up to 10 DNA constructs, often with efficiency on par with or superior to replicating plasmids. SAGE uses no replicating plasmids and thus lacks the host range limitations of other genome engineering technologies. We demonstrate the value of SAGE by characterizing genome integration efficiency in five bacteria that span multiple taxonomy groups and biotechnology applications and by identifying more than 95 heterologous promoters in each host with consistent transcription across environmental and genetic contexts. We anticipate that SAGE will rapidly expand the number of industrial and environmental bacteria compatible with high-throughput genetics and synthetic biology.

Diverse mechanisms of bioproduction heterogeneity in fermentation and their control strategies

Microbial bioproduction often faces challenges related to populational heterogeneity, where cells exhibit varying biosynthesis capabilities. Bioproduction heterogeneity can stem from genetic and non-genetic factors, resulting in decreased titer, yield, stability, and reproducibility. Consequently, understanding and controlling bioproduction heterogeneity are crucial for enhancing the economic competitiveness of large-scale biomanufacturing. In this review, we provide a comprehensive overview of current understandings of the various mechanisms underlying bioproduction heterogeneity. Additionally, we examine common strategies for controlling bioproduction heterogeneity based on these mechanisms. By implementing more robust measures to mitigate heterogeneity, we anticipate substantial enhancements in the scalability and stability of bioproduction processes.

ONE-SENTENCE SUMMARY

This review summarizes current understandings of different mechanisms of bioproduction heterogeneity and common control strategies based on these mechanisms.

Microbial lifelines in bioprocesses: From concept to application

Bioprocesses are scaled up for the production of large product quantities. With larger fermenter volumes, mixing becomes increasingly inefficient and environmental gradients get more prominent than in smaller scales. Environmental gradients have an impact on the microorganism's metabolism, which makes the prediction of large-scale performance difficult and can lead to scale-up failure. A promising approach for improved understanding and estimation of dynamics of microbial populations in large-scale bioprocesses is the analysis of microbial lifelines. The lifeline of a microbe in a bioprocess is the experience of environmental gradients from a cell's perspective, which can be described as a time series of position, environment and intracellular condition. Currently, lifelines are predominantly determined using models with computational fluid dynamics, but new technical developments in flow-following sensor particles and microfluidic single-cell cultivation open the door to a more interdisciplinary concept. We critically review the current concepts and challenges in lifeline determination and application of lifeline analysis, as well as strategies for the integration of these techniques into bioprocess development. Lifelines can contribute to a successful scale-up by guiding scale-down experiments and identifying strain engineering targets or bioreactor optimisations.

Substantial gradient mitigation in simulated large-scale bioreactors by optimally placed multiple feed points

The performance of large-scale stirred tank and bubble column bioreactors is often hindered by insufficient macromixing of feeds, leading to heterogeneities in pH, substrate, and oxygen, which complicates process scale-up. Appropriate feed placement or the use of multiple feed points could improve mixing. Here, theoretically optimal placement of feed points was derived using one-dimensional diffusion equations. The utility of optimal multipoint feeds was evaluated with mixing, pH control, and bioreaction simulations using three-dimensional compartment models of four industrially relevant bioreactors with working volumes ranging from 8 to 237 m3 . Dividing the vessel axially in equal-sized compartments and locating a feed point or multiple feed points symmetrically in each compartment reduced the mixing time substantially by more than a minute and mitigated gradients of pH, substrate, and oxygen. Performance of the large-scale bioreactors was consequently restored to ideal, homogeneous reactor performance: oxygen consumption and biomass yield were recovered and the phenotypical heterogeneity of the biomass population was diminished.

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.

Development of a miniature bioreactor model to study the impact of pH and DOT fluctuations on CHO cell culture performance as a tool to understanding heterogeneity effects at large-scale

Understanding the impact of spatial heterogeneities that are known to occur in large‐scale cell culture bioreactors remains a significant challenge. This work presents a novel methodology for mimicking the effects of pH and dissolved oxygen heterogeneities on Chinese hamster ovary (CHO) cell culture performance and antibody quality characteristics, using an automated miniature bioreactor system. Cultures of 4 different cell lines, expressing 3 IgG molecules and one fusion protein, were exposed to repeated pH and dissolved oxygen tension (DOT) fluctuations between pH 7.0–7.5 and DOT 10%–30%, respectively, for durations of 15, 30, and 60 min. Fluctuations in pH had a minimal impact on growth, productivity, and product quality although some changes in lactate metabolism were observed. DOT fluctuations were found to have a more significant impact; a 35% decrease in cell growth and product titre was observed in the fastest growing cell line tested, while all cell lines exhibited a significant increase in lactate accumulation. Product quality analysis yielded varied results; two cell lines showed an increase in the G0F glycan and decrease in G1F, G2F, and Man5; however, another line showed the opposite trend. The study suggests that the response of CHO cells to the effects of fluctuating culture conditions is cell line specific and that higher growing cell lines are most impacted. The miniature bioreactor system described in this work therefore provides a platform for use during early stage cell culture process development to identify cell lines that may be adversely impacted by the pH and DOT heterogeneities encountered on scale‐up. This experimental data can be combined with computational modeling approaches to predict overall cell culture performance in large‐scale bioreactors.

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.

Non-Idealities in Lab-Scale Kinetic Testing: A Theoretical Study of a Modular Temkin Reactor

The Temkin reactor can be applied for industrial relevant catalyst testing with unmodified catalyst particles. It was assumed in the literature that this reactor behaves as a cascade of continuously stirred tank reactors (CSTR). However, this assumption was based only on outlet gas composition or inert residence time distribution measurements. The present work theoretically investigates the catalytic CO2 methanation as a test case on different catalyst geometries, a sphere, and a ring, inside a single Temkin reaction chamber under isothermal conditions. Axial gas-phase species profiles from detailed computational fluid dynamics (CFD) are compared with a CSTR and 1D plug-flow reactor (PFR) model using a sophisticated microkinetic model. In addition, a 1D chemical reactor network (CRN) model was developed, and model parameters were adjusted based on the CFD simulations. Whereas the ideal reactor models overpredict the axial product concentrations, the CRN model results agree well with the CFD simulations, especially under low to medium flow rates. This study shows that complex flow patterns greatly influence species fields inside the Temkin reactor. Although residence time measurements suggest CSTR-like behavior, the reactive flow cannot be described by either a CSTR or PFR model but with the developed CRN model.

Design and construction of novel biocatalyst for bioprocessing: Recent advances and future outlook

Bioprocess, a biocatalysis-based technology, is becoming popular in many research fields and widely applied in industrial manufacturing. However, low bioconversion, low productivity, and high costs during industrial processes are usually the limitation in bioprocess. Therefore, many biocatalyst strategies have been developed to meet these challenges in recent years. In this review, we firstly discuss protein engineering strategies, which are emerged for improving the biocatalysis activity of biocatalysts. Then, we summarize metabolic engineering strategies that are promoting the development of microbial cell factories. Next, we illustrate the necessity of using the combining strategy of protein engineering and metabolic engineering for efficient biocatalysts. Lastly, future perspectives about the development and application of novel biocatalyst strategies are discussed. This review provides theoretical guidance for the development of efficient, sustainable, and economical bioprocesses mediated by novel biocatalysts.

The closure issue related to liquid-cell mass transfer and substrate uptake dynamics in biological systems

An original dynamic model for substrate uptake under transient conditions is established and used to simulate a variety of biological responses to external perturbations. The actual uptake and growth rates, treated as cell properties, are part of the model variables as well as the substrate concentration at the cell-liquid interface. Several regulatory loops inspired by the structure of the glycolytic chain are considered to establish a set of ordinary differential equations. The uptake rate evolves so as to reach an equilibrium between the cell demand and the environmental supply. This model does not contain any of the usual algebraic closure laws relating to the instantaneous uptake, growth rates, and the substrate concentration, nor does it enforce the continuity of mass fluxes at the liquid-cell interface. However, these relationships are found in the steady-state solution. Previously unexplained experimental observations are well reproduced by this model. Also, the model structure is suitable for further coupling with flux-based metabolic models and fluid-flow equations.

Euler-Lagrangian Simulations: A Proper Tool for Predicting Cellular Performance in Industrial Scale Bioreactors

Eulerian-Lagrangian approach to investigate cellular responses in a bioreactor has become the center of attention in recent years. It was introduced to biotechnological processes about two decades ago, but within the last few years, it proved itself as a powerful tool to address scale-up and -down topics of bioprocesses. It can capture the history of a cell and reveal invaluable information for, not only, bioprocess control and design but also strain engineering. This way it will be possible to shed light on the actual environment that cell experiences throughout its lifespan. Lifelines of a microorganism in a bioreactor can serve as the missing link that encompasses the biological timescales and the physical timescales. For this purpose digitalization of bioreactors provides us with new insights that are not achievable in industrial reactors easily if at all, namely, substrate and product gradients; high-shear regions are among the most interesting factors that can be reproduced adequately with help of a digital twin. In this chapter basic principles of this method will be introduced, and later on some practical aspects of particle tracking technique will be illustrated. In the final section, some of the advantages and challenges associated with this method will be discussed.

Overexpression of uracil permease and nucleoside transporter from Bacillus amyloliquefaciens improves cytidine production in Escherichia coli

Cytidine is an important raw material for nucleic acid health food and genetic engineering research. In recent years, it has shown irreplaceable effects in anti-virus, anti-tumor, and AIDS drugs. Its biosynthetic pathway is complex and highly regulated. In this study, overexpression of uracil permease and a nucleoside transporter from Bacillus amyloliquefaciens related to cell membrane transport in Escherichia coli strain BG-08 was found to increase cytidine production in shake flask cultivation by 1.3-fold (0.91 ± 0.03 g/L) and 1.8-fold (1.26 ± 0.03 g/L) relative to that of the original strain (0.70 ± 0.03 g/L), respectively. Co-overexpression of uracil permease and a nucleoside transporter further increased cytidine yield by 2.7-fold (1.59 ± 0.05 g/L) compared with that of the original strain. These results indicate that the overexpressed uracil permease and nucleoside transporter can promote the accumulation of cytidine, and the two proteins play a synergistic role in the secretion of cytidine in Escherichia coli.

Dynamic control in metabolic engineering: Theories, tools, and applications

Metabolic engineering has allowed the production of a diverse number of valuable chemicals using microbial organisms. Many biological challenges for improving bio-production exist which limit performance and slow the commercialization of metabolically engineered systems. Dynamic metabolic engineering is a rapidly developing field that seeks to address these challenges through the design of genetically encoded metabolic control systems which allow cells to autonomously adjust their flux in response to their external and internal metabolic state. This review first discusses theoretical works which provide mechanistic insights and design choices for dynamic control systems including two-stage, continuous, and population behavior control strategies. Next, we summarize molecular mechanisms for various sensors and actuators which enable dynamic metabolic control in microbial systems. Finally, important applications of dynamic control to the production of several metabolite products are highlighted, including fatty acids, aromatics, and terpene compounds. Altogether, this review provides a comprehensive overview of the progress, advances, and prospects in the design of dynamic control systems for improved titer, rate, and yield metrics in metabolic engineering.

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.

Predicting By-Product Gradients of Baker’s Yeast Production at Industrial Scale: A Practical Simulation Approach

Scaling up bioprocesses is one of the most crucial steps in the commercialization of bioproducts. While it is known that concentration and shear rate gradients occur at larger scales, it is often too risky, if feasible at all, to conduct validation experiments at such scales. Using computational fluid dynamics equipped with mechanistic biochemical engineering knowledge of the process, it is possible to simulate such gradients. In this work, concentration profiles for the by-products of baker’s yeast production are investigated. By applying a mechanistic black-box model, concentration heterogeneities for oxygen, glucose, ethanol, and carbon dioxide are evaluated. The results suggest that, although at low concentrations, ethanol is consumed in more than 90% of the tank volume, which prevents cell starvation, even when glucose is virtually depleted. Moreover, long exposure to high dissolved carbon dioxide levels is predicted. Two biomass concentrations, i.e., 10 and 25 g/L, are considered where, in the former, ethanol production is solely because of overflow metabolism while, in the latter, 10% of the ethanol formation is due to dissolved oxygen limitation. This method facilitates the prediction of the living conditions of the microorganism and its utilization to address the limitations via change of strain or bioreactor design or operation conditions. The outcome can also be of value to design a representative scale-down reactor to facilitate strain studies.

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.

Synergistically applying 1-D modeling and CFD for designing industrial scale bubble column syngas bioreactors

The reduction of greenhouse gas emissions and future perspectives of circular economy ask for new solutions to produce commodities and fine chemicals. Large-scale bubble columns operated by gaseous substrates such as CO, CO2, and H2 to feed acetogens for product formations could be promising approaches. Valid in silico predictions of large-scale performance are needed to dimension bioreactors properly taking into account biological constraints, too. This contribution deals with the trade-off between sophisticated spatiotemporally resolved large-scale simulations using computationally intensive Euler-Euler and Euler-Lagrange approaches and coarse-grained 1-D models enabling fast performance evaluations. It is shown that proper consideration of gas hold-up is key to predict biological performance. Intrinsic bias of 1-D models can be compensated by reconsideration of Sauter diameters derived from uniquely performed Euler-Lagrange computational fluid dynamics.

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.

A Novel Framework for Parameter and State Estimation of Multicellular Systems Using Gaussian Mixture Approximations

Multicellular systems play an important role in many biotechnological processes. Typically, these exhibit cell-to-cell variability, which has to be monitored closely for process control and optimization. However, some properties may not be measurable due to technical and financial restrictions. To improve the monitoring, model-based online estimators can be designed for their reconstruction. The multicellular dynamics is accounted for in the framework of population balance models (PBMs). These models are based on single cell kinetics, and each cellular state translates directly into an additional dimension of the obtained partial differential equations. As multicellular dynamics often require detailed single cell models and feature a high number of cellular components, the resulting population balance equations are often high-dimensional. Therefore, established state estimation concepts for PBMs based on discrete grids are not recommended due to the large computational effort. In this contribution a novel approach is proposed, which is based on the approximation of the underlying number density functions as the weighted sum of Gaussian distributions. Thus, the distribution is described by the characteristic properties of the individual Gaussians, like the mean and covariance. Thereby, the complex infinite dimensional estimation problem can be reduced to a finite dimension. The characteristic properties are estimated in a recursive approach. The method is evaluated for two academic benchmark examples, and the results indicate its potential for model-based online reconstruction for multicellular systems.

Heterogeneity in Pure Microbial Systems: Experimental Measurements and Modeling

Cellular heterogeneity influences bioprocess performance in ways that until date are not completely elucidated. In order to account for this phenomenon in the design and operation of bioprocesses, reliable analytical and mathematical descriptions are required. We present an overview of the single cell analysis, and the mathematical modeling frameworks that have potential to be used in bioprocess control and optimization, in particular for microbial processes. In order to be suitable for bioprocess monitoring, experimental methods need to be high throughput and to require relatively short processing time. One such method used successfully under dynamic conditions is flow cytometry. Population balance and individual based models are suitable modeling options, the latter one having in particular a good potential to integrate the various data collected through experimentation. This will be highly beneficial for appropriate process design and scale up as a more rigorous approach may prevent a priori unwanted performance losses. It will also help progressing synthetic biology applications to industrial scale.

Euler-Lagrange computational fluid dynamics for (bio)reactor scale down: An analysis of organism lifelines

The trajectories, referred to as lifelines, of individual microorganisms in an industrial scale fermentor under substrate limiting conditions were studied using an Euler‐Lagrange computational fluid dynamics approach. The metabolic response to substrate concentration variations along these lifelines provides deep insight in the dynamic environment inside a large‐scale fermentor, from the point of view of the microorganisms themselves. We present a novel methodology to evaluate this metabolic response, based on transitions between metabolic “regimes” that can provide a comprehensive statistical insight in the environmental fluctuations experienced by microorganisms inside an industrial bioreactor. These statistics provide the groundwork for the design of representative scale‐down simulators, mimicking substrate variations experimentally. To focus on the methodology we use an industrial fermentation of Penicillium chrysogenum in a simplified representation, dealing with only glucose gradients, single‐phase hydrodynamics, and assuming no limitation in oxygen supply, but reasonably capturing the relevant timescales. Nevertheless, the methodology provides useful insight in the relation between flow and component fluctuation timescales that are expected to hold in physically more thorough simulations. Microorganisms experience substrate fluctuations at timescales of seconds, in the order of magnitude of the global circulation time. Such rapid fluctuations should be replicated in truly industrially representative scale‐down simulators.