Next Article in Journal
The Wind Dynamics of Super-Eddington Sources in FRADO
Previous Article in Journal
Statistical Mechanics of Long Walks in Dynamic Complex Networks: Statistical Arguments for Diversifying Selection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiobjective Optimization of a Fed-Batch Bienzymatic Reactor for Mannitol Production

1
Department of Chemical and Biochemical Engineering, Politehnica University of Bucharest, Polizu Str. 1-7, 011061 Bucharest, Romania
2
Romanian Academy, Chemical Sciences Section, Calea Victoriei 125, 010071 Bucharest, Romania
3
National Institute for Research and Development in Environmental Protection, 060031 Bucharest, Romania
*
Author to whom correspondence should be addressed.
Submission received: 7 July 2022 / Revised: 10 August 2022 / Accepted: 17 August 2022 / Published: 19 August 2022

Abstract

:
Enzymatic reactions can successfully replace complex chemical syntheses using milder reaction conditions and generating less waste. The developed model-based numerical analysis turned out to be a beneficial tool to determine the optimal operating policies of complex multienzymatic reactors. As proved, for such cases, the determination of a Fed-Batch Reactor (FBR) optimal operating policy results in a difficult multiobjective optimization problem. Exemplification is made for the bienzymatic reduction of D-fructose to mannitol by using MDH (mannitol dehydrogenase) and nicotinamide adenine dinucleotide (NADH) cofactor with the in situ continuous regeneration of NADH at the expense of formate degradation in the presence of FDH (formate dehydrogenase). For such a coupled system, the model-based engineering evaluations must account for multiple competing (opposable) optimization objectives. Among the multiple novelty elements: i) an optimally operated FBR with a tightly controlled variable feeding (of the time stepwise type) during the batch can lead to higher performance; ii) the optimally operated FBR reported better performance compared to an optimally single or cyclic BR, or to optimally serial batch-to-batch reactors (SeqBR), when considering a multiobjective optimization; iii) the concomitant variable feeding with substrate, enzymes, and cofactor during the FBR “time-arcs” is an option seldom approached in the literature but which is proved here, leading to consistent economic benefits.

1. Introduction

“Remarkable progresses made in the development of new enzymes and in realizing complex coupled enzymatic systems, able to in-situ recover the reaction cofactor(s), reported important applications in the industrial biocatalysis, presenting important advantages. Examples include the large number of biosynthesis processes used to produce fine -chemicals, or organic compounds in food, pharma, or detergent industry, such as: the production of monosaccharides derivatives, organic acids, alcohols, amino -acids, etc., by using single- or multi-enzymatic reactors” [1,2].
“Mannitol is a natural hexitol with important applications in medicine and the food industry [3,4]. The present global market of mannitol is around $100 million in 2013, of an average price of $42–80 per kg, and with a production growth rate of 5% annually. Around 50,000 tons/year of mannitol are produced currently by the costly chemical hydrogenation alone [5], and the rest by the less expensive enzymatic routes” [6]. In short, the main routes to produce mannitol at a large scale are the following [4,7]:
(i) The chemical catalytic hydrogenation of fructose, sucrose (inverted sugar), or of the 50% glucose–50% fructose syrup (HFCS) coming from the enzymatic hydrolysis of starch [1,8]. The chemical route is very costly because it requires the use of high pressures (50–80 atm), high temperatures (120–160 °C), and a costly Raney nickel or another catalyst [4,5,7,8].
(ii) Biological routes. HFCS conversion (83% yield) uses a cell culture of Candida magnoliae [9,10] with high conversions of glucose or sucrose by using a culture of Streptococcus Mutans [11]. Reviews are given by [4,12].
As reviewed by Slatner et al. [13], these routes suffer a large number of disadvantages, which include: (a) the yields are much lower than those of the enzymatic alternatives; (b) a significant fraction (15%) of the substrate is converted into allergenic byproducts such as lactate, acetate, ethanol, and carbon dioxide, leading to costly purifying steps of the product; (c) the used cell cultures are too sensitive to the operating conditions; (d) a too-high yield sensitivity vs. temperature is reported in the range of 25–40 °C.
(iii) The bienzymatic process with the NADH cofactor in in situ continuous regeneration. This very effective technology proposed by Slatner et al. [13] uses two coupled enzymatic reactions: (1) “the enzymatic reduction of fructose to mannitol in the presence of mannitol dehydrogenase (MDH) and the cofactor NADH as a proton donor. (2) The resulted NAD+ is continuously regenerated in-situ by the expense of the enzymatic decomposition of ammonium formate in the presence of format dehydrogenase (FDH)”, according to Figure 1. The use of another cofactor, such as NADPH, is not recommended, being much more expensive [14], and very unstable [15]. This route presents multiple advantages: (i) It is less expensive while requiring mild operating conditions (normal pressure, pH = 7, 25 °C); (ii) the selectivity is practically 100%, with the separation of the mannitol at the batch being easy and less costly; (iii) as proved in this paper, the optimal FBR operation can lead to the much-reduced consumption of costly enzymes (FDH, MDH, and NADH) compared to optimally single or cyclic BR, or to an optimally serial SeqBR; (iv) the continuous in situ enzymatic regeneration of the NADH cofactor largely reduces the production costs [4,16,17].
A rough evaluation of the enzymatic process performances compared to those of the biological processes (classical fermentation), and those of the chemical catalytic processes, was performed by Moulijn et al. [18] by pointing out their large number of advantages compared to the classical syntheses.
Over the last decades, remarkable progress in the development of new enzymes has been made by integrating genetic and engineering methods [19,20], and by realizing coupled multienzymatic systems in the industrial biocatalysis.
Multienzymatic systems with parallel or sequential reactions are modern and effective alternatives to the complex chemical syntheses by covering several choices [21]: (i). For two enzymes on the same support, the substrate for the second enzyme is generated in situ by the first reaction [22]; (ii). The second enzymatic reaction regenerates the cofactor of the first enzymatic one (e.g., the regeneration of NADH in the present study); (iii). The second enzymatic reaction shifts the equilibrium of the main reaction by removing the intermediate or byproduct from the system [23]; (iv). The second enzyme prolongs the life of the first enzyme by a particular mechanism (e.g., catalase prolongs the life of an oxidase [24]).
“Even if the multi-enzymatic systems are advantageous, the engineering part for optimizing such a process is not an easy task because it must account for the interacting enzymatic reactions, enzymes deactivation kinetics (if significant), multiple and often opposed optimization objectives, technological constraints, and uncertainties coming from multiple sources (model/constraints inaccuracies, disturbances in the control variables), and a highly nonlinear process dynamics [25,26,27,28,29,30]. All these parametric/model/data uncertainties require to update (with a certain frequency) the enzymatic process model, the optimal operating policies of the reactor being determined by using rather deterministic (model-based) optimization rules [31]. Multi-objective criteria, including economic benefits, operating and materials costs, product quality, etc., are used to off-line, or to on-line derive feasible optimal operating/control policies for various bioreactor types” [30] by using specific numerical algorithms [26,32,33,34,35].
This is why enzymatic reactors have been developed in simple constructive/operating alternatives (see some examples in Figure 2, and [27,28]). “In spite of their larger volumes, continuously mixing tank reactors (CSTR), operated in a BR (simple, cyclic, parallel) mode, or in a FBR mode, are preferred for processes requiring a high mass transfer” (a homogeneous medium for quickly interacting multienzymatic reactions), and a rigorous temperature/pH control. As revealed by Maria and Peptănaru [27], the SeqBR also proved to be an effective alternative compared to simple, cyclic (repeated), or parallel optimal BR-s. Regardless if free or immobilized enzymes are used enzymes, other effective complex operating alternatives proved their effectiveness, in particular, “those in which substrate(s)/cofactor and/or enzyme(s) are intermittently (BRP, [28,36]), or continuously ([FXBR (fixed-bed solid–-liquid continuous reactor), MA(S)CR (mechanically agitated solid–-liquid (semi-)continuous reactor))] added during the batch by following a suitable optimally variable policy (usually off-line, model-based determined). A similar discussion was provided by Koller [37].
In spite of inherent difficulties, trials to improve the BR productivity have been reported in the literature, by using SeqBR, cyclic BR, or parallel BR. Examples includes a sequence of 2–5 FBR-s [38,39,40], or a 14-times cyclic BR [41], or 48 parallel BR-s [42], or 24 parallel BR-s [43]. As expected, the optimized FBR reported better performances compared to the simple BR, due to its higher operating flexibility [28,44,45,46]. Concerning the optimal SeqBR, by varying the initial condition of its every BR of the series” [27], it reproduces somewhat the FBR operation principle, namely the time stepwise variation of the control variables [28]. “However, FBR policy differs in principle from those of a SeqBR, as long as the set of control variables are kept constant over a ’time-arc’, which is not the case for the SeqBR when the variation of control variables concerns only the initial load of each BR of the series. On the other hand, the chemical reactor theory demonstrated that a series of a large number of CSTR-s tends to have the best performances of a plug-flow reactor [47]. Consequently, it is expected that an optimal FBR, or an optimal SeqBR will present better performances than a repeated operation of an optimal single BR with the same initial conditions, even if optimal.”
By using the kinetic model of Maria [36], built-up and validated based on the large amount of experimental kinetic data of Slatner et al. [13] for the approached bienzymatic process, the present study aims at deriving the optimal operating policy of a bienzymatic FBR to concomitantly fulfill multiple opposite objectives, including: the minimization of the costly enzymes (MDH, FDH) and cofactor (NADH) consumption, concomitantly with mannitol production maximization under some technological constraints, as formulated in Section 3.3. The study is accompanied by a model-based analysis of the optimized FBR performances compared to those of an optimally single or cyclic BR, or to an optimally serial SeqBR, from a multiobjective perspective. As proved, the FBR with both a variable feed flow rate and species inlet concentrations reported superior performances vs. the other operating alternatives. The main gain in the optimized FBR performance compared to a single simple BR repeatedly (cyclic) operated, or to an optimized series of BRs (SeqBR), comes from the simultaneous variation of a larger number (4) of control variables, not only of the reactor feed flow rate (FL), but also of the inlet fructose concentration [F]inlet,j, as well as the inlet concentrations of the two enzymes [MDH]inlet,j, and [FDH]inlet,j (see the results presented in Section 4).
The price paid for achieving the best performance of the FBR by using a time stepwise variable operating policy is the need to previously prepare different substrates/enzymes/cofactor stocks to be continuously fed during every time-arc (that is, a batch time division in which the feeding is constant as flow rate and composition). The feeding characteristics differ between time-arcs, being determined by the FBR optimization vs. several criteria. This FBR optimal operation disadvantage is offset by its net higher productivity compared to those of a single BR, cyclic BR, or SeqBR by consuming less costly enzymes and cofactors. Similar conclusions are also reported in the literature for both enzymatic and biological processes [25,26,32,44,45,48].
This paper presents a significant number of novel aspects, as follows: (i) “The in silico (model-based) demonstrative engineering analysis of a complex bi-enzymatic process, leading to obtain an optimal operating policy of the approached FBR.” (ii) This study is one of the first systematic comparative analyses in the literature concerning the efficiency of various batch operating modes (BR, FBR, and SeqBR) from a multiobjective perspective and for a concrete multienzymatic case study. (iii) To our knowledge, there are very few studies optimizing an FBR by concomitantly considering a larger (4) number of control variables simultaneously varied during the batch. (iv) This bienzymatic process is already known, but the associated engineering analyses related to its optimal development/control are missing in the literature. (v) “The scientific value of this paper is not virtual, as long as the numerical analysis is based on a kinetic model of Maria [36] constructed and validated by using the extensive experimental data sets of Slatner et al. [13]. (vi) The in-silico analysis suggests that an optimally operated FBR with a relatively small number of time-arcs (10 here) can lead to high performances.” (vii) The concomitant variable feeding with the substrate/enzymes/cofactor during the FBR process is an option seldom approached in the literature. (viii) As proved in this paper, the model-based offline optimal operation of the FBRs can lead to consistent economic benefits.
Figure 2. Simplified schemes of various operating modes of the BR to conduct enzymatic or biological processes. See [28] for the complete modeling hypotheses of each BR type. BR = “isothermal, iso-pH, and iso-DO; perfectly mixed liquid phase. Reactants/biocatalyst(s)/additives added at the beginning of the batch only. BRP = a BR with reactants and/or biocatalyst(s)/supplements added during the batch in a Pulse-like addition of equal/uneven solution volumes, with a certain frequency (to be determined) [28,49]. FBR = A BR with substrates/biocatalyst(s)/supplements added during the batch following a certain (optimal) policy (to be determined) [28]. SBR = in continuous mode, there is a CSTR with a continuous feeding and evacuation with equal flow -rates. Substrates/biomass/enzymes /supplements are added during the batch following a certain (optimal) policy (to be determined). SeqBR = a series of (usually identical) BR-s. The batch-to-batch operation implies that the every BR content is transferred to the next BR, with adjusting the reactant(s) and/or biocatalyst(s) amounts (concentrations) to reach optimal levels” (determined offline, see [27]).
Figure 2. Simplified schemes of various operating modes of the BR to conduct enzymatic or biological processes. See [28] for the complete modeling hypotheses of each BR type. BR = “isothermal, iso-pH, and iso-DO; perfectly mixed liquid phase. Reactants/biocatalyst(s)/additives added at the beginning of the batch only. BRP = a BR with reactants and/or biocatalyst(s)/supplements added during the batch in a Pulse-like addition of equal/uneven solution volumes, with a certain frequency (to be determined) [28,49]. FBR = A BR with substrates/biocatalyst(s)/supplements added during the batch following a certain (optimal) policy (to be determined) [28]. SBR = in continuous mode, there is a CSTR with a continuous feeding and evacuation with equal flow -rates. Substrates/biomass/enzymes /supplements are added during the batch following a certain (optimal) policy (to be determined). SeqBR = a series of (usually identical) BR-s. The batch-to-batch operation implies that the every BR content is transferred to the next BR, with adjusting the reactant(s) and/or biocatalyst(s) amounts (concentrations) to reach optimal levels” (determined offline, see [27]).
Dynamics 02 00014 g002aDynamics 02 00014 g002b

2. Dynamic Models for the Bienzymatic Process and FBR

The approached BR in this analysis is that of Slatner et al. [13]. Its main characteristics, presented in (Table 1), reveal a quite flexible operating domain, including large ranges for the initial concentrations of substrates (fructose [F]o and formate [HCOO]o), cofactor [NADH]o, and enzymes ([MDH]o, [FDH]o). Such an observation opens large number of options when optimizing the FBR operating mode. The BR or FBR constructive scheme is presented in (Figure 3).
For the same reason, Slatner et al. [13] used this BR to study the kinetics of the here approached bienzymatic process. The extended experiments were conducted at 25 °C, pH 7.0, under large ranges of initial conditions, that is: [F]o ∈ [0.1–3] M; [NADH]o ∈ [0.008–0.5] M; [HCOO]o = [F]o; [NAD]o = 0.0005 M; {[MDH]o, [FDH]o} ∈ 0.1–2 kU/L. The collected kinetic data allowed Maria [36] to build up a kinetic model of the bienzymatic process. Based on these experimental datasets and other qualitative observations, a simple “Michaelis-Menten kinetic model of Ping-Pong-Bi-Bi type was proposed for both reaction rates R1 and R2 (Table 2) by analogy with a similar process of pseudo 2-nd order kinetics [24]. For simplicity, this model includes a non-competitive inhibition with respect to reactants even if the mannitol inhibition might be significant [13]. Enzymes MDH and FDH inactivation during the reaction have been neglected due to lack of available data. The rate constants have been estimated by using an effective nonlinear least squares procedure [51,52], with adopting a simple dynamic model for the BR [36]. The adequacy of the resulted kinetic model of (Table 2) was proved to be very good vs. the experimental data (not presented here), for the all tested large number of initial conditions” [36].
As revealed by Slatner et al. [13], for this approached process (Table 2), the pH is maintained constant at 7 ± 0.1 (optimum for MDH) by an “automatic addition of HCl, due to the shift in pH towards the alkaline region observed when FDH-catalyzes oxidation of ammonium formate. Besides, the reactor initial load contains a medium of Tris-HCI buffer of pH 7.0. Some other compounds (Dithiothreitol) are added to stabilize the MDH.”
To simulate the FBR dynamics in the present engineering analysis, the simple mathematical model of (Table 3) was adopted. The model includes the differential mass balances of the process key species, that is: F, HCOO, M, MDH, FDH, NADH, NAD, CO2, besides the liquid volume dynamics (VL) in the reactor. The FBR dynamic model is presented in the below Equation (2i) and Table 3.
This classic ideal model of the FBR assumes the following hypotheses [53]: (a) “isothermal, iso-pH; (b) additives (for the pH-control) are added initially and during the batch operation in recommended quantities; (c) perfectly mixed liquid phase (with no concentration gradients)” [28] and Table 3.
Table 1. Nominal reaction conditions of Slatner et al. [13] “for the enzymatic reduction of D-fructose to mannitol using MDH, and NADH cofactor in an experimental BR, with the in-situ continuous regeneration of the cofactor at the expense of formate degradation in the presence of FDH. The used FDH (EC 1.2.1.2) from Candida boidinii has a specific NAD+-dependent activity of 2.4 U/mg, measured at 25 °C and pH 7.0. The MDH (EC 1.1.1.67) from Pseudomonas fluorescens DSM 50106 was over-expressed in E. coli JM 109.” The values in parentheses are used in the present engineering calculations.
Table 1. Nominal reaction conditions of Slatner et al. [13] “for the enzymatic reduction of D-fructose to mannitol using MDH, and NADH cofactor in an experimental BR, with the in-situ continuous regeneration of the cofactor at the expense of formate degradation in the presence of FDH. The used FDH (EC 1.2.1.2) from Candida boidinii has a specific NAD+-dependent activity of 2.4 U/mg, measured at 25 °C and pH 7.0. The MDH (EC 1.1.1.67) from Pseudomonas fluorescens DSM 50106 was over-expressed in E. coli JM 109.” The values in parentheses are used in the present engineering calculations.
ParameterValues and Remarks [13]
Temperature/Pressure/pH (buffer solution)25 °C/Normal (1 atm.)/7
Molar initial concentrations
Fructose, [F]o0.1–1 M ([13])
(0.1–4 M) (this paper, in the range of kinetic model validity); (a)
[NADH]o0.008–0.5 M ([13])
(0.1–0.5 M) (this paper); 0.1 M (initial guess for optimal FBR)
[NAD+]o0.0005 M [13]
Formate initial [HCOO]oIdentical to initial [F]o
Others: [M]o = [CO2]o = 0None
[ C O 2 ] * = CO2 saturation level at 25 °C and pH = 7 ± 0.10.0313 M [54,55]
pH control by an “automatic addition of HCl” [13]
Batch   time   ( t f )48 h (this paper), and for an optimal single BR [27]
40 h for each BR in the SeqBR series of NBR = 10 BR [27]
Bioreactor   liquid   initial   volume   ( V L , 0 )0.5 L (adopted); variable due to the continuous feeding of the FBR.
Max. 3 L
The   feed   flow   rate   ( F L )0.02 L h−1 (initial). Its time stepwise dynamics are to be optimized in the adopted range of 0.01–0.04 L h−1 (limited by the reactor volume)
FDH enzyme (referred to the reactor liquid)Initial value = 0. Its time stepwise additions are to be optimized in the range of 0.1–2 kU/L.
MDH enzyme (referred to the reactor liquid)Initial value = 0. Its time stepwise additions are to be optimized in the range of 0.1–2 kU/L.
Some other compounds (dithiothreitol) are added by [13] to stabilize the MDH.
Footnotes: (a) The work at higher fructose concentrations is not recommended due to several reasons: lack of kinetic model validity; a higher reaction medium viscosity; limited solubility at 25 °C of fructose (4000 g/L) and of mannitol (ca. 200 g/L); a higher substrate inhibition effect.
Table 2. The kinetic model proposed by Maria [36] for the approached bienzymatic process. “The two coupled enzymatic reactions are: (R1) reduction of D-fructose to mannitol using MDH and NADH cofactor and, (R2) in-situ continuous regeneration of the cofactor NADH by the expense of formate degradation in the presence of FDH (Figure 1).” Rate constants have been estimated to match the experimental kinetic data of [13].
Table 2. The kinetic model proposed by Maria [36] for the approached bienzymatic process. “The two coupled enzymatic reactions are: (R1) reduction of D-fructose to mannitol using MDH and NADH cofactor and, (R2) in-situ continuous regeneration of the cofactor NADH by the expense of formate degradation in the presence of FDH (Figure 1).” Rate constants have been estimated to match the experimental kinetic data of [13].
Reactions:Rate Expressions:
F + N A D H   ( + H + ) M D H M + N A D + R 1 = k c 1 C M D H C F C N A D H K M 1 + K F C F + K N H C N A D H
H C O O + N A D + F D H C O 2 + N A D H R 2 = k c 2 C F D H C H C O O C N A D K M 2 + K H C C H C O O + K N A D C N A D
Rate constants k c 1 = 2 × 10−3; k c 2 = 8.3259 × 10−3; 1/h//(U/L); K M 1 = 7.2367 × 10−2 M; K M 2 = 8.8047 × 10−2 M;
K F = 1; K N H = 1; K H C = 5.0061 × 10−2; K N A D = 90.181 [27]
Process stoichiometry: d c F d t =   R 1   ;   d c N A D H d t =   R 1 + R 2 ; d c N A D d t =     d c N A D H d t ;
d c H C O O d t =   R 2 ; d c M d t = + R 1   ;   d c C O 2 d t = + R 2
Table 3. “Key-species mass balances in the FBR model, by including the enzymatic process kinetic model together with the associated rate constants of (Table 2).” The ideal model hypotheses of Maria [28] assumes a homogeneous liquid composition with negligible mass transport resistance in the bulk phase.
Table 3. “Key-species mass balances in the FBR model, by including the enzymatic process kinetic model together with the associated rate constants of (Table 2).” The ideal model hypotheses of Maria [28] assumes a homogeneous liquid composition with negligible mass transport resistance in the bulk phase.
SpeciesRemarks
Substrate (fructose, F) and cofactor (NADH) mass balances:
d c F d t = F L , j V L c F , i n l e t , j c F     R 1
c F , i n l e t , j = control   variable ;   j = 1 , . . , N d i v time stepwise unknown values to be determined from the FBR optimization;
c F , 0 = c F , i n l e t , 1
d c N A D H d t = F L , j V L c N A D H , i n l e t , j c N A D H                           R 1 + R 2
c N A D H , i n l e t , j = control   variable ;   j = 1 , . . , N d i v time stepwise unknown values to be determined from the FBR optimization;
c N A D H , 0 = c N A D H , i n l e t , 1
For the optimal FBR with adopted Ndiv = 10, the feeding policy is (a):
c F , i n l e t c N A D H , i n l e t = c F , i n l e t , 1 c N A D H , i n l e t , 1 if 0 t T 1 c F , i n l e t , 2 c N A D H , i n l e t , 2 if T 1 t T 2 c F , i n l e t , 3 c N A D H , i n l e t , 3 if T 2 t T 3 c F , i n l e t , 4 c N A D H , i n l e t , 4 if T 3 t T 4 c F , i n l e t , 5 c N A D H , i n l e t , 5 if T 4 t T 5 c F , i n l e t , 6 c N A D H , i n l e t , 6 if T 5 t T 6 c F , i n l e t , 7 c N A D H , i n l e t , 7 if T 6 t T 7 c F , i n l e t , 8 c N A D H , i n l e t , 8 if T 7 t T 8 c F , i n l e t , 9 c N A D H , i n l e t , 9 if T 8 t T 9 c F , i n l e t , 10 c N A D H , i n l e t , 10 if T 9 t T 10 = t f
Formate (HCOO) and exhausted cofactor (NAD+) mass balances:
d c H C O O d t = F L , j V L c H C O O , i n l e t , j c H C O O     R 2
d c N A D d t = F L , j V L c N A D , i n l e t , j c N A D + R 1     R 2
j = 1,.., N d i v time stepwise values
c H C O O , i n l e t , j = 0   ,     j
c H C O O , 0 = c H C O O   ( t = 0 ) = c F , 0

c N A D , i n l e t , j = 0   ,     j
c N A D , 0 = c N A D   ( t = 0 ) = 5   ·   10 4 M [13]
Mannitol (M) and CO2 mass balances:
d c M d t = F L , j V L c M , i n l e t , j c M + R 1
d c C O 2 d t = F L , j V L c C O 2 , i n l e t , j c C O 2 + R 2
j = 1,.., N d i v time stepwise values
c M , i n l e t , j = 0   ,     j
c M , 0 = c M (t = 0) = 0.

c C O 2 , i n l e t , j = 0   ,     j
c C O 2 , 0 = c C O 2 (t = 0) = 0.
If   c C O 2     c C O 2 * ,   then   c C O 2 = c C O 2 * ; (b)
Enzymes MDH and FDH mass balances (c):
d c M D H d t = F L , j V L c M D H , i n l e t , j c M D H  
c M D H , i n l e t , j = control   variable ;   j = 1 , . . , N d i v time stepwise unknown values to be determined from the FBR optimization;
c M D H , 0 = 0.
d c F D H d t = F L , j V L c F D H , i n l e t , j c F D H  
c F D H , i n l e t , j = control   variable ;   j = 1 , . . , N d i v time stepwise unknown values to be determined from the FBR optimization;
c F D H , 0 = 0.
For the optimal FBR with adopted Ndiv = 10, the feeding policy is (a):
c M D H , i n l e t c F D H , i n l e t = c M D H , i n l e t , 1 c F D H , i n l e t , 1 if 0 t T 1 c M D H , i n l e t , 2 c F D H , i n l e t , 2 if T 1 t T 2 c M D H , i n l e t , 3 c F D H , i n l e t , 3 if T 2 t T 3 c M D H , i n l e t , 4 c F D H , i n l e t , 4 if T 3 t T 4 c M D H , i n l e t , 5 c F D H , i n l e t , 5 if T 4 t T 5 c M D H , i n l e t , 6 c F D H , i n l e t , 6 if T 5 t T 6 c M D H , i n l e t , 7 c F D H , i n l e t , 7 if T 6 t T 7 c M D H , i n l e t , 8 c F D H , i n l e t , 8 if T 7 t T 8 c M D H , i n l e t , 9 c F D H , i n l e t , 9 if T 8 t T 9 c M D H , i n l e t , 10 c F D H , i n l e t , 10 if T 9 t T 10 = t f
Liquid volume in the reactor (d):
d V L d t = F L , j
F L , j = control   variable ;   j = 1 , . . , N d i v time stepwise unknown values to be determined from the FBR optimization. The unknown
F L , 0 = F L   ( t = 0 ) = F L , 1
is   determined   together   with   the   all   F L , j values from the FBR optimization.
For the optimal FBR with adopted Ndiv =10, the feeding policy is (a): F L = F L , 1   if   0   t   <   T 1   F L , 2   if   T 1   t   <   T 2 F L , 3   if   T 2   t   <   T 3 F L , 4   if   T 3   t   <   T 4 F L , 5   if   T 4   t   <   T 5 F L , 6   if   T 5   t   <   T 6 F L , 7   if   T 6   t   <   T 7 F L , 8   if   T 7   t   <   T 8 F L , 9   if   T 8   t   <   T 9 F L , 10   if   T 9   t   <   T 10 = t f  
Footnotes: (a) For the adopted Ndiv = 10, the j = 1,…, N d i v time-arcs approx. switching points are: T1 = 4.8 h.; T2 = 9.6 h.; T3 = 14.4 h.; T4 = 19.2 h.; T5 = 24 h.; T6 = 28.8 h.; T7 = 33.6 h.; T8 = 38.4 h.; T9 = 43.2 h.; T10 = t f = 48 h. (b) The saturation level of CO2 is given in Table 1. (c) Enzymes deactivation is neglected. (d) The F L , j time stepwise feed flow rates are determined simultaneously with the other control variables (Section 3.1) to ensure the FBR optimal operation (see Section 3.2).

3. Optimization Problem Formulation

By using a step-by-step approach, this section presents the numerical methodology used to solve the FBR optimization problem, that is: the math translation of the formulated objectives, the problem constraints, and the used numerical rules. All these aspects are discussed in detail by intentionally using Section 3.1 to Section 3.6, and a regular, but simple, (bio)chemical engineering language, while keeping the necessary math rigor.

3.1. The Control Variable Choice

Due to its improvements, the approached bienzymatic alternative proposed by Slatner et al. [13] is the most promising and cheaper technology to produce mannitol at an industrial scale [6,7,13]. However, due to the costly enzymes, reflected in the product cost, one of the sensitive issues to be solved is an engineering one, related to the chosen reactor operating strategy able to minimize the enzymes (MDH, FDH) consumption but keeping the mannitol high yields ([M]final) for the imposed technological limits for the F/HCOO substrates and NADH cofactor.
Recent advances try to bypass this problem by using the biological route, which couples the two reactions of (Figure 1) “in the same genetic modified micro-organism (Bacillus megaterium) used as host for both enzymes synthesis, and cofactor regeneration” [12]. However, this last route is not yet available at an industrial scale due to the numerous disadvantages related to the control of the biological reactor.
In the present case, such problems related to cell cultures are inexistent. Additionally, to greatly increase the FBR flexibility and efficiency, a larger number of control variables can be used. By analyzing the species relative importance, the following five control variables have been chosen:
-
The fed enzymes ( c M D H , i n l e t , j ; c F D H , i n l e t , j ),
-
The fed substrate ( c F , i n l e t , j ), and the fed cofactor ( c N A D H , i n l e t , j );
-
The feed flow rate ( F L , j ) of the solution containing the substrates, enzymes, and the cofactor.
In the above notations, the index “j = 1, Ndiv” denotes the time stepwise value to be determined by optimization for every of Ndiv “time-arc” in which the FBR batch-time is divided.
As an important observation, the optimal FBR operating strategy to be determined in this study is more complex than the simple operation of a BR or of other reactors from Figure 2. Mainly, the substrate (F), enzymes (MDH, FDH), and the cofactor (NADH) concentrations in the feeding solution of the variable feed flow rate (FL) are no longer kept constant as well. Thus, in the FBR case, the most convenient optimal operation is the following: (i) “the batch-time is divided in Ndiv (equal time-arcs) of equal lengths, and (ii) the control variables are kept constant only over every time-arc at optimal values determined from solving an optimization problem (e.g., maximization of the mannitol (M) production in this case). The time-intervals of equal lengths Δt = tf /Ndiv are obtained by dividing the batch time tf into Ndiv parts tj−1 ≤ t ≤ tj, where tj = j·Δt are switching points (where the reactor input is continuous and differentiable)”. Time-intervals for the present case study with an adopted Ndiv = 10 are shown in footnote (a) of Table 3. Consequently, there are 5 × Ndiv = 50 unknown in total to be determined by applying an optimization procedure.
For this case, initial concentrations of the other compounds are null (that is, [M]o, and [CO2]o), or at values recommended by Slatner et al. [13] (that is [NAD]o, and [HCOO]o = [F]o, see Table 3).
Here, it is to observe that there are several other operating variables that could be considered when optimizing the FBR (e.g., [HCOO]o ≠ [F]o, or even the batch time ( t f ) imposed by Slatner et al. [13]). However, as the number of control variables increases, the associated nonlinear optimization problem, subjected to nonlinear constraints (NLP), becomes more multimodal, and even more difficult to solve, by greatly increasing the computation time to find the global feasible optimum of the NLP problem.
The practical implementation of such a variable feeding strategy (as the feed flow rate and composition) over the batch assumes the previous preparation of the stock solutions with different concentrations for the substrate (F), enzymes (MDH, FDH), and cofactor (NADH), at levels which resulted from the optimization calculations to be fed for every time-arc (that is a batch-time division in which the feeding is constant as flow-rate and composition). Of course, every stock solution presents species (F, NDH, MDH, FDH) concentrations that are different from each other.
In the present case of the FBR operation adopting Ndiv = 10 (equal time-arcs), it results in the requirement to prepare in advance 10 stock solutions, each with the concentrations of F, NADH, MDH, and FDH to be determined from the optimization calculations (Section 4). Then, each of these stock solutions are fed into the FBR according to the optimal feed flow rate determined (by optimization, Section 4) for each time-arc. Of course, the feedings of the 10 stock solutions during the 10 time-arcs are different, not only in species concentrations (F, NADH, MDH, FDH), but also in flow rates.

3.2. Optimization Objectives and the Composite Objective Function (W) Choice

To explore the efficiency limits of the studied FBR of Table 1 by testing a larger number of economical objectives, the following goals are selected:
(i)
maximum FBR productivity in mannitol;
(ii)
minimum consumption of the costly enzymes MDH and FDH;
(iii)
minimum consumption of the cofactor NADH—even if NADH is not very expensive [14], its separation over the product purification is relatively costly [56];
(iv)
an adjustable substrate c F , i n l e t , j consumption (within defined limits) to ensure its maximum conversion to mannitol.
As revealed by Slatner et al. [13], the enzymatic reduction of D-fructose to mannitol by using MDH and NADH as a cofactor takes place very selectively and without byproducts. So, the fructose reduction occurs quantitatively. Thus, here the fructose conversion is the most influencing indicator for FBR productivity (Fobj1 in Equation (1B)).
In mathematical terms, the above multiobjectives (i–iv) problem translates in the following NLP problem:
Given   [ NAD ] o ,   Find   the   control   variables : F L , j ;   c M D H , i n l e t , j ;   c F D H , i n l e t , j ;   c F , i n l e t , j ;   c N A D H , i n l e t , j ,   For   j   =   1 , ,   N div ,   with   the   adopted   N div   =   10   time - arcs   and   the FBR   initial   condition   of   Table   3 ,   so   as   to   obtain :   Min   W ( c ,   c o ,   k ) , with   the   following   composite   objective   function : W ( c ,   c o ,   k ) = Min ( 1 / Fobj 1 )     Min ( Fobj 2 )     Min ( Fobj 3 )     Min ( Fobj 4 )
“The time-arcs are of equal length Δt = tf /Ndiv, being obtained by dividing the batch time tf into Ndiv parts tj−1 ≤ t ≤ tj, where tj = j·Δt are the switching points between arcs. Notations: Λ = ‘and’ (simultaneously), in the math sense; c = species concentration vector; co = initial value of c; k = vector of the kinetic model rate constants”; W = the overall (composite) objective function. The components Fobj1–Fobj4 of the composite objective function are the following:
Fobj 1 = M t f ,   with   [ M ]   in   M   units . Fobj 2 = j = 1 N d i v c M D H , i n l e t , j   F L , j   Δ t j ,   with   MDH   conc .   in   kU / L   units . Fobj 3 = j = 1 N d i v c F D H , i n l e t , j   F L , j   Δ t j ,   with   FDH   conc .   in   kU / L   units . Fobj 4 = j = 1 N d i v c N A D H , i n l e t , j   F L , j   Δ t j ,   with   NADH   conc .   in   M   units . Δ t j = t f   /   N d i v ,   for   all   j ( that   is ,   equal   time - arcs   in   the   present   case   study )
Equation (1A,B) can be re-written in a form much more accessible to the numerical calculus, as follows:
Given   [ NAD ] o ,   Find   the   control   variables : F L , j ;   c M D H , i n l e t , j ;   c F D H , i n l e t , j ;   c F , i n l e t , j ;   c N A D H , i n l e t , j , For   j   =   1 , ,   N div ,   with   the   adopted   N div   =   10   time - arcs   and   the FBR   initial   condition   of   Table   3 ,   so   as   to   obtain : Min   W ( c ,   c o ,   k ) , with   the   composite   objective   function : W ( c ,   c o ,   k )   =   ( Fobj 2   +   Fobj 3   +   Fobj 4 ) / Fobj 1 Where   the   components   Fobj 1 - Fobj 4   are   defined   in   the   Equation   ( 1 B )
As an observation, the use of the overall objective function ’W’ in Equation (1C) corresponds to the so-called ‘inverted utility function method’ [57]. Because the four components Fobj1–Fobj4 of W are opposed (that is, the realization of one objective is done to the detriment of another objective), other optimization procedures can be applied as well, for instance the Pareto front method [26]. However, the adopted multiobjective optimization rule with this composite objective function presents the advantage of simplicity, with an easy application and interpretation.
Concerning the adopted overall objective function ‘W’ in (1C) some observations are necessary: (i) “The chosen units for M, MDH, FDH, and NADH allow the direct comparison of Fobj1–Fobj4, and their concomitant use in ‘W’, because they present the same order of magnitude; (ii) The way by which “W was elaborated, implicitly ensures the simultaneous realization of the all objectives formulated in Equation (1A), that is consumption minimization of the two enzymes (MDH, FDH) and of the cofactor (NADH), together with the mannitol M production maximization. (iii) Of course, other formulations of the ‘W’ function are also possible, depending on the weight (relative importance) given to each component objective (Fobj1–Fobj4) (see the example [58])”.

3.3. Optimization Problem Constraints

The NLP optimization problem Equation (1C) is highly nonconvex and nonlinear, being subjected to the following technological/physical meaning/model constraints:
Nonlinear process and reactor model:
d c j d t = F L , u V L   c j , i n l e t , u c j + i = 1 n r ν i j r i   ( FBR   nonlinear   dynamic   model   of   Table   3 ) ;   J = species   index   ( F ,   M ,     HCOO ,   NADH ,   NAD + ,   CO 2 ,   MDH ,   FDH ) ;   ν i j = the   stoichiometric   coefficient   of   the   species   j   in   the   i - th   reaction ; r i = the   i - th   reaction   rate ;   u = 1 , ,   N div   denotes   the   time - arc   number   of   the   FBR   optimal   feeding   policy ;   n r = no .   of   reactions . d V L d t = F L , u ,   the   liquid   volume   dynamics   of   ( Table   3 ) , With   the   initial   conditions   of   ( Table 3 ) ,   that   is : c j , o = c j t = 0 ;   where   j = ( F ,   NADH ,   MDH ,   FDH )   are   to   be   optimized ; c j , o = 0 ,   for   j = ( M ,   CO 2 ) ,   and   for   j = ( NAD + )   are   given   in   ( Table   3 ) , c j , o   for   j = ( HCOO )   is   [ HCOO ] o = [ F ] o ,
Physical significance constraints:
c j ( t ) 0 ,   for   all   the   species   j ,   and   for   all   t     [ 0 - t f ] ,  
Searching ranges for the control variables are those suggested and experimentally validated by Slatner et al. [13], as below indicated. It no to note that [HCOO]o = [F]o, as recommended by [13].
[ c M D H , i n l e t , u ;   c F D H , i n l e t , u ]     [ 0.1 2 ]   kU / L ;   c F , i n l e t , u     [ 0.1 4 ]   M ;   c N A D H , i n l e t , u     [ 0.1 0.5 ]   M ;   F L , u     [ 0.01 0.04 ]   L / h .   ( due   to   the   FBR   capacity ) .  
u= 1,…, Ndiv denotes the time-arc number of the FBR optimal feeding policy.
V L = variable   ( see   2 i ) ;   ( FBR   liquid   volume ) ;
One excludes the unfeasible trivial solution:
W = Fobj1 = Fobj2 = Fobj3 = Fobj4 = 0

3.4. Time Stepwise Divisions Choice

The simplest FBR operating mode is to keep constant during the batch the control variables selected in the Section 3.1. It is understood that such a policy will not necessarily ensure the optimal FBR performances. That is because it lacks of adapting the feeding dynamics to the needs of the reactions progress.
Actually, the optimal variable control policy “implies a time step-wise variable feeding of the bioreactor for all the selected control variables, over an adopted (Ndiv = 10 here) equal time-arcs that covers the whole batch time. In the present study, the time-arcs are of equal length Δt = tf /Ndiv, being obtained by dividing the batch time tf into Ndiv parts tj−1 ≤ t ≤ tj, where tj = j·Δt are the switching points between arcs. Each time-arc “j”, of length Δtj (j = 1, Ndiv) is characterized by optimal levels (to be determined) of the control variables,” that is: the feed flow rate F L , j , the substrate, and the main cofactor ( c F , i n l e t , j ; c N A D H , i n l e t , j ), and the two enzymes ( c M D H , i n l e t , j ; c F D H , i n l e t , j ), (see Equation (1C)).
Besides seeking the optimal operation of Equation (1C), this type of FBR variable feeding operation “includes enough degrees of freedom to offer a wide range of FBR operating facilities that, in principle, might be investigated (see also the discussion given by” [28,59], as, for instance:
(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 (Ndiv) to obtain a more refined and adaptable FBR operating policy, but keeping the same non-uniform feeding policy for the chosen control variables.
(d)
by extend the search min/max limits of the control variables, or even by considering them as unknown (to be determined).
(e)
by feeding the bioreactor with a variable feed flow-rate, but with a solution of an uniform concentration for the all other control variables, over a small or large number (Ndiv) of time-arcs.
All the alternatives (a–e) are not approached here from the following reasons:
  • Alternatives (a–c) are not good options, because as (Ndiv) increases, the necessary computational effort to solve the NLP problem {Equation (1C) + Equation (2i–v)} grows significantly due to considerable increase in the number of searching variables. Thus, the quick (real-time) implementation of the derived FBR operating policy becomes questionable. Additionally, multiple optimal operating policies can exist for the resulted over-parameterized constrained NLP optimization problem, thus increasing the difficulty to quickly locate a feasible globally optimal solution of the FBR optimization problem.”
The alternative (c) is also not advisable because, “as the (Ndiv) increases, the operating policy is more difficult and more costly to be implemented, since the optimal feeding policy requires a larger number of stocks with feeding substrate/NADH/enzymes solutions of different concentrations, separately prepared in advance to be fed for every time-arc Δtj (j = 1,…, Ndiv) of the FBR operation. Also, the NLP optimization problem is more difficult to solve because of the increased number of operating parameters to be determined from locating the feasible minimum of the multi-modal composite objective function—a challenging problem due to the multiple solutions difficult to discriminate and evaluate. Besides, FBR operation with using a too large number (Ndiv) of small time-arcs (Δtj) can raise special operating problems when including PAT (Process Analytical Technology) tools [60].
A brief survey of the FBR optimization literature [27,28,46,61] reveals that a small number (Ndiv) ≤ 10 is commonly used for such FBR due to the above-mentioned reasons. In fact, the present numerical analysis does not intend to exhaust all the possibilities of the approached FBR optimization. Thus, an extended analysis of the operating Alternatives (a)–(d) of the FBR operation, or the influence of the parametric uncertainty deserves a separate investigation, being beyond the scope of this paper. To not complicate the computational analysis, a value of (Ndiv) = 10 equal time-arcs have been adopted here”, leading to equal time-arcs of length t f /(Ndiv) = 4.8 h.
  • The alternative (d) is unlikely because it cannot limit the excessive consumption of substrates and cannot prevent the hydrodynamic stress or the excessive product dilution. Moreover, such an alternative might indicate unrealistic results outside the validity of the used kinetic model [36], whose plausibility is strictly restricted to the investigated experimental domain (Table 1). “In our numerical analysis, carefully documented upper bounds of control variables Equation (2iii) were tested to ensure realistic results for the FBR optimal operating policy.”
  • “The alternative (e) is also not feasible, even if a larger (Ndiv) will be used (see (c)). That is because, it is well-known that the variability of the FBR feeding over the batch time-arcs is the main degree-of-freedom used to obtain FBR optimal operating policies of superior quality” [25,28,45,59]. By giving up to the variable feed flow rate, the enzymes, substrate, and cofactor concentrations, suboptimal FBR operating policies will be obtained and of low performance.

3.5. The Used Solver and the Problem Solution Particularities

During the FBR optimization, when simulating the FBR dynamics, the time evolution of the species concentrations in the bulk phase and of the liquid–volume are obtained by solving the FBR dynamic model of Table 3 “with the initial condition c j , o (t = 0) stipulated in Equation (2i). The imposed batch time t f , and the optimal medium conditions are those of (Table 1). The dynamic model solution was obtained with a high precision, by using the variable-order stiff integrator (ode15s) of the MATLAB™.” All calculations were performed in the framework of the MATLAB™ computational platform with an original code developed by the authors.
Because the FBR differential model of Table 3, the optimization objective Equation (1C), and the problem constraints Equation (2i–v) are all highly nonlinear, the formulated problem of Equation (1C) with the explicit/implicit nonlinear constraints Equation (2i–v) “translates into a nonlinear optimization problem (NLP) with a multi-modal composite objective function, and a non-convex searching domain. To obtain the global feasible solution with enough precision, by avoiding local sub-optimal solutions, the multi-modal optimization solver MMA of Maria [51,52,62] has been used, as being proved in previous works to be more effective compared to the common (commercial) algorithms. The computational time was reasonably short (minutes) by using a common PC, thus offering a quick implementation of the obtained FBR optimal operating policy.”

3.6. Optimized BR and SEQBR Policies

To further compare the FBR optimal operating policy with equivalent batch operations, two alternatives have been chosen from the literature, those who analyzed the same bienzymatic process of mannitol production, respectively, namely:
(a)
The optimization of the initial load of a simple single BR of Maria [27,36]. The optimal BR policies have been derived by using an equivalent composite multiobjective function seeking the minimum enzyme (MDH,FDH) consumption simultaneously with realizing a maximum mannitol production under the multiple constraints similar to Equation (2i–v). To realize the equivalence in the direct comparison terms of the present FBR with (Ndiv) = 10 equal time-arcs, a N = 10× times repeated (cyclic) BR (of the same size) and an optimized series SeqBR of NBR = 10 equal BR-s (of the same size) are considered below in the discussion (Section 4).
(b)
The optimization of a SeqBR with a series of NBR = 10 BR-s (Figure 2) [27]. At the end of each batch, the content of the BR from the series is passed to the next BR but adjusting its initial load (that is [F]o, [NADH]o, [MDH]o, [FDH]o) according to the optimal operating strategy derived by using an equivalent composite multiobjective function realizing the minimum enzyme (MDH,FDH) consumption simultaneously with the maximum mannitol production under multiple constraints similar to Equation (2i–v). For a direct comparison to the analyzed FBR with (Ndiv) = 10 equal time-arcs, NBR = 10 equal BR-s (of the same size) are considered below in the discussion (Section 4).
“To determine the optimal BR policy (a), and the optimal SeqBR policy (b), the same numerical rules of (Section 3.5) have been used, that is the variable-order stiff integrator (“ode15s”) of the MATLAB™, and the multi-modal optimization solver MMA of Maria [51,52,62], together with the authors’ orginal code, all being implemented on the MATLAB™ computational platform. To confront the BR optimization result, a simple exhaustive numerical algorithm has also been used by Maria [36]. As revealed by Maria and Peptănaru [27], optimal SeqBR is superior to the optimal (cyclic) BR, and to the experimental BR trials of Slatner et al. [13] in terms of enzymes consumption (3×–12× less for FDH, and MDH),” to obtain the same high fructose conversions.

4. FBR Optimization Results and Discussion

4.1. The Resulted FBR Optimal Operating Policy

The optimal operating policy of the FBR, obtained by solving the optimization problem Equation (1C) with the control variables of Section 3.1, and the constraints Equation (2i–v), for an adopted Ndiv =10 (Section 3.4), is presented as follows:
(a)
In Figure 4A, the simulated key species concentrations dynamics, that is c j ( t ) in Equation (2a).
(b)
In Figure 4B, the variable policy of the feed flow rate ( F L , j , j = 1,…, Ndiv) and the liquid volume dynamics.
(c)
In Figure 4C, the variable feeding policy with the substrate ( c F , i n l e t , j , j = 1,…,Ndiv), and with NADH ( c N A D H , i n l e t , j , j = 1,…, Ndiv),
(d)
In Figure 4D, the variable feeding policy with MDH enzyme ( c M D H , i n l e t , j , j=1,…,Ndiv).
(e)
In Figure 4E, the variable feeding policy with FDH enzyme ( c F D H , i n l e t , j , j=1,…,Ndiv).
Figure 4. (A) Simulated dynamic trajectories of species concentrations in the optimally operated FBR with the time stepwise feeding of figures (BE). (B) The optimal feed flow rate (FL) policy and the liquid volume (VL) increase during the batch. (C) The optimal feeding policy with the substrate (F) and cofactor NADH. (D) The optimal feeding policy with the MDH enzyme and its accumulation in the reactor. (E) The optimal feeding policy with the FDH enzyme and its accumulation in the reactor. Initial conditions and the control variables ranges are given in Table 1 and Table 3.
Figure 4. (A) Simulated dynamic trajectories of species concentrations in the optimally operated FBR with the time stepwise feeding of figures (BE). (B) The optimal feed flow rate (FL) policy and the liquid volume (VL) increase during the batch. (C) The optimal feeding policy with the substrate (F) and cofactor NADH. (D) The optimal feeding policy with the MDH enzyme and its accumulation in the reactor. (E) The optimal feeding policy with the FDH enzyme and its accumulation in the reactor. Initial conditions and the control variables ranges are given in Table 1 and Table 3.
Dynamics 02 00014 g004aDynamics 02 00014 g004b
“It is to observe that, due to the above formulated optimization problem, the FBR optimal operating policy is given for each of the time-intervals (of equal length) uniformly distributed over the batch-time. It is also to underline that, as revealed by (Figure 4A) the simulated key-species concentrations time-trajectories are monotonous functions, despite their time step-wise variable addition in the FBR, thus revealing the reaction efficiency.
Such an optimal time step-wise variable feeding of the bioreactor presents advantages and inherent disadvantages. The advantages are related to the higher flexibility of the FBR operation, leading to a higher productivity in mannitol as proved in the below Section 4.2. Beside, the imposed limits of the control variables prevent excessive substrate consumption, or an excessive reactor content dilution, or violation of validity limits of the process kinetic model.
As a disadvantage, the FBR-s with such a time-variable control is more difficult to be operated compared to the simple BR, as long as the time step-wise optimal feeding policy requires different stocks of feeding substrates/cofactor/enzymes solutions of different concentrations prepared to be used over the batch. This is the price paid for achieving FBR best performances. This need to previously prepare different solutions 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 below discussed, and pointed-out in the literature [26,27,45,46,58,59]. In fact, the best operating alternative (FBR vs. BR) is related to many others economic factors (operating policy implementation costs, product cost compared to its production costs, product price fluctuation, etc.)”, not discussed here.
The present biochemical engineering numerical analysis is of an in silico type, being based on the dynamic model of the bienzymatic process and of the enzymatic reactor. Even if we have no experimental data with which to compare our optimization results, their credentials (i.e., a valid optimal operating policy for the studied FBR) are high enough, as long as the kinetic model is enough adequate, because:
(i)
The used kinetic model was built-up by Maria [36] by using the large number of experimental kinetic data sets of Slatner et al [13] (Section 2).
(ii)
The experimental program conducted by Slatner et al. [13] covers a large domain for the control variables (mentioned in Section 2).
However, if after implementing such an optimal operating policy, “significant inconsistencies are observed between the model-predicted reactor dynamics and the experimental evidence, the numerical optimization stage is applied again, after an intermediate step (between batches) necessary to off-line improve the model adequacy (by updating the model parameters)” [28].

4.2. Comparative Discussion of the Optimization Results

By analyzing the resulted optimal FBR feeding policy, and the species dynamics in the reactor bulk presented in (Figure 4A–E), several observations are to be made:
(i)
“There is a close connection between the coupling reactions, enzyme concentrations, and the supra-unitary NADH/NAD+ ratio over the batch (Figure 4A). This means that the two reactions are well coupled, [NADH] being higher than [NAD+] for most of the time, that is a NADH/NAD+ ratio higher than 1 most of the time, thus maintaining a process efficiency [36].
(ii)
The cofactor NADH regeneration is very efficient, the formate decomposition being quasi-complete and leading to saturation [CO2]* = 0.0313 M (Table 1) in short time (earlier than 10 h in Figure 4A), with removal of the CO2 excess from the system over the rest of the batch-time.”
(iii)
“As revealed by repeated simulations of Maria [36] (not presented here), the BR performances are more sensitive to [MDH] and [NADH] than to [FDH].”
(iv)
Bulk [MDH] grows continuously due to its continuous feeding (Figure 4D). By contrast, [FDH] evolution in the bulk phase is more special (Figure 4E). It grows continuously due to the continuous feeding with FDH until a certain time (10 h). After this time, [FDH] in bulk declines because its reduced feeding rate cannot compensate the reactor content dilution rate.
On the other hand, optimal FBR performances can be compared with those of the 10× repeated optimal BR, and with those of an optimized SeqBR with NBR = 10 BR in series, as reported by Maria and Peptănaru [27] for the same process (see Section 3.6). Such a direct comparison is presented in Table 4. Species consumption was evaluated by using the following relationship:
species   consumption = j = 1 N d i v F L , j × [ s p e c i e s ] i n l e t , j × Δ t
In the Equation (3), the resulted consumption units are the following: “moles” for j = F, NADH, and “kU” for j = MDH, FDH. The hourly consumptions in Table 4 refer to the overall batch time stipulated in the column “remarks”.
By comparatively analyzing the optimal operating policies resulted from the optimal operation of the repeated BR, FBR, and SeqBR over a significant range of initial/fed [F]o, and [NADH]o, presented in Table 4 and by Maria and Peptănaru [27], several conclusions can be derived, as follows:
(v)
As a general conclusion, the optimally operated FBR reported better performance compared to the optimal SeqBR, which, in turn, reported better performance than the repeated optimally operated BR. Thus, for a comparable mannitol (M) production: (a) the MDH enzyme consumption of FBR is much smaller, that is 10x (vs. SeqBR) and 130x (vs. BR); (b) the FDH enzyme consumption in FBR is also smaller at 9x (vs. SeqBR) and 15–50x (vs. BR); (c) the NADH consumption in FBR is smaller by 2–9x (vs. SeqBR) and 1.5–7x (vs. BR); (d) compared to the experimental (nonoptimal) BR trials of Slatner et al. [13] (within the range of [F]o = 0.1–1 M, [NADH]o = 0.25–1 mM; [NAD+]o = 0.0005 M; [HCOO]o = [F]o; [ c M D H , 0 ; c F D H , 0 ] ∈ [0.1–2] kU/L), the optimal FBR reported much smaller enzyme consumptions, that is, 10x for MDH and 4.5x for FDH.
(vi)
As revealed by Figure 4D,E, it is not necessary to keep the enzymes feeding at a high level over the whole batch. On the contrary, an initial high addition over the first 1–2 time-arcs, followed by a continuous addition of a low level, are enough to ensure the necessary enzyme concentrations so to result in the required high reaction rates. The same observation is also valid for the NADH feeding policy.
(vii)
One of the better FBR performance costs is the final 2× dilution of the reactor content (Figure 4B). However, balancing the product dilution with the considerable savings in terms of the expensive enzymes and the regenerated cofactor, the use of FBR is worth the investment in its (model-based) offline optimal control.
(viii)
The relatively large time-arcs (of ca. 5 h) of the FBR allow an easy operation control over each time-arc, but also when switching between them.
(ix)
As proved by the present numerical analysis, “solving this engineering problem by using an experimental procedure, as tried by Slatner et al. [13], may not be the best choice because it involves high costs due to the required large number of batch tests”, even if at a small lab-scale.
(x)
In the FBR variable operating case, conversion is more difficult to interpret because the substrate is continuously added to the initial one. Due to this FBR operation particularity, the substrate consumption cannot be referred solely to its initial concentration but rather to the added substrate during the batch, following the variable concentration of Figure 4C. This is why the substrate F consumption is computed with the Equation (3) formula. In the present case, the formal conversion is high, that is: (transformed F in M)/(totally fed F) = 72.5%. The lack of substrate near total conversion is not a critical issue as long as there are lot of solutions to solve the product separations problems, that is, the use: (i) of immobilized enzymes and (ii) of efficient technologies to separate saccharides, as is well-presented in the literature [63,64,65].

5. Conclusions

“A relevant case study analyzed in this paper proves that, for the coupled multi-enzymatic systems, derivation of the multi-objective optimal operating policies (minimum enzyme/cofactor consumption, with maximum reactor productivity) for an enzymatic reactor is not a trivial engineering problem” [36,66,67,68].
Even if it is case dependent, the present model-based numerical analysis indicates that, for multienzymatic systems, the performance of an optimally operated FBR is much better in terms of enzyme and cofactor consumption (10x–130x less MDH, 5x–50x less FDH, and 7–9x less NADH) compared to the optimized repeated (cyclic) BR, or to optimal serial batch-to-batch reactors (SeqBR) from the same multiobjective perspective. The optimized FBR also presents much better performance compared to the experimental BR’s large number of costly and exhaustive trials, as performed by Slatner et al. [13], to obtain a high fructose conversion.
The use of an adequate process model can offer an approximate, if not exact, solution to the reactor multiobjective optimization problem with a moderate computational effort. Such an in silico engineering approach allows saving considerable experimental effort at the expense of a computationally intensive step.
The numerical methodology presented here can be repeated with any supplementary problems if the kinetic models for enzymatic inactivation are available. The same rule can also be applied to derive multiobjective operating policies to other enzymatic reactor types (BRP, FXBR, MA(S)CR, [28]) with free or immobilized enzymes by using suitable reactor dynamic models.
This numerical analysis can also be extended to the cases where the control variables are not deterministic (as here) but which present a random (known) distribution centered on the set point (see [26] for an example).

Author Contributions

Conceptualization, G.M.; Investigation, G.M., L.R. and C.M.; Methodology, G.M.; Project administration, G.M.; Software, G.M.; Supervision, G.M.; Validation, G.M.; Data curation, L.R. and C.M.; Formal analysis, L.R. and C.M.; Writing—original draft, G.M., L.R. and C.M.; Writing—review & editing, G.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no funding of any kind.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

c j “Species j concentration
c j * Species j saturation level
FLFeed flow rate of the liquid
kj, k, Kj, k c 1 , k c 2 , KMj Rate constants
n r Number of reactions
NdivNumber of time-arcs of equal lengths in which the batch time is uniformly divided; on every time-arc, the control variables are kept constant
r j , r i , R1, R2Species j reaction rate, or the i-th reaction rate
t Time
t f Final batch time
VLThe BR liquid volume
WThe overall (composite) objective function of the optimization problem Equation (1A,C)
Greek Symbols
ΔtTime interval
ν i j stoichiometric coefficient of species j in the reaction i
Index
o, 0Initial
inletReferring to the FBR feeding
F, finalFinal
Abbreviations
BR“batch reactor
BRPBR with intermittent addition of enzymes/substrate(s)/additives
Eenzyme
CSTRContinuous stirred tank reactor
FD-Fructose
FBRFed-batch reactor
FDHFormate dehydrogenase
FXBRFixed-bed solid-liquid continuous reactor
Fobj1-Fobj4Objective functions defined in eqn.(1B)
HFCSfructose/glucose syrup
HCOOFormate ion
MMannitol, or molar concentration
MA(S)CRMechanically agitated solid-liquid (semi-)continuous reactor
MDHMannitol dehydrogenase
MMAThe numerical algorithm of Maria [51,52] to solve optimization problems
Min/Maxminimum/maximum
NAD(P)Hnicotinamide adenine dinucleotide (phosphate)
NAD, NAD+Nicotinamide adenine dinucleotide (oxidized form) »
NLPNonlinear optimization problem, subjected to nonlinear constraints
NBRNumber of equal BR included in the SeqBR series
SeqBRSequential batch-to-batch reactor
SBRSemi-batch reactor”
T1 …T10Switching times between time-arcs
[X]Concentration of the species X
Λ« and » (simultaneously), in the math sense.

References

  1. Liese, A.; Seelbach, K.; Wandrey, C. Industrial Biotransformations; Wiley-VCH: Weinheim, Germany, 2006. [Google Scholar]
  2. Straathof, A.J.J.; Adlercreutz, P. Applied Biocatalysis; Harwood Academic Publ.: Amsterdam, The Netherlands, 2005. [Google Scholar]
  3. Wikipedia, “Mannitol”. 2021. Available online: https://en.wikipedia.org/wiki/Mannitol (accessed on 7 August 2021).
  4. Saha, B.C.; Racine, F.M. Biotechnological production of mannitol and its application. Appl. Microbiol. Biotechnol. 2011, 89, 879–891. [Google Scholar] [CrossRef] [PubMed]
  5. Bhatt, S.M.; Mohan, A.; Srinivastava, S.K. Challenges in enzymatic route of mannitol production. ISRN Biotechnol. 2013, 2013, 914187. [Google Scholar] [CrossRef] [PubMed]
  6. Von Weymarn, N. Process Development for Mannitol Production by Lactic acid Bacteria. Ph.D. Thesis, Helsinki University of Technology, Laboratory of Bioprocess Engineering, Beijing, China, 2002. Available online: http://lib.tkk.fi/Diss/2002/isbn9512258854/ (accessed on 7 August 2021).
  7. Ghoreishi, S.M.; Shahrestani, R.G. Innovative strategies for engineering mannitol production. Trends Food Sci. Technol. 2009, 20, 263–268. [Google Scholar] [CrossRef]
  8. Makkee, M.; Kieboom, A.P.G.; Van Bekkum, H. Production methods of D-mannitol. Starch/Staerke 1985, 37, 136–141. [Google Scholar] [CrossRef]
  9. Khan, A.; Bhide, A.; Gadre, R. Mannitol production from glycerol by resting cells of Candida magnoliae. Bioresour. Technol. 2009, 100, 4911–4913. [Google Scholar] [CrossRef] [PubMed]
  10. Song, K.H.; Lee, J.K.; Song, J.Y.; Hong, S.G.; Baek, H.; Kim, S.Y.; Hyun, H.H. Production of mannitol by a novel strain of Candida magnoliae. Biotechnol. Lett. 2002, 24, 9–12. [Google Scholar] [CrossRef]
  11. Loesche, W.J.; Kornman, K.S. Production of mannitol by Streptococcus mutans. Arch. Oral Biol. 1976, 21, 551–553. [Google Scholar] [CrossRef]
  12. Bäumchen, C.; Roth, A.H.F.J.; Biedendieck, R.; Malten, M.; Follmann, M.; Sahm, H.; Bringer-Meyer, S.; Jahn, D. D-Mannitol production by resting state whole cell biotransformation of D-fructose by heterologous mannitol and formate dehydrogenase gene expression in Bacillus megentarium. Biotechnol. J. 2007, 2, 1408–1416. [Google Scholar] [CrossRef]
  13. Slatner, M.; Nagl, G.; Haltrich, D.; Kulbe, K.D.; Nidetzky, B. Enzymatic production of pure D-mannitol at high productivity. Biocatal. Biotransform. 1998, 16, 351–363. [Google Scholar] [CrossRef]
  14. Chenault, H.K.; Whitesides, G.M. Regeneration of nicotinamide cofactors for use in organic synthesis. Appl. Biochem. Biotechnol. 1987, 14, 147–197. [Google Scholar] [CrossRef]
  15. Maria, G.; Ene, M.D. Modelling enzymatic reduction of 2-keto-D-glucose by suspended aldose reductase. Chem. Biochem. Eng. Q. 2013, 27, 385–395. Available online: http://silverstripe.fkit.hr/cabeq/assets/Uploads/Cabeq-2013-04-01.pdf (accessed on 30 December 2021).
  16. Peacock, D.; Boulter, D. Kinetic studies of formate dehydrogenase. Biochem. J. 1970, 120, 763–769. [Google Scholar] [CrossRef] [PubMed]
  17. Schmidt, T.; Michalik, C.; Zavrel, M.; Spiess, A.; Marquardt, W.; Ansorge-Schumacher, M.B. Mechanistic model for prediction of formate dehydrogenase kinetics under industrially relevant conditions. Biotechnol. Prog. 2010, 26, 73–78. [Google Scholar] [CrossRef] [PubMed]
  18. Moulijn, J.A.; Makkee, M.; van Diepen, A. Chemical Process Technology; Wiley: New York, NY, USA, 2001. [Google Scholar]
  19. Buchholz, K.; Hempel, D.C. From Gene to Product (Editorial). Eng. Life Sci. 2006, 6, 437. [Google Scholar] [CrossRef]
  20. Hempel, D.C. Development of biotechnological processes by integrating genetic and engineering methods. Eng. Life Sci. 2006, 6, 443–447. [Google Scholar] [CrossRef]
  21. Xue, R.; Woodley, J.M. Process technology for multi-enzymatic reaction systems. Bioresour. Technol. 2012, 115, 183–195. [Google Scholar] [CrossRef] [PubMed]
  22. Blanch, H.W.; Clark, D.S. Biochemical Engineering; Marcel Dekker: New York, NY, USA, 1997. [Google Scholar]
  23. Hames, B.D.; Hooper, N.M.; Houghton, J.D. Instant Notes in Biochemistry; Springer Verlag: New York, NY, USA, 1997. [Google Scholar]
  24. Maria, G.; Ene, M.D.; Jipa, I. Modelling enzymatic oxidation of D-glucose with pyranose 2-oxidase in the presence of catalase. J. Mol. Catal. B Enzym. 2012, 74, 209–218. [Google Scholar] [CrossRef]
  25. 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/Heidelberg, Germany, 2008; pp. 299–317. [Google Scholar]
  26. Maria, G.; Crisan, 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]
  27. Maria, G.; Peptănaru, I.M. Model-based optimization of mannitol production by using a sequence of batch reactors for a coupled bi-enzymatic process–A dynamic approach. Dyn.-Basel 2021, 1, 134–154. [Google Scholar] [CrossRef]
  28. Maria, G. Model-based optimization of a fed-batch bioreactor for mAb production using a hybridoma cell culture. Mol.-Basel–Org. Chem. 2020, 25, 5648. [Google Scholar] [CrossRef]
  29. Bonvin, D. Optimal operation of batch reactors—A personal view. J. Process Control. 1998, 8, 355–368. [Google Scholar] [CrossRef]
  30. 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]
  31. Srinivasan, B.; Bonvin, D.; Visser, E.; Palanki, S. Dynamic optimization of batch processes: II. Role of measurements in handling uncertainty. Comput. Chem. Eng. 2003, 27, 27–44. [Google Scholar] [CrossRef]
  32. 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, 14–18 August 2005; p. 20. [Google Scholar]
  33. 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. [Google Scholar]
  34. 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]
  35. 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]
  36. Maria, G. Model-based optimisation of a batch reactor with a coupled bi-enzymatic process for mannitol production. Comput. Chem. Eng. 2020, 133, 106628–106635. [Google Scholar] [CrossRef]
  37. Koller, M. A review on established and emerging fermentation schemes for microbial production of polyhydroxyalkanoate (PHA) biopolyesters. Fermentation 2018, 4, 30. [Google Scholar] [CrossRef]
  38. Irvine, R.L.; Busch, A.W. Sequencing batch biological reactors–An overview. J. Water Pollut. Control. Fed. 1979, 51, 235–243. Available online: http://0-www-jstor-org.brum.beds.ac.uk/stable/25039819 (accessed on 20 January 2022).
  39. Shacham, M. Semicontinuous fed-batch and cyclic-fed batch operation. In Proceedings of the Workshop “Modern Problem Solving Techniques in Engineering with POLYMATH, Excel and MATLAB”, Tel-Aviv, Israel, 23 September 2008; pp. 1439–1441, Chapter 14.13. Available online: http://www.eng.tau.ac.il/~brauner/Workshop08Participant/Example-7/Prob-14-13.pdf (accessed on 23 February 2021).
  40. US EPA. Wastewater Technology Fact Sheet: Sequencing Batch Reactors; OWM, Water Permits Division, Municipal Branch: Washington, WC, USA, 2012. Available online: https://www3.epa.gov/npdes/pubs/sbr_new.pdf (accessed on 23 February 2021).
  41. Rätze, K.H.G.; Jokiel, M.; Sundmacher, K. Cyclic operation of a semibatch reactor for the hydroformylation of long-chain olefins and integration in a continuous production process. Chem. Eng. J. 2019, 377, 120453. [Google Scholar] [CrossRef]
  42. 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]
  43. Hans, S.; Ulmer, C.; Narayanan, H.; Brautaset, T.; Krausch, N.; Neubauer, P.; Schäffl, I.; Sokolov, M.; Bournazou, M.N.C. Monitoring parallel robotic cultivations with online multivariate analysis. Processes 2020, 8, 582. [Google Scholar] [CrossRef]
  44. Maria, G.; Crisan, M. Evaluation of optimal operation alternatives of reactors used for D-glucose oxidation in a bi-enzymatic system with a complex deactivation kinetics. Asia-Pac. J. Chem. Eng. 2014, 10, 22–44. [Google Scholar] [CrossRef]
  45. 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. Jl. Biotechnol. 2012, 10, 263–269. Available online: https://www.researchgate.net/publication/288966207_Comparison_between_batch_and_fed-batch_production_of_rhamnolipid_by_Pseudomonas_aeruginosa (accessed on 9 January 2022).
  46. 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]
  47. Froment, G.F.; Bischoff, K.B. Chemical Reactor Analysis and Design; Wiley: New York, NY, USA, 1990. [Google Scholar]
  48. Franco-Lara, E.; Weuster-Botz, D. Estimation of optimal feeding strategies for fed-batch bioprocesses. Bioprocess Biosyst. Eng. 2005, 28, 71. [Google Scholar] [CrossRef]
  49. Maria, G. Model-based heuristic optimised 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]
  50. Wikipedia, “Bioreactor”. 2021. Available online: https://en.wikipedia.org/wiki/Bioreactor (accessed on 30 December 2021).
  51. Hamza, M.H. (Ed.) Modelling, Identification and Control; IASTED/ACTA Press: Anaheim, CA, USA, 2003; pp. 112–118. Available online: https://www.actapress.com/Content_of_Proceeding.aspx?proceedingID=213 (accessed on 7 March 2021).
  52. 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] [CrossRef]
  53. Moser, A. Bioprocess Technology: Kinetics and Reactors; Springer: New York, NY, USA, 1988. [Google Scholar]
  54. Carroll, J.J.; Mather, A.E. The system carbon dioxide-water and the Krichevsky-Kasarnovsky equation. J. Solut. Chem. 1992, 21, 607–621. [Google Scholar] [CrossRef]
  55. Reid, R.C.; Prausnitz, J.M.; Poling, B.E. The Properties of Gases and Liquids, 4th ed.; McGraw-Hill: Boston, MA, USA, 1987. [Google Scholar]
  56. Dutta, R. Fundamentals of Biochemical Engineering; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  57. Rao, S.S. Engineering Optimization—Theory and Practice; Wiley: New York, NY, USA, 2009; Chapter 14.10. [Google Scholar]
  58. Scoban, A.G.; Maria, G. Model-based optimization of the feeding policy of a fluidized bed bioreactor for mercury uptake by immobilized P. putida cells. Asia-Pac. J. Chem. Eng. 2016, 11, 721–734. [Google Scholar] [CrossRef]
  59. Maria, G.; Renea, L. Tryptophan production maximization in a fed-batch bioreactor with modified E. coli cells, by optimizing its operating policy based on an extended structured cell kinetic model. Bioeng.-Basel 2021, 8, 210. [Google Scholar] [CrossRef]
  60. 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 21 October 2021).
  61. Loeblein, C.; Perkins, J.; Srinivasan, B.; Bonvin, D. Performance analysis of on-line batch optimization systems. Comput. Chem. Eng. 1997, 21, S867–S872. [Google Scholar] [CrossRef]
  62. Maria, G. Adaptive random search and short-cut techniques for process model identification and monitoring. In Proceedings of the FOCAPO98 International Conference on Foundations of Computer Aided Process Operations, Snowbird, UT, USA, 5–10 July 1998; pp. 351–359. [Google Scholar]
  63. Gillarová, S.; Henke, S.; Svoboda, T.; Kadlec, P.; Hinková, A.; Bubník, Z.; Pour, V.; Sluková, M. Chromatographic separation of mannitol from mixtures of other carbohydrates in aqueous solutions. Czech J. Food Sci. 2021, 39, 281–288. [Google Scholar] [CrossRef]
  64. Blignaut, J.; Albataineh, K.; Banat, Y.; Abu El-Rub, Z. Comparison of Adsorbent Behavior in Glucose/Fructose Separation by Simulated Moving Bed (SMB) Chromatography, Res. Report; University of Twente: Enschede, The Netherlands, 2001. [Google Scholar]
  65. Lee, K.N. Continuous Separation of Glucose and Fructose at High Concentration Using Two-Section Simulated Moving Bed Process. Korean J. Chem. Eng. 2003, 20, 532–537. [Google Scholar] [CrossRef]
  66. Rapaport, A.; Dochain, D. Minimal time control of fed-batch processes with growth functions having several maxima. IEEE Trans. Autom. Control. 2011, 56, 2671–2676. Available online: https://hal.archives-ouvertes.fr/hal-00529001/document (accessed on 15 January 2022). [CrossRef]
  67. Gajardo, P.; Ramirez, H.C.; Rapaport, A.E. Minimal time sequential batch reactors with bounded and impulse controls for one or more species. SIAM J. Control. Optim. 2008, 47, 2827–2856. [Google Scholar] [CrossRef]
  68. Flores, V.H.R. Optimal Control Problems for Bioremediation of Water Resources. Ph.D. Thesis, Universidad de Chile, Santiago, Chile, 2016. Available online: https://repositorio.uchile.cl/bitstream/handle/2250/141336/Problemas-de-control-optimo-para-la-biorremediacion-de-recursos-acuiferos.pdf?sequence=1&isAllowed=y (accessed on 15 January 2022).
Figure 1. Simplified “reaction scheme of the two coupled enzymatic reactions: (top) D-fructose (F) reduction to mannitol (M) by using suspended MDH (mannitol dehydrogenase), and the cofactor NADH (Nicotinamide adenine dinucleotide). (Down) NADH cofactor continuous regeneration by the expense of formate (HCOO) degradation in the presence of suspended FDH (Formate dehydrogenase) [13]. The use of NADPH cofactor is not recommended, being much more expensive [14], and very unstable” [15].
Figure 1. Simplified “reaction scheme of the two coupled enzymatic reactions: (top) D-fructose (F) reduction to mannitol (M) by using suspended MDH (mannitol dehydrogenase), and the cofactor NADH (Nicotinamide adenine dinucleotide). (Down) NADH cofactor continuous regeneration by the expense of formate (HCOO) degradation in the presence of suspended FDH (Formate dehydrogenase) [13]. The use of NADPH cofactor is not recommended, being much more expensive [14], and very unstable” [15].
Dynamics 02 00014 g001
Figure 3. “Simplified scheme of a BR or of a FBR used to conduct enzymatic or biological processes. In the BR operating mode, substrate(s), biocatalyst, and additives are initially loaded in the recommended amounts (concentrations). In the FBR operating mode, the substrate(s)/biocatalyst(s) and additives (nutrients, pH-control substances) are continuously fed, following a certain (optimal) policy.” Adapted from [50].
Figure 3. “Simplified scheme of a BR or of a FBR used to conduct enzymatic or biological processes. In the BR operating mode, substrate(s), biocatalyst, and additives are initially loaded in the recommended amounts (concentrations). In the FBR operating mode, the substrate(s)/biocatalyst(s) and additives (nutrients, pH-control substances) are continuously fed, following a certain (optimal) policy.” Adapted from [50].
Dynamics 02 00014 g003
Table 4. Comparison of the operating alternatives [27,36] of the studied “bi-enzymatic reduction of D-fructose to mannitol with the continuous in-situ regeneration of the NADH cofactor”.
Table 4. Comparison of the operating alternatives [27,36] of the studied “bi-enzymatic reduction of D-fructose to mannitol with the continuous in-situ regeneration of the NADH cofactor”.
Reactor TypeF consumption
(Moles/h)
M
Production (M/h)
MDH Consumption
(kU/L)
FDH Consumption
(kU/L)
NADH Consumption
(moles/h)
Remarks
Optimal BR (a,c)0.021–0.062 (c)0.017–0.054 (c)1–13.88 (c)3–10 (c)0.002–0.010 (c)10 BR cycles of
48 h each
Optimal SeqBR (b)0.025–0.075 (b)0.012–0.044 (b)1.1121.7550.0025–0.0125 (b)Series of 10 BR,
40 h run for each BR
Optimal FBR (this paper)0.03620.02620.0882 kU/(0.5–1.2 L FBR)0.1686 kU/(0.5–1.2 L FBR)0.0014Single batch of 48 h
(a) BR with a batch time of t f = 48 h [27,36]; (b) optimal SeqBR with search intervals corresponding to initial load of [F]o = 1–4 M and [NADH]o = 0.1–0.5 M for every BR of the series [27]. Search intervals are [ c M D H , 0 ; c F D H , 0 ] ∈ [0.1–2] kU/L. The interval corresponds to an initial load of [F]o = 1–3 M and [NADH]o = 0.1–0.5 M. All other initial loads (that is [NAD+]o, [HCOO]o = [F]o) are the same as those of Table 1. (c) The interval corresponds to an initial load of [F]o = 0.1–3 M, [NADH]o = 0.008–0.5 M; [NAD+]o = 0.0005 M; [HCOO]o = [F]o. Search intervals are [ c M D H , 0 ; c F D H , 0 ] ∈ [0.1–2] kU/L [27,36].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maria, G.; Renea, L.; Maria, C. Multiobjective Optimization of a Fed-Batch Bienzymatic Reactor for Mannitol Production. Dynamics 2022, 2, 270-294. https://0-doi-org.brum.beds.ac.uk/10.3390/dynamics2030014

AMA Style

Maria G, Renea L, Maria C. Multiobjective Optimization of a Fed-Batch Bienzymatic Reactor for Mannitol Production. Dynamics. 2022; 2(3):270-294. https://0-doi-org.brum.beds.ac.uk/10.3390/dynamics2030014

Chicago/Turabian Style

Maria, Gheorghe, Laura Renea, and Cristina Maria. 2022. "Multiobjective Optimization of a Fed-Batch Bienzymatic Reactor for Mannitol Production" Dynamics 2, no. 3: 270-294. https://0-doi-org.brum.beds.ac.uk/10.3390/dynamics2030014

Article Metrics

Back to TopTop