Model-Based Optimization of a Fed-Batch Bioreactor for mAb Production Using a Hybridoma Cell Culture

Maria, Gheorghe

doi:10.3390/molecules25235648

Open AccessArticle

Model-Based Optimization of a Fed-Batch Bioreactor for mAb Production Using a Hybridoma Cell Culture

by

Gheorghe Maria

^1,2

¹

Department of Chemical and Biochemical Engineering, University Politehnica of Bucharest, Polizu Str. 1-7, P.O. 35-107, 011061 Bucharest, Romania

²

Romanian Academy, Calea Victoriei, 125, 010071 Bucharest, Romania

Molecules 2020, 25(23), 5648; https://doi.org/10.3390/molecules25235648

Submission received: 28 September 2020 / Revised: 7 November 2020 / Accepted: 26 November 2020 / Published: 30 November 2020

(This article belongs to the Special Issue New Insight of Bio-Mimetic Emulations)

Download

Browse Figures

Review Reports Versions Notes

Abstract

:

Production of monoclonal antibodies (mAbs) is a well-known method used to synthesize a large number of identical antibodies, which are molecules of huge importance in medicine. Due to such reasons, intense efforts have been invested to maximize the mAbs production in bioreactors with hybridoma cell cultures. However, the optimal control of such sensitive bioreactors is an engineering problem difficult to solve due to the large number of state-variables with highly nonlinear dynamics, which often translates into a non-convex optimization problem that involves a significant number of decision (control) variables. Based on an adequate kinetic model adopted from the literature, this paper focuses on developing an in-silico (model-based, offline) numerical analysis of a fed-batch bioreactor (FBR) with an immobilized hybridoma culture to determine its optimal feeding policy by considering a small number of control variables, thus ensuring maximization of mAbs production. The obtained time stepwise optimal feeding policies of FBR were proven to obtain better performances than those of simple batch operation (BR) for all the verified alternatives in terms of raw material consumption and mAbs productivity. Several elements of novelty (i–iv) are pointed out in the “conclusions” section (e.g., considering the continuously added biomass as a control variable during FBR).

Keywords:

monoclonal antibodies (mAbs) maximization; hybridoma cell culture; fed-batch bioreactor dynamics optimization; raw material consumption; time stepwise operating policy

1. Introduction

Over the last decades, there has been a continuous trend in developing increasing numbers of effective cell or enzymatic reactors [1] to industrialize various biosynthesis technologies for producing industrial fine chemicals (including the bio-mimetic emulations, eventually by using Systems Biology tools, among others [2,3]) by using free-suspended or immobilized cell cultures/enzymes [4]. Effective and selective biosyntheses replaced energetically intensive chemical processes [5,6].

Bioreactors with microbial/animal cell cultures used for the production of a lot of products have been developed in simple constructive/operating alternatives as reviewed in Table 1 together with some examples. Complex alternatives are extensively discussed by [7,8]. In spite of their larger volumes, continuously mixing aerated tank reactors, operated in BR (batch), or FBR (semi-batch) modes, are preferred for processes requiring high oxygen transfers, and rigorous temperature/pH control, as the case of the approached mAbs production.

In addition to production capacity optimization, a crucial engineering problem to be solved concerns the development of optimal operating policies based on an available process model derived from on-/offline measurements. The model-based optimal operation of a bioreactor can be applied in two ways: (a) offline, in which an optimal operating policy is determined by using an adequate kinetic model of the previously identified process, and (b) online, using a simplified (often empirical) model and a classic state-parameter estimator, based on online recorded data [9,13,26,27].

Even if the bioprocess kinetics and biomass characteristics are known, solving this engineering problem is not an easy task due to the presence of multiple (often competing) objectives, technological constraints, and a significant degree of uncertainty originating from multiple sources [10,14,15]. All these parametric/model/data uncertainties require a frequent update of the bioprocess model, the optimal operating policies being determined by using both heuristic and deterministic (model-based) optimization rules [28]. Multi-objective criteria, including economic benefits, operating and materials costs, product quality, etc., are used to derive feasible optimal operating/control policies for various bioreactor types [15] by using specific numerical algorithms [11,13,14,29,30]

Due to such multiple reasons, for a given bioprocess, it is difficult to a priori indicate the least costly optimal operating policy of the chosen bioreactor. This problem is case-dependent.

To be consistent, the in-silico derived optimal operation policy of a bioreactor must be based on a simulation model that must include most of the key variables (particularly the biomass viability) in an adequate bioprocess kinetic model. The use of a deterministic (process mechanism-based) model, as the case here, is often preferred due to the physical significance of the terms/parameters, which make their verification possible vs. observations and the literature, even if repeated model updating is often necessary due to the high variability of the bioprocess. Typical objective functions were reviewed by [5,16].

On the other hand, this preliminary in-silico (model-based) analysis is necessary because: (a) New complex bioprocesses that are developed every year, requiring an engineering analysis to correctly determine the most profitable bioreactor operating policy; (b) There are multiple constructive/operating alternatives (Table 1), which need to be comparatively analyzed.

The offline (model-based) numerical derivation of an optimal operating policy of a bioreactor (FBR in this work) is very commonly used in the engineering practice and achieves satisfactory results [13,14,16,17,18,31,32,33], being supported by several arguments: (i) The obtained problem solution is susceptible to an immediate implementation, as long as the obtained FBR operating policy is based on an adequate (validated) structured kinetic model; (ii) Optimal operation of bioreactors is a crucial engineering issue because it can lead to consistent economic benefits. Our results presented in the Section 5 prove it; (iii) Most of the kinetic models of moderate complexity are very flexible. Thus if, at a certain batch, significant inconsistencies are observed between the model-predicted bioreactor dynamics and the experimental evidence, the optimization stage is applied again, after an intermediate step (between batches) necessary to offline improve the model adequacy (model updating).

In the simplest BR operating alternative, biomass/substrates and other nutrients are initially added. Even if this bioreactor type has a low production capacity, it is commonly used for slow bioprocesses. Despite several drawbacks, BRs are highly flexible and easy to operate [31]. Due to these reasons, various BRs alternatives are often model-based verified, namely: (i) multi-objective BR optimization (that is of the batch time, substrate/biocatalyst initial load) [9,10,22,28,34]; (ii) a SeqBR optimization by including a model updating step based on acquired information from past batches (“tendency modeling”, not approached here) [4,11,12,17,19].

In contrast, other operating variants are more effective, in particular, those in which substrates/nutrients and/or biomass/enzymes are intermittently (BRP, [23]), or continuously (FBR, cyclic FBR, or MASCR) added during the batch by following a certain optimal policy (usually in-silico determined). A comparative discussion of the all mentioned bioreactor types and operating alternatives is provided by [35]. Recent advances reported better performances obtained by using multiple (48) fully automated mini-FBR operated in parallel [36].

A simpler but faster bioreactor optimization route was adopted in this paper, by comparing different operating alternatives of the FBR vs. BR.

However, the flexibility of FBR involves a more considerable modeling and computational effort to achieve a time-varying optimal operating policy (that is, an adjustable feeding policy with substrates/biomass) to obtain an efficient control of the biological process (see some examples in Table 1). The optimization problem often translates into a constrained multi-objective FBR optimization [14,16,19,20]. The superiority of FBR vs. BR should be model-based verified for every bioprocess case. Additionally, FBRs are more difficult to operate than BRs, as long as the time stepwise optimal feeding policy requires different feeding substrate solutions of different concentrations, and separate different cell cultures to be fed over the batch. This is the price paid for achieving FBR best performances. This need to previously prepare different cell cultures, and substrate stocks to be fed for every “time-arc” (that is a batch-time division in which the feeding is constant) is offset by the net higher productivity of FBR compared to those of BR as discussed in the below Section 5, and proved in the literature [14,16,23,37,38].

One very important application of BRs/FBRs is the industrial production of mAbs (of quasi-uniform characteristics), that is a molecule of major importance in medicine, by the so-called “Hybridoma technology” [39,40] by using antibody-secreting hybridoma cell cultures. However, the large-scale production of mAbs by using mammalian cells in optimized BRs, or FBRs (that is, with the continuous feeding of glucose (GLC) and glutamine (GLN) substrates and nutrients with variable concentrations/feeding rates; [33,40]) is limited due to the engineering problems associated with the model-based optimization, i.e., (i) the presence of multiple/contradictory objectives; (ii) the significant degree of uncertainty originating from multiple sources, and (iii) the unwanted decline in the cell viability that may occur during the batch. A short review of how some such engineering problems are approached is made by [5,9,17,21,41,42]. Also, process inhibition by several by-products (such as ammonium, and lactic acid, LAC, etc.) and/or hyper-osmotic stress related to nutrient feeds and base additions to control pH, all raise serious issues regarding the optimal operation of the batch-bioreactor. To reduce the LAC production, the use of adapted CHO (Chinese Hamster Ovary) cells was suggested by [43].

Other optimization routes of the bioreactor, not approached here, are focused on (i) the bioprocess model extension to account for the intracellular factors related to the cellular metabolic fluxes, or (ii) by improving the FBR operation by also accounting for the uncertainty in the model structure and its parameters, that is, the so-called “robust/stochastic” optimization, (see [44] denoted LG17, or [28,45,46]). Thus, several FBR optimization alternatives have been proposed in the literature, as follows: (a) offline one-time optimization not considering (nominal) or by considering (robust) the parameter uncertainty, with an offline derived process model, see LG17; or [47]; (b) optimization of the expected value of the objective function by considering the parametric uncertainty [46]; and (c) online FBR optimization, in which the feeding policy is adapted based on the online acquired information, and used to improve the model adequacy and the operating policy for the rest of the batch [21,29]. These formulations are solved using various nonlinear programming (NLP), mixed-integer (MINLP), or other solvers [48,49], due to the presence of a large number of constraints [11,20,21,30].

This paper aims to in-silico (offline) determine the optimal (time stepwise) feeding policy of a FBR with substrates, and with an immobilized hybridoma culture (on porous alginate beads) to maximize the production of mAbs over the batch by considering a small number of key control variables and an experimentally validated bioprocess dynamic model. The performance of the FBR in all checked alternatives is compared to those of a BR, in terms of raw material consumption and mAbs productivity. The study will illustrate: (i) in which cases the FBR is more effective than the BR, by using a simple operating policy given on every time-interval of equal length (i.e., over N_div “time-arcs” in which the batch-time is divided; [37]), even if small N_div values are checked (2 or 5 here, explained in Section 3.6), and (ii) the major influence of the control/decision variable searching ranges (set to be narrower or wider) on the optimal FBR efficiency. Multiple elements of novelty (i–iv) of the paper are underlined in Section 6 “Conclusions”. The in-silico analysis of the paper proved how such an optimal FBR operation is simpler to be implemented and more flexible (by presenting a larger number of degree of freedom, due to multiple control variables, in spite of a moderate small N_div) compared to some optimal policies reported in the literature (e.g., [18]), that uses an exponential trajectory of the feeding liquid flow rate, and only the inlet levels of [GLC] and [GLN] as control variables, all being obtained by means of a hybrid deterministic (differential, intrinsic)–empiric (macroscopic) model.

As long as classical optimization of hybridoma FBR most of the time leads to some singular arc trajectories for biomass optimization, the advantage of the optimized FBR considering multiple control variables and a certain number of feed arcs that are not necessarily constant (that is, more subdivisions in the feed policy) becomes obvious, as this paper proves by using a detailed numerical analysis.

2. SBR Culture and Bioprocess Dynamics

A typical BR or FBR reactor is equipped with temperature, pressure (mechanical agitation, air sparger, pH, and dissolved oxygen (DO) control systems; see a simplified scheme in Figure 1). In the present study, the BR studied by LG17 (that is [44]) was adapted to investigate a semi-continuous fed-batch bioreactor (FBR) operation mode, allowing for the addition of the substrate solution (with a variable feed flow rate), and of viable biomass (immobilized inside alginate porous beads [50,51,52]) during the batch according to optimal policies to be determined. This in-silico (model-based) analysis will allow comparing the FBR vs. BR performances under various operating alternatives.

To in-silico solve this problem, the use of both bioprocess, and bioreactor dynamic models is necessary. To support engineering calculations, several attempts have been reported in the literature to obtain adequate dynamic models to predict the key-species behavior in either BR or FBR systems as a function of the extracellular nutrient/metabolite concentrations, such as: (i) Models in reasonably reduced formulations but including the dominant factors involved in the optimization of the mAbs production [12,17,19,33], or LG17. (ii) More sophisticated bioprocess kinetic models that explicitly include intracellular factors related to the cellular metabolic fluxes [33,34,41,53], which are associated with the central carbon metabolism (CCM) [3]. However, such complex structured models of (ii) type are very difficult to identify, and their use for bioreactor optimization is very limited.

In the present study, the bioprocess kinetic model (denoted by LGM) adopted for mAbs production is that proposed by LG17 ([44]). LGM is presented in Table 2, together with the associated rate constants. The hypotheses used to develop the kinetic expressions of LGM have been presented by LG17 and are not further discussed in this paper. LGM was verified over extensive experiments performed by the authors. Besides, it is worth mentioning that the LGM is the reduced form of a previously published extended kinetic model (EKM), also experimentally validated by [33,54]. In the LGM, the essential terms account for the inhibition/limitation effects on the reaction rates due to the presence of various intra-/extracellular metabolites, substrates, or by-products in a similar way done by other reported models [34,41,55,56,57,58].

As presented in Table 2, the bioprocess LGM was included in the dynamic model of the FBR to describe the key-species dynamics during the batch. The bioreactor initial conditions and time stepwise values of the control variables will be further explored.

To not complicate the computational step, the adopted ideal model of FBR is a classical one [1], developed with the following simplifying hypotheses: (i) the operation is isothermal, iso-pH, and iso-DO; (ii) it is self-understood that nutrients (that is, compounds playing roles of sources of carbon, nitrogen, and phosphorus; [55]) are added initially and during the FBR operation, in recommended quantities, and of a C/N/P ratio of ca. 100/5/1 wt., together with an excess of aeration (pure oxygen, if necessary) for ensuring an optimally biomass maintenance, and any growth limitation due to such factors; (iii) the volume of the perfectly mixed liquid phase (with no concentration gradients) increases according to the liquid feed flow rate time-varying policy; (iv) the limits of the volumetric liquid feed flow rate (F_L,j in Table 2) are adjusted to ensure a maximum reactor content dilution of 10–25% relative to the initial liquid volume (V_L,0) to not increase separation costs significantly; (v) there is negligible mass resistance for the transport of nutrients/substrates/products/oxygen into the liquid and in porous alginate beads; (vi) the substrates and solid carrier (less than one millimeter size), including the immobilized biomass, are initially added to the bioreactor, and then added during the batch according to an optimal feeding policy of the FBR to be determined; (vii) uniform solid particles are considered uniformly distributed in the homogeneous liquid phase, due perfect mixing conditions. In practice, perfect mixing conditions are never met. However, this is a usual hypothesis adopted for most of FBR models to not complicate them with hydrodynamic parameters difficult to be experimentally determined. A vigorous air sparging, and the use of a suitable mechanical mixing can support such a hypothesis.

The FBR is considered here with using immobilized biomass, for several reasons: (i) Its reported higher stability [50,51,52], and (ii) FBR operation with an online addition of the viable biomass (of variable concentration, coming from different cell cultures stocks prepared separately to be fed for every “time-arc”) should eventually be considered in our analysis (an option seldom described in the literature; refer to the optimal MASCR operation reported by [16]). To simplify the numerical analysis, the cell cultures stocks are assumed quasi-homogeneous and included as the Xv lump into the kinetic model.

The higher biomass fed concentration is used, the more vigorous aeration is employed (pure oxygen, if necessary). (iii) The used immobilized biomass is best suitable alternative for a variable (optimal) feeding of the FBR with biomass, because it is easy to dose [50,51,52].

From a mathematical point of view, the FBR dynamic model translates to a set of differential mass balances written for every considered species in the following general form:

\frac{{dC}_{i}}{dt} = \frac{F_{L}}{V_{L}} (C_{inlet, i} - C_{i}) \pm r_{i} (C (t), C_{0}, k); C_{i, 0} = C_{i} (t = 0)

(1)

where species index i relates to X_v, X_t, GLC, GLN, LAC, AMM, and mAbs from the abbreviation list.

The reaction rate r_i expressions together with the associated rate constants and other details are given in Table 2. In Equation (1), C = the vector of the species concentrations; C_o = initial vector C (at time t = 0); k = the model rate constant vector. The reactor content dilution (determined by the increasing V_L) is due to the variable F_L(t) term.

In the model (1), the viable biomass (X_v) dynamics is of the following generic form.

\frac{{dX}_{v}}{dt} = X_{v} f (ξ (t)) - X_{v} g (ξ (t)),

(2)

The rate expressions f(ξ(t)) and g(ξ(t)) in (2) proposed in the literature are of a Monod-type and include terms accounting for the growth and death inhibition, respectively, caused by the concentration of the extracellular species, such as [GLC], [GLN], [LAC], [AMM], etc., see LG17 and [33,41,56,57,58], or by intracellular metabolites related to the CCM [34]. The (X_v) mass balance in the LGM of Table 2 is made fairly adequate by including the main influential species and terms, as pointed out and experimentally proved by LG17. Other authors [55,57] have proposed empirical kinetic forms for (2). Regardless of the approach, the initial/fed cell density, its quality, immobilization type, and medium characteristics clearly play the central role on the biomass dynamics, and mAb production yield.

To determine the species dynamics over the batch time (t_f), the model (1) and (2) is solved with a proposed initial condition of C_i,0 = C_i (t = 0), and using the best medium conditions shown in Table 3. Except for the control variables F_L,0, [GLC]₀, [GLN]₀, and X_V,0, (X_t,0, = X_V,0) whose initial values are to be determined by the optimization.

There are several reasons why the adopted LGM is preferred in our optimization analysis: (i) Its adequacy was experimentally validated and was reportedly fair; (ii) The LGM is the slightly reduced form of moderate complexity (19 rate constants in Table 2) of an EKM (of 31 rate constants, developed and validated by [33,54]). (iii) Even if the EKM is more complex by including differential mass balances of additional 8 intermediates, the dynamics of the 8 key state-variables (X_v, X_t, GLC, GLN, V_L, LAC, AMM, and mAbs.) are fairly represented by both models. (iv) Being a reduced form of the EKM, the LGM is easier to be used, as it has fewer parameters (parsimonious principle). (v) Being simpler but adequate, the LGM is expected to offer interpretable results with less computational effort.

3. Formulation of the Bioreactor Optimization Problem

The optimal FBR operation will be compared with the reference BR used by LG17 with a nominal SPBR (SetPoint of the BR) presented in Table 3, corresponding to the species dynamics plotted in Figure 2. The realized BR modest performance is of Max[mAb] = 1254.6 (mg/L) over a batch of t_f = 100 h, as presented in Table 4 and Table 5 (SPBR line) In this BR simple operation, the substrates and immobilized biomass are initially loaded, and the mAb product is separated at the end of the batch. FBR operation is more complex, as below described.

3.1. Control Variables Selection

By analyzing the process model of Table 2, the natural option is to choose as control variables those that are related to reactor feeding with raw materials (GLC, GLN) and biomass (X_v), whose concentrations play the major role in mAbs production. Additionally, the liquid feed flow rate F_L will also be considered, being responsible for the reactor content dilution. Consequently, the selected control variables are as follows, including the inoculum size [34], for each time-arc (index “j”):

(i).: the continuously added liquid flow rate F_L,j (j = 1, …, N_div);
(ii).: the time stepwise added [GLC]_inlet,j; [GLN]_inlet,j; [X_v]_inlet,j (j = 1, …, N_div);
(iii).: the FBR initial condition, that is, the initial liquid flow rate F_L,0, and the initial substrates (as shown in Table 2), that is:

[GLC]₀ = [GLC](t = 0) = [GLC]_inlet,1;

(3)

[GLN]₀ = [GLN](t = 0) = [GLN]_inlet_,1; [X_v]₀ = [X_v](t = 0) = [X_v]_inlet_,1.

(4)

3.2. Objective Function (Ω) Choice

By considering the mentioned control variables, the FBR optimization consists of determining the initial conditions, and the optimal feeding policy for every time-interval during the batch that leads to the maximization of the [mAb] produced during the batch, that is:

Max Ω,

(5)

where: Ω = Max [mAb (t)].

The mAb(C(t), C_o, k, F_L) (t) in (5) is model-based evaluated over the whole batch time

(t) ∈ [0, t_f ],

(6)

According to Equations (5) and (6), the objective function (Ω) consists in maximization of the [mAb] produced during the whole batch, and not only of the final [mAb]. That is because, once [mAb] reaches its maximum, the batch process can be stopped. In the present analysis, it happened that, for the derived optimum policies SP1–SP3 (SetPoints of the FBR defined in Table 4), and for the SPBR (SetPoint of the BR presented in Table 3), the maximum [mAb] is reached at the batch end (Figure 2, Figure 3, Figure 4 and Figure 5).

3.3. Problem Constraints

(a): The FBR model (1)–(4) including the bioprocess kinetic model (Table 2);
(b): The FBR initial condition, that is: [GLC]₀, [GLN]₀; F_L,0; [X_v]₀ = [X_t]₀ (adopted);
(c): The initial [mAb]₀, [AMM]₀, [LAC]₀, adopted at values given in Table 3;
(d): To limit the excessive consumption of raw materials, feasible searching ranges are imposed to the control/decision variable (with limits specified in Table 4), that is:

[GLC]_inlet_,min ≤ [GLC]_inlet_,j ≤ [GLC]_inlet_,max; F_L,min ≤ F_L,j; F_L,0 ≤ F_L,max

(7)

[GLN]_inlet,min ≤ [GLN]_inlet,j ≤ [GLN]_inlet,max;
[X_v]_inlet,min ≤ [X_v]_inlet,j ≤ [X_v]_inlet,max

(8)

The imposed ranges for the control variables (tested wider or narrower) are related to not only the implementation possibilities, discussed by [12,17,52] and LG17, but also to economic reasons (minimum substrate consumption; effective control).

3.4. Searching (Control) Variables and Problem Formulation

To conclude, the FBR optimization problem consists in finding the optimal values of the initial and of the input levels of the four selected control variables, that is {[GLC]_inlet_,j; [GLN]_inlet_,j; [X_v]_inlet_,j; F_L_,0; F_L_,_j; j = 1, …, N_div}, over N_div time-intervals (“arcs”) of equal lengths Δt = t_f/N_div and under the specified operating constraints Equations (7) and (8), that maximize the chosen objective function Equations (5) and (6). In total, there are 4 × N_div searching variables (i.e., 20 if N_div = 5, or 8 if N_div = 2). The time-intervals of equal lengths Δt = t_f/N_div are obtained by dividing the batch time into N_div parts t_j−1 ≤ t ≤ t_j, where t_j = jΔt are switching points (where the reactor input is continuous and differentiable). Time-intervals are shown in the “Liquid volume dynamics:” row of Table 2.

3.5. The Problem Solution

The problem solution will indicate the optimal running conditions of the FBR using the best medium conditions shown in Table 3, which maximize the objective Equations (5) and (6) in the presence of constraints Equations (7) and (8). As above formulated, the optimal operating policy will be given for every of time-intervals (of equal lengths) distributed throughout the batch-time.

3.6. (N_div) and Operating Alternative Choice

The adopted LGM and the bioreactor dynamic model includes enough degrees of freedom to offer a wide range of FBR optimal operating alternatives that might be investigated (a–e) [not all being approached here], as follows:

(a): by choosing unequal time-arcs, of lengths to be determined by the optimization rule;
(b): by considering the whole batch time as an optimization variable;
(c): by increasing the number of equal time-arcs (N_div) to obtain a more “refined” and versatile FBR operating policy;
(d): by considering the search min/max limits of the control variables as unknown (to be determined);
(e): by feeding the bioreactor with solutions of uniform concentrations over a small/large number (N_div) of time-arcs.

The alternative (d) is unlikely because it might indicate unrealistic results, such as unreachable [X_v]_inlet_,j high levels, or substrate high levels which can inhibit the bioprocess; or even [X_v]_inlet_,j/substrates of too low levels that can stop the bioprocess. In our numerical analysis, carefully documented [X_v]_inlet_,j, and substrates upper bounds were tested instead to ensure the practical implementation possibility of the optimal results. The alternative (e) is also not feasible, even if a larger (N_div) will be used. That is because it is well-known that the variability in the feeding solutions over the batch time-arcs is the main degree-of-freedom used to obtain FBR optimal policies of superior quality [23,24]. By giving up to the variable concentrations in the feeding solutions, sub-optimally policies will be obtained.

Concerning the alternative (c), it is worth noting that, as (N_div) increases, and the above alternatives (a,b) are considered as well, the necessary computational effort grows significantly (due to considerable increase in the number of searching variables), thus hindering the quick (real-time) implementation of the derived operating policy. Additionally, multiple optimal operating policies can exist for an over-parameterized constrained optimization problem with a high nonlinearity, increasing the difficulty to quickly locate a feasible globally optimal operating policy.

Additionally, as the (N_div) increases, the operating policy is more difficult to implement since the optimal feeding policy requires a larger number of stocks with feeding substrate solutions of different concentrations, and different cell cultures stocks separately prepared to be fed for every time-arc of the FBR operation (a too expensive alternative). Also, the NLP problem is more difficult to solve because multiple optimal solutions may exist difficult to be discriminated and implemented. This is the case, for instance of an obtained optimal policy with a very high [X_v]_inlet_,j difficult to be ensured due to oxygenation limitations. Besides, FBR operation with using a larger number of small time-arcs (N_div) can raise special operating problems when including PAT (Process Analytical Technology) tools [59].

A brief survey of the FBR (SBR) optimization literature [37,60] reveals that a small number (N_div) <10 is commonly used due to the above-mentioned reasons. In fact, the present numerical analysis does not intend to exhaust all the possibilities of FBR optimization. Thus, an extended analysis of alternatives (a–e) of the FBR operation, or the influence of the parametric uncertainty deserves a separate investigation, which is beyond the scope of this paper.

To not complicate the computational analysis, only equal time-arcs have been tested, with the batch time t_f = 100 h. Two alternatives were adopted: (i) (N_div) = 5, in which equal time-arc-lengths of t_f/(N_div) = 20 h were used, and ii) (N_div) = 2, with equal time-arc-lengths of t_f/(N_div) = 50 h. However, the search min/max limits of the control variables have been varied in each tested case.

3.7. The Used Solvers

The mAb time-evolution in Equation (5) is determined by solving the bioreactor dynamic model Equations (1)–(4) with a tested initial condition of C_j,0 = C_j (t = 0), an imposed batch time t_f, and the optimal medium conditions of Table 3. The dynamic model solution was obtained with enough precision by using the low-order stiff integrator (“ode23s”) of the MATLAB™.

Because the bioreactor/bioprocess nonlinear model and optimization objective Equations (1)–(8) are constrained, the problem translates into a nonlinear optimization problem (NLP). To obtain a global solution with enough precision, the multi-modal optimization solver MMA of [48,61] has been used, as being proved in previous works to be more effective than common/commercial algorithms. The computational time was reasonably short (minutes) using a common PC (Core-I7 processor), thus offering a quick implementation of the FBR optimal policy.

4. In-Silico Optimization Results

Further FBR optimization calculations will be made by considering the biomass with the characteristics of LG17. It is also assumed that the biomass retains its characteristics after immobilization, reflected by the same LGM rate constants of Table 2.

The offline (model-based) FBR optimization problem Equations (5) and (6) is realized gradually, by using various searching ranges to limit the raw materials consumption, and the liquid volume excessive increase. Three optimization alternatives have been checked (SP denotes a set-point), as follows:

(1): FBR-SP1. For the adopted N_div = 5, with equal time-arcs, and by using narrow search intervals for the control variables [GLC]_inlet_,j; [GLN]_inlet_,j; [X_v]_inlet_,j; and F_L_,j specified in Table 4 (in the SP1 row), the obtained optimal operating policy SP1 (for every time-arc) is presented in Table 4, together with the key-species dynamics in Figure 3. Final liquid volume is 1.27 × V_L_,0
(2): FBR-SP2. For N_div = 5 and equal time-arcs, but using wider search intervals for the (above mentioned) control variables, as specified in Table 4 (in the SP2 row), the obtained optimal operating policy SP2 of the control variables is presented for every time-arc in Table 4, together with the key-species dynamics in Figure 4. The final liquid volume is 1.1 × V_L_,0.
(3): FBR-SP3. For an adopted smaller N_div = 2 with equal time-arcs, but using the same wide search intervals for the control variables as for the SP2 case (SP3 row in Table 4), the obtained optimal operating policy SP3 of the control variables is presented for every time-arc in Table 4, together with the key-species dynamics during the batch in Figure 5. The final liquid volume is 1.1 × V_L_,0.

A comparison of substrate and biomass consumption among the derived optimal operating alternatives of BR and FBR is presented in Table 5 together with the realized performances, that is the follows setpoints:

—: BR (with nominal SPBR of Figure 2, and Table 4) [44], or
—: BR of [54], or
—: FBR (with set-points SP1, SP2, SP3 of Figure 3, Figure 4 and Figure 5, Table 4),
—: FBR of [33].

The substrate consumption for FBR was evaluated with the formula:

\sum_{j = 1}^{N_{div}} F_{L} {[conc . species]}_{inlet, j} Δ t_{j}

. Alternatively, for the BR case, raw materials consumption is based on the only initial load. The optimization alternative (d) (mentioned in the above Section 3.6), that is the choice of even wider search intervals for the control variables, was not approached, being not feasible due to the previously explained reasons. Carefully documented [X_v]_inlet_,j upper bounds were tested instead to ensure the practical implementation of the obtained optimal operating policy.

5. Results and Discussion

(I).

In all the simulated alternatives, the FBR performance (in terms of produced mg mAb/L) is better than that realized by the BR (see the results summarized in Table 5), even if the overall batch time is the same (100 h), and a simple operating policy with equal time-arcs, in a small number (2–5) is used. The FBR productivity is up to 6× higher than that of the BR, while FBR is using fewer raw materials (Table 5). In Table 5, the mAb productivity is expressed in the absolute terms of Max [mAb] [mg/L]. Other indices, such as Max [mAb][mg/cells.h] can be used as well, by combining the data of Table 4 and Table 5. Being an intensive index, according to “Max [mAb][mg/cells.h]”, the BR appears more favorable because it uses less total [X_v]. However, due to the large value of the product vs. the used biomass, such a poor advantage of the BR becomes negligible.

(II).

The study points out the major influence of the control variable setting ranges (narrower or wider) used by the optimization rule, on the obtained efficiency of the FBR optimal policy. More specifically, according to the results of Table 5, it turns out that:

(IIa).: The GLC consumption during FBR is ca. 1/2 for SP1, or 1/3 for SP3 than that of the BR case. Similarly, fewer GLN was also consumed. The biomass (X_v) consumption is roughly the same because of a smaller N_div (SP3), or of narrower search ranges of control variables (SP1).
(IIb).: Non-uniform adding policy of (X_v) and substrates is better in the FBR case (SP2 and SP3; (Figure 4 and Figure 5d–g) than in the other cases, as the compensation of the biomass death is attempted (Figure 4 and Figure 5a,g), while maintaining a continuous increase in the produced mAb.
(IIc).: The price paid by the FBR-SP2 to achieve the best performances compared to BR-SPBR, and FBR-SP1, or FBR-SP3 is a higher consumption of raw materials, i.e., (vs. BR) of ca. 1.5× more GLC and GLN, and 10× more biomass.
(IId).: The FBR-SP3 (X_v_,0, and mAb net productivity) policy appears to be somehow intermediate between the FBR-SP2 and BR-SPBR. Compared to the SPBR, the raw material consumptions are smaller, but the realized FBR performances are better, because an operation with N_div = 2 is more versatile than that of the BR with the initial load being the only optimization option.
(IIe).: The used biomass is generally higher in the FBR cases compared to those of BR (Table 5). Thus, the used (X_v) is roughly the same for SP3 (N_div = 2), but 2× for FBR-SP1, or 10× for FBR-SP2.

(III).

The GLC consumption for FBR operating case depends on the inlet GLC policy (the [GLC]_inlet_,j term in Table 2), and on the used control variable dynamics (X_v, GLN).

(IV).

A comparison indicates that the GLC dynamics of the SP1 vs. SP2 of the FBR (with N_div = 5) is depending on not only the searching interval chosen for the control variable [GLC]_inlet_,j, but also on the other species dynamics. Thus, if one compares the (X_v) plots of SP1 in Figure 3 to that of FBR-SP2 in Figure 4, it is easy to observe and explain that when (X_v) is high, GLC consumption is also high, in spite of a larger inlet [GLC]. This clearly shows that the FBR optimization must consider all variables simultaneously.

(V).

In the BR-SPBR case of LG17, a less flexible feeding explains its modest mAb productivity. The BR species dynamics (in Figure 2) is comparable to those of FBR-SP3 (N_div = 2; in Figure 5), that is, GLC is quickly consumed during the first half of the batch, and the biomass displays a pronounced peak in the first half of the batch.

In the end, it is worth mentioning that the present numerical engineering analysis presents multiple elements of novelty as briefly mentioned in Section 6 “Conclusions”.

(VI).

Our results prove the multiple advantages obtained when using FBR operated with multiple control variables following very versatile optimal feeding policies consisting in time stepwise variable of: (i) the feeding liquid flow rate, (ii) the added [GLC], (iii) the added [GLN], and (iv) the added [X_v] over the batch.

The in-silico analysis of the paper proved how such an optimal FBR operation is leading to quick results, easy to interpret and to implement, being more flexible and effective due to a larger number of degree of freedom (coming from the multiple control variables, and from their variable time stepwise policy), in spite of an economically advantageous small number of employed time-arcs (N_div) compared to some optimal policies of a similar FBR reported in the literature. For instance, [18] uses only an exponential trajectory of the feeding liquid flow rate, and only the inlet levels of [GLC] and [GLN] as control variables, all being obtained by using a hybrid deterministic (differential, intrinsic)–empiric (macroscopic) model.

(VII).

The present in-silico (model-based) analysis have not been experimentally validated. However, as long as various forms of the used LGM were experimentally validated in a multiple and independent manner by LG17, and by [33,54] (Section 2), the results obtained by our numerical analysis shows sufficient credibility from the engineering point of view, from the following reasons:

(a): Even if an experimental validation of the derived optimal policy FBR-SP2 policy is missing, our paper presents a very strong engineering value by exemplifying, in a relatively simple manner, a numerical procedure (process model-based) that can be used to solve similar complex optimization problems of FBR.
(b): Such an FBR optimal control rule is possible because most kinetic models of moderate complexity are very flexible. Thus, if significant inconsistencies are observed between the model-predicted bioreactor dynamics (e.g., optimal policy SP2 in the present case) and the experimental data, then the optimization stage is applied again by using the same rule, but after performing an intermediate numerical-analysis step (between batches) necessary to improve the model adequacy (the so-called “model updating” based on the online measurements).

(VIII).

As displayed in Figure 2, Figure 3, Figure 4 and Figure 5 inhibition given by the increasing LAC, and AMM by-product concentration cannot be diminished by simple manipulations of the chosen control variables, even if the derived operating policy is an optimal one. However, the adopted kinetic model is able to fairly predict the dynamics of these inhibitory species. But the adverse side effects, such as a low pH, or a hyper-osmotic stress (due to the nutrient feeds and base additions to control pH) cannot be avoided by the above-used engineering (model-based) rules. As revealed in the literature, “biological” solutions are used instead to cope with such a problem. For instance, to reduce the LAC production, the use of adapted CHO (Chinese Hamster Ovary) cells can be an alternative [43].

(IX).

The used time-arcs of constant control variables are of 20 h (for SP1, SP2), and of 50 h (for SP3). Such an operation cannot raise special operating problems for a FBR, with also including PAT (Process analytical technology) tools. This is an additional argument not using a larger number of small time-arcs (N_div, see Section 3.6).

(X).

The comparison in Table 5 of the obtained optimal policies FBR-SP2 and FBR-SP3 with those from the literature for a BR [44,54], or a FBR [33] operated with the same cell culture, indicates better performances despite a longer batch time, and a larger number of substrate feeding solutions of those in the literature.

6. Conclusions

To conclude, the optimized FBR operation with a time stepwise control of the feeding policy reported better performances than the simple BR operation due to its higher flexibility in using biomass and substrates, even if a small number of equal time-arcs is used. The FBR major drawback is coming from its difficult operation, as long as the time stepwise optimal feeding policy requires different feeding substrate solution stocks of different concentrations, and separate different cell cultures stocks to be fed over the batch. This is the price paid for achieving improved FBR performances. However, its reported excellent performances fully justify the extra investment in implementing the FBR optimal operating policy. An economic global evaluation accounting for the product/raw material value can give a more accurate answer to such a sensitive issue.

The present optimization analysis proves its worth by including multiple elements of novelty, as follows: (i) An optimally operated FBR with a small number of time-arcs (below 10) and using wider but feasible ranges for setting the control variables can lead to high performances of the bioreactor. (ii) The major role played by the variable feeding with the viable biomass, leading to consider (X_v) as a control variable during FBR optimization (an option seldom discussed in the literature). (iii) The major influence of the control variable searching ranges (set to be narrower or wider) on the resulted efficiency of the FBR optimal policy. (iv) The model-based optimal operation of bioreactors is a very important engineering issue because it can lead to consistent economic benefits, as proved by the results presented in Table 5.

The present in-silico engineering analysis preserves its value by offering a useful example and comparative results for further engineering applications seeking for optimization of a FBR used to conduct complex bioprocesses.

The paper proves in a simple yet suggestive way how a lumped but adequately detailed dynamic model of a bioprocess can successfully support in-silico engineering evaluations by aiming to optimize the FBR operation if the nano-scale cell metabolism (including metabolic by-products) is somehow reflected by the bioreactor macro-scale dynamic model. In such a way, the considerable experimental and computational effort required to develop hybrid cell/bioreactor models of good quality is fully justified through the benefits of subsequent model-based engineering analyses that assist the FBR optimal operation.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

C_i	Species i concentration
$f (ξ (t)), g (ξ (t))$	Complex functions in Equation (2) accounting for the biomass growth inhibition or death inhibition respectively
F_L	Liquid feed flow rate
$k$	Model rate constant vector
$k_{j}$ , $K_{j}$ , $Q_{j}$ , $Y_{j}$ , $m_{j}$ , etc.	Rate constants given in Table 2, the index “j” relates to the all forms of these constants.
N_div	Number of equal divisions (“time-arcs”) of the batch time t_f
r_i	Species ”i” reaction rate
t	Time
t_f	Batch time
V_L	Liquid volume in the bioreactor
X_V	Viable cell density
X_t	Total cell density
Index
f	Final
inlet	Value in the bioreactor inlet
max	Maximum
min	Minimum
0,o	Initial
Greeks
$Δ$	Discrete interval
$α_{1}$ , γ, λ, μ_max, μ_d_,max	Rate constants given in Table 2
μ	Specific growth rate of viable biomass in Table 2
μ_d	Specific cell death rate in Table 2
Ω	The optimization objective function in eqn.(1)
AMM	Ammonia
BR	Batch reactor
BRP	Batch reactor with intermittent addition of biomass/raw materials
CCM	Central carbon metabolism
DO	Dissolved oxygen
EKM	Extended kinetic model of [33,54]
FBR	Fed-batch bioreactor
GLC	Glucose
GLN	Glutamine
LAC	Lactate
LDH	Lactate dehydrogenase
LGM	The kinetic model of [44]
LG17	Notation for the reference [44]
MA(S)CR	Mechanically agitated (semi-)continuous reactor
mAb	Monoclonal antibody
Max	Maximum
MINLP	Mixed-integer NLP
MMA	The multi-modal optimization solver of [48,61]
NADH	Nicotinamide adenine dinucleotide (reduced form)
NLP	Nonlinear programming (numerical rules for solving optimization problems)
PAT	Process Analytical Technology tools.
SP	Bioreactor setpoint
SPBR	BR nominal setpoint of [44] (Table 3)
SBR	Semi-batch reactor
SeqBR	Sequential batch reactor
TPFB	Three-phase fluidized-bed bioreactor
[x]	x species concentration

References

Moser, A. Bioprocess Technology—Kinetics and Reactors; Springer: Berlin, Germany, 1988. [Google Scholar]
Heinemann, M.; Panke, S. Synthetic Biology—Putting engineering into biology. Bioinformatics 2006, 22, 2790–2799. [Google Scholar] [CrossRef] [Green Version]
Maria, G. In-Silico Design of Genetic Modified Micro-Organisms (GMO) of Industrial Use, by Using Systems Biology and (Bio)chemical Engineering Tools, Simi Valley; Juniper Publ.: Irvine, CA, USA, 2018; ISBN 978-1-946628-12-1. Available online: https://juniperpublishers.com/ebook-info.php (accessed on 28 November 2020).
Lübbert, A.; Jørgensen, S.B. Bioreactor performance: A more scientific approach for practice. J. Biotechnol. 2001, 85, 187–212. [Google Scholar] [CrossRef]
Engasser, J.-M. Bioreactor engineering: The design and optimization of reactors with living cells. Chem. Eng. Sci. 1988, 43, 1739–1748. [Google Scholar] [CrossRef]
Gavrilescu, M.; Chisti, Y. Biotechnology—A sustainable alternative for chemical industry. Biotechnol. Adv. 2005, 23, 471–499. [Google Scholar] [CrossRef] [PubMed]
Chaudhuri, J. Bioreactors for Tissue Engineering Principles, Design and Operation; Al-Rubeai, M., Ed.; Springer: Berlin, Germany, 2005. [Google Scholar]
Nedovic, V.; Willaert, R. Applications of Cell Immobilisation Technology; Springer: Amsterdam, The Netherlands, 2005. [Google Scholar]
Abel, O.; Marquardt, W. Scenario-integrated online optimisation of batch reactors. J. Process Control 2003, 13, 703–715. [Google Scholar] [CrossRef]
Bonvin, D. Optimal operation of batch reactors—A personal view. J. Process Control 1998, 8, 355–368. [Google Scholar] [CrossRef]
Martinez, E. Batch-to-batch optimization of batch processes using the STATSIMPLEX search method. In Proceedings of the 2nd Mercosur Congress on Chemical Engineering, Rio de Janeiro, Costa Verde, Brasil, 31 March 2005. [Google Scholar]
Dewasme, L.; Côte, F.; Filee, P.; Hantson, A.-L.; Wouwer, A.V. Macroscopic Dynamic Modeling of Sequential Batch Cultures of Hybridoma Cells: An Experimental Validation. Bioengineering 2017, 4, 17. [Google Scholar] [CrossRef] [Green Version]
DiBiasio, D. Introduction to the control of biological reactors. In Chemical Engineering Problems in Biotechnology; Shuler, M.L., Ed.; American Institute of Chemical Engineers: New York, NY, USA, 1989; pp. 351–391. ISBN 9780816904693. [Google Scholar]
Maria, G.; Crişan, M. Operation of a mechanically agitated semi-continuous multi-enzymatic reactor by using the Pareto-optimal multiple front method. J. Process Control 2017, 53, 95–105. [Google Scholar] [CrossRef]
Smets, I.Y.; Claes, J.E.; November, E.J.; Bastin, G.P.; Van Impe, J.F. Optimal adaptive control of (bio)chemical reactors: Past, present and future. J. Process Control 2004, 14, 795–805. [Google Scholar] [CrossRef]
Şcoban, A.G.; Maria, G. Model-based optimization of the feeding policy of a fluidized bed bioreactor for mercury uptake by immobilizedPseudomonas putidacells. Asia Pac. J. Chem. Eng. 2016, 11, 721–734. [Google Scholar] [CrossRef]
Dewasme, L.L.; Amribt, Z.; Santos, L.O.; Hantson, A.-L.; Bogaerts, P.; Wouwer, A.V. Hybridoma cell culture optimization using nonlinear model predictive control. In Proceedings of the 12th IFAC Symposium on Computer Applications in Biotechnology, Mumbai, India, 16–18 December 2013; International Federation of Automatic Control: Laxenburg, Austria, 2013; Volume 46, pp. 60–65. [Google Scholar]
Amribt, Z.; Dewasme, L.; Wouwer, A.V.; Bogaerts, P. Optimization and robustness analysis of hybridoma cell fed-batch cultures using the overflow metabolism model. Bioprocess Biosyst. Eng. 2014, 37, 1637–1652. [Google Scholar] [CrossRef] [PubMed]
Binette, J.-C.; Srinivasan, B. On the Use of Nonlinear Model Predictive Control without Parameter Adaptation for Batch Processes. Processes 2016, 4, 27. [Google Scholar] [CrossRef] [Green Version]
Franco-Lara, E.; Weuster-Botz, D. Estimation of optimal feeding strategies for fed-batch bioprocesses. Bioprocess Biosyst. Eng. 2005, 28, 71. [Google Scholar] [CrossRef] [Green Version]
Mendes, R.; Rocha, I.; Pinto, J.P.; Ferreira, E.C.; Rocha, M. Differential Evolution for the Offline and Online Optimization of Fed-Batch Fermentation Processes. In Advances in Differential Evolution. Studies in Computational Intelligence, Chakraborty; Chakraborty, U.K., Ed.; Springer: Berlin, Germany, 2008; pp. 299–317. [Google Scholar]
Crisan, M.; Maria, G. Modular Simulation to Determine the Optimal Operating Policy of a Batch Reactor for the Enzymatic Fructose Reduction to Mannitol with the in situ Continuous Enzymatic Regeneration of the NAD Cofactor. Rev. Chim. 2017, 68, 2196–2203. [Google Scholar] [CrossRef]
Maria, G. Enzymatic reactor selection and derivation of the optimal operation policy, by using a model-based modular simulation platform. Comput. Chem. Eng. 2012, 36, 325–341. [Google Scholar] [CrossRef]
Maria, G.; Crisan, M. Evaluation of optimal operation alternatives of reactors used ford-glucose oxidation in a bi-enzymatic system with a complex deactivation kinetics. Asia Pac. J. Chem. Eng. 2014, 10, 22–44. [Google Scholar] [CrossRef]
Maria, G. Model-based heuristic optimized operating policies for d-glucose oxidation in a batch reactor with pulsate addition of enzyme. Comput. Chem. Eng. 2007, 31, 1231–1241. [Google Scholar] [CrossRef]
Lee, J.; Lee, K.S.; Lee, J.H.; Park, S. An online batch span minimization and quality control strategy for batch and semi-batch processes. Control Eng. Pract. 2001, 9, 901–909. [Google Scholar] [CrossRef]
Ruppen, D.; Bonvin, D.; Rippin, D. Implementation of adaptive optimal operation for a semi-batch reaction system. Comput. Chem. Eng. 1998, 22, 185–199. [Google Scholar] [CrossRef]
Srinivasan, B.; Bonvin, D.; Visser, E.; Palanki, S. Dynamic optimization of batch processes. Comput. Chem. Eng. 2003, 27, 27–44. [Google Scholar] [CrossRef]
Rao, M.; Qiu, H. Process Control Engineering: A Textbook for Chemical, Mechanical and Electrical Engineers; Gordon and Breach Science Publ.: Amsterdam, The Netherlands, 1993. [Google Scholar]
Hartig, F.; Keil, F.J.; Luus, R. Comparison of Optimization Methods for a Fed-Batch Reactor. Hung. J. Ind. Chem. 1995, 23, 81–160. [Google Scholar]
Koutinas, M.; Kiparissides, A.; Pistikopoulos, E.N.; Mantalaris, A. Bioprocess systems engineering: Transferring traditional process engineering principles to industrial biotechnology. Comput. Struct. Biotechnol. J. 2012, 3, e201210022. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Tsangaris, D.; Baltzis, B. Evaluation of batch and semi-batch reactor operation for enzymatic reactions with inhibitory kinetics. Chem. Eng. Sci. 1996, 51, 2757–2762. [Google Scholar] [CrossRef]
Kontoravdi, C.; Pistikopoulos, E.N.; Mantalaris, A. Systematic development of predictive mathematical models for animal cell cultures. Comput. Chem. Eng. 2010, 34, 1192–1198. [Google Scholar] [CrossRef] [Green Version]
Ozturk, S.S.; Palsson, B.Ø. Effect of initial cell density on hybridoma growth, metabolism, and monoclonal antibody production. J. Biotechnol. 1990, 16, 259–278. [Google Scholar] [CrossRef] [Green Version]
Koller, M. A Review on Established and Emerging Fermentation Schemes for Microbial Production of Polyhydroxyalkanoate (PHA) Biopolyesters. Fermentation 2018, 4, 30. [Google Scholar] [CrossRef] [Green Version]
Sawatzki, A.; Hans, S.; Narayanan, H.; Haby, B.; Krausch, N.; Sokolov, M.; Glauche, F.; Riedel, S.L.; Neubauer, P.; Bournazou, M.N.C. Accelerated Bioprocess Development of Endopolygalacturonase-Production with Saccharomyces cerevisiae Using Multivariate Prediction in a 48 Mini-Bioreactor Automated Platform. Bioengineering 2018, 5, 101. [Google Scholar] [CrossRef] [Green Version]
Maria, G.; Dan, A. Derivation of optimal operating policies under safety and technological constraints for the acetoacetylation of pyrrole in a semi-batch catalytic reactor. Comput. Chem. Eng. 2011, 35, 177–189. [Google Scholar] [CrossRef]
Avili, M.G.; Fazaelipoor, M.H.; Jafari, S.A.; Ataei, S.A. Comparison between Batch and Fed-Batch Production of Rhamnolipid by Pseudomonas aeruginosa. Iran. J. Biotechnol. 2012, 10, 263–269. Available online: https://www.researchgate.net/publication/288966207 (accessed on 3 June 2004).
Lee, G.M.; Varma, A.H.; Palsson, B.O. Production of monoclonal antibody using free-suspended and immobilized hybridoma cells: Effect of serum. Biotechnol. Bioeng. 1991, 38, 821–830. [Google Scholar] [CrossRef] [Green Version]
Hauser, H.; Wagner, R. Mammalian Cell Biotechnology in Protein Production; Walter de Gruyter Publ.: New York, NY, USA, 1997. [Google Scholar]
Dorka, P. Modelling Batch and Fed-Batch Mammalian Cell Cultures for Optimizing MAb Productivity. Master’s Thesis, University of Waterloo, Waterloo, ON, Canada, 2007. [Google Scholar]
Li, F.; Vijayasankaran, N.; Shen, A.; Kiss, R.; Amanullah, A. Cell culture processes for monoclonal antibody production. mAbs 2010, 2, 466–479. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Freund, N.W.; Croughan, M.S. A Simple Method to Reduce both Lactic Acid and Ammonium Production in Industrial Animal Cell Culture. Int. J. Mol. Sci. 2018, 19, 385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Liu, Y.; Gunawan, R. Bioprocess optimization under uncertainty using ensemble modeling. J. Biotechnol. 2017, 244, 34–44. [Google Scholar] [CrossRef] [PubMed]
Kadam, J.; Schlegel, M.; Srinivasan, B.; Bonvin, D.; Marquardt, W. Dynamic optimization in the presence of uncertainty: From offline nominal solution to measurement-based implementation. J. Process Control 2007, 17, 389–398. [Google Scholar] [CrossRef]
Terwiesch, P.; Agarwal, M.; Rippin, D.W. Batch unit optimization with imperfect modelling: A survey. J. Process Control 1994, 4, 238–258. [Google Scholar] [CrossRef]
Watanabe, N.; Nishimura, Y.; Matsubara, M. Optimal design of chemical processes involving parameter uncertainty. Chem. Eng. Sci. 1973, 28, 905–913. [Google Scholar] [CrossRef]
Maria, G. A review of algorithms and trends in kinetic model identification for Chemical and Biochemical Systems. Chem. Biochem. Eng. Q. 2004, 18, 195–222. [Google Scholar]
Maria, G. Adaptive random search and shortcut techniques for process model identification and monitoring. In Foundations of Computer Aided Process Operations; Pekny, J.F., Blau, G.E., Carnahan, B., Eds.; AIChE symposium series; American Institute of Chemical Engineers: New York, NY, USA, 1998; Volume 94, pp. 351–359. [Google Scholar]
Selimoglu, S.M.; Ayyildiz-Tamis, D.; Gurhan, I.D.; Elibol, M. Purification of alginate and feasible production of monoclonal antibodies by the alginate-immobilized hybridoma cells. J. Biosci. Bioeng. 2012, 113, 233–238. [Google Scholar] [CrossRef]
Selimoglu, S.M.; Elibol, M. Alginate as an immobilization material for MAb production via encapsulated hybridoma cells. Crit. Rev. Biotechnol. 2010, 30, 145–159. [Google Scholar] [CrossRef]
Lee, G.M.; Kim, S.J.; Palsson, B.O. Enhanced specific antibody productivity of calcium alginate-entrapped hybridoma is cell line-specific. Cytotechnology 1994, 16, 1–15. [Google Scholar] [CrossRef] [Green Version]
Quiroga-Campano, A.L.; Panoskaltsis, N.; Mantalaris, A. Energy-based culture medium design for biomanufacturing optimization: A case study in monoclonal antibody production by GS-NS0 cells. Metab. Eng. 2018, 47, 21–30. [Google Scholar] [CrossRef] [PubMed]
Kiparissides, A.; Koutinas, M.; Kontoravdi, C.; Mantalaris, A.; Pistikopoulos, E.N. ‘Closing the loop’ in biological systems modeling—From the in silico to the in vitro. Automatica 2011, 47, 1147–1155. [Google Scholar] [CrossRef]
Venables, D. Two Stage Chemostat Studies of Hybridoma Growth, Nutrient Utilization, and Monoclonal Antibody Production. Ph.D. Thesis, University of Surrey (MI), Stag Hill, Guildford, UK, 1994. [Google Scholar]
Bree, M.A.; Dhurjati, P.; Geoghegan, R.F.; Robnett, B. Kinetic modelling of hybridoma cell growth and immunoglobulin production in a large-scale suspension culture. Biotechnol. Bioeng. 1988, 32, 1067–1072. [Google Scholar] [CrossRef]
Zeng, A.P.; Deckwer, W.D.; Hu, W.S. Determinants and rate laws of growth and death of hybridoma cells in continuous culture. Biotechnol. Bioeng. 1998, 57, 642–654. [Google Scholar] [CrossRef]
Lee, Y.-K.; Yap, P.-K.; Teoh, A.-P. Correlation between steady-state cell concentration and cell death of hybridoma cultures in chemostat. Biotechnol. Bioeng. 1995, 45, 18–26. [Google Scholar] [CrossRef] [PubMed]
Bharat, A. Process Analytical Technology (PAT), Msc-Diss., P.D.V.V.P.F’S. College of Pharmacy, AhmedNagar (India). 2013. Available online: https://www.slideshare.net/anjalibharat19/process-analytical-tchnology (accessed on 28 November 2020).
Loeblein, C.; Perkins, J.; Srinivasan, B.; Bonvin, D. Performance analysis of online batch optimization systems. Comput. Chem. Eng. 1997, 21, S867–S872. [Google Scholar] [CrossRef] [Green Version]
Maria, G. ARS Combination with an Evolutionary Algorithm for Solving MINLP Optimization Problems. In Modelling, Identification and Control; Hamza, M.H., Ed.; IASTED/ACTA Press: Anaheim, CA, USA, 2003; pp. 112–118. Available online: http://www.sigmod.org/sigmod/dblp/db/conf/mic/mic2003.html (accessed on 3 June 2004).

Sample Availability: Samples of the compounds are not available from the authors.

Figure 1. The fed-batch bioreactor (FBR) simplified scheme. The bioreactor is operated in a fed-batch mode, with time stepwise continuous feeding addition of substrates, nutrients, and immobilized biomass (under one millimeter size alginate beads) at levels to be determined by optimization for each “time-arc”.

Figure 2. The simulated SPBR (SetPoint of the BR given in Table 3). Key species dynamics (that is time-trajectories) are generated by using the LGM of LG17 [44].

Figure 3. The simulated SP1 (the optimum policy SetPoint no.1 of the FBR, defined in Table 4). The plots refer to the dynamics of the key species (a,b), and of the liquid volume (c). The plots (d–g) refer to the time stepwise optimal policy of the control variables [GLC]_inlet(t) (d); F_L(t) (e); [GLN]_inlet(t) (f), and [X_v]_inlet(t) (g) for the approached FBR (the running details are given in Table 4). Species trajectories are generated by using the LGM model (Table 2).

Figure 4. The simulated SP2 (the derived optimum policy SetPoint no.2 of the FBR, defined in Table 4). The plots refer to the dynamics of the key species (a,b), and of the liquid volume (c). The plots (d–g) refer to the time stepwise optimal policy of the control variables [GLC]_inlet(t) (d); F_L(t) (e); [GLN]_inlet(t) (f), and [X_v]_inlet(t) (g) for the FBR approached in this paper (running details are given in Table 4). Species trajectories are generated by using the of LGM model (Table 2).

Figure 5. The simulated SP3 (the derived optimum policy SetPoint no.3 of the FBR, defined in Table 4). The plots refer to the dynamics of the key species (a,b), and of the liquid volume (c). The plots (d,g) refer to the time stepwise optimal policy of the control variables [GLC]_inlet(t) (d); F_L(t) (e); [GLN]_inlet(t) (f), and [X_v]_inlet(t) (g) for the FBR approached in this paper (running details are given in Table 4). Species trajectories are generated by using the of LGM model (Table 2).

Table 1. The main constructive and operating alternatives of bioreactors [1,4,5].

Reactor Type	Notation [Examples]	Operation; Modeling Hypotheses
	(i) simple Batch Reactor (BR) Examples: [9,10]	isothermal, iso-pH, and iso-DO (air sparger); perfectly mixed liquid phase (with no concentration gradients, by using mechanical agitation). Reactants/biomass added at the beginning of the batch only.
	(ii) SeqBR (Sequential Batch-to-batch Reactor) Examples: [11,12]	Ibidem. Reactants and/or biomass added at the beginning of each batch, in optimized amounts (to be determined)
	(iii) Semi-Batch (fed-batch) Reactor (SBR or FBR). Examples: [13,14,15,16,17,18,19,20,21]	Ibidem. Substrates/biomass/supplements added during the batch by following a certain (optimal) policy (to be determined)
	(iv) BRP (Batch Reactor with intermittent/Pulse-like additions of biocatalyst/substrates). Examples: [22,23,24,25]	Ibidem. Reactants and/or biomass added during the batch in a Pulse-like additions of equal/uneven solution volumes, with a certain frequency (to be determined)
	(v) continuously operated packed-bed columns, FXBR (FiXed-Bed continuous bioReactor) Examples: [23]	immobilized enzyme on a porous support packed in columns; continuous fed of the substrate/nutrient solution; continuous solution output; various aeration alternatives. Model hypotheses: isothermal, ideal plug-flow reactor of constant volume, with model dynamic terms allowing simulating transient operating conditions and the continuous enzyme/biomass deactivation
	(vi) MA(S)CR (Mechanically Agitated (Semi-)Continuous Reactor) (three-phases). Examples: [14]	immobilized enzyme on porous support suspended in the mechanically agitated bioreactor, with sparged gas (air); continuous fed of the substrate/nutrient solution, with/without continuous evacuation; Model hypotheses: isothermal, ideal perfectly mixed liquid phase (with no concentration gradients, by using mechanical agitation, aeration), with model dynamic terms allowing simulating transient operating conditions and the continuous enzyme/biomass deactivation. Substrates/biomass can be added with a constant/variable feed flow rate (to be determined).

Table 2. Key-species mass balances in the fed-batch bioreactor FBR model, including the bioprocess kinetic model LGM of LG17, together with the associated rate constants. Note: the ideal model below (of homogeneous liquid composition), is neglecting the mass transport resistance in the porous beads. Rate constants have been estimated by [44] from experiments that use the mammalian hybridoma cell culture of [54].

Species	Parameters	Remarks
Biomass balance: Viable biomass balance: $\frac{d X_{v}}{d t} = \frac{F_{L}}{V_{L}} (X_{v, i n l e t, j} - X_{v}) + (μ - μ_{d}) X_{v}$ $X_{v, i n l e t, j}$ = control variable; j = 1, …, $N_{d i v}$ time stepwise values to be optimized; $X_{v, 0}$ = $X_{v, i n l e t, 1}$ ; $X_{v, i n l e t} (t)$ to be optimized; Total biomass balance: $\frac{d X_{t}}{d t} = \frac{F_{L}}{V_{L}} (X_{t, i n l e t} - X_{v}) + μ X_{v} - K_{l y s i s} (X_{t} - X_{v})$ $X_{t, i n l e t} = 0$ (adopted); $X_{t, 0}$ = $X_{v, 0}$ where: $μ = μ_{\max} (\frac{[G L C]}{K_{g l c} + [G L C]}) (\frac{[G L N]}{K_{g l n} + [G L N]})$ $(\frac{K I_{l a c}}{K I_{l a c} + [L A C]}) (\frac{K I_{a m m}}{K I_{a m m} + [A M M]})$ $μ_{d} = \frac{μ_{d, \max}}{1 + {(K_{d, a m m} / [A M M])}^{2}}$	$K_{l y s i s}$ = 0.0551 h⁻¹ $μ_{\max}$ = 0.058 h⁻¹ $K_{g l c}$ = 0.75 mM $K_{g l n}$ = 0.075 mM $K I_{l a c}$ = 172 mM $K I_{a m m}$ = 28.5 mM $μ_{d, \max}$ = 0.03 h⁻¹ $K_{d, a m m}$ = 1.76 mM	LGM
Balance of other species: $\frac{d [G L C]}{d t} = \frac{F_{L}}{V_{L}} ({[G L C]}_{i n l e t, j} - [G L C]) - Q_{g l c} X_{v}$ , ${[G L C]}_{i n l e t, j}$ = control variable; j = 1, …, $N_{d i v}$ time stepwise values to be optimized; ${[G L C]}_{0}$ = ${[G L C]}_{i n l e t, 1}$ ; ${[G L C]}_{i n l e t} (t)$ to be optimized where: $Q_{g l c} = \frac{μ}{Y_{x, g l c}} + m_{g l c}$	$Y_{x, g l c}$ = 1.06 × 10⁸ cell/mmol $m_{g l c}$ = 4.85 × 10⁻¹⁴ mmol/cell.h	LGM
$\frac{d [G L N]}{d t} = \frac{F_{L}}{V_{L}} ({[G L N]}_{i n l e t, j} - [G L N]) - Q_{g_{l n}} X_{v} - K_{d, g_{l n}} [G L N]$ , ${[G L N]}_{i n l e t, j}$ = control variable; j = 1, …, $N_{d i v}$ time stepwise values to be optimized; ${[G L N]}_{0}$ = ${[G L N]}_{i n l e t, 1}$ ; ${[G L N]}_{i n l e t} (t)$ to be optimized, where: $Q_{g_{l n}} = \frac{μ}{Y_{x, g_{l n}}} + m_{g_{l n}}$ $K_{d, g_{l n}}$ = 0.0096 h⁻¹; $Y_{x, g_{l n}}$ = 5.57 × 10⁸ cell/mmol $m_{g_{l n}}$ = $α_{1}$ = −0.00067 mmol/cell.h		LGM
Liquid volume dynamics: $\frac{d V_{L}}{d t} = F_{L, j}; F_{L} = {\begin{matrix} F_{L, 0} if 0 \leq t < T 1 \\ F_{L, 1} if T 1 \leq t < T 2 \\ F_{L, 2} if T 2 \leq t < T 3 \\ F_{L, 3} if T 3 \leq t < T 4 \\ F_{L, 4} if T 4 \leq t < t_{f} \end{matrix}$ , (a) For the adopted $N_{d i v}$ = 5, the j = 1, …, $N_{d i v}$ time-arcs switching points are: T1 = 20 h.; T2 = 40 h.; T3 = 60 h.; T4 = 80 h.; $t_{f}$ = 100 h., where $F_{L, 0} - F_{L, 4}$ time stepwise values are to be determined together with the other control variables to ensure an optimal FBR operation; (b) For the adopted $N_{d i v}$ = 2, the j = 1, …, $N_{d i v}$ time-arcs switching points are: T1 = 50 h.; $t_{f}$ = 100 h., with $F_{L} = {\begin{matrix} F_{L, 0} if 0 \leq t < T 1 \\ F_{L, 1} if T 1 \leq t < t_{f} \end{matrix}$ where $F_{L, 0} - F_{L, 1}$ time stepwise values are to be determined together with the other control variables to ensure an optimal FBR operation;		This paper
$\frac{d [L A C]}{d t} = \frac{F_{L}}{V_{L}} ({[L A C]}_{i n l e t} - [L A C]) + Q_{l a c} X_{v}$ , ${[L A C]}_{i n l e t} = {[L A C]}_{0}$ = 0 where: $Q_{l a c} = Y_{l a c, g l c} Q_{g l c}$	$Y_{l a c, g l c}$ = 1.4 L	LGM
$\frac{d [A M M]}{d t} = \frac{F_{L}}{V_{L}} ({[A M M]}_{i n l e t} - [A M M]) + Q_{a m m} X_{v} + K_{d, g l n} [G L N]$ ${[A M M]}_{i n l e t}$ = 0; ${[A M M]}_{0}$ = 0.31 mM where: $Q_{a m m} = Y_{a m m, g_{l n}} Q_{g_{l n}}$ $K_{d, g_{l n}}$ = 0.0096 h⁻¹; $Y_{a m m, g_{l n}}$ = 0.427 L		LGM
$\frac{d [m A b]}{d t} = \frac{F_{L}}{V_{L}} ({[m A b]}_{i n l e t} - [m A b]) + (2 - γ μ) λ X_{v}$ ${[m A b]}_{i n l e t}$ = 0; ${[m A b]}_{0}$ = 80.6 mg/L	γ = 0.1 h λ = 7.21 × 10⁻⁹ mg/(cell·h)	LGM

Table 3. The nominal operating conditions (SPBR) of LG17 for the batch bioreactor BR with suspended mammalian hybridoma cell culture.

Parameter	Nominal Value	Remarks (*)
Total cell initial density (X_t,0)	2 × 10⁸ Cell/L	Ref. to reactor-lq.
Viable cell initial density (X_V,0)	2 × 10⁸ Cell/L	Ref. to reactor-lq.
Glucose initial concentration, [GLC]₀	29.1, mM
Glutamine initial concentration, [GLN]₀	4.9, mM
Lactate initial concentration, [LAC]₀	0, mM
Ammonia initial concentration, [AMM]₀	0.31, mM
Monoclonal antibody initial concentration, [mAb]₀	80.6, mg/L	Ref. to reactor-lq.
Temperature	35–37 °C	[42]
pH (buffer, using CO₂ injection)	7	See an optimal policy given by [42]
Aeration in excess, nutrients in sufficient amounts		[42,55]
Initial volume of the liquid in the bioreactor (V_L,0)	1 L	LGM
Batch time (t_f)	approx. 100 h.	LGM

(*) Ref. to reactor-lq. = Value relative to the liquid volume of the reactor.

Table 4. Model-based derived optimal operating policies for the approached FBR comparatively to performances of LG17 batch bioreactor (BR) with immobilized mammalian hybridoma cell culture. Biomass and medium characteristics are those given by LG17 in Table 3. The larger number of displayed digits comes from the numerical simulations.

Reactor SP	Searching Policy	Control Variables					Obs.
SPBR (BR) Optimal Values	Sensitivity Analysis (Exhaustive)	Initial Values of the BR Content
		[GLC], mM		[GLN], mM	X_v,0 = X_t,0, Cell/L	Max [mAb](t), (mg/L)
		29.1		4.9	2 × 10⁸	1254.6	LGM
SP1^(a) (FBR) Optimal values ^(d)	Searching variables	F_L, L/h.	[GLC]_inlet, mM	[GLN]_inlet, mM	X_v,inlet Cell/L		This paper
	Searching ranges	(10⁻⁴–10⁻²)	(25–100)	(5–25)	(2 × 10⁸–2 × 10⁹)
	Multi-dimensional optimization	Inlet optimal values of the FBR control variables
	Multi-dimensional optimization	F_L, ^(b,c), L/h.	[GLC]_inlet mM	[GLN]_inlet mM	X_v,inlet Cell/L	Max [mAb](t), (mg/L)
	Time interval (0, 20) h.	10⁻³	96.62	17.75	2 × 10⁸		This paper
	Time interval (20, 40) h.	9.55 × 10⁻³	45.13	9.52	1.7 × 10⁹
	Time interval (40, 60) h.	10⁻³	26.99	16.58	1.62 × 10⁹
	Time interval (60, 80) h.	10⁻³	87.80	21.77	1.33 × 10⁹
	Time interval (80, 100) h	10⁻³	68.42	14.43	5.74 × 10⁸
		Optimal value of Max [mAb](t), (mg/L)				1351.3	This paper
SP2^(a)(FBR) Optimal values ^(d)	Searching variables	F_L, L/h.	[GLC]_inlet, mM	[GLN]_inlet, mM	X_v,inlet, Cell/L		This paper
	Searching ranges	(10⁻⁴–5 × 10⁻²)	(25–150)	(5–25)	(2 × 10⁸–5 × 10⁹)
	Multi-dimensional optimization	Inlet optimal values of the FBR control variables
	Multi-dimensional optimization	F_L ^(b,c), L/h.	[GLC]_inlet mM	[GLN]_inlet mM	X_v,inlet, Cell/L	Max [mAb](t), (mg/L)
	Time interval (0, 20) h.	10⁻³	141.63	17.76	4.38 × 10⁹		This paper
	Time interval (20, 40) h.	10⁻³	55.81	9.52	4.20 × 10⁹
	Time interval (40, 60) h.	10⁻³	25.60	16.58	3.98 × 10⁹
	Time interval (60, 80) h.	10⁻³	126.92	21.77	3.21 × 10⁹
	Time interval (80, 100) h	10⁻³	94.62	14.43	1.20 × 10⁹
		Optimal value of Max [mAb](t), (mg/L)				6098.4	This paper
SP3^(a) (FBR) Optimal values ^(d)	Searching variables	F_L, L/h.	[GLC]_inlet, mM	[GLN]_inlet, mM	X_v,inlet, Cell/L		This paper
	Searching ranges	(10⁻⁴–5 × 10⁻²)	(25–150)	(5–25)	(2 × 10⁸–5 × 10⁹)
	Multi-dimensional optimization	Inlet optimal values of the FBR control variables
	Multi-dimensional optimization	F_L ^(b,c) L/h.	[GLC]_inlet mM	[GLN]_inlet mM	X_v,inlet, Cell/L	Max [mAb](t), (mg/L)
	Time interval (0, 50) h.	10⁻³	88.65	21.58	3.21 × 10⁹		This paper
	Time interval (50, 100) h.	10⁻³	137.97	20.65	1.2 × 10⁹		This paper
		Optimal value of Max [mAb](t), (mg/L)				5700.1	This paper

^(a) Time step -wise optimal operating policy of the FBR (plotted in Figure 3, Figure 4 and Figure 5). ^(b) The minimum feed flow rate of the inlet liquid F_L was set to be around 10%V_L,0/t_f = 10⁻³ L/h., or even below, to avoid excessive dilution of the bioreactor content [23]. The F_L(t) policy during the batch is to be optimized, being adjusted so that the final dilution of reactor content doesn’t exceed 10–25% of the initial liquid volume. ^(c) The initial liquid volume in the FBR (V_L,0), was adopted to 1 L, as for the BR case, that is 5× larger than that of [33]. ^(d) The optimal values refer to the inlet levels of the control variables for every time-arc, including the adjustable immobilized inlet X_t,inlet. In the X_t case, the initial value is X_t,0 = X_v,0. During the batch, the inlet X_t,inlet is taken 0.

Table 5. Substrate and biomass consumption, and realized performances by BR, and by the optimally operated FBR. The BR initial load, and the FBR optimal feeding policy are presented in Table 4. The equal time-arcs of the FBR are of 20 h. each for SP1 and SP2, and of 50 h. each for SP3. The BR, and FBR initial volume is of 1 L. The batch time is 100 h. in all cases, excepting for the last two lines.

Bioreactor Operation			Raw Material Consumption ^(b)			Reactor Performance Max [mAb](t), ^(b)		FBR Dilution
Type	N_div	Set-Point	Consumed GLC (mmoles)	Consumed GLN (mmoles)	X_v,0 (cells) ^(c)	(mg/L)	(mg/cells/h)	(%) ^(a)
BR	1	Nominal [44] SPBR ^(d)	29.1	4.9	2 × 10⁸	1254	6.3 × 10⁻⁸	0
FBR	5	Optimal SP1	14.22	3.23	ca.4 × 10⁸	1351	3.4 × 10⁻⁸	27
FBR	5	Optimal SP2	44.46	8.00	1.7 × 10⁹	6098	3.6 × 10⁻⁸	10
FBR	2	Optimal SP3	11.33	2.11	2.2 × 10⁸	5700	2.6 × 10⁻⁷	10
BR	1	[54] ^(e)			1–5 × 10⁸	~1100		0
FBR	7–13	[33] ^(f)			2 × (10⁸–10⁹)	~2400		?

^(a) Referring to the reactor liquid initial volume. ^(b) The larger number of displayed digits comes from the numerical simulations. ^(c) Referred to the FBR initial volume of 1 L (Table 3). ^(d) The BR nominal set-point (Table 3) of [44]. ^(e) The same cell culture, 120 h batch time (results from experimental plots). ^(f) The same cell culture, 168 h batch time (results from experimental plots); Initial [GLC] = 5.5–25 mM; Initial [GLN] = 3.74 mM; F_L,j = 0–12.5 mL/h; V₀ = 0.2 L.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

© 2020 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Maria, G. Model-Based Optimization of a Fed-Batch Bioreactor for mAb Production Using a Hybridoma Cell Culture. Molecules 2020, 25, 5648. https://doi.org/10.3390/molecules25235648

AMA Style

Maria G. Model-Based Optimization of a Fed-Batch Bioreactor for mAb Production Using a Hybridoma Cell Culture. Molecules. 2020; 25(23):5648. https://doi.org/10.3390/molecules25235648

Chicago/Turabian Style

Maria, Gheorghe. 2020. "Model-Based Optimization of a Fed-Batch Bioreactor for mAb Production Using a Hybridoma Cell Culture" Molecules 25, no. 23: 5648. https://doi.org/10.3390/molecules25235648

Article Menu

Model-Based Optimization of a Fed-Batch Bioreactor for mAb Production Using a Hybridoma Cell Culture

Abstract

1. Introduction

2. SBR Culture and Bioprocess Dynamics